Study of Differential Genes Expression of Uterine tissue Related to EggShell Using Transcriptomics data

پژوهشهای علوم دامی ایران

Article 9, Volume 15, Issue 1 - Serial Number 53, March 2023, Pages 107-121   

Document Type: Research Articles

DOI: 10.22067/ijasr.2022.73912.1055

Authors

Milad Gholami Tahoone^* ; Hossein Moradi Shahrebabak

Department of Animal Science, College of Agriculture & Natural Resources, University of Tehran, Karaj, Iran.

Abstract

Introduction: Due to the high importance of egg shell quality in terms of hardness, extensive transportation with high volume of this product between different countries, the importance of shell quality in terms of protecting the internal content from extra-shell contamination and also one of the ways to improve traits, relying on The genetic structure of a trait is the application of genetic selection. Therefore, the aim this study is to use transcriptome data obtained from uterine tissue of two groups of laying hens with hard and weak eggshell, in each of the three stages of the calcification period (initiation, growth and termination). Gene expression profile was determined and differential expression of genes is analyzed so that different index genes can be expressed and their ontology can be studied to help the results. As a result, the genetic structure of the trait and the list of index genes have been regularly completed, and with the information obtained from the present study, the correct and practical genetic selection can be made in laying hen breeds.
Materials and Methods: The aim of this study was to investigate the gene expression profile between two groups of chickens with hard and weak eggshells in three stages of calcification period (initiation, growth and termination), to identify significant index genes with different expression and analysis. Ontology and the paths involved are related to them. Therefore, in this study, transcripts (total mRNA sequences) of 18 samples of uterine tissue including three replicates in each stage of each group were prepared. The process of converting the initial format of readings was done using SRAtoolkit software, measuring the quality of readings with fastQC software, editing readings with Trimmomatic software and re-measuring the quality of readings after editing with fastQC software; Also, the alignment and mapping of the readings were done using Hisat2 software and format conversion and sorting of results was performed using Samtools software. Finally, the differential expression analysis of genes was performed using Cufflinks software package (Cufflink, Cuffmerge and Cuffdiff software) and Then, considering p-value <0.00025 and the criterion for comparing the value of differential gene expression, log2 (fold_change) less than -4 and greater than +4 were identified as the most significant genes with different and significant expression. Histogram of the results was performed using CummeRbund package in R software environment and ontological analysis was performed using DAVID and Ensembl databases.
Results and Discussion: According to the results of gene expression profile, 28833 genes were identified on the transcript of the samples. Finally, by analyzing the differential expression of genes, 1338 genes at the initiation stage, 81 genes at the growth stage and 190 genes at the termination stage were identified with significantly different expressions, that Considering log2(fold_change) > +4 and log2(fold_change) < -4, 51 genes at the initiation stage, two genes at the growth stage and four genes at the termination stage had different and significant expression. The largest difference in gene expression between the two groups at the initiation stage, ENSGALG0000000044418 gene with log2 (fold_change) equal to 8/63 and at the growth stage, ENSGALG0000000049618 gene with log2 (fold_change) equal to 6/82 and at the termination stage, ENSGALG0000000049618 gene with log2 (fold_change) equal to 5/14 was observed. Ontological analysis of index genes showed that they are mainly involved in protein binding activities (ENSGAL000000444418, NMRAL1, SLITRK4, CHL1, ENSGAL00000009041, PHACTR3 and CORIN genes), DNA transcription regulation (LHX1, TFAP2A, RUNX3, ETV1, NR4A3 and RUNX1 genes), immune response (ENSGAL0000000043996, TMEM117, AvBD9, GCNT3, ENSGAL00000046947 and ENSGAL00000051617 genes), fat metabolism (PNPLA3, NR4A3, CD36, FFAR4 and DGKB genes) and calcium ion binding (ANXA10, CAPN8, CDH18 and DGKB genes), respectively.
Conclusion: In the differential expression analysis of genes, considering log2 (fold_change) more than +4 and less than -4, in the initial stage, 51 genes, in the growth stage, two genes and in the termination stage, four genes were identified. The highest differences in gene expression between the two groups were observed in the initial stage, ENSGALG00000044418 and PDZK1IP1 genes, in the growth stage, ENSGALG00000049618, ENSGALG00000048945 and TMEM63C genes and in the termination stage, ENSGALG000000004400 and ENSG26 genes; The results of ontological analysis of index genes showed that they are mainly involved in protein binding, DNA transcription regulation, immune response, lipid metabolism and calcium ion binding, respectively.

