Molecular characterization, expression patterns and polymorphism analysis of porcine Six1 gene

Chronological Expression of PITX2 and SIX1 Genes and the Association between Their Polymorphisms and Chicken Meat Quality Traits

Meat quality is closely related to the development of skeletal muscle, in which PITX2 and SIX1 genes play important regulatory roles. The present study firstly provided the data of chronological expression files of PITX2 and SIX1 genes in the post-hatching pectoral muscle and analyzed the association of their polymorphisms with the meat quality traits of Wuliang Mountain Black-bone (WLMB) chickens. The results showed that both PITX2 and SIX1 genes were weakly expressed in the second and third weeks, and then increased significantly from the third week to the fourth week. Furthermore, there was a significant positive correlation between the expression levels of the two genes. Twelve and one SNPs were detected in the chicken PITX2 and SIX1 genes, respectively, of which four SNPs (g.9830C > T, g.10073C > T, g.13335G > A, g.13726A > G) of the PITX2 gene and one SNP (g.564G > A) of the SIX1 gene were significantly associated with chicken meat quality traits. For the PITX2 gene, chickens with the CT genotype of g.9830C > T showed the highest meat color L*, shear force (SF), pH, and the lowest electrical conductivity (EC), and drip loss (DL) (p < 0.05 or p < 0.01); chickens with the CC genotype of g.10073C > T had the lowest L*, pH, and the highest DL (p < 0.01). For the SIX1 gene, chickens with the GG genotype of g.564G > A had the highest (p < 0.05) SF and pH. Furthermore, pH had a significant correlation with all the other meat quality traits. The current study could contribute to the research of regulatory mechanisms of meat quality and lay the foundation for improving meat quality based on marker-assisted selection in chickens.

DNA methylation across the genome in aged human skeletal muscle tissue and muscle-derived cells: the role of HOX genes and physical activity

Skeletal muscle tissue demonstrates global hypermethylation with age. However, methylome changes across the time-course of differentiation in aged human muscle derived cells, and larger coverage arrays in aged muscle tissue have not been undertaken. Using 850K DNA methylation arrays we compared the methylomes of young (27 ± 4.4 years) and aged (83 ± 4 years) human skeletal muscle and that of young/aged heterogenous muscle-derived human primary cells (HDMCs) over several time points of differentiation (0, 72 h, 7, 10 days). Aged muscle tissue was hypermethylated compared with young tissue, enriched for; pathways-in-cancer (including; focal adhesion, MAPK signaling, PI3K-Akt-mTOR signaling, p53 signaling, Jak-STAT signaling, TGF-beta and notch signaling), rap1-signaling, axon-guidance and hippo-signalling. Aged cells also demonstrated a hypermethylated profile in pathways; axon-guidance, adherens-junction and calcium-signaling, particularly at later timepoints of myotube formation, corresponding with reduced morphological differentiation and reductions in MyoD/Myogenin gene expression compared with young cells. While young cells showed little alterations in DNA methylation during differentiation, aged cells demonstrated extensive and significantly altered DNA methylation, particularly at 7 days of differentiation and most notably in focal adhesion and PI3K-AKT signalling pathways. While the methylomes were vastly different between muscle tissue and HDMCs, we identified a small number of CpG sites showing a hypermethylated state with age, in both muscle tissue and cells on genes KIF15, DYRK2, FHL2, MRPS33, ABCA17P. Most notably, differential methylation analysis of chromosomal regions identified three locations containing enrichment of 6–8 CpGs in the HOX family of genes altered with age. With HOXD10, HOXD9, HOXD8, HOXA3, HOXC9, HOXB1, HOXB3, HOXC-AS2 and HOXC10 all hypermethylated in aged tissue. In aged cells the same HOX genes (and additionally HOXC-AS3) displayed the most variable methylation at 7 days of differentiation versus young cells, with HOXD8, HOXC9, HOXB1 and HOXC-AS3 hypermethylated and HOXC10 and HOXC-AS2 hypomethylated. We also determined that there was an inverse relationship between DNA methylation and gene expression for HOXB1, HOXA3 and HOXC-AS3. Finally, increased physical activity in young adults was associated with oppositely regulating HOXB1 and HOXA3 methylation compared with age. Overall, we demonstrate that a considerable number of HOX genes are differentially epigenetically regulated in aged human skeletal muscle and HDMCs and increased physical activity may help prevent age-related epigenetic changes in these HOX genes.

