Single nucleotide polymorphisms, haplotypes and combined genotypes in MYH₃ gene and their associations with growth and carcass traits in Qinchuan cattle

Effect of bovine myosin heavy chain 3 on proliferation and differentiation of myoblast

The myosin heavy chain 3 (MYH3) gene is an essential gene that affects muscle development. This study aimed to discuss the expression characteristics of the MYH3 gene and its effect on the proliferation and differentiation of bovine myoblasts. Quantitative real time-PCR results display that the expression level of MYH3 was higher in muscle tissue, and the expression increased in the early stage of myoblast differentiation. Interfering with the MYH3 gene in myoblasts resulted in fewer EDU-positive cells and decreased expression of proliferation marker genes. Interference with MYH3 can also affect the differentiation process of myoblasts. Regarding phenotype, myotube differentiation in the interference group was slowed or even stopped. Interference with the expression of MYH3 could significantly reduce the expression of myogenic differentiation marker genes. The above results show that MYH3 is mainly expressed in muscle tissue and is highly expressed in the early stage of differentiation of bovine myoblasts, and interfering with the MYH3 can promote the proliferation and inhibit the differentiation of bovine myoblasts. This study provides a theoretical basis for revealing the regulatory process of bovine myoblast proliferation and differentiation and bovine molecular breeding.

Transcriptome sequencing to explore the effect of miR-214 on chicken primary myoblasts

MicroRNAs are involved in a series of biological processes, such as proliferation, differentiation and apoptosis of primary myoblasts. The research group found that miR-214 is highly expressed in chicken primary myoblasts (CPMs), so we used miR-214 as a starting point to explore the biological function of miR-214 in skeletal muscle growth and development. In this experiment, CPMs were cultured in vitro; miR-214 was overexpressed in CPMs; and cell samples were collected for subsequent transcriptome sequencing (RNA-seq). After miR-214 overexpression, we identified 97 differentially expressed genes (DEGs), of which 21 DEGs were up-regulated and 76 DEGs were down-regulated. After bioinformatics analysis, these DEGs were found to be significantly enriched in myofibrils, muscle system processes, myofibril assembly and other biological processes related to muscle development. The significantly enriched KEGGs include focal adhesion and type II diabetes mellitus. The protein network of DEGs was drawn by STRING and Cytoscape software, and 5 DEGs were randomly selected to verify the sequencing results by real-time fluorescence quantification. CAV3 is not only an important node protein in the protein network but also a member of the focal adhesion signaling pathway. It is speculated that miR-214 may regulate muscle development through the focal adhesion signaling pathway.

Transcriptomic Profile of Primary Culture of Skeletal Muscle Cells Isolated from Semitendinosus Muscle of Beef and Dairy Bulls

The aim of the study was to identify differences in the transcriptomic profiles of primary muscle cell cultures derived from the semitendinosus muscle of bulls of beef breeds (Limousin (LIM) and Hereford (HER)) and a dairy breed (Holstein-Friesian (HF)) (n = 4 for each breed). Finding a common expression pattern for proliferating cells may point to such an early orientation of the cattle beef phenotype at the transcriptome level of unfused myogenic cells. To check this hypothesis, microarray analyses were performed. The analysis revealed 825 upregulated and 1300 downregulated transcripts similar in both beef breeds (LIM and HER) and significantly different when compared with the dairy breed (HF) used as a reference. Ontological analyses showed that the largest group of genes were involved in muscle organ development. Muscle cells of beef breeds showed higher expression of genes involved in myogenesis (including erbb-3, myf5, myog, des, igf-1, tgfb2) and those encoding proteins comprising the contractile apparatus (acta1, actc1, myh3, myh11, myl1, myl2, myl4, tpm1, tnnt2, tnnc1). The obtained results confirmed our hypothesis that the expression profile of several groups of genes is common in beef breeds at the level of proliferating satellite cells but differs from that observed in typical dairy breeds.

Transcriptome profile analysis of leg muscle tissues between slow- and fast-growing chickens

