Genome-wide analysis of DNA methylation patterns in horse.

BACKGROUND: DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known. RESULTS: We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density. CONCLUSIONS: We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.

Novel variable number of tandem repeats of gibbon MAOA gene and its evolutionary significance.

Variable number of tandem repeats (VNTRs) are scattered throughout the primate genome, and genetic variation of these VNTRs have been accumulated during primate radiation. Here, we analyzed VNTRs upstream of the monoamine oxidase A (MAOA) gene in 11 different gibbon species. An abundance of truncated VNTR sequences and copy number differences were observed compared to those of human VNTR sequences. To better understand the biological role of these VNTRs, a luciferase activity assay was conducted and results indicated that selected VNTR sequences of the MAOA gene from human and three different gibbon species (Hylobates klossii, Hylobates lar, and Nomascus concolor) showed silencing ability. Together, these data could be useful for understanding the evolutionary history and functional significance of MAOA VNTR sequences in gibbon species.

Epigenetic regulation of miR-29a/miR-30c/DNMT3A axis controls SOD2 and mitochondrial oxidative stress in human mesenchymal stem cells.

The use of human mesenchymal stem cells (hMSCs) in clinical applications requires large-scale cell expansion prior to administration. However, the prolonged culture of hMSCs results in cellular senescence, impairing their proliferation and therapeutic potentials. To understand the role of microRNAs (miRNAs) in regulating cellular senescence in hMSCs, we globally depleted miRNAs by silencing the DiGeorge syndrome critical region 8 (DGCR8) gene, an essential component of miRNA biogenesis. DGCR8 knockdown hMSCs exhibited severe proliferation defects and senescence-associated alterations, including increased levels of reactive oxygen species (ROS). Transcriptomic analysis revealed that the antioxidant gene superoxide dismutase 2 (SOD2) was significantly downregulated in DGCR8 knockdown hMSCs. Moreover, we found that DGCR8 silencing in hMSCs resulted in hypermethylation in CpG islands upstream of SOD2. 5-aza-2'-deoxycytidine treatment restored SOD2 expression and ROS levels. We also found that these effects were dependent on the epigenetic regulator DNA methyltransferase 3 alpha (DNMT3A). Using computational and experimental approaches, we demonstrated that DNMT3A expression was regulated by miR-29a-3p and miR-30c-5p. Overexpression of miR-29a-3p and/or miR-30c-5p reduced ROS levels in DGCR8 knockdown hMSCs and rescued proliferation defects, mitochondrial dysfunction, and premature senescence. Our findings provide novel insights into hMSCs senescence regulation by the miR-29a-3p/miR-30c-5p/DNMT3A/SOD2 axis.

Lactate Activates the E2F Pathway to Promote Cell Motility by Up-Regulating Microtubule Modulating Genes.

Excess lactate production due to enhanced aerobic glycolysis is characteristic of malignant cancers, which is also intimately associated with poor cancer prognoses. Although tumor-associated lactate contributes to all major steps in carcinogenesis, its action mechanism remains obscure. To understand the molecular mechanism of the lactate-induced tumor metastatic process, we identified an array of lactate-responsive genes via transcriptome analysis of a metformin-induced hyper-glycolytic liver cancer model. Gene set enrichment analysis suggested E2F-RB pathway as the dominant regulator of the lactate-induced gene expression. We experimentally verified that lactate indeed activates E2F-mediated transcription by promoting E2F1 protein accumulation through a posttranscriptional mechanism. Literature-based analysis of target pathways potentially modulated by 136 top-ranked genes indicated that genes functioning in cell-cell or cell-matrix communications dominate the lactate-induced gene expression. Especially, those regulating microtubule functions, including a group of kinesin family members, were significantly up-regulated in lactate- and E2F1-dependent manners. Depletion of E2F1 or kinesins (KIF2C, KIF18B, KIF20A) led to deformation of microtubule structures, impairing cell motility as much as the deficit in lactate production. These results indicate that E2F pathway activation by tumor-associated lactate and subsequent transcriptional activation of microtubule functions play crucial roles in tumor metastasis, providing mechanistic clues to cell motility-directed anti-cancer strategies.

