Promoter variant-dependent mRNA expression of the MEF2A in longissimus dorsi muscle in cattle

Chromatin accessibility and regulatory vocabulary across indicine cattle tissues

BACKGROUND

Spatiotemporal changes in the chromatin accessibility landscape are essential to cell differentiation, development, health, and disease. The quest of identifying regulatory elements in open chromatin regions across different tissues and developmental stages is led by large international collaborative efforts mostly focusing on model organisms, such as ENCODE. Recently, the Functional Annotation of Animal Genomes (FAANG) has been established to unravel the regulatory elements in non-model organisms, including cattle. Now, we can transition from prediction to validation by experimentally identifying the regulatory elements in tropical indicine cattle. The identification of regulatory elements, their annotation and comparison with the taurine counterpart, holds high promise to link regulatory regions to adaptability traits and improve animal productivity and welfare.

RESULTS

We generate open chromatin profiles for liver, muscle, and hypothalamus of indicine cattle through ATAC-seq. Using robust methods for motif discovery, motif enrichment and transcription factor binding sites, we identify potential master regulators of the epigenomic profile in these three tissues, namely HNF4, MEF2, and SOX factors, respectively. Integration with transcriptomic data allows us to confirm some of their target genes. Finally, by comparing our results with Bos taurus data we identify potential indicine-specific open chromatin regions and overlaps with indicine selective sweeps.

CONCLUSIONS

Our findings provide insights into the identification and analysis of regulatory elements in non-model organisms, the evolution of regulatory elements within two cattle subspecies as well as having an immediate impact on the animal genetics community in particular for a relevant productive species such as tropical cattle.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.

Study on the transcriptional regulatory mechanism of the MyoD1 gene in Guanling bovine

The MyoD1 gene plays a key role in regulating the myoblast differentiation process in the early stage of skeletal muscle development. To understand the functional elements of the promoter region and transcriptional regulation of the bovine MyoD1 gene, we cloned eight fragments from the sequence region of the MyoD1 gene promoter and inserted them into eukaryotic expression vectors for cotransfection with the mouse myoblast cell line C2C12 and Madin-Darby bovine kidney (MDBK) line. A variety of transcription factor binding sites in the longest 5'-flanking fragment from Guanling cattle MyoD1-P1 were predicted by using the online software TFSEARCH and ALGGEN PROMO as well as validated by the promoter-binding TF profiling assay II and yeast one-hybrid (Y1H) assay, including MyoD, VDR, MEF1, MEF2, SF1, and Myf6. Myf6 strongly activated the MyoD1 promoter, while MyoD1 was also capable of efficiently activating the expression of its own promoter. The transcription factors MEF2A, SF1, and VDR were further confirmed to be capable of binding to MyoD1 by Y1H system experiments. The effects of the Guanling cattle MyoD1 gene on the mRNA expression of the MEF2A, SF1, and VDR genes were determined by using a lentivirus-mediated overexpression technique, confirming that overexpression of the MyoD1 gene upregulated the mRNA expression of MEF2A as well as downregulated the expression of SF1 and VDR in the process of muscle myogenesis. Our study revealed the effects of transcription factors including MEF2A, SF1 and VDR on regulatory aspects of MyoD1, providing abundant information for transcriptional regulation of MyoD1 in muscle differentiation.

Study of bovine gene: the temporal-spatial expression patterns, polymorphism and association analysis with meat production traits

