Chuzhao Lei

Genome-wide DNA Methylation Profiles and Their Relationships with mRNA and the microRNA Transcriptome in Bovine Muscle Tissue (Bos taurine)

DNA methylation is a key epigenetic modification in mammals and plays important roles in muscle development. We sampled longissimus dorsi muscle (LDM) from a well-known elite native breed of Chinese Qinchuan cattle living within the same environment but displaying distinct skeletal muscle at the fetal and adult stages. We generated and provided a genome-wide landscape of DNA methylomes and their relationship with mRNA and miRNA for fetal and adult muscle studies. Integration analysis revealed a total of 77 and 1,054 negatively correlated genes with methylation in the promoter and gene body regions, respectively, in both the fetal and adult bovine libraries. Furthermore, we identified expression patterns of high-read genes that exhibit a negative correlation between methylation and expression from nine different tissues at multiple developmental stages of bovine muscle-related tissue or organs. In addition, we validated the MeDIP-Seq results by bisulfite sequencing PCR (BSP) in some of the differentially methylated promoters. Together, these results provide valuable data for future biomedical research and genomic and epigenomic studies of bovine skeletal muscle that may help uncover the molecular basis underlying economically valuable traits in cattle. This comprehensive map also provides a solid basis for exploring the epigenetic mechanisms of muscle growth and development.

Morphological and genetic evidence for early Holocene cattle management in northeastern China

The domestication of cattle is generally accepted to have taken place in two independent centres: around 10,500 years ago in the Near East, giving rise to modern taurine cattle, and two millennia later in southern Asia, giving rise to zebu cattle. Here we provide firmly dated morphological and genetic evidence for early Holocene management of taurine cattle in northeastern China. We describe conjoining mandibles from this region that show evidence of oral stereotypy, dated to the early Holocene by two independent ¹⁴C dates. Using Illumina high-throughput sequencing coupled with DNA hybridization capture, we characterize 15,406 bp of the mitogenome with on average 16.7-fold coverage. Phylogenetic analyses reveal a hitherto unknown mitochondrial haplogroup that falls outside the known taurine diversity. Our data suggest that the first attempts to manage cattle in northern China predate the introduction of domestic cattle that gave rise to the current stock by several thousand years.

Multiple maternal origins of native modern and ancient horse populations in China

To obtain more knowledge of the origin and genetic diversity of domestic horses in China, this study provides a comprehensive analysis of mitochondrial DNA (mtDNA) D-loop sequence diversity from nine horse breeds in China in conjunction with ancient DNA data and evidence from archaeological and historical records. A 247-bp mitochondrial D-loop sequence from 182 modern samples revealed a total of 70 haplotypes with a high level of genetic diversity. Seven major mtDNA haplogroups (A-G) and 16 clusters were identified for the 182 Chinese modern horses. In the present study, nine 247-bp mitochondrial D-loop sequences of ancient remains of Bronze Age horse from the Chifeng region of Inner Mongolia in China (c. 4000-2000a bp) were used to explore the origin and diversity of Chinese modern horses and the phylogenetic relationship between ancient and modern horses. The nine ancient horses carried seven haplotypes with rich genetic diversity, which were clustered together with modern individuals among haplogroups A, E and F. Modern domestic horse and ancient horse data support the multiple origins of domestic horses in China. This study supports the argument that multiple successful events of horse domestication, including separate introductions of wild mares into the domestic herds, may have occurred in antiquity, and that China cannot be excluded from these events. Indeed, the association of Far Eastern mtDNA types to haplogroup F was highly significant using Fishers exact test of independence (P = 0.00002), lending support for Chinese domestication of this haplogroup. High diversity and all seven mtDNA haplogroups (A-G) with 16 clusters also suggest that further work is necessary to shed more light on horse domestication in China.

Identification and profiling of conserved and novel microRNAs from Chinese Qinchuan bovine longissimus thoracis

BackgroundMicroRNAs (miRNAs) are a family of ~22 nucleotide small RNA molecules that regulate gene expression by fully or partially binding to their complementary sequences. Recently, a large number of miRNAs and their expression patterns have been identified in various species. However, to date no miRNAs have been reported to modulate muscle development in beef cattle.ResultsTotal RNAs from the Chinese Qinchuan bovine longissimus thoracis at fetal and adult stages were used to construct small RNA libraries for Solexa SBS technology sequencing. A total of 15,454,182 clean reads were obtained from the fetal bovine library and 13,558,164 clean reads from the adult bovine library. In total, 521 miRNAs including 104 novel miRNA candidates were identified. Furthermore, the nucleotide bias, base edit and family of the known miRNAs were also analyzed. Based on stem-loop qPCR, 25 high-read miRNAs were detected, and the results showed that bta-miRNA-206, miRNA-1, miRNA-133, miRNAn12, and miRNAn17 were highly expressed in muscle-related tissue or organs, suggesting that these miRNAs may play a role in the development of bovine muscle tissues.ConclusionsThis study confirmed the authenticity of 417 known miRNAs, discovered 104 novel miRNAs in bos taurus, and identified five muscle-specific miRNAs. The identification of novel miRNAs significantly expanded the repertoire of bovine miRNAs and could contribute to further studies on the muscle development of cattle.

Long non-coding RNA ADNCR suppresses adipogenic differentiation by targeting miR-204

