Target-Specific Profiling of RNA m5C Methylation Level Using Amplicon Sequencing.

Mapping the position and quantifying the level of 5-methylcytosine (m5C) as a modification in different types of cellular RNA is an important objective in the field of epitranscriptomics. Bisulfite conversion has long been the gold standard for the detection of m5C in DNA, but it can also be applied to RNA. Here, we detail methods for bisulfite treatment of RNA, locus-specific PCR amplification, and detection of candidate sites by sequencing on the Illumina MiSeq platform.

A cis-regulatory-directed pipeline for the identification of genes involved in cardiac development and disease.

BACKGROUND: Congenital heart diseases are the major cause of death in newborns, but the genetic etiology of this developmental disorder is not fully known. The conventional approach to identify the disease-causing genes focuses on screening genes that display heart-specific expression during development. However, this approach would have discounted genes that are expressed widely in other tissues but may play critical roles in heart development. RESULTS: We report an efficient pipeline of genome-wide gene discovery based on the identification of a cardiac-specific cis-regulatory element signature that points to candidate genes involved in heart development and congenital heart disease. With this pipeline, we retrieve 76% of the known cardiac developmental genes and predict 35 novel genes that previously had no known connectivity to heart development. Functional validation of these novel cardiac genes by RNAi-mediated knockdown of the conserved orthologs in Drosophila cardiac tissue reveals that disrupting the activity of 71% of these genes leads to adult mortality. Among these genes, RpL14, RpS24, and Rpn8 are associated with heart phenotypes. CONCLUSIONS: Our pipeline has enabled the discovery of novel genes with roles in heart development. This workflow, which relies on screening for non-coding cis-regulatory signatures, is amenable for identifying developmental and disease genes for an organ without constraining to genes that are expressed exclusively in the organ of interest.

Loss of Foxd4 Impacts Neurulation and Cranial Neural Crest Specification During Early Head Development.

The specification of anterior head tissue in the late gastrulation mouse embryo relies on signaling cues from the visceral endoderm and anterior mesendoderm (AME). Genetic loss-of-function studies have pinpointed a critical requirement of LIM homeobox 1 (LHX1) transcription factor in these tissues for the formation of the embryonic head. Transcriptome analysis of embryos with gain-of-function LHX1 activity identified the forkhead box gene, Foxd4, as one downstream target of LHX1 in late-gastrulation E7.75 embryos. Our analysis of single-cell RNA-seq data show Foxd4 is co-expressed with Lhx1 and Foxa2 in the anterior midline tissue of E7.75 mouse embryos, and in the anterior neuroectoderm (ANE) at E8.25 alongside head organizer genes Otx2 and Hesx1. To study the role of Foxd4 during early development we used CRISPR-Cas9 gene editing in mouse embryonic stem cells (mESCs) to generate bi-allelic frameshift mutations in the coding sequence of Foxd4. In an in vitro model of the anterior neural tissues derived from Foxd4-loss of function (LOF) mESCs and extraembryonic endoderm cells, expression of head organizer genes as well as Zic1 and Zic2 was reduced, pointing to a need for FOXD4 in regulating early neuroectoderm development. Mid-gestation mouse chimeras harbouring Foxd4-LOF mESCs displayed craniofacial malformations and neural tube closure defects. Furthermore, our in vitro data showed a loss of FOXD4 impacts the expression of cranial neural crest markers Twist1 and Sox9. Our findings have demonstrated that FOXD4 is essential in the AME and later in the ANE for rostral neural tube closure and neural crest specification during head development.

Multiple links between 5-methylcytosine content of mRNA and translation.