The gene () encodes a transcription factor belonging to the MEF2 family that plays an important role in myogenesis by transcriptional regulation of genes involved in skeletal muscle growth and development. Despite the established importance of the factors in the muscular growth and development, the temporal-spatial expression and biological function of have not been reported in cattle. The aim of this study was to analyze the level of expression in the developing longissimus dorsi muscle (LM) of 4 cattle breeds (Polish Holstein-Friesian [HF], Limousine [LIM], Hereford [HER], Polish Red [PR]), differing in terms of meat production and utility type, at 6, 9, and 12 mo of age. The genetic polymorphism and expression patterns in 6 tissues (heart, spleen, liver, semitendinosus muscle [ST], gluteus medius muscle [GM], and LM) were also investigated. The results showed that mRNA was expressed at a high level in adult skeletal and cardiac muscles. Moreover, expression was markedly greater in the GM than in the LM ( 0.05) and ST ( 0.01). An age-dependent and breed-specific comparison of mRNA level in skeletal muscle of HF, LIM, HER, and PR bulls showed that age was significant differentiating factor of transcript/protein abundance in the LM of HER and LIM ( 0.001) compared to HF and PR, for which the differences in mRNA level were not significant ( > 0.05). Regarding the breed effect on the expression, significantly greater mRNA/protein level was noticed in the LM of 9 and 12 mo-old HER than of LIM ( 0.01), HF ( 0.001), and PR ( 0.001). Four novel SNP, namely, (promoter), (exon 7), (exon 8), and (3'UTR), were identified. We found that 3'UTR variant, situated within the seed region of the miR-5187-3p and miR-6931-5p binding sites, was associated with the level of mRNA/protein in LM of 12-mo-old HF bulls. In addition, we observed a significant association between some carcass quality traits, including meat and carcass fatness quality traits, and various 3'UTR genotypes in the investigated population of HF cattle. Our finding provides new evidence of the significant role in the postnatal muscle growth and development in cattle, and indicates that can be a promising molecular marker for carcass quality-related traits in adult cattle.

The molecular characterization and temporal-spatial expression of myocyte enhancer factor 2 genes in the goat and their association with myofiber traits

The myocyte enhancer factor-2 (MEF2) gene family in vertebrates includes MEF2A, MEF2B, MEF2C, and MEF2D, which have important functions in the regulation of muscular growth and development. To investigate their temporal-spatial expression and functions in the goat, these genes were cloned (accession nos. JN967621-24) and their expression patterns characterized at five postnatal stages (3, 30, 60, 90, and 120days). Association analysis was then applied regarding MEF2 expression levels and myofiber diameter and density. MEF2B was shown to be weakly homologous with other species, the distant branches with other members and the lowest expression levels, suggesting that it is distinct from other family members. Expression of the other three MEF2 genes was widely distributed, but this was largely accumulated in the skeletal muscle and myocardium compared with the viscera at all developmental stages. MEF2A and MEF2D expression levels were higher overall than MEF2B and MEF2C in six tissues, and were significantly positively correlated with the myofiber diameter of the longissimus dorsi. These findings suggest that goat MEF2 genes mainly function in the skeletal muscle and myocardium, and that MEF2A and MEF2D are likely to effectively promote muscular growth and development during postnatal stages. MEF2A expression was highest in the myocardium, where MEF2C expression increased with age, implying that both gene products are related to the growth and development of postnatal myocardium.

Postnatal expression patterns and polymorphism analysis of the bovine myocyte enhancer factor 2C (Mef2C) gene

The aim of this study was to analyze the level of expression of the Mef2C gene in the developing bovine longissimus dorsi (LD) muscle (at 6, 9 and 12months of age) and to evaluate differences in expression among Polish Holstein-Friesian (HO) and Limousine (LM) bulls. Moreover, the expression patterns of Mef2C in different tissues were determined. The results showed that Mef2C mRNA was expressed at a high level in adult skeletal and cardiac muscles. Moreover, Mef2C expression was markedly lower in the semitendinosus (ST) than in the gluteus medius (GM) and LD muscles. A relatively higher Mef2C mRNA and MEF2C protein level was estimated in the muscles of HO bulls at the age of 12months in comparison with its lower expression in LM bulls. Furthermore, we found that the Mef2C promoter variant (GU211004:g.-1606C>T) does not affect the level of mRNA in the LD and ST muscles of 12-month-old HO bulls.

Effects of new polymorphisms in the bovine myocyte enhancer factor 2D (MEF2D) gene on the expression rates of the longissimus dorsi muscle

Myocyte enhancer factor 2D (MEF2D), a product of the MEF2D gene, belongs to the myocyte enhancer factor 2 (MEF2) protein family which is involved in vertebrate skeletal muscle development and differentiation during myogenesis. The aim of the present study was to search for polymorphisms in the bovine MEF2D gene and to analyze their effect on MEF2D mRNA and on protein expression levels in the longissimus dorsi muscle of Polish Holstein–Friesian cattle. Overall, three novel variations, namely, insertion/deletion g.−818_−814AGCCG and g.−211C<A transversion in the promoter region as well as g.7C<T transition in the 5′untranslated region (5′UTR), were identified by DNA sequencing. A total, 375 unrelated bulls belonging to six different cattle breeds were genotyped, and three combined genotypes (Ins-C-C/Ins-C-C, Del-A-T/Del-A-T and Ins-C-C/Del-A-T) were determined. The frequency of the combined genotype Ins-C-C/Ins-C-C and Del-A-T/Del-A-T was varied between the breeds and the average frequency was 0.521 and 0.037, respectively. Expression analysis showed that the MEF2D variants were highly correlated with MEF2D mRNA and protein levels in the longissimus dorsi muscle of Polish Holstein–Friesian bulls carrying the three different combined genotypes. The highest MEF2D mRNA and protein levels were estimated in the muscle of bulls with the Ins-C-C/Ins-C-C homozygous genotype as compared to the Del-A-T/Del-A-T homozygotes (P < 0.01) and Ins-C-C/Del-A-T heterozygotes (P < 0.05). A preliminary association study showed no significant differences in the carcass quality traits between bulls with various MEF2D combined genotypes in the investigated population of Polish Holstein–Friesian cattle.