Adipogenesis is a complex and precisely orchestrated process mediated by a network of adipogenic regulatory factors. Several studies have highlighted the relevance of lncRNAs in adipocyte differentiation, but the precise molecular mechanism has largely remained elusive. In the present study, we performed Ribo-Zero RNA-Seq to investigate both the poly(A)+and poly(A)-lncRNAs of in vitro cultured bovine preadipocytes and differentiated adipocytes. A stringent set of 2882 lncRNAs was finally identified. A comparison of the lncRNAs expression profiles revealed that 16 lncRNAs are differentially expressed during adipocyte differentiation. We focused on the most downregulated lncRNA, which we named adipocyte differentiation-associated long noncoding RNA (ADNCR). Mechanistically, ADNCR inhibited adipocyte differentiation by functioning as a competing endogenous RNA (ceRNA) for miR-204, thereby augmenting the expression of the miR-204 target gene, SIRT1, which is known to inhibit adipocyte differentiation and adipogenic gene expression by docking with NCoR and SMART to repress PPARγ activity. Our data not only provide a valuable genomic resource for the identification of lncRNAs with functional roles in adipocyte differentiation but also reveal new insights into understanding the mechanisms of adipogenic differentiation.

Haplotype distribution in the GLI3 gene and their associations with growth traits in cattle

The glioma-associated oncogene family zinc finger 3 gene (GLI3) mediates in all vertebrates hedgehog (Hh) signaling that plays an essential role in the induction and patterning of numerous cell types during invertebrate and vertebrate development. In this study, a total of 6 single nucleotide polymorphisms (SNPs: 1-6) were identified by polymerase chain reaction-single stranded conformational polymorphism (PCR-SSCP) and DNA pool sequencing, including all 13 exons and 12 exon-intron boundaries within the bovine GLI3 gene. 16 haplotypes and 13 combined genotypes were revealed and the linkage disequilibrium was assessed in 708 individuals representing three main cattle breeds from China. The statistical analyses indicated that the SNP2, 3 and 4 are associated with the body weight at birth and 6 months in Nanyang cattle population (P<0.05). No significant association was detected between 11 combined genotypes and body weight at five different ages. Our results provide evidence that polymorphisms in the GLI3 gene are associated with growth traits, and may be used for marker-assisted selection in beef cattle breeding program.

Novel 12-bp deletion in the coding region of the bovine NPM1 gene affects growth traits.

The nucleophosmin 1 gene (NPM1) encodes a multifunctional nucleolar phosphoprotein that plays a crucial role in the control of various aspects of cell growth and homeostasis. In this study, the coding region of theNPM1 gene was screened in 1035 individuals of 4 Chinese cattle breeds by DNA sequencing and polyacrylamide gel electrophoresis. A novel 12-bp deletion mutation was identified in the coding region of theNPM1 gene. The PCR products of primerNPM1-P2 exhibited 3 genotypes and 2 alleles: 178 bp (denoted asW) and 166 bp (denoted asD). GenotypeDD and alleleD were predominant in the studied populations. Association analysis with growth traits in the Nanyang breed (N = 265) showed that the animals with genotypeDD had significantly greater birth weight, body weight, body length, and heart girth than those with genotypeWD (P < 0.01 orP < 0.05) at birth and after 6 months and 12 months, but not at 18 and 24 months of age. Results of this study suggest that theNPM1 gene is a candidate gene for growth traits in cattle.

Detection of copy number variations and their effects in Chinese bulls

BackgroundCopy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.ResultsWe identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.ConclusionsThe results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Sequence variants in the bovine nucleophosmin 1 gene, their linkage and their associations with body weight in native cattle breeds in China

The nucleophosmin (nucleolar phosphoprotein B23, numatrin) gene (NPM1, previously known as nucleophosmin/nucleoplasmin family, member 1) encodes a multifunctional nucleolar phosphoprotein that plays a crucial role in cell growth and homeostasis. Seven sequence variants (SVs) were identified in the coding region of bovine NPM1, five of which were in complete linkage disequilibrium. Eight different haplotypes were identified, of which two major haplotypes have a frequency of 23.2% and 20.4%. Three SVs were significantly associated with body weight in the Nanyang population as analysed at different ages. No significant association was detected between 18 combined genotypes and body weight at five different ages. Our results suggest that some polymorphisms in NPM1 are associated with body weight at some ages and may be used as candidates for marker-assisted selection and management in beef cattle breeding programmes.

Genetic variation of mitochondrial D-loop region and evolution analysis in some Chinese cattle breeds.

The complete mitochondrial D-loop region from 123 individuals in 12 Chinese cattle breeds and two individuals in Germany Yellow cattle breed was sequenced and analyzed. The results were shown as follows: 93 variations and 57 haplotypes were detected, and the average number of nucleotide difference was 22.708, nucleotide diversity (d) was 0.0251 +/- 0.00479, and haplotype diversity (Hd) was 0.888 +/- 0.026, indicating very high genetic diversity in Chinese cattle breeds. In the Neighbor-Joining tree, 13 cattle breeds were divided into two main clades, Bos taurus and Bos indicus; new Clade ? had only one sequence from Apeijiaza cattle breed in Tibet, which was similar to that of yak at a higher level than other cattle breeds, proving the introgression of genes from the yak. The proportions of Bos taurus and Bos indicus were 64.3% and 35.7% respectively in Xigazê Humped cattle breed, and 50.0% and 50.0% respectively in Apeijiaza cattle breed, which revealed that Tibet cattle included Bos indicus haplotypes. The importance of Yunnan cattle in the origin of Chinese cattle was also confirmed based on their abundant haplotypes. Then, a very special haplotype i1 discovered in 27 Chinese cattle breeds, including 11 breeds in this study and 16 breeds in the GenBank, played the role of a nucleus in Chinese zebu and was further discussed. At the same time, the construction of Chinese zebu core group based on haplotype i1 validated the distinct origin of Bos indicus in Tibet, which was different from that of the other cattle breeds with zebu haplotypes in China.