1
|
Cortellari M, Barbato M, Talenti A, Bionda A, Carta A, Ciampolini R, Ciani E, Crisà A, Frattini S, Lasagna E, Marletta D, Mastrangelo S, Negro A, Randi E, Sarti FM, Sartore S, Soglia D, Liotta L, Stella A, Ajmone-Marsan P, Pilla F, Colli L, Crepaldi P. The climatic and genetic heritage of Italian goat breeds with genomic SNP data. Sci Rep 2021; 11:10986. [PMID: 34040003 PMCID: PMC8154919 DOI: 10.1038/s41598-021-89900-2] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2020] [Accepted: 04/29/2021] [Indexed: 02/04/2023] Open
Abstract
Local adaptation of animals to the environment can abruptly become a burden when faced with rapid climatic changes such as those foreseen for the Italian peninsula over the next 70 years. Our study investigates the genetic structure of the Italian goat populations and links it with the environment and how genetics might evolve over the next 50 years. We used one of the largest national datasets including > 1000 goats from 33 populations across the Italian peninsula collected by the Italian Goat Consortium and genotyped with over 50 k markers. Our results showed that Italian goats can be discriminated in three groups reflective of the Italian geography and its geo-political situation preceding the country unification around two centuries ago. We leveraged the remarkable genetic and geographical diversity of the Italian goat populations and performed landscape genomics analysis to disentangle the relationship between genotype and environment, finding 64 SNPs intercepting genomic regions linked to growth, circadian rhythm, fertility, and inflammatory response. Lastly, we calculated the hypothetical future genotypic frequencies of the most relevant SNPs identified through landscape genomics to evaluate their long-term effect on the genetic structure of the Italian goat populations. Our results provide an insight into the past and the future of the Italian local goat populations, helping the institutions in defining new conservation strategy plans that could preserve their diversity and their link to local realities challenged by climate change.
Collapse
Affiliation(s)
- Matteo Cortellari
- Dipartimento di Scienze Agrarie e Ambientali - Produzione, Territorio, Agroenergia, Università degli Studi di Milano, Via Celoria 2, 20133, Milan, Italy
| | - Mario Barbato
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti and BioDNA Centro di ricerca sulla Biodiversità e sul DNA Antico, Università Cattolica del Sacro Cuore, Via Emilia Parmense 84, 29122, Piacenza, Italy
| | - Andrea Talenti
- Dipartimento di Scienze Agrarie e Ambientali - Produzione, Territorio, Agroenergia, Università degli Studi di Milano, Via Celoria 2, 20133, Milan, Italy.
- The Roslin Institute, University of Edinburgh, Easter Bush Campus, Midlothian, EH25 9RG, UK.
| | - Arianna Bionda
- Dipartimento di Scienze Agrarie e Ambientali - Produzione, Territorio, Agroenergia, Università degli Studi di Milano, Via Celoria 2, 20133, Milan, Italy
| | - Antonello Carta
- Unità di Ricerca di Genetica e Biotecnologie, Agris Sardegna, 07100, Sassari, Italy
| | - Roberta Ciampolini
- Dipartimento di Scienze Veterinarie, Università di Pisa, Viale delle Piagge 2, 56124, Pisa, Italy
| | - Elena Ciani
- Dipartimento di Bioscienze Biotecnologie e Biofarmaceutica, Università degli Studi di Bari, Via Orabona 4, 70126, Bari, Italy
| | - Alessandra Crisà
- Consiglio per la ricerca in agricoltura e l'analisi dell'economia agraria (CREA) - Research Centre for Animal Production and Acquaculture, 00015, Monterotondo, Rome, Italy
| | - Stefano Frattini
- Dipartimento di Scienze Agrarie e Ambientali - Produzione, Territorio, Agroenergia, Università degli Studi di Milano, Via Celoria 2, 20133, Milan, Italy
| | - Emiliano Lasagna
- Department of Agricultural, Food and Environmental Sciences, University of Perugia, 06121, Perugia, Italy
| | - Donata Marletta
- Department of Agriculture, Food and Environment, University of Catania, Via Valdisavoia 5, 95123, Catania, Italy
| | - Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128, Palermo, Italy
| | - Alessio Negro
- Dipartimento di Scienze Agrarie e Ambientali - Produzione, Territorio, Agroenergia, Università degli Studi di Milano, Via Celoria 2, 20133, Milan, Italy
| | - Ettore Randi
- Department of Chemistry and Bioscience, Faculty of Engineering and Science, University of Aalborg, Aalborg, Denmark
| | - Francesca M Sarti
- Department of Agricultural, Food and Environmental Sciences, University of Perugia, 06121, Perugia, Italy
| | - Stefano Sartore
- Dipartimento di Scienze Veterinarie, Università degli Studi di Torino, largo Braccini 2, 10095, Grugliasco, Italy
| | - Dominga Soglia
- Dipartimento di Scienze Veterinarie, Università degli Studi di Torino, largo Braccini 2, 10095, Grugliasco, Italy
| | - Luigi Liotta
- Dipartimento di Scienze Veterinarie, University of Messina, Messina, Italy
| | - Alessandra Stella
- Institute of Biology and Biotechnology in Agriculture, National Research Council (CNR), Milan, Italy
| | - Paolo Ajmone-Marsan
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti and BioDNA Centro di ricerca sulla Biodiversità e sul DNA Antico, Università Cattolica del Sacro Cuore, Via Emilia Parmense 84, 29122, Piacenza, Italy
| | - Fabio Pilla
- Dipartimento Agricoltura, Ambiente e Alimenti Universitá degli Studi del Molise, 86100, Campobasso, Italy
| | - Licia Colli
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti and BioDNA Centro di ricerca sulla Biodiversità e sul DNA Antico, Università Cattolica del Sacro Cuore, Via Emilia Parmense 84, 29122, Piacenza, Italy
| | - Paola Crepaldi
- Dipartimento di Scienze Agrarie e Ambientali - Produzione, Territorio, Agroenergia, Università degli Studi di Milano, Via Celoria 2, 20133, Milan, Italy
| |
Collapse
|
2
|
Jiao Q, Wang YR, Zhao JY, Wang ZY, Guo D, Bai WL. Identification and molecular analysis of cashmere goat lncRNAs reveal their integrated regulatory network and potential roles in secondary hair follicle. Anim Biotechnol 2020; 32:719-732. [PMID: 32233965 DOI: 10.1080/10495398.2020.1747477] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
Long non-coding RNAs (lncRNAs) is a class of eukaryotic transcripts with length of more than 200 bp. They contribute to the regulation of gene expressions involved in multiple processes including the skin cell proliferation, differentiation, and reconstruction of the secondary hair follicles (SHFs). In this study, firstly, we identified 16 putative lncRNAs from SHFs of cashmere goat based on the EST sequences from NCBI database. Secondly, we investigated their transcriptional pattern in SHFs of cashmere goat, and constructed their ceRNA regulatory networks. The RT-qPCR results showed four lncRNAs (lncRNA-475074, -052149, -052140, and -051789) were significantly up-regulated, and nine lncRNAs (lncRNA-711032, -475083, -475070, -052139, -052127, -052037, -051903, -051847, and -051804) were significantly down-regulatd in anagen SHFs of cashmere goat. CeRNA networks analysis revealed complex interactional relationship among lncRNAs, miRNAs and mRNAs. Further, the KEGG pathway enrichment was performed for the potential target genes of the identified lncRNAs based on bioinformatics technique, and the results indicated that differentially expressed lncRNAs directly or indirectly might regulate potential target genes. Our results from this study will provide a significant information for further exploring the functions and possible mechanisms of the identified lncRNAs in SHFs of cashmere goat.
Collapse
Affiliation(s)
- Qian Jiao
- College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, P.R. China
| | - Yan R Wang
- College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, P.R. China
| | - Jun Y Zhao
- College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, P.R. China
| | - Ze Y Wang
- College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, P.R. China
| | - Dan Guo
- Academy of Animal Husbandry Science of Liaoning Province, Liaoyang, P.R. China
| | - Wen L Bai
- College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, P.R. China
| |
Collapse
|
3
|
Characterization of Embryonic Skin Transcriptome in Anser cygnoides at Three Feather Follicles Developmental Stages. G3-GENES GENOMES GENETICS 2020; 10:443-454. [PMID: 31792007 PMCID: PMC7003092 DOI: 10.1534/g3.119.400875] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In order to enrich the Anser cygnoides genome and identify the gene expression profiles of primary and secondary feather follicles development, de novo transcriptome assembly of skin tissues was established by analyzing three developmental stages at embryonic day 14, 18, and 28 (E14, E18, E28). Sequencing output generated 436,730,608 clean reads from nine libraries and de novo assembled into 56,301 unigenes. There were 2,298, 9,423 and 12,559 unigenes showing differential expression in three stages respectively. Furthermore, differentially expressed genes (DEGs) were functionally classified according to genes ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and series-cluster analysis. Relevant specific GO terms such as epithelium development, regulation of keratinocyte proliferation, morphogenesis of an epithelium were identified. In all, 15,144 DEGs were clustered into eight profiles with distinct expression patterns and 2,424 DEGs were assigned to 198 KEGG pathways. Skin development related pathways (mitogen-activated protein kinase signaling pathway, extra-cellular matrix -receptor interaction, Wingless-type signaling pathway) and genes (delta like canonical Notch ligand 1, fibroblast growth factor 2, Snail family transcriptional repressor 2, bone morphogenetic protein 6, polo like kinase 1) were identified, and eight DEGs were selected to verify the reliability of transcriptome results by real-time quantitative PCR. The findings of this study will provide the key insights into the complicated molecular mechanism and breeding techniques underlying the developmental characteristics of skin and feather follicles in Anser cygnoides.
Collapse
|
4
|
Comprehensive analysis of circRNAs from cashmere goat skin by next generation RNA sequencing (RNA-seq). Sci Rep 2020; 10:516. [PMID: 31949277 PMCID: PMC6965140 DOI: 10.1038/s41598-019-57404-9] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Accepted: 12/23/2019] [Indexed: 02/07/2023] Open
Abstract
Circular RNA (circRNA) is endogenous non-coding RNA (ncRNA) with a covalently closed circular structure. It is mainly generated through RNA alternative splicing or back-splicing. CircRNA is known in the majority of eukaryotes and very stable. However, knowledge of the circRNA involved in regulating cashmere fineness is limited. Skin samples were collected from Liaoning cashmere goats (LCG) and Inner Mongolia cashmere goats (MCG) during the anagen period. For differentially expressed circRNAs, RNA sequencing was performed, and the analysis led to an identification of 17 up-regulated circRNAs and 15 down-regulated circRNAs in LCG compared with MCG skin samples. In order to find the differentially expressed circRNAs in LCG, we carried out qPCRs on 10 candidate circRNAs in coarse type skin of LCG (CT-LCG) and fine type skin of LCG (FT-LCG). Four circRNAs: ciRNA128, circRNA6854, circRNA4154 and circRNA3620 were confirmed to be significantly differential expression in LCG. Also, a regulatory network of circRNAs-miRNAs was bioinformatically deduced and may help to understand molecular mechanisms of potential circRNA involvement in regulating cashmere fineness.
Collapse
|
5
|
Tang H, Xiao X, Xu Y, Li C, Cheng K, Pan X, Li W. Utilization of carbon sources in the rice rhizosphere and nonrhizosphere soils with different long-term fertilization management. J Basic Microbiol 2019; 59:621-631. [DOI: 10.1002/jobm.201800736] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2019] [Revised: 03/07/2019] [Accepted: 03/31/2019] [Indexed: 01/26/2023]
Affiliation(s)
- Haiming Tang
- Farming Ecology, Hunan Academy of Agricultural Sciences; Institute of Soil and Fertilizer; Changsha China
| | - Xiaoping Xiao
- Farming Ecology, Hunan Academy of Agricultural Sciences; Institute of Soil and Fertilizer; Changsha China
| | - Yilan Xu
- Plant Science and Technology; Hunan Biological and Electromechanical Polytechnic; Changsha China
| | - Chao Li
- Farming Ecology, Hunan Academy of Agricultural Sciences; Institute of Soil and Fertilizer; Changsha China
| | - Kaikai Cheng
- Farming Ecology, Hunan Academy of Agricultural Sciences; Institute of Soil and Fertilizer; Changsha China
| | - Xiaochen Pan
- Farming Ecology, Hunan Academy of Agricultural Sciences; Institute of Soil and Fertilizer; Changsha China
| | - Weiyan Li
- Farming Ecology, Hunan Academy of Agricultural Sciences; Institute of Soil and Fertilizer; Changsha China
| |
Collapse
|
6
|
Seasonal differences in the testicular transcriptome profile of free-living European beavers (Castor fiber L.) determined by the RNA-Seq method. PLoS One 2017; 12:e0180323. [PMID: 28678806 PMCID: PMC5498055 DOI: 10.1371/journal.pone.0180323] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2017] [Accepted: 06/14/2017] [Indexed: 11/20/2022] Open
Abstract
The European beaver (Castor fiber L.) is an important free-living rodent that inhabits Eurasian temperate forests. Beavers are often referred to as ecosystem engineers because they create or change existing habitats, enhance biodiversity and prepare the environment for diverse plant and animal species. Beavers are protected in most European Union countries, but their genomic background remains unknown. In this study, gene expression patterns in beaver testes and the variations in genetic expression in breeding and non-breeding seasons were determined by high-throughput transcriptome sequencing. Paired-end sequencing in the Illumina HiSeq 2000 sequencer produced a total of 373.06 million of high-quality reads. De novo assembly of contigs yielded 130,741 unigenes with an average length of 1,369.3 nt, N50 value of 1,734, and average GC content of 46.51%. A comprehensive analysis of the testicular transcriptome revealed more than 26,000 highly expressed unigenes which exhibited the highest homology with Rattus norvegicus and Ictidomys tridecemlineatus genomes. More than 8,000 highly expressed genes were found to be involved in fundamental biological processes, cellular components or molecular pathways. The study also revealed 42 genes whose regulation differed between breeding and non-breeding seasons. During the non-breeding period, the expression of 37 genes was up-regulated, and the expression of 5 genes was down-regulated relative to the breeding season. The identified genes encode molecules which are involved in signaling transduction, DNA repair, stress responses, inflammatory processes, metabolism and steroidogenesis. Our results pave the way for further research into season-dependent variations in beaver testes.
Collapse
|
7
|
Li Z, Chen F, Huang C, Zheng W, Yu C, Cheng H, Zhou R. Genome-wide mapping and characterization of microsatellites in the swamp eel genome. Sci Rep 2017; 7:3157. [PMID: 28600492 PMCID: PMC5466649 DOI: 10.1038/s41598-017-03330-7] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2017] [Accepted: 04/26/2017] [Indexed: 11/09/2022] Open
Abstract
We described genome-wide screening and characterization of microsatellites in the swamp eel genome. A total of 99,293 microsatellite loci were identified in the genome with an overall density of 179 microsatellites per megabase of genomic sequences. The dinucleotide microsatellites were the most abundant type representing 71% of the total microsatellite loci and the AC-rich motifs were the most recurrent in all repeat types. Microsatellite frequency decreased as numbers of repeat units increased, which was more obvious in long than short microsatellite motifs. Most of microsatellites were located in non-coding regions, whereas only approximately 1% of the microsatellites were detected in coding regions. Trinucleotide repeats were most abundant microsatellites in the coding regions, which represented amino acid repeats in proteins. There was a chromosome-biased distribution of microsatellites in non-coding regions, with the highest density of 203.95/Mb on chromosome 8 and the least on chromosome 7 (164.06/Mb). The most abundant dinucleotides (AC)n was mainly located on chromosome 8. Notably, genomic mapping showed that there was a chromosome-biased association of genomic distributions between microsatellites and transposon elements. Thus, the novel dataset of microsatellites in swamp eel provides a valuable resource for further studies on QTL-based selection breeding, genetic resource conservation and evolutionary genetics.
Collapse
Affiliation(s)
- Zhigang Li
- Hubei Key Laboratory of Cell Homeostasis, Laboratory of Molecular and Developmental Genetics, College of Life Sciences, Wuhan University, Wuhan, 430072, P. R. China
| | - Feng Chen
- Hubei Key Laboratory of Cell Homeostasis, Laboratory of Molecular and Developmental Genetics, College of Life Sciences, Wuhan University, Wuhan, 430072, P. R. China
| | - Chunhua Huang
- Hubei Key Laboratory of Cell Homeostasis, Laboratory of Molecular and Developmental Genetics, College of Life Sciences, Wuhan University, Wuhan, 430072, P. R. China
| | - Weixin Zheng
- Hubei Key Laboratory of Cell Homeostasis, Laboratory of Molecular and Developmental Genetics, College of Life Sciences, Wuhan University, Wuhan, 430072, P. R. China
| | - Chunlai Yu
- Hubei Key Laboratory of Cell Homeostasis, Laboratory of Molecular and Developmental Genetics, College of Life Sciences, Wuhan University, Wuhan, 430072, P. R. China
| | - Hanhua Cheng
- Hubei Key Laboratory of Cell Homeostasis, Laboratory of Molecular and Developmental Genetics, College of Life Sciences, Wuhan University, Wuhan, 430072, P. R. China.
| | - Rongjia Zhou
- Hubei Key Laboratory of Cell Homeostasis, Laboratory of Molecular and Developmental Genetics, College of Life Sciences, Wuhan University, Wuhan, 430072, P. R. China.
| |
Collapse
|
8
|
Comparative and temporal transcriptome analysis of peste des petits ruminants virus infected goat peripheral blood mononuclear cells. Virus Res 2017; 229:28-40. [DOI: 10.1016/j.virusres.2016.12.014] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2016] [Revised: 12/20/2016] [Accepted: 12/20/2016] [Indexed: 11/22/2022]
|
9
|
Bai WL, Wang JJ, Yin RH, Dang YL, Wang ZY, Zhu YB, Cong YY, Deng L, Guo D, Wang SQ, Yang SH, Xue HL. Molecular characterization of HOXC8 gene and methylation status analysis of its exon 1 associated with the length of cashmere fiber in Liaoning cashmere goat. Genetica 2017; 145:115-126. [DOI: 10.1007/s10709-017-9950-5] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2016] [Accepted: 01/07/2017] [Indexed: 11/29/2022]
|
10
|
Jia BY, Ba HX, Wang GW, Yang Y, Cui XZ, Peng YH, Zheng JJ, Xing XM, Yang FH. Transcriptome analysis of sika deer in China. Mol Genet Genomics 2016; 291:1941-53. [PMID: 27423230 DOI: 10.1007/s00438-016-1231-y] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2016] [Accepted: 07/11/2016] [Indexed: 12/17/2022]
Abstract
Sika deer is of great commercial value because their antlers are used in tonics and alternative medicine and their meat is healthy and delicious. The goal of this study was to generate transcript sequences from sika deer for functional genomic analyses and to identify the transcripts that demonstrate tissue-specific, age-dependent differential expression patterns. These sequences could enhance our understanding of the molecular mechanisms underlying sika deer growth and development. In the present study, we performed de novo transcriptome assembly and profiling analysis across ten tissue types and four developmental stages (juvenile, adolescent, adult, and aged) of sika deer, using Illumina paired-end tag (PET) sequencing technology. A total of 1,752,253 contigs with an average length of 799 bp were generated, from which 1,348,618 unigenes with an average length of 590 bp were defined. Approximately 33.2 % of these (447,931 unigenes) were then annotated in public protein databases. Many sika deer tissue-specific, age-dependent unigenes were identified. The testes have the largest number of tissue-enriched unigenes, and some of them were prone to develop new functions for other tissues. Additionally, our transcriptome revealed that the juvenile-adolescent transition was the most complex and important stage of the sika deer life cycle. The present work represents the first multiple tissue transcriptome analysis of sika deer across four developmental stages. The generated data not only provide a functional genomics resource for future biological research on sika deer but also guide the selection and manipulation of genes controlling growth and development.
Collapse
Affiliation(s)
- Bo-Yin Jia
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Heng-Xing Ba
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Gui-Wu Wang
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Ying Yang
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Xue-Zhe Cui
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Ying-Hua Peng
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Jun-Jun Zheng
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Xiu-Mei Xing
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China
| | - Fu-He Yang
- State Key Laboratory for Molecular Biology of Special Economical Animals, Institute of Special Economic Animals and Plants, Chinese Academy of Agricultural Sciences, 4899 Juye Street, Changchun, 130112, China.
| |
Collapse
|
11
|
Gao Y, Wang X, Yan H, Zeng J, Ma S, Niu Y, Zhou G, Jiang Y, Chen Y. Comparative Transcriptome Analysis of Fetal Skin Reveals Key Genes Related to Hair Follicle Morphogenesis in Cashmere Goats. PLoS One 2016; 11:e0151118. [PMID: 26959817 PMCID: PMC4784850 DOI: 10.1371/journal.pone.0151118] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2015] [Accepted: 02/22/2016] [Indexed: 11/18/2022] Open
Abstract
Cashmere goat skin contains two types of hair follicles (HF): primary hair follicles (PHF) and secondary hair follicles (SHF). Although multiple genetic determinants associated with HF formation have been identified, the molecules that determine the independent morphogenesis of HF in cashmere goats remain elusive. The growth and development of SHF directly influence the quantity and quality of cashmere production. Here, we report the transcriptome profiling analysis of nine skin samples from cashmere goats using 60- and 120-day-old embryos (E60 and E120, respectively), as well as newborns (NB), through RNA-sequencing (RNA-seq). HF morphological changes indicated that PHF were initiated at E60, with maturation from E120, while differentiation of SHF was identified at E120 until formation of cashmere occurred after birth (NB). The RNA-sequencing analysis generated over 20.6 million clean reads from each mRNA library. The number of differentially expressed genes (DEGs) in E60 vs. E120, E120 vs. NB, and E60 vs. NB were 1,024, 0 and 1,801, respectively, indicating that no significant differences were found at transcriptomic levels between E120 and NB. Key genes including B4GALT4, TNC, a-integrin, and FGFR1, were up-regulated and expressed in HF initiation from E60 to E120, while regulatory genes such as GPRC5D, PAD3, HOXC13, PRR9, VSIG8, LRRC15, LHX2, MSX-2, and FOXN1 were up-regulated and expressed in HF keratinisation and hair shaft differentiation from E120 and NB to E60. Several genes belonging to the KRT and KRTAP gene families were detected throughout the three HF developmental stages. The transcriptional trajectory analyses of all DEGs indicated that immune privilege, glycosaminoglycan biosynthesis, extracellular matrix receptor interaction, and growth factor receptors all played dominant roles in the epithelial-mesenchymal interface and HF formation. We found that the Wnt, transforming growth factor-beta/bone morphogenetic protein, and Notch family members played vital roles in HF differentiation and maturation. The DEGs we found could be attributed to the generation and development of HF, and thus will be critically important for improving the quantity and quality of fleece production in animals for fibres.
Collapse
Affiliation(s)
- Ye Gao
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
- College of Life Science, Yulin University, Yulin, People’s Republic of China
| | - Xiaolong Wang
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
| | - Hailong Yan
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
- College of Life Science, Yulin University, Yulin, People’s Republic of China
| | - Jie Zeng
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
| | - Sen Ma
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
| | - Yiyuan Niu
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
| | - Guangxian Zhou
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
| | - Yu Jiang
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
| | - Yulin Chen
- College of Animal Science and Technology, Northwest A&F University, Yangling, People’s Republic of China
- * E-mail:
| |
Collapse
|
12
|
Deng T, Pang C, Lu X, Zhu P, Duan A, Tan Z, Huang J, Li H, Chen M, Liang X. De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery. PLoS One 2016; 11:e0147132. [PMID: 26766209 PMCID: PMC4713091 DOI: 10.1371/journal.pone.0147132] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2015] [Accepted: 12/29/2015] [Indexed: 01/11/2023] Open
Abstract
The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.
Collapse
Affiliation(s)
- Tingxian Deng
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Chunying Pang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Xingrong Lu
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Peng Zhu
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Anqin Duan
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Zhengzhun Tan
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Jian Huang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Hui Li
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Mingtan Chen
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Xianwei Liang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
- * E-mail:
| |
Collapse
|
13
|
De Novo Assembly and Transcriptome Characterization of Canine Retina Using High-Throughput Sequencing. GENETICS RESEARCH INTERNATIONAL 2015; 2015:638679. [PMID: 26788372 PMCID: PMC4695645 DOI: 10.1155/2015/638679] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/31/2015] [Revised: 11/10/2015] [Accepted: 11/11/2015] [Indexed: 01/29/2023]
Abstract
We performed transcriptome sequencing of canine retinal tissue by 454 GS-FLX and Ion Torrent PGM platforms. RNA-Seq analysis by CLC Genomics Workbench mapped expression of 10,360 genes. Gene ontology analysis of retinal transcriptome revealed abundance of transcripts known to be involved in vision associated processes. The de novo assembly of the sequences using CAP3 generated 29,683 contigs with mean length of 560.9 and N50 of 619 bases. Further analysis of contigs predicted 3,827 full-length cDNAs and 29,481 (99%) open reading frames (ORFs). In addition, 3,782 contigs were assigned to 316 KEGG pathways which included melanogenesis, phototransduction, and retinol metabolism with 33, 15, and 11 contigs, respectively. Among the identified microsatellites, dinucleotide repeats were 68.84%, followed by trinucleotides, tetranucleotides, pentanucleotides, and hexanucleotides in proportions of 25.76, 9.40, 2.52, and 0.96%, respectively. This study will serve as a valuable resource for understanding the biology and function of canine retina.
Collapse
|
14
|
Manjunath S, Kumar GR, Mishra BP, Mishra B, Sahoo AP, Joshi CG, Tiwari AK, Rajak KK, Janga SC. Genomic analysis of host - Peste des petits ruminants vaccine viral transcriptome uncovers transcription factors modulating immune regulatory pathways. Vet Res 2015; 46:15. [PMID: 25827022 PMCID: PMC4337102 DOI: 10.1186/s13567-015-0153-8] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2014] [Accepted: 01/16/2015] [Indexed: 12/21/2022] Open
Abstract
Peste des petits ruminants (PPR), is an acute transboundary viral disease of economic importance, affecting goats and sheep. Mass vaccination programs around the world resulted in the decline of PPR outbreaks. Sungri 96 is a live attenuated vaccine, widely used in Northern India against PPR. This vaccine virus, isolated from goat works efficiently both in sheep and goat. Global gene expression changes under PPR vaccine virus infection are not yet well defined. Therefore, in this study we investigated the host-vaccine virus interactions by infecting the peripheral blood mononuclear cells isolated from goat with PPRV (Sungri 96 vaccine virus), to quantify the global changes in the transcriptomic signature by RNA-sequencing. Viral genome of Sungri 96 vaccine virus was assembled from the PPRV infected transcriptome confirming the infection and demonstrating the feasibility of building a complete non-host genome from the blood transcriptome. Comparison of infected transcriptome with control transcriptome revealed 985 differentially expressed genes. Functional analysis showed enrichment of immune regulatory pathways under PPRV infection. Key genes involved in immune system regulation, spliceosomal and apoptotic pathways were identified to be dysregulated. Network analysis revealed that the protein - protein interaction network among differentially expressed genes is significantly disrupted in infected state. Several genes encoding TFs that govern immune regulatory pathways were identified to co-regulate the differentially expressed genes. These data provide insights into the host - PPRV vaccine virus interactome for the first time. Our findings suggested dysregulation of immune regulatory pathways and genes encoding Transcription Factors (TFs) that govern these pathways in response to viral infection.
Collapse
|
15
|
Tchitchek N, Safronetz D, Rasmussen AL, Martens C, Virtaneva K, Porcella SF, Feldmann H, Ebihara H, Katze MG. Sequencing, annotation and analysis of the Syrian hamster (Mesocricetus auratus) transcriptome. PLoS One 2014; 9:e112617. [PMID: 25398096 PMCID: PMC4232415 DOI: 10.1371/journal.pone.0112617] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2014] [Accepted: 10/06/2014] [Indexed: 11/18/2022] Open
Abstract
Background The Syrian hamster (golden hamster, Mesocricetus auratus) is gaining importance as a new experimental animal model for multiple pathogens, including emerging zoonotic diseases such as Ebola. Nevertheless there are currently no publicly available transcriptome reference sequences or genome for this species. Results A cDNA library derived from mRNA and snRNA isolated and pooled from the brains, lungs, spleens, kidneys, livers, and hearts of three adult female Syrian hamsters was sequenced. Sequence reads were assembled into 62,482 contigs and 111,796 reads remained unassembled (singletons). This combined contig/singleton dataset, designated as the Syrian hamster transcriptome, represents a total of 60,117,204 nucleotides. Our Mesocricetus auratus Syrian hamster transcriptome mapped to 11,648 mouse transcripts representing 9,562 distinct genes, and mapped to a similar number of transcripts and genes in the rat. We identified 214 quasi-complete transcripts based on mouse annotations. Canonical pathways involved in a broad spectrum of fundamental biological processes were significantly represented in the library. The Syrian hamster transcriptome was aligned to the current release of the Chinese hamster ovary (CHO) cell transcriptome and genome to improve the genomic annotation of this species. Finally, our Syrian hamster transcriptome was aligned against 14 other rodents, primate and laurasiatheria species to gain insights about the genetic relatedness and placement of this species. Conclusions This Syrian hamster transcriptome dataset significantly improves our knowledge of the Syrian hamster's transcriptome, especially towards its future use in infectious disease research. Moreover, this library is an important resource for the wider scientific community to help improve genome annotation of the Syrian hamster and other closely related species. Furthermore, these data provide the basis for development of expression microarrays that can be used in functional genomics studies.
Collapse
Affiliation(s)
- Nicolas Tchitchek
- Department of Microbiology, University of Washington, Seattle, Washington, United States of America
| | - David Safronetz
- Laboratory of Virology, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rocky Mountain Laboratories, Hamilton, Montana, United States of America
| | - Angela L Rasmussen
- Department of Microbiology, University of Washington, Seattle, Washington, United States of America
| | - Craig Martens
- Genomics Unit, Research Technologies Section, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rocky Mountain Laboratories, Hamilton, Montana, United States of America
| | - Kimmo Virtaneva
- Genomics Unit, Research Technologies Section, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rocky Mountain Laboratories, Hamilton, Montana, United States of America
| | - Stephen F Porcella
- Genomics Unit, Research Technologies Section, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rocky Mountain Laboratories, Hamilton, Montana, United States of America
| | - Heinz Feldmann
- Laboratory of Virology, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rocky Mountain Laboratories, Hamilton, Montana, United States of America
| | - Hideki Ebihara
- Laboratory of Virology, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rocky Mountain Laboratories, Hamilton, Montana, United States of America
| | - Michael G Katze
- Department of Microbiology, University of Washington, Seattle, Washington, United States of America; Washington National Primate Research Center, University of Washington, Seattle, Washington, United States of America
| |
Collapse
|