Activity analysis of LTR12C as an effective regulatory element of the RAE1 gene.

Ribonucleic acid export 1 (RAE1) plays an important role in the export of mature mRNAs from the nucleus to the cytoplasm. Long terminal repeats (LTRs) became integrated into the human genome during primate evolution. One such repeat element, LTR12C, lies within a predicted regulatory region located upstream of the RAE1 gene. We examined the transcriptional activity of LTR12C by using the luciferase assay, and showed that the tandem repeat region (TRR) located within LTR12C was required for its regulatory function. A bioinformatics analysis revealed that the LTR12C element had multiple transcription factor binding sites specific for nuclear transcription factor Y (NF-Y), and the promoter activity of LTR12C was significantly decreased after NF-Y knockdown. Additionally, we discovered novel data indicating that LTR12C was initially inserted into the gorilla genome. Taken together, our results reveal that the TRR of LTR12C has powerful regulatory activity due to its NF-Y binding sites, and the integration of the LTR12C element into the primate genome during evolution may have affected RAE1 transcription.

Crystal Structure of Streptococcus pyogenes Cas1 and Its Interaction with Csn2 in the Type II CRISPR-Cas System.

CRISPRs and Cas proteins constitute an RNA-guided microbial immune system against invading nucleic acids. Cas1 is a universal Cas protein found in all three types of CRISPR-Cas systems, and its role is implicated in new spacer acquisition during CRISPR-mediated adaptive immunity. Here, we report the crystal structure of Streptococcus pyogenes Cas1 (SpCas1) in a type II CRISPR-Cas system and characterize its interaction with S. pyogenes Csn2 (SpCsn2). The SpCas1 structure reveals a unique conformational state distinct from type I Cas1 structures, resulting in a more extensive dimerization interface, a more globular overall structure, and a disruption of potential metal-binding sites for catalysis. We demonstrate that SpCas1 directly interacts with SpCsn2, and identify the binding interface and key residues for Cas complex formation. These results provide structural information for a type II Cas1 protein, and lay a foundation for studying multiprotein Cas complexes functioning in type II CRISPR-Cas systems.

Downregulation of UHRF1 promotes EMT via inducing CXCR4 in human cancer cells.

Activation of epithelial-mesenchymal transition (EMT) is important for malignant tumor progression exhibiting migratory and invasive properties. UHRF1 (ubiquitin-like, with PHD and RING finger domains 1), as an epigenetic regulator, plays a crucial role in DNA CpG methylation, chromatin remodeling and gene expression. Many studies demonstrated that UHRF1 is aberrantly expressed in various types of human cancer. However, the precise role of UHRF1 in human cancers remains highly controversial. In the present study, we found that downregulation of UHRF1 enhances the migratory and invasive properties of human cancer cells by inducing EMT, and that the CXCR4 signaling pathway is strictly necessary for UHRF1 deficiency-mediated induction of EMT. Downregulation of UHRF1 induced the expression of the EMT-regulating transcription factors, Zeb1, Slug and Snail and then led to decreased protein level of E-cadherin, and increased protein level of N-cadherin and vimentin, including increased migratory and invasive properties of human cancer cells. In addition, siRNA targeting of Zeb1 or Snail effectively attenuated UHRF1 deficiency-induced EMT, but siRNA targeting of Slug did not, indicating that Zeb1 and Snail play key roles in this event. Moreover, downregulation of UHRF1 induced the expression of CXCR4 in HepG2 cells. siRNA targeting of CXCR4 greatly suppressed the UHRF1 deficiency-induced EMT, as evidenced by a reversal of expression patterns of Snail and Zeb1, and by reduced migratory and invasive properties of HepG2 cells. In conclusion, our results demonstrate that downregulation of UHRF1 contributes to the induction of EMT in human cancer cells via the activation of CXCR4 signaling pathway. Our observation also suggests that UHRF1 may play a pivotal role in suppressing the malignant alteration of cancer cells.