Chicken is widely favored by consumers because of some unique features. The leg muscles occupy an important position in the market. However, the specific mechanism for regulating muscle growth speed is not clear. In this experiment, we used Jinghai yellow chickens with different body weights at 300 days as research subjects. The chickens were divided into fast- and slow-growing groups, and we collected leg muscles after slaughtering for use in RNA-seq. After comparing the two groups, 87 differentially expressed genes (DEGs) were identified (fold change ≥ 2 and FDR < 0.05). The fast-growing group had 42 up-regulated genes and 45 down-regulated genes among these DEGs compared to the slow-growing group. Six items were significantly enriched in the biological process: embryo development ending in birth or egg hatching, chordate embryonic development, embryonic skeletal system development, and embryo development as well as responses to ketones and the sulfur compound biosynthetic process. Two significantly enriched pathways were found in the KEGG pathway analysis (P-value < 0.05): the insulin signaling pathway and the adipocytokine signaling pathway. This study provides a theoretical basis for the molecular mechanism of chicken growth and for improving the production of Jinghai yellow chicken.

Bidirectional adaptive introgression between two ecologically divergent sparrow species

Natural hybrid zones can be used to dissect the mechanisms driving key evolutionary processes by allowing us to identify genomic regions important for establishing reproductive isolation and that allow for transfer of adaptive variation. We leverage whole-genome data in a system where two bird species, the saltmarsh (Ammospiza caudacuta) and Nelson's (A. nelsoni) sparrow, hybridize despite their relatively high background genetic differentiation and past ecological divergence. Adaptive introgression is plausible in this system because Nelson's sparrows are recent colonists of saltwater marshes, compared to the specialized saltmarsh sparrow that has a longer history of saltmarsh adaptation. Comparisons among whole-genome sequences of 34 individuals from allopatric and sympatric populations show that ongoing gene flow is shaping the genomic landscape, with allopatric populations exhibiting genome-wide F_ST estimates close to double of that observed in sympatry. We characterized patterns of introgression across the genome and identify regions that exhibit biased introgression into hybrids from one parental species. These regions offer compelling candidates for genes related to tidal marsh adaptations suggesting that adaptive introgression may be an important consequence of hybridization. These findings highlight the value of considering the landscapes of both genome-wide introgression and divergence when characterizing the evolutionary forces that drive speciation.

MiR-2425-5p targets RAD9A and MYOG to regulate the proliferation and differentiation of bovine skeletal muscle-derived satellite cells

Our group previously identified miR-2425-5p, a unique bovine miRNA; however, its biological function and regulation in muscle-derived satellite cells (MDSCs) remain unclear. Herein, stem-loop RT-PCR results showed that miR-2425-5p increased during MDSCs proliferation, but decreased during differentiation. Cell proliferation was examined using EdU assays, cyclin B1 (CCNB1) and proliferating cell nuclear antigen (PCNA) western blot (WB) and flow cytometry analysis. These results showed that miR-2425-5p mimics (miR-2425-M) enhanced MDSCs proliferation, whereas, miR-2425-5p inhibitor (miR-2425-I) had opposite effect. Conversely, cell differentiation studies by desmin (DES) immunofluorescence, myotubes formation, and myosin heavy chain 3 (MYH3) WB analyses revealed that miR-2425-M and miR-2425-I blocked and promoted MDSCs differentiation, respectively. Moreover, luciferase reporter, RT-PCR, and WB assays showed that miR-2425-5p directly targeted the 3′-UTR of RAD9 homolog A (RAD9A) and myogenin (MYOG) to regulate their expression. Rescue experiment showed RAD9A inhibited the proliferation of MDSCs through miR-2425-5p. In addition, we found that miR-2425-5p expression was regulated by its host gene NCK associated protein 5-like (NCKAP5L) rather than being transcribed independently as a separate small RNA. Collectively, these data indicate that miR-2425-5p is a novel regulator of bovine MDSCs proliferation and differentiation and provides further insight into the biological functions of miRNA in this species.

Detection of genomic signatures of recent selection in commercial broiler chickens

Background