BACKGROUND: 5-Methylcytosine (m5C) is a prevalent base modification in tRNA and rRNA but it also occurs more broadly in the transcriptome, including in mRNA, where it serves incompletely understood molecular functions. In pursuit of potential links of m5C with mRNA translation, we performed polysome profiling of human HeLa cell lysates and subjected RNA from resultant fractions to efficient bisulfite conversion followed by RNA sequencing (bsRNA-seq). Bioinformatic filters for rigorous site calling were devised to reduce technical noise. RESULTS: We obtained ~ 1000 candidate m5C sites in the wider transcriptome, most of which were found in mRNA. Multiple novel sites were validated by amplicon-specific bsRNA-seq in independent samples of either human HeLa, LNCaP and PrEC cells. Furthermore, RNAi-mediated depletion of either the NSUN2 or TRDMT1 m5C:RNA methyltransferases showed a clear dependence on NSUN2 for the majority of tested sites in both mRNAs and noncoding RNAs. Candidate m5C sites in mRNAs are enriched in 5'UTRs and near start codons and are embedded in a local context reminiscent of the NSUN2-dependent m5C sites found in the variable loop of tRNA. Analysing mRNA sites across the polysome profile revealed that modification levels, at bulk and for many individual sites, were inversely correlated with ribosome association. CONCLUSIONS: Our findings emphasise the major role of NSUN2 in placing the m5C mark transcriptome-wide. We further present evidence that substantiates a functional interdependence of cytosine methylation level with mRNA translation. Additionally, we identify several compelling candidate sites for future mechanistic analysis.

Mechanistic insights from the LHX1-driven molecular network in building the embryonic head.

Development of an embryo is driven by a series of molecular instructions that control the differentiation of tissue precursor cells and shape the tissues into major body parts. LIM homeobox 1 (LHX1) is a transcription factor that plays a major role in the development of the embryonic head of the mouse. Loss of LHX1 function disrupts the morphogenetic movement of head tissue precursors and impacts on the function of molecular factors in modulating the activity of the WNT signaling pathway. LHX1 acts with a transcription factor complex to regulate the transcription of target genes in multiple phases of development and in a range of embryonic tissues of the mouse and Xenopus. Determining the interacting factors and transcriptional targets of LHX1 will be key to unraveling the ensemble of factors involved in head development and building a head gene regulatory network.

Gene Editing of Mouse Embryonic and Epiblast Stem Cells.

Efficient and reliable methods for gene editing are critical for the generation of loss-of-gene function stem cells and genetically modified mice. Here, we outline the application of CRISPR-Cas9 technology for gene editing in mouse embryonic stem cells (mESCs) to generate knockout ESC chimeras for the fast-tracked analysis of gene function. Furthermore, we describe the application of gene editing directly to mouse epiblast stem cells (mEpiSCs) for modelling germ layer differentiation in vitro.

A gene regulatory network anchored by LIM homeobox 1 for embryonic head development.

Development of the embryonic head is driven by the activity of gene regulatory networks of transcription factors. LHX1 is a homeobox transcription factor that plays an essential role in the formation of the embryonic head. The loss of LHX1 function results in anterior truncation of the embryo caused by the disruption of morphogenetic movement of tissue precursors and the dysregulation of WNT signaling activity. Profiling the gene expression pattern in the Lhx1 mutant embryo revealed that tissues in anterior germ layers acquire posterior tissue characteristics, suggesting LHX1 activity is required for the allocation and patterning of head precursor tissues. Here, we used LHX1 as an entry point to delineate its transcriptional targets and interactors and construct a LHX1-anchored gene regulatory network. Using a gain-of-function approach, we identified genes that immediately respond to Lhx1 activation. Meta-analysis of the datasets of LHX1-responsive genes and genes expressed in the anterior tissues of mouse embryos at head-fold stage, in conjunction with published Xenopus embryonic LHX1 (Xlim1) ChIP-seq data, has pinpointed the putative transcriptional targets of LHX1 and an array of genetic determinants functioning together in the formation of the mouse embryonic head.

Transcriptome-Wide Mapping of RNA 5-Methylcytosine in Arabidopsis mRNAs and Noncoding RNAs.