Demethylation and derepression of genomic retroelements in the skeletal muscles of aged mice

Changes in DNA methylation influence the aging process and contribute to aging phenotypes, but few studies have been conducted on DNA methylation changes in conjunction with skeletal muscle aging. We explored the DNA methylation changes in a variety of retroelement families throughout aging (at 2, 20, and 28 months of age) in murine skeletal muscles by methyl‐binding domain sequencing (MBD‐seq). The two following contrasting patterns were observed among the members of each repeat family in superaged mice: (a) hypermethylation in weakly methylated retroelement copies and (b) hypomethylation in copies with relatively stronger methylation levels, representing a pattern of “regression toward the mean” within a single retroelement family. Interestingly, these patterns depended on the sizes of the copies. While the majority of the elements showed a slight increase in methylation, the larger copies (>5 kb) displayed evident demethylation. All these changes were not observed in T cells. RNA sequencing revealed a global derepression of retroelements during the late phase of aging (between 20 and 28 months of age), which temporally coincided with retroelement demethylation. Following this methylation drift trend of “regression toward the mean,” aging tended to progressively lose the preexisting methylation differences and local patterns in the genomic regions that had been elaborately established during the early period of development.

Identification of an (AC)n microsatellite in the Six1 gene promoter and its effect on production traits in Pietrain × Duroc × Landrace × Yorkshire pigs

The Sine oculis homeobox 1 (Six1) gene is important for skeletal muscle growth and fiber specification; therefore, it is considered as a promising candidate gene that may influence porcine growth and meat quality traits. Nevertheless, the association of Six1 with these processes and the mechanisms regulating its expression remain unclear. The objectives of this study were to identify variant sites of Six1 in different pig breeds, conduct association analysis to evaluate the relationship between polymorphisms of these variants and porcine production traits in Pietrain × Duroc × Landrace × Yorkshire commercial pigs, and explore the potential regulatory mechanisms of Six1 affecting production traits. A total of 12 variants were identified, including 10 single- nucleotide variations (SNVs), 1 insertion- deletion (Indel), and 1 (AC)n microsatellite. Association analysis demonstrated that the SNV, g.1595A>G, was significantly associated with meat color (redness, a*); individuals with the G allele had greater a* values (P < 0.05). Moreover, our results demonstrated that the (AC)n polymorphism in the Six1 promoter was significantly associated with weaning weight (P < 0.05), carcass weight (P < 0.05), and thoracic and lumbar back fat (P < 0.01).In addition, we found that the (AC)n variant was closely related with Six1 expression levels and demonstrated this polymorphism on promoter activity by in vitro experiments. Overall, this study provides novel evidence for elucidating the effects of Six1 on porcine production traits as promising candidate and describes two variants with these traits, which are potential reference markers for pig molecular breeding. In addition, our data on the relationship between porcine Six1 expression and the polymorphic (AC)n microsatellite in its promoter may facilitate similar studies in other species.

Characterization of the promoter region of the bovine SIX1 gene: Roles of MyoD, PAX7, CREB and MyoG