Keywords

Calcification; Differential gene expression; Eggshell; Gene ontology; Transcriptome

References

1. Abdelmanova, A. S., Dotsev, A. V., Romanov, M. N., Stanishevskaya, O. I., Gladyr, E. A., Rodionov, A. N., Vetokh, A. N., Volkova, N. A., Fedorova, E. S., Gusev, I. V., Griffin, D. K., Brem, G., & Zinovieva, N. A. (2021). Unveiling comparative genomic trajectories of selection and key candidate genes in egg-type russian white and meat-type white cornish chickens. Biology, 10(9), 876. DOI: 10.3390/biology10090876
2. Andrews, S. (2010). FastQC: A quality control tool for high throughput sequence data. From: http://www.bioinformatics.babraham.ac.uk/projects/fastqc
3. Bahrami, A. (2020). Which aligner software is the best for our study?. Genetics Research, 7, 048. DOI: 10.23937/2378-3648/1410048
4. Banabazi, M. H., Naserkheil, M., & Miraei-Ashtiani, R. (2012). Network regulatory gene expreesion of cerevisiae cell cycle with correlation weight method. The third national conference on agricultural biotechnology iran. Ferdowsi University of Mashhad, Iran. (In Persian).
5. Banabazi, M. H., Naserkheil, M., & Miraei-Ashtiani, R. (2012). Algorithm to identify genes expressed differently in the yeast Saccharomyces cerevisiae microarray data using software packages R. The third national conference on agricultural biotechnology Iran. Ferdowsi University of Mashhad, Iran. (In Persian).
6. Bolger, A. M., Lohse, M., & Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics, 30, 2114-2120. DOI: 1093/bioinformatics/btu170
7. Ekblom, R., & Galindo, J. (2010). Applications of next generation sequencing in molecular ecology of non-model organisms. Heredity, 107, 1-15. DOI:1038/hdy.2010.152
8. (2020). FAOSTAT. https://www.fao.org/faostat/en/#data/QCL/visualize/
9. Gan, Q., Chepelev, I., Wei, G., Tarayrah, L., Cui, K., Zhao, K., & Chen, X. (2010). Dynamic regulation of alternative splicing and chromatin structure in Drosophila gonads revealed by RNA-seq. Cell Research, 20 (7), 763-783. DOI: 1038/cr.2010.64
10. Github. (2017). Cufflinks. http://cole-trapnell-lab.github.io/cufflinks/
11. Haas, B. J., & Zody, M. C. (2010). Advancing RNA-seq analysis. Nature Biotechnology, 28, 421. DOI: 1038/nbt0510-421
12. Jonchere, V., Rehault-Godbert, S., Hennequet-Antier, C., Cabau, C., Sibut, V., Cogburn, L. A., Nys, Y., & Gautron, J. (2010). Gene expression profiling to identify eggshell proteins involved in physical defense of the chicken egg. BMC Genomic, 11, 57. DOI: 1186/1471-2164-11-57
13. Kim, D., Langmead, B., & Salzberg, S. L. (2015). HISAT: A fast spliced aligner with low memory requirements. Nature Methods, 12, 357-360. DOI: 1038/nmeth.3317
14. Marguerat, S., & Bahler, J. (2010). RNA-seq: From technology to biology. Cellular and Molecular Life Sciences, 67, 569-579. DOI: 1007/s00018-009-0180-6
15. Mortazavi, A., Williams, B., McCue, K., Schaeffer, L., & Wold, B. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods, 5(7), 621-628. DOI: 1038/nmeth.1226
16. Sultan, M., Schulz, M. H., Richard, H., Magen, A., Klingenhoff, A., Scherf, M., & Schmidt, D. (2008). A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science, 321(5891), 956-960. DOI: 1126/science.1160342
17. Wilhelm, B. T., & Landry, J. R. (2009). RNA-seq quantitative measurement of expression through massively parallel RNA-sequencing. Methods, 48, 249-257. DOI: 1016/j.ymeth.2009.03.016
18. Wolc, A., Drobik-Czwarno, W., Jankowski, T., Arango, J., Settar, P., Fulton, J. E., Fernando, R. L., Garrick, D. J., & Dekkers, J. C. M. (2020). Accuracy of genomic prediction of shell quality in a white leghorn line. Poultry Science, 99, 2833-2840. DOI: 10.1016/j.psj.2020.01.019