Genome-wide analysis of DNA methylation before-and after exercise in the thoroughbred horse with MeDIP-Seq.

Athletic performance is an important criteria used for the selection of superior horses. However, little is known about exercise-related epigenetic processes in the horse. DNA methylation is a key mechanism for regulating gene expression in response to environmental changes. We carried out comparative genomic analysis of genome-wide DNA methylation profiles in the blood samples of two different thoroughbred horses before and after exercise by methylated-DNA immunoprecipitation sequencing (MeDIP-Seq). Differentially methylated regions (DMRs) in the pre-and post-exercise blood samples of superior and inferior horses were identified. Exercise altered the methylation patterns. After 30 min of exercise, 596 genes were hypomethylated and 715 genes were hypermethylated in the superior horse, whereas in the inferior horse, 868 genes were hypomethylated and 794 genes were hypermethylated. These genes were analyzed based on gene ontology (GO) annotations and the exercise-related pathway patterns in the two horses were compared. After exercise, gene regions related to cell division and adhesion were hypermethylated in the superior horse, whereas regions related to cell signaling and transport were hypermethylated in the inferior horse. Analysis of the distribution of methylated CpG islands confirmed the hypomethylation in the gene-body methylation regions after exercise. The methylation patterns of transposable elements also changed after exercise. Long interspersed nuclear elements (LINEs) showed abundance of DMRs. Collectively, our results serve as a basis to study exercise-based reprogramming of epigenetic traits.

No expression of porcine endogenous retrovirus after pig to monkey xenotransplantation.

This study was performed to investigate the expression of two porcine endogenous retrovirus (PERV) elements, PERV gag and full-length conserved PERV, in blood cells collected periodically from organ-recipient monkeys that underwent pig to non-human primate xenotransplantation. The heart and kidney-respectively acquired from α-1,3-galactosyltransferase knockout (GT-KO) pigs that survived for24 and 25 days-were xenografted into cynomolgus monkeys. The two PERV elements expressed in the xenografted GT-KO pig organs were not present in the blood cells of the recipient monkeys. In the present study, we deduced that PERVs are not transmitted during GT-KO pig to monkey xenotransplantation.

A New Exon Derived from a Mammalian Apparent LTR Retrotransposon of the SUPT16H Gene.

The SUPT16H gene known as FACTP140 is required for the transcription of other genes. For transcription, genes need to be complexed with accessory factors, including transcription factors and RNA polymerase II. One such factor, FACT, interacts with histones H2A/H2B for nucleosome disassembly and transcription elongation. The SUPT16H gene has a transcript and many expressed sequence tags (ESTs). We were especially interested in an MaLR-derived transcript (EST, BX333035) that included a new exon introduced by a transposable element, a mammalian apparent LTR retrotransposon (MaLR). The MaLR was detected ranging from humans to galagos, indicating the MaLR in the SUPT16H gene is integrated into the primate ancestor genome. A new exon was created by alternative donor site provided by the MaLR. The original transcript and the MaLR-derived transcript were expressed in various human, rhesus monkey, and other primate tissues. Additionally, we identified a new alternative transcript that included the MaLR, but there was no significant difference in the expression of the original transcript and the MaLR-derived transcript. Interestingly, the new alternative transcript and the MaLR-derived transcript had the MaLR sequence in the new exon, but they had different structures by adopting different 3' splice sites. From this study, we verified transposable elements that contributed to transcriptome diversity.

Alu-related transcript of TJP2 gene as a marker for colorectal cancer.