The SIX1 gene belongs to the family of six homeodomain transcription factors (TFs), that regulates the extracellular signal-regulated kinase 1/2 (ERK1/2) pathway and mediate skeletal muscle growth and regeneration. Previous studies have demonstrated that SIX1 is positively correlated with body measurement traits (BMTs). However, the transcriptional regulation of SIX1 remains unclear. In the present study, we determined that bovine SIX1 was highly expressed in the longissimus thoracis. To elucidate the molecular mechanisms involved in bovine SIX1 regulation, 2-kb of the 5' regulatory region were obtained. Sequence analysis identified neither a consensus TATA box nor a CCAAT box in the 5' flanking region of bovine SIX1. However, a CpG island was predicted in the region -235 to +658 relative to the transcriptional start site (TSS). An electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP) assay in combination with serial deletion constructs of the 5' flanking region, site-directed mutation and siRNA interference demonstrated that MyoD, PAX7 and CREB binding occur in region -689/-40 and play important roles in bovine SIX1 transcription. In addition, MyoG drives SIX1 transcription indirectly via the MEF3 motif. Taken together these interactions suggest a key functional role for SIX1 in mediating skeletal muscle growth in cattle.

NRF1 and ZSCAN10 bind to the promoter region of the SIX1 gene and their effects body measurements in Qinchuan cattle

The SIX1 homeobox gene belongs to the six homeodomain family and is widely thought to play a principal role in mediating of skeletal muscle development. In the present study, we determined that the bovine SIX1 gene was highly expressed in the longissimus thoracis and physiologically immature individuals. DNA sequencing of 428 individual Qinchuan cattle identified nine single nucleotide polymorphisms (SNPs) in the promoter region of the SIX1 gene. Using a series of 5′ deletion promoter plasmid luciferase reporter assays and 5′-rapid amplification of cDNA end analysis (RACE), two of these SNPs were found to be located in the proximal minimal promoter region −216/−28 relative to the transcriptional start site (TSS). Correlation analysis showed the combined haplotypes H₁-H₂ (-GG-GA-) was significantly greater in the body measurement traits (BMTs) than the others, which was consistent with the results showing that the transcriptional activity of Hap2 was higher than the others in Qinchuan cattle myoblast cells. Furthermore, the electrophoretic mobility shift assays (EMSA) and chromatin immunoprecipitation assay (ChIP) demonstrated that NRF1 and ZSCAN10 binding occurred in the promoter region of diplotypes H₁-H₂ to regulate SIX1 transcriptional activity. This information may be useful for molecular marker-assisted selection (MAS) in cattle breeding.

LncRNA-Six1 Encodes a Micropeptide to Activate Six1 in Cis and Is Involved in Cell Proliferation and Muscle Growth

Long non-coding RNAs (lncRNAs) play important roles in epigenetic regulation of skeletal muscle development. In our previous RNA-seq study (accession number GSE58755), we found that lncRNA-Six1 is an lncRNA that is differentially expressed between White Recessive Rock (WRR) and Xinghua (XH) chicken. In this study, we have further demonstrated that lncRNA-Six1 is located 432 bp upstream of the gene encoding the protein Six homeobox 1 (Six1). A dual-luciferase reporter assay identified that lncRNA-Six1 overlaps the Six1 proximal promoter. In lncRNA-Six1, a micropeptide of about 7.26 kDa was found to play an important role in the lncRNA-Six1 in cis activity. Overexpression of lncRNA-Six1 promoted the mRNA and protein expression level of the Six1 gene, while knockdown of lncRNA-Six1 inhibited Six1 expression. Moreover, tissue expression profiles showed that both the lncRNA-Six1 and the Six1 mRNA were highly expressed in chicken breast tissue. LncRNA-Six1 overexpression promoted cell proliferation and induced cell division. Conversely, its loss of function inhibited cell proliferation and reduced cell viability. Similar effects were observed after overexpression or knockdown of the Six1 gene. In addition, overexpression or knockdown of Six1 promoted or inhibited, respectively, the expression levels of muscle-growth-related genes, such as MYOG, MYHC, MYOD, IGF1R, and INSR. Taken together, these data demonstrate that lncRNA-Six1 carries out cis-acting regulation of the protein-encoding Six1 gene, and encodes a micropeptide to activate Six1 gene, thus promoting cell proliferation and being involved in muscle growth.