Posttranscriptional methylation of RNA cytosine residues to 5-methylcytosine (m5C) is an important modification with diverse roles, such as regulating stress responses, stem cell proliferation, and RNA metabolism. Here, we used RNA bisulfite sequencing for transcriptome-wide quantitative mapping of m5C in the model plant Arabidopsis thaliana We discovered more than a thousand m5C sites in Arabidopsis mRNAs, long noncoding RNAs, and other noncoding RNAs across three tissue types (siliques, seedling shoots, and roots) and validated a number of these sites. Quantitative differences in methylated sites between these three tissues suggest tissue-specific regulation of m5C. Perturbing the RNA m5C methyltransferase TRM4B resulted in the loss of m5C sites on mRNAs and noncoding RNAs and reduced the stability of tRNAAsp(GTC) We also demonstrate the importance of m5C in plant development, as trm4b mutants have shorter primary roots than the wild type due to reduced cell division in the root apical meristem. In addition, trm4b mutants show increased sensitivity to oxidative stress. Finally, we provide insights into the targeting mechanism of TRM4B by demonstrating that a 50-nucleotide sequence flanking m5C C3349 in MAIGO5 mRNA is sufficient to confer methylation of a transgene reporter in Nicotiana benthamiana.

The emerging epitranscriptomics of long noncoding RNAs.

The pervasive transcription of genomes into long noncoding RNAs has been amply demonstrated in recent years and garnered much attention. Similarly, emerging 'epitranscriptomics' research has shown that chemically modified nucleosides, thought to be largely the domain of tRNAs and other infrastructural RNAs, are far more widespread and can exert unexpected influence on RNA utilization. Both areas are characterized by the often-ephemeral nature of the subject matter in that few individual examples have been fully assessed for their molecular or cellular function, and effects might often be subtle and cumulative. Here we review available information at the intersection of these two exciting areas of biology, by focusing on four RNA modifications that have been mapped transcriptome-wide: 5-methylcytidine, N6-methyladenosine, pseudouridine as well as adenosine to inosine (A-to-I) editing, and their incidence and function in long noncoding RNAs. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

Nucleotide-Level Profiling of m5C RNA Methylation.

Mapping the position and quantifying the level of 5-methylcytosine (m(5)C) as a modification in different types of cellular RNA is an important objective in the emerging field of epitranscriptomics. Bisulfite conversion has long been the gold standard for detection of m(5)C in DNA but it can also be applied to RNA. Here, we detail methods for bisulfite treatment of RNA, locus-specific PCR amplification and detection of candidate sites by sequencing on the Illumina MiSeq platform.

Mapping and significance of the mRNA methylome.

Internal methylation of eukaryotic mRNAs in the form of N6-methyladenosine (m(6)A) and 5-methylcytidine (m(5)C) has long been known to exist, but progress in understanding its role was hampered by difficulties in identifying individual sites. This was recently overcome by high-throughput sequencing-based methods that mapped thousands of sites for both modifications throughout mammalian transcriptomes, with most sites found in mRNAs. The topology of m(6)A in mouse and human revealed both conserved and variable sites as well as plasticity in response to extracellular cues. Within mRNAs, m(5)C and m(6)A sites were relatively depleted in coding sequences and enriched in untranslated regions, suggesting functional interactions with post-transcriptional gene control. Finer distribution analyses and preexisting literature point toward roles in the regulation of mRNA splicing, translation, or decay, through an interplay with RNA-binding proteins and microRNAs. The methyltransferase (MTase) METTL3 'writes' m(6)A marks on mRNA, whereas the demethylase FTO can 'erase' them. The RNA:m(5)C MTases NSUN2 and TRDMT1 have roles in tRNA methylation but they also act on mRNA. Proper functioning of these enzymes is important in development and there are clear links to human disease. For instance, a common variant of FTO is a risk allele for obesity carried by 1 billion people worldwide and mutations cause a lethal syndrome with growth retardation and brain deficits. NSUN2 is linked to cancer and stem cell biology and mutations cause intellectual disability. In this review, we summarize the advances, open questions, and intriguing possibilities in this emerging field that might be called RNA modomics or epitranscriptomics.

Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA.

