1
|
Li L, Xu X, Xiao M, Huang C, Cao J, Zhan S, Guo J, Zhong T, Wang L, Yang L, Zhang H. The Profiles and Functions of RNA Editing Sites Associated with High-Altitude Adaptation in Goats. Int J Mol Sci 2023; 24:3115. [PMID: 36834526 PMCID: PMC9964554 DOI: 10.3390/ijms24043115] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Revised: 01/27/2023] [Accepted: 01/28/2023] [Indexed: 02/08/2023] Open
Abstract
High-altitude environments dramatically influenced the genetic evolution of vertebrates. However, little is known about the role of RNA editing on high-altitude adaptation in non-model species. Here, we profiled the RNA editing sites (RESs) of heart, lung, kidney, and longissimus dorsi muscle from Tibetan cashmere goats (TBG, 4500 m) and Inner Mongolia cashmere goats (IMG, 1200 m) to reveal RNA editing-related functions of high-altitude adaptation in goats. We identified 84,132 high-quality RESs that were unevenly distributed across the autosomes in TBG and IMG, and more than half of the 10,842 non-redundant editing sites were clustered. The majority (62.61%) were adenosine-to-inosine (A-to-I) sites, followed by cytidine-to-uridine (C-to-U) sites (19.26%), and 32.5% of them had a significant correlation with the expression of catalytic genes. Moreover, A-to-I and C-to-U RNA editing sites had different flanking sequences, amino acid mutations, and alternative splicing activity. TBG had higher editing levels of A-to-I and C-to-U than IMG in the kidney, whereas a lower level was found in the longissimus dorsi muscle. Furthermore, we identified 29 IMG and 41 TBG population-specific editing sites (pSESs) and 53 population-differential editing sites (pDESs) that were functionally involved in altering RNA splicing or recoding protein products. It is worth noting that 73.3% population-differential, 73.2% TBG-specific, and 80% IMG-specific A-to-I sites were nonsynonymous sites. Moreover, the pSESs and pDESs editing-related genes play critical functions in energy metabolisms such as ATP binding molecular function, translation, and adaptive immune response, which may be linked to goat high-altitude adaptation. Our results provide valuable information for understanding the adaptive evolution of goats and studying plateau-related diseases.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | | | - Liu Yang
- Farm Animal Genetic Resources Exploration Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Hongping Zhang
- Farm Animal Genetic Resources Exploration Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| |
Collapse
|
2
|
Wang YM, Ye LQ, Wang MS, Zhang JJ, Khederzadeh S, Irwin DM, Ren XD, Zhang YP, Wu DD. Unveiling the functional and evolutionary landscape of RNA editing in chicken using genomics and transcriptomics. Zool Res 2022; 43:1011-1022. [PMID: 36266925 PMCID: PMC9700494 DOI: 10.24272/j.issn.2095-8137.2022.331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Accepted: 10/17/2022] [Indexed: 09/10/2024] Open
Abstract
The evolutionary and functional features of RNA editing are well studied in mammals, cephalopods, and insects, but not in birds. Here, we integrated transcriptomic and whole-genomic analyses to exhaustively characterize the expansive repertoire of adenosine-to-inosine (A-to-I) RNA editing sites (RESs) in the chicken. In addition, we investigated the evolutionary status of the chicken editome as a potential mechanism of domestication. We detected the lowest editing level in the liver of chickens, compared to muscles in humans, and found higher editing activity and specificity in the brain than in non-neural tissues, consistent with the brain's functional complexity. To a certain extent, specific editing activity may account for the specific functions of tissues. Our results also revealed that sequences critical to RES secondary structures remained conserved within avian evolution. Furthermore, the RNA editome was shaped by purifying selection during chicken domestication and most RESs may have served as a selection pool for a few functional RESs involved in chicken domestication, including evolution of nervous and immune systems. Regulation of RNA editing in chickens by adenosine deaminase acting on RNA (ADAR) enzymes may be affected by non-ADAR factors whose expression levels changed widely after ADAR knockdown. Collectively, we provide comprehensive lists of candidate RESs and non-ADAR-editing regulators in the chicken, thus contributing to our current understanding of the functions and evolution of RNA editing in animals.
Collapse
Affiliation(s)
- Yun-Mei Wang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
| | - Ling-Qun Ye
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
| | - Ming-Shan Wang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Department of Ecology and Evolutionary Biology, Howard Hughes Medical Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Jin-Jin Zhang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
| | - Saber Khederzadeh
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
| | - David M Irwin
- Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Ontario M5S 1A8, Canada
| | - Xiao-Die Ren
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
| | - Ya-Ping Zhang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650201, China. E-mail:
| | - Dong-Dong Wu
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- Kunming College of Life Science, University of the Chinese Academy of Sciences, Kunming, Yunnan, 650223, China
- KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650201, China
- Kunming Natural History Museum of Zoology, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650223, China. E-mail:
| |
Collapse
|
3
|
Zhou X, Bao P, Zhang X, Guo X, Liang C, Chu M, Wu X, Yan P. Genome-wide detection of RNA editing events during the hair follicles cycle of Tianzhu white yak. BMC Genomics 2022; 23:737. [PMID: 36316632 PMCID: PMC9624038 DOI: 10.1186/s12864-022-08951-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2022] [Accepted: 10/18/2022] [Indexed: 11/06/2022] Open
Abstract
BACKGROUND The hair coat is available for the yak to live in the harsh environment of the plateau. Besides, improving the hair production of yak is necessary for its textile industry development. Hair grows from hair follicles (HFs). The HFs undergo periodic growth after birth and are regulated by the complex gene regulatory network. However, the molecular mechanism of HFs regeneration in the Tianzhu white yak remains unclear. RNA editing is a post-transcriptional mechanism that regulates gene expression and produces new transcripts. Hence, we investigated the influence of the A-to-I RNA editing events on the HFs cycle of the Tianzhu white yak. RESULTS We finally identified 54,707 adenosine-to-inosine (A-to-I) RNA editing sites (RESs) from RNA sequencing data of the HFs cycle in the Tianzhu white yak. Annotation results showed RESs caused missense amino acid changes in 7 known genes. And 202 A-to-I editing sites altered 23 target genes of 140 microRNAs. A total of 1,722 differential RESs were identified during the HFs cycle of Tianzhu white yak. GO and KEGG enrichment analysis revealed several signaling pathways and GO terms involved skin development, hair growth, and HFs cycle. Such as genes with differential RNA editing levels were significantly enriched in the peroxisome, metabolic pathways, Notch signaling pathway, and PPAR signaling pathway. Besides, the editing sites in HFs development-related genes FAS, APCDD1, WWOX, MPZL3, RUNX1, KANK2, DCN, DSC2, LEPR, HEPHL1, and PTK2B were suggested as the potential RESs involving HFs development. CONCLUSION This study investigated the global A-to-I RNA editing events during the HFs cycle of yak skin tissue and expanded the knowledge of A-to-I RNA editing on the HFs cycle. Furthermore, this study revealed that RNA editing-influenced genes may regulate the HFs cycle by participating in the HFs development-related pathways. The findings might provide new insight into the regulation of RNA editing in hair growth.
Collapse
Affiliation(s)
- Xuelan Zhou
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China
| | - Pengjia Bao
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China
| | - Xiaolan Zhang
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China
| | - Xian Guo
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China
| | - Chunnian Liang
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China
| | - Min Chu
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China
| | - Xiaoyun Wu
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China.
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China.
| | - Ping Yan
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, 730050, Lanzhou, P.R. China.
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, 730050, Lanzhou, P.R. China.
| |
Collapse
|
4
|
Zhang Y, Liu X, Zhang L, Wang L, He J, Ma H, Wang L. Preliminary identification and analysis of differential RNA editing between higher and lower backfat thickness pigs using DNA-seq and RNA-seq data. Anim Genet 2022; 53:327-339. [PMID: 35342974 DOI: 10.1111/age.13193] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Revised: 12/24/2021] [Accepted: 03/03/2022] [Indexed: 12/31/2022]
Abstract
RNA editing is an essential post-transcriptional regulatory mechanism. However, few studies have investigated the functional RNA edits in the economic traits of livestock on a genome-wide scale. Pigs are one of the most important livestock species and their fat is the principal organ involved in the regulation of adipose deposition. Here, we used three full-sibling pairs, with each pair comprising a pig with higher backfat (BF) thickness and lower backfat thickness, to identify RNA editing events based on whole-genome and transcriptome sequencing data. A total of 60,903 non-redundant RNA editing sites with 59,472 (97.7%) A-to-G edits were detected using a revised bioinformatics pipeline. A specific sequence context with G preference was found one base downstream of the edited site, and the editing level was associated with the distribution of nucleotides across nearly sites. Moreover, the A-to-G editing sites mostly occurred in the pig-special short interspersed nuclear elements, Pre0_SS. Comparing the difference between pigs with higher BF and lower BF, we found 211 differentially edited sites (DESites). Functional enrichment analyses revealed a significant enrichment of genes containing DESites in terms of adipose deposition. The DESites located in the six adipose-related genes (SKP1, GSK3B, COL5A3, MDM4, NT5C2, and DENND2A) were selected as candidate RNA editing sites associated with adipose deposition, and thus require further evaluation. This study mined the potentially functional RNA editing sites in pig adipose tissue and indicated that RNA editing may play an important role in adipose deposition, which provides a new insight into the post-transcriptionally mediated regulation mechanism of fat development.
Collapse
Affiliation(s)
- Yuebo Zhang
- College of Animal Science and Technology, Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, Hunan Agricultural University, Changsha, China
| | - Xin Liu
- Institute of Animal Science, State Key Laboratory of Animal Nutrition, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Longchao Zhang
- Institute of Animal Science, State Key Laboratory of Animal Nutrition, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Ligang Wang
- Institute of Animal Science, State Key Laboratory of Animal Nutrition, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Jun He
- College of Animal Science and Technology, Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, Hunan Agricultural University, Changsha, China
| | - Haiming Ma
- College of Animal Science and Technology, Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, Hunan Agricultural University, Changsha, China
| | - Lixian Wang
- Institute of Animal Science, State Key Laboratory of Animal Nutrition, Chinese Academy of Agricultural Sciences, Beijing, China
| |
Collapse
|
5
|
Wu X, Ayalew W, Chu M, Pei J, Liang C, Bao P, Guo X, Yan P. Characterization of RNA Editome in the Mammary Gland of Yaks during the Lactation and Dry Periods. Animals (Basel) 2022; 12:ani12020207. [PMID: 35049829 PMCID: PMC8773173 DOI: 10.3390/ani12020207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Revised: 01/01/2022] [Accepted: 01/14/2022] [Indexed: 12/03/2022] Open
Abstract
Simple Summary In order to study the influence of RNA editing sites on lactation and mammary gland development process in yaks, we comprehensively characterized the RNA editome of the yak mammary gland during the lactation period and dry period by using the transcriptome and genome sequencing data. The results revealed 82,872 nonredundant RNA editing sites, 14,159 of which were differentially edited between the lactation period and dry period. Enrichment analysis showed that the genes harboring differential editing sites were mainly associated with mammary gland development-related pathways, such as MAPK pathway, PI3K-Akt pathway, FoxO signaling pathway, GnRH signaling pathway, and focal adhesion pathway. Our findings offer some novel insights into the RNA editing function in the mammary gland of yaks. Abstract The mammary gland is a complicated organ comprising several types of cells, and it undergoes extensive morphogenetic and metabolic changes during the female reproductive cycle. RNA editing is a posttranscriptional modification event occurring at the RNA nucleotide level, and it drives transcriptomic and proteomic diversities, with potential functional consequences. RNA editing in the mammary gland of yaks, however, remains poorly understood. Here, we used REDItools to identify RNA editing sites in mammary gland tissues in yaks during the lactation period (LP, n = 2) and dry period (DP, n = 3). Totally, 82,872 unique RNA editing sites were identified, most of which were detected in the noncoding regions with a low editing degree. In the coding regions (CDS), we detected 5235 editing sites, among which 1884 caused nonsynonymous amino acid changes. Of these RNA editing sites, 486 were found to generate novel possible miRNA target sites or interfere with the initial miRNA binding sites, indicating that RNA editing was related to gene regulation mediated by miRNA. A total of 14,159 RNA editing sites (involving 3238 common genes) showed a significant differential editing level in the LP when compared with that in the DP through Tukey’s Honest Significant Difference method (p < 0.05). According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, genes that showed different RNA editing levels mainly participated in pathways highly related to mammary gland development, including MAPK, PI3K-Akt, FoxO, and GnRH signaling pathways. Collectively, this work demonstrated for the first time the dynamic RNA editome profiles in the mammary gland of yaks and shed more light on the mechanism that regulates lactation together with mammary gland development.
Collapse
Affiliation(s)
- Xiaoyun Wu
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
| | - Wondossen Ayalew
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
- Department of Animal Production and Technology, Wolkite University, Wolkite P.O. Box 07, Ethiopia
| | - Min Chu
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
| | - Jie Pei
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
| | - Chunnian Liang
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
| | - Pengjia Bao
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
| | - Xian Guo
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
- Correspondence: (X.G.); (P.Y.)
| | - Ping Yan
- Key Laboratory of Yak Breeding Engineering, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China; (X.W.); (W.A.); (M.C.); (J.P.); (C.L.); (P.B.)
- Correspondence: (X.G.); (P.Y.)
| |
Collapse
|
6
|
Adetula AA, Fan X, Zhang Y, Yao Y, Yan J, Chen M, Tang Y, Liu Y, Yi G, Li K, Tang Z. Landscape of tissue-specific RNA Editome provides insight into co-regulated and altered gene expression in pigs ( Sus-scrofa). RNA Biol 2021; 18:439-450. [PMID: 34314293 PMCID: PMC8677025 DOI: 10.1080/15476286.2021.1954380] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Revised: 07/02/2021] [Accepted: 07/07/2021] [Indexed: 11/08/2022] Open
Abstract
RNA editing generates genetic diversity in mammals by altering amino acid sequences, miRNA targeting site sequences, influencing the stability of targeted RNAs, and causing changes in gene expression. However, the extent to which RNA editing affect gene expression via modifying miRNA binding site remains unexplored. Here, we first profiled the dynamic A-to-I RNA editome across tissues of Duroc and Luchuan pigs. The RNA editing events at the miRNA binding sites were generated. The biological function of the differentially edited gene in skeletal muscle was further characterized in pig muscle-derived satellite cells. RNA editome analysis revealed a total of 171,909 A-to-I RNA editing sites (RESs), and examination of its features showed that these A-to-I editing sites were mainly located in SINE retrotransposons PRE-1/Pre0_SS element. Analysis of differentially edited sites (DESs) revealed a total of 4,552 DESs across tissues between Duroc and Luchuan pigs, and functional category enrichment analysis of differentially edited gene (DEG) sets highlighted a significant association and enrichment of tissue-developmental pathways including TGF-beta, PI3K-Akt, AMPK, and Wnt signaling pathways. Moreover, we found that RNA editing events at the miRNA binding sites in the 3'-UTR of HSPA12B mRNA could prevent the miRNA-mediated mRNA downregulation of HSPA12B in the muscle-derived satellite (MDS) cell, consistent with the results obtained from the Luchuan skeletal muscle. This study represents the most systematic attempt to characterize the significance of RNA editing in regulating gene expression, particularly in skeletal muscle, constituting a new layer of regulation to understand the genetic mechanisms behind phenotype variance in animals.Abbreviations: A-to-I: Adenosine-to-inosine; ADAR: Adenosine deaminase acting on RNA; RES: RNA editing site; DEG: Differentially edited gene; DES: Differentially edited site; FDR: False discovery rate; GO: Gene Ontology; KEGG: Kyoto Encyclopaedia of Genes and Genomes; MDS cell: musclederived satellite cell; RPKM: Reads per kilobase of exon model in a gene per million mapped reads; UTR: Untranslated coding regions.
Collapse
Affiliation(s)
- Adeyinka A. Adetula
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Xinhao Fan
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Yongsheng Zhang
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Yilong Yao
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Junyu Yan
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Muya Chen
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Yijie Tang
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Yuwen Liu
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Kunpeng Institute of Modern Agriculture at Foshan, Foshan, China
- GuangXi Engineering Centre for Resource Development of Bama Xiang Pig, Bama, China
| | - Guoqiang Yi
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Kunpeng Institute of Modern Agriculture at Foshan, Foshan, China
- GuangXi Engineering Centre for Resource Development of Bama Xiang Pig, Bama, China
| | - Kui Li
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Zhonglin Tang
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
- Group of Pig Genome and Design Breeding, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- Kunpeng Institute of Modern Agriculture at Foshan, Foshan, China
- GuangXi Engineering Centre for Resource Development of Bama Xiang Pig, Bama, China
| |
Collapse
|
7
|
Jehl F, Degalez F, Bernard M, Lecerf F, Lagoutte L, Désert C, Coulée M, Bouchez O, Leroux S, Abasht B, Tixier-Boichard M, Bed'hom B, Burlot T, Gourichon D, Bardou P, Acloque H, Foissac S, Djebali S, Giuffra E, Zerjal T, Pitel F, Klopp C, Lagarrigue S. RNA-Seq Data for Reliable SNP Detection and Genotype Calling: Interest for Coding Variant Characterization and Cis-Regulation Analysis by Allele-Specific Expression in Livestock Species. Front Genet 2021; 12:655707. [PMID: 34262593 PMCID: PMC8273700 DOI: 10.3389/fgene.2021.655707] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Accepted: 06/01/2021] [Indexed: 12/19/2022] Open
Abstract
In addition to their common usages to study gene expression, RNA-seq data accumulated over the last 10 years are a yet-unexploited resource of SNPs in numerous individuals from different populations. SNP detection by RNA-seq is particularly interesting for livestock species since whole genome sequencing is expensive and exome sequencing tools are unavailable. These SNPs detected in expressed regions can be used to characterize variants affecting protein functions, and to study cis-regulated genes by analyzing allele-specific expression (ASE) in the tissue of interest. However, gene expression can be highly variable, and filters for SNP detection using the popular GATK toolkit are not yet standardized, making SNP detection and genotype calling by RNA-seq a challenging endeavor. We compared SNP calling results using GATK suggested filters, on two chicken populations for which both RNA-seq and DNA-seq data were available for the same samples of the same tissue. We showed, in expressed regions, a RNA-seq precision of 91% (SNPs detected by RNA-seq and shared by DNA-seq) and we characterized the remaining 9% of SNPs. We then studied the genotype (GT) obtained by RNA-seq and the impact of two factors (GT call-rate and read number per GT) on the concordance of GT with DNA-seq; we proposed thresholds for them leading to a 95% concordance. Applying these thresholds to 767 multi-tissue RNA-seq of 382 birds of 11 chicken populations, we found 9.5 M SNPs in total, of which ∼550,000 SNPs per tissue and population with a reliable GT (call rate ≥ 50%) and among them, ∼340,000 with a MAF ≥ 10%. We showed that such RNA-seq data from one tissue can be used to (i) detect SNPs with a strong predicted impact on proteins, despite their scarcity in each population (16,307 SIFT deleterious missenses and 590 stop-gained), (ii) study, on a large scale, cis-regulations of gene expression, with ∼81% of protein-coding and 68% of long non-coding genes (TPM ≥ 1) that can be analyzed for ASE, and with ∼29% of them that were cis-regulated, and (iii) analyze population genetic using such SNPs located in expressed regions. This work shows that RNA-seq data can be used with good confidence to detect SNPs and associated GT within various populations and used them for different analyses as GTEx studies.
Collapse
Affiliation(s)
- Frédéric Jehl
- INRAE, INSTITUT AGRO, PEGASE UMR 1348, Saint-Gilles, France
| | - Fabien Degalez
- INRAE, INSTITUT AGRO, PEGASE UMR 1348, Saint-Gilles, France
| | - Maria Bernard
- INRAE, SIGENAE, Genotoul Bioinfo MIAT, Castanet-Tolosan, France.,INRAE, AgroParisTech, Université Paris-Saclay, GABI UMR 1313, Jouy-en-Josas, France
| | | | | | - Colette Désert
- INRAE, INSTITUT AGRO, PEGASE UMR 1348, Saint-Gilles, France
| | - Manon Coulée
- INRAE, INSTITUT AGRO, PEGASE UMR 1348, Saint-Gilles, France
| | - Olivier Bouchez
- INRAE, US 1426, GeT-PlaGe, Genotoul, Castanet-Tolosan, France
| | - Sophie Leroux
- INRAE, INPT, ENVT, Université de Toulouse, GenPhySE UMR 1388, Castanet-Tolosan, France
| | - Behnam Abasht
- Department of Animal and Food Sciences, University of Delaware, Newark, DE, United States
| | | | - Bertrand Bed'hom
- INRAE, AgroParisTech, Université Paris-Saclay, GABI UMR 1313, Jouy-en-Josas, France
| | | | | | - Philippe Bardou
- INRAE, SIGENAE, Genotoul Bioinfo MIAT, Castanet-Tolosan, France
| | - Hervé Acloque
- INRAE, AgroParisTech, Université Paris-Saclay, GABI UMR 1313, Jouy-en-Josas, France
| | - Sylvain Foissac
- INRAE, INPT, ENVT, Université de Toulouse, GenPhySE UMR 1388, Castanet-Tolosan, France
| | - Sarah Djebali
- INRAE, INPT, ENVT, Université de Toulouse, GenPhySE UMR 1388, Castanet-Tolosan, France
| | - Elisabetta Giuffra
- INRAE, AgroParisTech, Université Paris-Saclay, GABI UMR 1313, Jouy-en-Josas, France
| | - Tatiana Zerjal
- INRAE, AgroParisTech, Université Paris-Saclay, GABI UMR 1313, Jouy-en-Josas, France
| | - Frédérique Pitel
- INRAE, INPT, ENVT, Université de Toulouse, GenPhySE UMR 1388, Castanet-Tolosan, France
| | | | | |
Collapse
|
8
|
Genome-Wide Characterization of RNA Editing Sites in Primary Gastric Adenocarcinoma through RNA-seq Data Analysis. Int J Genomics 2020; 2020:6493963. [PMID: 33415135 PMCID: PMC7768588 DOI: 10.1155/2020/6493963] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 07/28/2020] [Accepted: 12/07/2020] [Indexed: 12/15/2022] Open
Abstract
RNA editing is a posttranscriptional nucleotide modification in humans. Of the various types of RNA editing, the adenosine to inosine substitution is the most widespread in higher eukaryotes, which is mediated by the ADAR family enzymes. Inosine is recognized by the biological machinery as guanosine; therefore, editing could have substantial functional effects throughout the genome. RNA editing could contribute to cancer either by exclusive editing of tumor suppressor/promoting genes or by introducing transcriptomic diversity to promote cancer progression. Here, we provided a comprehensive overview of the RNA editing sites in gastric adenocarcinoma and highlighted some of their possible contributions to gastric cancer. RNA-seq data corresponding to 8 gastric adenocarcinoma and their paired nontumor counterparts were retrieved from the GEO database. After preprocessing and variant calling steps, a stringent filtering pipeline was employed to distinguish potential RNA editing sites from SNPs. The identified potential editing sites were annotated and compared with those in the DARNED database. Totally, 12362 high-confidence adenosine to inosine RNA editing sites were detected across all samples. Of these, 12105 and 257 were known and novel editing events, respectively. These editing sites were unevenly distributed across genomic regions, and nearly half of them were located in 3′UTR. Our results revealed that 4868 editing sites were common in both normal and cancer tissues. From the remaining sites, 3985 and 3509 were exclusive to normal and cancer tissues, respectively. Further analysis revealed a significant number of differentially edited events among these sites, which were located in protein coding genes and microRNAs. Given the distinct pattern of RNA editing in gastric adenocarcinoma and adjacent normal tissue, edited sites have the potential to serve as the diagnostic biomarkers and therapeutic targets in gastric cancer.
Collapse
|
9
|
Wang L, Li J, Hou X, Yan H, Zhang L, Liu X, Gao H, Zhao F, Wang L. Genome-Wide Identification of RNA Editing Sites Affecting Intramuscular Fat in Pigs. Animals (Basel) 2020; 10:E1616. [PMID: 32927662 PMCID: PMC7552122 DOI: 10.3390/ani10091616] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 09/03/2020] [Accepted: 09/07/2020] [Indexed: 12/19/2022] Open
Abstract
Intramuscular fat (IMF) is essential for improving the palatability and flavor of meat, and it is strongly associated with human insulin resistance. RNA editing is a widespread regulating event in different tissues. Here, we investigated the global RNA editing difference of two groups of pig with different IMF contents to find the potential editing sites affecting IMF. In this research, RES-Scanner and REDItools were used to identify RNA editing sites based on the whole genome and transcriptome sequencing data of the high and low groups composed of three full-sib pairs with opposite IMF phenotypes. A total of 295 RNA editing sites were investigated in the longissimus dorsi muscle, and 90.17% of these sites caused A to G conversion. After annotation, most editing sites were located in noncoding regions (including five sites located on the 3' UTR regions). Five editing sites (including two sites that could lead to nonsynonymous amino acid changes) were located in the exons of genes. A total of 36 intergroup (high and low IMF) differential RNA editing sites were found in 33 genes. Some candidate editing sites, such as sites in acyl-coenzymeA: cholesterol acyltransferase 1 (ACAT1), coatomer protein, subunit alpha (COPA), and nuclear receptor coactivator 3 (NCOA3), were selected as candidate RNA editing sites associated with IMF. One site located on the 3' UTR region of growth hormone secretagogue receptor (GHSR) may regulate GHSR expression by affecting the interaction of miRNA and mRNA. In conclusion, we identified a total of 36 nonredundant RNA editing sites in the longissimus dorsi muscle, which may reveal the potential importance of RNA editing in IMF. Four were selected as candidate sites associated with IMF. Our findings provide some new insights of RNA editing function in pig longissimus dorsi muscle which useful for pig IMF breeding or human insulin resistances research.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | - Lixian Wang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing 100193, China; (L.W.); (J.L.); (X.H.); (H.Y.); (L.Z.); (X.L.); (H.G.); (F.Z.)
| |
Collapse
|
10
|
Bakhtiarizadeh MR, Alamouti AA. RNA-Seq based genetic variant discovery provides new insights into controlling fat deposition in the tail of sheep. Sci Rep 2020; 10:13525. [PMID: 32782325 PMCID: PMC7419499 DOI: 10.1038/s41598-020-70527-8] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 07/27/2020] [Indexed: 01/09/2023] Open
Abstract
Genetic basis of fat deposition in sheep tail have not been completely elucidated yet. Understanding the genetic mechanisms controlling fat-tail size can improve breeding strategies to modulate fat deposition. RNA sequencing has made it possible to discover genetic variants that may underlie various phenotypic differences. Hence, to identify genetic variants that are important for describing different fat-tail phenotypes in sheep, RNA sequencing was used for single nucleotide polymorphism (SNP) calling in two Iranian sheep breeds (Lori-Bakhtiari, fat-tailed; n = 4, vs Zel, thin-tailed; n = 4). Using a stringent pipeline, a total of 112,344 known SNPs were genotyped, of which 30,550 and 42,906 SNPs were shared by at least two Lori-Bakhtiari and Zel, respectively. Comparing these SNPs showed 2,774 (including 209 missense and 25 deleterious SNPs) and 10,470 (including 1,054 missense and 116 deleterious SNPs) breed-specific SNPs in Lori-Bakhtiari and Zel sheep, respectively. Potential breed-specific SNPs were detected by considering those located in QTL regions associated with fatness or reported as important candidates in previous similar studies. Of the breed-specific SNPs, 724 and 2,905 were located in the QTL regions. Functional enrichment analysis of the affected genes revealed several enriched gene ontologies and KEGG pathways related to fat metabolism. Based on the results, several affected genes were proposed to be strongly linked with fat deposition such as DGAT2, ACSL1, ACACA, ADIPOQ, ACLY, FASN, CPT2, SCD, ADCY6, PER3, CSF1R, SLC22A4, GFPT1, CDS2, BMP6, ACSS2, ELOVL6, HOXA10 and FABP4. Moreover, several SNPs were found in the candidate genes related to fatty acid oxidation introducing them as promising candidates responsible for lower fat content in tail of Zel. Our findings provided new insights into the genetic mechanisms of fat deposition in sheep, which can serve to designing appropriate breeding programs.
Collapse
Affiliation(s)
| | - Ali A Alamouti
- Department of Animal and Poultry Science, College of Aburaihan, University of Tehran, Tehran, Iran
| |
Collapse
|
11
|
Adetunji MO, Lamont SJ, Abasht B, Schmidt CJ. Variant analysis pipeline for accurate detection of genomic variants from transcriptome sequencing data. PLoS One 2019; 14:e0216838. [PMID: 31545812 PMCID: PMC6756534 DOI: 10.1371/journal.pone.0216838] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2019] [Accepted: 09/10/2019] [Indexed: 12/27/2022] Open
Abstract
The wealth of information deliverable from transcriptome sequencing (RNA-seq) is significant, however current applications for variant detection still remain a challenge due to the complexity of the transcriptome. Given the ability of RNA-seq to reveal active regions of the genome, detection of RNA-seq SNPs can prove valuable in understanding the phenotypic diversity between populations. Thus, we present a novel computational workflow named VAP (Variant Analysis Pipeline) that takes advantage of multiple RNA-seq splice aware aligners to call SNPs in non-human models using RNA-seq data only. We applied VAP to RNA-seq from a highly inbred chicken line and achieved high accuracy when compared with the matching whole genome sequencing (WGS) data. Over 65% of WGS coding variants were identified from RNA-seq. Further, our results discovered SNPs resulting from post transcriptional modifications, such as RNA editing, which may reveal potentially functional variation that would have otherwise been missed in genomic data. Even with the limitation in detecting variants in expressed regions only, our method proves to be a reliable alternative for SNP identification using RNA-seq data. The source code and user manuals are available at https://modupeore.github.io/VAP/.
Collapse
Affiliation(s)
- Modupeore O. Adetunji
- Department of Animal and Food Sciences, University of Delaware, Newark, Delaware, United States of America
- * E-mail:
| | - Susan J. Lamont
- Department of Animal Science, Iowa State University, Ames, Iowa, United States of America
| | - Behnam Abasht
- Department of Animal and Food Sciences, University of Delaware, Newark, Delaware, United States of America
| | - Carl J. Schmidt
- Department of Animal and Food Sciences, University of Delaware, Newark, Delaware, United States of America
| |
Collapse
|