Identification of the genomic signatures of recent selection may help uncover causal polymorphisms controlling traits relevant to recent decades of selective breeding in livestock. In this study, we aimed at detecting signatures of recent selection in commercial broiler chickens using genotype information from single nucleotide polymorphisms (SNPs). A total of 565 chickens from five commercial purebred lines, including three broiler sire (male) lines and two broiler dam (female) lines, were genotyped using the 60K SNP Illumina iSelect chicken array. To detect genomic signatures of recent selection, we applied two methods based on population comparison, cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio (XP-CLR), and further analyzed the results to find genomic regions under recent selection in multiple purebred lines.

Results

A total of 321 candidate selection regions spanning approximately 1.45 % of the chicken genome in each line were detected by consensus of results of both XP-EHH and XP-CLR methods. To minimize false discovery due to genetic drift, only 42 of the candidate selection regions that were shared by 2 or more purebred lines were considered as high-confidence selection regions in the study. Of these 42 regions, 20 were 50 kb or less while 4 regions were larger than 0.5 Mb. In total, 91 genes could be found in the 42 regions, among which 19 regions contained only 1 or 2 genes, and 9 regions were located at gene deserts.

Conclusions

Our results provide a genome-wide scan of recent selection signatures in five purebred lines of commercial broiler chickens. We found several candidate genes for recent selection in multiple lines, such as SOX6 (Sex Determining Region Y-Box 6) and cTR (Thyroid hormone receptor beta). These genes may have been under recent selection due to their essential roles in growth, development and reproduction in chickens. Furthermore, our results suggest that in some candidate regions, the same or opposite alleles have been under recent selection in multiple lines. Most of the candidate genes in the selection regions are novel, and as such they should be of great interest for future research into the genetic architecture of traits relevant to modern broiler breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0430-1) contains supplementary material, which is available to authorized users.

An integrated in silico approach for functional and structural impact of non- synonymous SNPs in the MYH1 gene in Jeju Native Pigs

BACKGROUND

This study was performed to identify the non- synonymous polymorphisms in the myosin heavy chain 1 gene (MYH1) association with skeletal muscle development in economically important Jeju Native Pig (JNP) and Berkshire breeds. Herein, we present an in silico analysis, with a focus on (a) in silico approaches to predict the functional effect of non-synonymous SNP (nsSNP) in MYH1 on growth, and (b) molecular docking and dynamic simulation of MYH1 to predict the effects of those nsSNP on protein-protein association.

RESULTS

The NextGENe (V 2.3.4.) tool was used to identify the variants in MYH1 from JNP and Berkshire using RNA seq. Gene ontology analysis of MYH1 revealed significant association with muscle contraction and muscle organ development. The 95 % confidence intervals clearly indicate that the mRNA expression of MYH1 is significantly higher in the Berkshire longissimus dorsi muscle samples than JNP breed. Concordant in silico analysis of MYH1, the open-source software tools identified 4 potential nsSNP (L884T, K972C, N981G, and Q1285C) in JNP and 1 nsSNP (H973G) in Berkshire pigs. Moreover, protein-protein interactions were studied to investigate the effect of MYH1 mutations on association with hub proteins, and MYH1 was found to be closely associated with the protein myosin light chain, phosphorylatable, fast skeletal muscle MYLPF. The results of molecular docking studies on MYH1 (native and 4 mutants) and MYLFP demonstrated that the native complex showed higher electrostatic energy (-466.5 Kcal mol(-1)), van der Walls energy (-87.3 Kcal mol(-1)), and interaction energy (-835.7 Kcal mol(-1)) than the mutant complexes. Furthermore, the molecular dynamic simulation revealed that the native complex yielded a higher root-mean-square deviation (0.2-0.55 nm) and lower root-mean-square fluctuation (approximately 0.08-0.3 nm) as compared to the mutant complexes.

CONCLUSIONS

The results suggest that the variants at L884T, K972C, N981G, and Q1285C in MYH1 in JNP might represent a cause for the poor growth performance for this breed. This study is a pioneering in-depth in silico analysis of polymorphic MYH1 and will serve as a valuable resource for further targeted molecular diagnosis and population-based studies conducted for improving the growth performance of JNP.

Associations of MYH3 gene copy number variations with transcriptional expression and growth traits in Chinese cattle

