Yongzhen Huang

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.

Comparative Transcriptome Profiling of Dairy Goat MicroRNAs from Dry Period and Peak Lactation Mammary Gland Tissues

Background MicroRNAs (miRNAs) are small noncoding RNA molecules that serve as important post-transcriptional gene expression regulators by targeting messenger RNAs for post-transcriptional endonucleolytic cleavage or translational inhibition. miRNAs play important roles in many biological processes. Extensive high-throughput sequencing studies of miRNAs have been performed in several animal models. However, little is known about the diversity of these regulatory RNAs in goat (Capra hircus), which is one of the most important agricultural animals and the oldest domesticated species raised worldwide. Goats have long been used for their milk, meat, hair (including cashmere), and skins throughout much of the world. Results In this study, two small RNA libraries were constructed based on dry period and peak lactation dairy goat mammary gland tissues and sequenced using the Illumina-Solexa high-throughput sequencing technology. A total of 346 conserved and 95 novel miRNAs were identified in the dairy goat. miRNAs expression was confirmed by qRT-PCR in nine tissues and in the mammary gland during different stages of lactation. In addition, several candidate miRNAs that may be involved in mammary gland development and lactation were found by comparing the miRNA expression profiles in different tissues and developmental stages of the mammary gland. Conclusions This study reveals the first miRNAs profile related to the biology of the mammary gland in the dairy goat. The characterization of these miRNAs could contribute to a better understanding of the molecular mechanisms of lactation physiology and mammary gland development in the dairy goat.

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.

Two Novel Coding SNPs of SREBP1c Gene are Associated with Body Weight and Average Daily Gain in Bovine

It is known that the SREBP1c gene is an important gene responsible for adipogenesis and regulation of the expression of genes controlling fatty acid biosynthesis. Its expression levels increase in parallel with obesity. Therefore, the present study focused on screening the genetic variation within bovine SREBP1c gene and analyzing its effect on growth traits in 1035 individuals belonging to four Chinese cattle breeds (QC, NY, JX, CH) using PCR-SSCP, DNA sequencing, and forced PCR-RFLP methods. The results revealed two novel mutations: NC_007317: g. 10781 C > A (457aa), 10914 G > A (502aa). Association analysis with growth traits in the Nangyang breed indicated that: The SNPs in the bovine SREBP1c gene had significant effects on body weight and average daily gain at birth, 6 and 12 months old (P < 0.05 or P < 0.01). Therefore, these results suggest that the SREBP1c gene is a strong candidate gene that affects growth traits in cattle.

The developmental transcriptome sequencing of bovine skeletal muscle reveals a long noncoding RNA, lncMD, promotes muscle differentiation by sponging miR-125b.

Pervasive transcription of the mammalian genome generates numerous long noncoding RNAs (lncRNAs), which are of crucial importance in diverse biological processes. Recent advances in high throughput sequencing technology have helped to accelerate the pace of lncRNA discovery. However, no study on the overall expression patterns of lncRNAs during muscle development has been conducted. We reported here the first analysis of lncRNA landscape in bovine embryonic, neonatal and adult skeletal muscle using Ribo-Zero RNA-Seq, a technology which can capture both poly(A)+ and poly(A)- transcripts. We finally defined 7692 high-confidence lncRNAs and uncovered 401 lncRNAs differentially expressed among three developmental stages, including lncMD, a novel muscle-specific lncRNA which is gradually up-regulated during myoblast differentiation. lncMD overexpression upregulated, whereas lncMD silencing decreased the expression of two well-established myogenic markers, myosin heavy chain (MHC) and myogenin (MyoG). In-depth analyses showed that lncMD acts as a molecular sponge for miR-125b and that insulin-like growth factor 2 (IGF2) is a direct target of miR-125b in cattle. Moreover, lncMD level was positively correlated with IGF2 mRNA level in bovine muscle tissues, a vital corollary to ceRNA function. Altogether, our research showed that lncMD acts as a ceRNA to sequester miR-125b, leading to heightened IGF2 expression and thus promotes muscle differentiation. Our findings also complement the reference genome annotation of cattle, which will likely be useful for further functional lncRNA cloning and more comprehensive studies on lncRNA regulation in muscle development.

Intragenic DNA methylation status down-regulates bovine IGF2 gene expression in different developmental stages

DNA methylation is a key epigenetic modification in mammals and has an essential and important role in muscle development. Insulin-like growth factor 2 (IGF2) is a fetal growth and differentiation factor that plays an important role in muscle growth and in myoblast proliferation and differentiation. The aim of this study was to evaluate the expression of IGF2 and the methylation pattern on the differentially methylated region (DMR) of the last exon of IGF2 in six tissues with two different developmental stages. The DNA methylation pattern was compared using bisulfite sequencing polymerase chain reaction (BSP) and combined bisulfite restriction analysis (COBRA). The quantitative real-time PCR (qPCR) analysis indicated that IGF2 has a broad tissue distribution and the adult bovine group showed significant lower mRNA expression levels than that in the fetal bovine group (P<0.05 or P<0.01). Moreover, the DNA methylation level analysis showed that the adult bovine group exhibited a significantly higher DNA methylation levels than that in the fetal bovine group (P<0.05 or P<0.01). These results indicate that IGF2 expression levels were negatively associated with the methylation status of the IGF2 DMR during the two developmental stages. Our results suggest that the methylation pattern in this DMR may be a useful parameter to investigate as a marker-assisted selection for muscle developmental in beef cattle breeding program and as a model for studies in other species.

Circular RNA profiling reveals an abundant circLMO7 that regulates myoblasts differentiation and survival by sponging miR-378a-3p

Circular RNAs (circRNAs) have been identified from various tissues and species, but their regulatory functions during developmental processes are not well understood. We examined circRNA expression profiles of two developmental stages of bovine skeletal muscle (embryonic and adult musculus longissimus) to provide first insights into their potential involvement in bovine myogenesis. We identified 12 981 circRNAs and annotated them to the Bos taurus reference genome, including 530 circular intronic RNAs (ciRNAs). One parental gene could generate multiple circRNA isoforms, with only one or two isoforms being expressed at higher expression levels. Also, several host genes produced different isoforms when comparing development stages. Most circRNA candidates contained two to seven exons, and genomic distances to back-splicing sites were usually less than 50 kb. The length of upstream or downstream flanking introns was usually less than 105 nt (mean≈11 000 nt). Several circRNAs differed in abundance between developmental stages, and real-time quantitative PCR (qPCR) analysis largely confirmed differential expression of the 17 circRNAs included in this analysis. The second part of our study characterized the role of circLMO7—one of the most down-regulated circRNAs when comparing adult to embryonic muscle tissue—in bovine muscle development. Overexpression of circLMO7 inhibited the differentiation of primary bovine myoblasts, and it appears to function as a competing endogenous RNA for miR-378a-3p, whose involvement in bovine muscle development has been characterized beforehand. Congruent with our interpretation, circLMO7 increased the number of myoblasts in the S-phase of the cell cycle and decreased the proportion of cells in the G0/G1 phase. Moreover, it promoted the proliferation of myoblasts and protected them from apoptosis. Our study provides novel insights into the regulatory mechanisms underlying skeletal muscle development and identifies a number of circRNAs whose regulatory potential will need to be explored in the future.