The modified base 5-methylcytosine (m(5)C) is well studied in DNA, but investigations of its prevalence in cellular RNA have been largely confined to tRNA and rRNA. In animals, the two m(5)C methyltransferases NSUN2 and TRDMT1 are known to modify specific tRNAs and have roles in the control of cell growth and differentiation. To map modified cytosine sites across a human transcriptome, we coupled bisulfite conversion of cellular RNA with next-generation sequencing. We confirmed 21 of the 28 previously known m(5)C sites in human tRNAs and identified 234 novel tRNA candidate sites, mostly in anticipated structural positions. Surprisingly, we discovered 10,275 sites in mRNAs and other non-coding RNAs. We observed that distribution of modified cytosines between RNA types was not random; within mRNAs they were enriched in the untranslated regions and near Argonaute binding regions. We also identified five new sites modified by NSUN2, broadening its known substrate range to another tRNA, the RPPH1 subunit of RNase P and two mRNAs. Our data demonstrates the widespread presence of modified cytosines throughout coding and non-coding sequences in a transcriptome, suggesting a broader role of this modification in the post-transcriptional control of cellular RNA function.

MicroRNA alterations and associated aberrant DNA methylation patterns across multiple sample types in oral squamous cell carcinoma.

BACKGROUND: MicroRNA (miRNA) expression is broadly altered in cancer, but few studies have investigated miRNA deregulation in oral squamous cell carcinoma (OSCC). Epigenetic mechanisms are involved in the regulation of >30 miRNA genes in a range of tissues, and we aimed to investigate this further in OSCC. METHODS: TaqMan® qRT-PCR arrays and individual assays were used to profile miRNA expression in a panel of 25 tumors with matched adjacent tissues from patients with OSCC, and 8 control paired oral stroma and epithelium from healthy volunteers. Associated DNA methylation changes of candidate epigenetically deregulated miRNA genes were measured in the same samples using the MassArray® mass spectrometry platform. MiRNA expression and DNA methylation changes were also investigated in FACS sorted CD44(high) oral cancer stem cells from primary tumor samples (CSCs), and in oral rinse and saliva from 15 OSCC patients and 7 healthy volunteers. RESULTS: MiRNA expression patterns were consistent in healthy oral epithelium and stroma, but broadly altered in both tumor and adjacent tissue from OSCC patients. MiR-375 is repressed and miR-127 activated in OSCC, and we confirm previous reports of miR-137 hypermethylation in oral cancer. The miR-200 s/miR-205 were epigenetically activated in tumors vs normal tissues, but repressed in the absence of DNA hypermethylation specifically in CD44(high) oral CSCs. Aberrant miR-375 and miR-200a expression and miR-200c-141 methylation could be detected in and distinguish OSCC patient oral rinse and saliva from healthy volunteers, suggesting a potential clinical application for OSCC specific miRNA signatures in oral fluids. CONCLUSIONS: MiRNA expression and DNA methylation changes are a common event in OSCC, and we suggest miR-375, miR-127, miR-137, the miR-200 family and miR-205 as promising candidates for future investigations. Although overall activated in OSCC, miR-200/miR-205 suppression in oral CSCs indicate that cell specific silencing of these miRNAs may drive tumor expansion and progression.

Discovery pipeline for epigenetically deregulated miRNAs in cancer: integration of primary miRNA transcription.

BACKGROUND: Cancer is commonly associated with widespread disruption of DNA methylation, chromatin modification and miRNA expression. In this study, we established a robust discovery pipeline to identify epigenetically deregulated miRNAs in cancer. RESULTS: Using an integrative approach that combines primary transcription, genome-wide DNA methylation and H3K9Ac marks with microRNA (miRNA) expression, we identified miRNA genes that were epigenetically modified in cancer. We find miR-205, miR-21, and miR-196b to be epigenetically repressed, and miR-615 epigenetically activated in prostate cancer cells. CONCLUSIONS: We show that detecting changes in primary miRNA transcription levels is a valuable method for detection of local epigenetic modifications that are associated with changes in mature miRNA expression.