Copy number variations (CNVs) have been recently recognized as another important genetic variability complementary to single nucleotide polymorphisms (SNPs). Compelling evidence has indicated that CNVs are responsible for phenotypic traits by changing the copy numbers of functional genes. Myosin heavy chain 3 (MYH3) gene is a critical regulatory factor in skeletal muscle development, and has been detected in the CNVs region by comparative genomic hybridization (CGH) array. This study was conducted to validate and detect the distribution of MYH3 copy numbers (relative to Angus cattle) in four Chinese cattle breeds (NY, QC, LX, and CY), and further to investigate the associations of the copy number changes with its transcriptional expression and cattle growth traits. Substantial genetic differences of MYH3 copy numbers were identified between NY and the other three breeds (P<0.01). The copy numbers of MYH3 gene presented the positive correlations with the transcript level of MYH3 gene in both fetal and adult skeletal muscles (P<0.05). Statistical analysis revealed that CNVs of MYH3 gene were significantly associated with growth traits of NY cattle, and the individuals with copy number gain showed better phenotypes than the loss and/or median groups (P<0.05). This study firstly attempted to establish the correlations between CNVs of candidate genes and growth traits, and our results suggested that the CNVs of MYH3 gene may be utilized as the potential markers for economic traits in selection breeding programs of Chinese cattle.

Gene expression profiles in engineered cardiac tissues respond to mechanical loading and inhibition of tyrosine kinases

Engineered cardiac tissues (ECTs) are platforms to investigate cardiomyocyte maturation and functional integration, the feasibility of generating tissues for cardiac repair, and as models for pharmacology and toxicology bioassays. ECTs rapidly mature in vitro to acquire the features of functional cardiac muscle and respond to mechanical load with increased proliferation and maturation. ECTs are now being investigated as platforms for in vitro models for human diseases and for pharmacologic screening for drug toxicities. We tested the hypothesis that global ECT gene expression patterns are complex and sensitive to mechanical loading and tyrosine kinase inhibitors similar to the maturing myocardium. We generated ECTs from day 14.5 rat embryo ventricular cells, as previously published, and then conditioned constructs after 5 days in culture for 48 h with mechanical stretch (5%, 0.5 Hz) and/or the p38 MAPK (p38 mitogen-activated protein kinase) inhibitor BIRB796. RNA was isolated from individual ECTs and assayed using a standard Agilent rat 4 × 44k V3 microarray and Pathway Analysis software for transcript expression fold changes and changes in regulatory molecules and networks. Changes in expression were confirmed by quantitative-polymerase chain reaction (q-PCR) for selected regulatory molecules. At the threshold of a 1.5-fold change in expression, stretch altered 1559 transcripts, versus 1411 for BIRB796, and 1846 for stretch plus BIRB796. As anticipated, top pathways altered in response to these stimuli include cellular development, cellular growth and proliferation; tissue development; cell death, cell signaling, and small molecule biochemistry as well as numerous other pathways. Thus, ECTs display a broad spectrum of altered gene expression in response to mechanical load and/or tyrosine kinase inhibition, reflecting a complex regulation of proliferation, differentiation, and architectural alignment of cardiomyocytes and noncardiomyocytes within ECT.

Genetic diversity of MYH3 gene associated with growth and carcass traits in Chinese Qinchuan cattle

MYH3, whose function is to convert chemical energy to mechanical energy through ATP hydrolysis, is mainly expressed in skeletal muscle at various stages and is indispensable in the procedure of development of skeletal muscle and heart. In the study, genetic variations and genotypes of MYH 3 gene in a total of 365 Qinchuan cattles were analyzed by polymerase chain reaction-restriction fragment length polymorphism, as well as verified the effect on growth and carcass traits. After PCR products were digested by restriction enzymes, eight SNPs were identified and individuals were genotyped. It showed that the SNPs at nucleotides were all in low linkage disequilibrium, therefore no dominated haplotype was found in the population. The result of statistic analysis indicated seven SNPs were significantly associated with growth and carcass traits (P < 0.05, N = 365) except locus G13791A. To sum up, the result of the study proved that polymorphisms in MYH3 gene are associated with the growth performance of Chinese Qinchuan cattle, so the variations of the gene could be used as possible molecular assisted-makers in the beef cattle breeding program and management.