Identifying diagnostic markers is very important in cancer biology. Transposable elements (TEs) account for 45% of the human genome and are usually silent in normal tissues. Contrary to this, in cancer tissues, TEs have been found to be expressed specifically with genomic instability. Here, we investigated the transcripts related to TEs using bioinformatic tools and RT-PCR analysis in tumor and adjacent normal tissues as well as in cancer cell lines. From this analysis, we identified the Alu-related transcript of TJP2 gene (TJP2-Alu transcript) that was differentially expressed between tumor and normal tissues. Its expression was higher in colon cancer cell lines and tumor tissues of colorectal cancer patients. These results could contribute to a better understanding of the involvement of TEs in human colorectal cancer, and to evaluating their potential as diagnostic markers for colorectal cancer.

Retroelements: molecular features and implications for disease.

Eukaryotic genomes comprise numerous retroelements that have a major impact on the structure and regulation of gene function. Retroelements are regulated by epigenetic controls, and they generate multiple miRNAs that are involved in the induction and progression of genomic instability. Elucidation of the biological roles of retroelements deserves continuous investigation to better understand their evolutionary features and implications for disease.

Structure and Expression Analyses of SVA Elements in Relation to Functional Genes.

SINE-VNTR-Alu (SVA) elements are present in hominoid primates and are divided into 6 subfamilies (SVA-A to SVA-F) and active in the human population. Using a bioinformatic tool, 22 SVA element-associated genes are identified in the human genome. In an analysis of genomic structure, SVA elements are detected in the 5' untranslated region (UTR) of HGSNAT (SVA-B), MRGPRX3 (SVA-D), HYAL1 (SVA-F), TCHH (SVA-F), and ATXN2L (SVA-F) genes, while some elements are observed in the 3'UTR of SPICE1 (SVA-B), TDRKH (SVA-C), GOSR1 (SVA-D), BBS5 (SVA-D), NEK5 (SVA-D), ABHD2 (SVA-F), C1QTNF7 (SVA-F), ORC6L (SVA-F), TMEM69 (SVA-F), and CCDC137 (SVA-F) genes. They could contribute to exon extension or supplying poly A signals. LEPR (SVA-C), ALOX5 (SVA-D), PDS5B (SVA-D), and ABCA10 (SVA-F) genes also showed alternative transcripts by SVA exonization events. Dominant expression of HYAL1_SVA appeared in lung tissues, while HYAL1_noSVA showed ubiquitous expression in various human tissues. Expression of both transcripts (TDRKH_SVA and TDRKH_noSVA) of the TDRKH gene appeared to be ubiquitous. Taken together, these data suggest that SVA elements cause transcript isoforms that contribute to modulation of gene regulation in various human tissues.

Identification and promoter analysis of PERV LTR subtypes in NIH-miniature pig.

Porcine endogenous retroviruses (PERVs) are integrated into the genomes of all pigs. Since some PERVs can also infect human cells, they represent a potential risk for xenotransplantation involving pig cells or organs. The long terminal repeat (LTR) elements of PERVs show promoter activity that can affect human functional genes; therefore, we examined these elements in this study. We detected several expressed LTRs in the NIH-miniature pig liver, among which we identified 9 different subtypes. When these LTRs were compared, distinct structures that contained several insertion and deletion (INDEL) events and tandem repeats were identified in the U3 region. The transcriptional activity of the 9 LTR subtypes was analyzed in the PK15 porcine cell line and in the HepG2 and Hep3B human liver cell lines, and transcriptional regulation was found to be different in the 3 cell lines. The D LTR subtype was found to have stronger promoter activity than all other types in 4 different human cell lines (HepG2, Hep3B, U251, and 293). Using computational approaches, the D type was shown to contain 4 transcription factor-binding sites distinct from those in the U3 regions of the other subtypes. Therefore, deletion mutants were constructed and examined by a transient transfection luciferase assay. The results of this analysis indicated that the binding site for the Hand1:E47 transcription factor might play a positive role in the transcriptional regulation of PERV LTR subtype D in human liver cell lines.

Promoter activity analysis and methylation characterization of LTR elements of PERVs in NIH miniature pig.