DNA methylation dynamics in muscle development and disease

DNA methylation is an essential epigenetic modification for mammalian development and is crucial for the establishment and maintenance of cellular identity. Traditionally, DNA methylation has been considered as a permanent repressive epigenetic mark. However, the application of genome-wide approaches has allowed the analysis of DNA methylation in different genomic contexts revealing a more dynamic regulation than originally thought, since active DNA methylation and demethylation occur during cellular differentiation and tissue specification. Satellite cells are the primary stem cells in adult skeletal muscle and are responsible for postnatal muscle growth, hypertrophy, and muscle regeneration. This review outlines the published data regarding DNA methylation changes along the skeletal muscle program, in both physiological and pathological conditions, to better understand the epigenetic mechanisms that control myogenesis.

The role of Six1 in the genesis of muscle cell and skeletal muscle development

The sine oculis homeobox 1 (Six1) gene encodes an evolutionarily conserved transcription factor. In the past two decades, much research has indicated that Six1 is a powerful regulator participating in skeletal muscle development. In this review, we summarized the discovery and structural characteristics of Six1 gene, and discussed the functional roles and molecular mechanisms of Six1 in myogenesis and in the formation of skeletal muscle fibers. Finally, we proposed areas of future interest for understanding Six1 gene function.

Overexpression of Six1 gene suppresses proliferation and enhances expression of fast-type muscle genes in C2C12 myoblasts

Sine oculis homeobox 1 (Six1) homeodomain transcription factor is implicated in the genesis of muscle fiber type diversity, but its regulatory mechanisms on the formation of muscle fiber type are still poorly understood. To elucidate the biological roles of Six1 gene in muscle fiber formation, we established C2C12 cell line overexpressing Six1 and determined the effects of forced Six1 expression on muscle-specific genes expression, cell proliferation, and cell cycles. Our results indicated that Six1 overexpression could significantly promote the expression of fast-type muscle genes Atp2a1, Srl, and Mylpf. Furthermore, Six1 overexpressing C2C12 cells displayed a relative lower proliferative potential, and cell cycle analysis showed that Six1 exerted its role in cell cycle primarily through the regulation of G1/S and G2/M phases. In conclusion, Six1 plays an essential role in modulation of the fast-twitch muscle fiber phenotype through up-regulating fast-type muscle genes expression, and it could suppress the proliferation of muscle cells.

Selection of reference genes for quantitative real-time RT-PCR analysis in citrus

Quantitative real-time reverse transcription polymerase chain reaction (qPCR) has become the preferred method for studying low-abundant mRNA expression. Appropriate application of qPCR in such studies requires the use of reference gene(s) as an internal control in order to normalize the mRNA levels between different samples for an exact comparison of gene expression levels. Expression of the reference gene should be independent from development stage, cell/tissue types, treatments and environmental conditions. Recognizing the importance of reference gene(s) in normalization of qPCR data, various reference genes have been evaluated for stable expression under specific conditions in various organisms. In plants, only a few of them have been investigated, and very few reports about such reference genes in citrus. In the present study, seven candidate reference genes (18SrRNA, ACTB, rpII, UBQI, UBQ10, GAPDH and TUB) were tested, and three of them (18SrRNA, ACTB and rpII) proved to be the most stable ones among six leaf samples of different citrus genotypes. The three candidate reference genes were further analyzed for their stability of expression in five different tissues, and the results indicated that they were not completely stable. It is commonly accepted that gene expression studies should be normalized using more than one reference gene. Based on our results, we propose the use of the mean result rendered by18SrRNA, ACTB and rpII as reference genes to normalize mRNA levels in qPCR analysis of diverse cultivars and tissues of citrus. These results may provide a guideline for future works on gene expression in citrus by using qPCR.