The potential risk of porcine endogenous retrovirus (PERV) transmission is an important issue in xenotransplantation (pig-to-human transplantation). Long terminal repeats (LTRs) in PERV elements show promoter activity that could affect neighboring functional genes. The methylation status and promoter activities of 3 LTR structures (PERV-LTR1, LTR2, and LTR3 elements) belonging to the PERV-A family were examined using luciferase reporter genes in human liver cell lines (HepG2 and Hep3B). The PERV LTR3 element exhibited hypomethylation and stronger promoter activity than the other LTR elements in human liver cells. We also performed comparative sequences analysis of the PERV LTR elements by using bioinformatics tools. Our findings showed that several transcription factors such as Nkx2-2 and Elk-1 positively influenced the high transcriptional activity of the PERV LTR3 element.

Radiation-induced human endogenous retrovirus (HERV)-R env gene expression by epigenetic control.

It is commonly accepted that ionizing radiation induces genomic instability by changes in genomic structure, epigenetic regulation and gene expression. Human endogenous retroviruses (HERV)-R also are often differentially expressed between normal and disease tissues under unstable genomic conditions and are implicated in the pathogenesis of several human diseases. To understand the influence of ionizing radiation on HERV-R expression, we performed quantitative reverse transcription-polymerase chain reaction (RT-PCR) analyses using γ-irradiated normal human cells. Compared to nonirradiated cells, HERV-R expression was up-regulated in γ-irradiated cells. The regulatory mechanism of HERV-R expression in irradiated cells was investigated by methylation analyses of HERV-R 5'LTRs and treatment with garcinol. These data indicated that the up-regulated transcription of HERV-R may be regulated by radiation-induced epigenetic changes induced by histone modification, and thus could be of great importance for understanding the relationship between radiation-induced biological effects and transposable elements.

Quantitative analysis of transcript variants of CHM gene containing LTR12C element in humans.

Choroideremia (CHM) is essential for the posttranslational activation of retina-specific Rab protein. Transcript variants (a and b) of the CHM gene were detected in human cancer cells and tissues. Sequence analysis of the both variants found that isoform b is caused by an LTR12C element offering an alternative splicing site within the CHM gene. We performed real-time RT-PCR analysis to study expression levels of the CHM transcript variants in tumor and adjacent normal tissue samples. Our results showed that CHM isoform b was highly expressed in tumor tissues but its expression levels were relatively low in those of adjacent normal tissues including colon, testis, and lung. In addition, high expression levels of CHM isoform b were detected in the cancer cell lines of colon and lung, and colon cancer patient tissues. Thus, we suggest that expression levels of alternative transcripts of the CHM gene could be used as a molecular marker system to identify human cancers.

Analysis of the molecular and regulatory properties of active porcine endogenous retrovirus gamma-1 long terminal repeats in kidney tissues of the NIH-Miniature pig.

The pig genome contains the gamma 1 family of porcine endogenous retroviruses (PERVs), which are a major obstacle to the development of successful xenotransplantation from pig to human. Long terminal repeats (LTRs) found in PERVs are known to be essential elements for the control of the transcriptional activity of single virus by different transcription factors (TFs). To identify transcribed PERV LTR elements, RT-PCR and DNA sequencing analyses were performed. Twenty-nine actively transcribed LTR elements were identified in the kidney tissues of the NIH-Miniature pig. These elements were divided into two major groups (I and II), and four minor groups (I-1, I-2, I-3, and II-1), by the presence of insertion and deletion (INDEL) sequences. Group I elements showed strong transcriptional activity compared to group II elements. Four different LTR elements (PL1, PL2, PL3, and PL4) as representative of the groups were analyzed by using a transient transfection assay. The regulation of their promoter activity was investigated by treatment with M.SssI (CpG DNA methyltransferase) and garcinol (histone acetyltransferase inhibitor). The transcriptional activity of PERV LTR elements was significantly reduced by treatment with M.SssI. These data indicate that transcribed PERV LTR elements harbor sufficient promoter activity to regulate the transcription of a single virus, and the transcriptional activity of PERV LTRs may be controlled by DNA methylation events.