Genome-Wide Profiling of DNA Methylation Reveals Transposon Targets of CHROMOMETHYLASE3

DNA methylation has been implicated in a variety of epigenetic processes, and abnormal methylation patterns have been seen in tumors 1, 2, 3. Analysis of methylation patterns has traditionally been conducted either by using Southern analysis after cleavage with methyl-sensitive restriction endonucleases or by bisulfite sequencing [4]. However, neither method is practical for analyzing more than a few genes. Here, we describe a simple technique for genome-wide mapping of DNA methylation patterns. Fragmentation by a methyl-sensitive restriction endonuclease is followed by size fractionation and hybridization to microarrays. We demonstrate the utility of this method by characterizing methylation patterns in Arabidopsis methylation mutants. This analysis reveals that CHROMOMETHYLASE3 (CMT3) [5], which was previously shown to maintain CpXpG methylation 6, 7, preferentially methylates transposons, even when they are present as single copies within the genome. Methylation profiling has potential applications in disease research and diagnostic screening.     

Genome-Wide Sequencing Reveals Small Nucleolar RNAs Downregulated in Cerebral Cavernous Malformations

Cerebral cavernous malformations (CCM) are vascular malformations associated with abnormally dilated blood vessels and leaky capillaries that often result in hemorrhages. Despite recent advances, precise understanding of the cellular and molecular mechanism leading to the pathogenesis of CCM remains elusive. Emerging evidence indicates that small nucleolar RNAs (snoRNAs), belonging to the class of non-coding RNAs, may play a significant role as diagnostic markers in human diseases. However, there is no report till date that studied the role of snoRNAs in CCM biology. The objective of the current study was to identify snoRNAs associated with CCM pathogenesis. Using genome-wide small RNA sequencing, we identified a total of 271 snoRNAs reliably expressed in CCM. By applying additional statistical stringency, three snoRNAs (SNORD115-32, SNORD114-22, and SNORD113-3) were found to be significantly downregulated in CCM patient tissue samples (n?=?3) as compared to healthy brains (n?=?3). Deregulation of the selected snoRNAs was further validated by qRT-PCR. Further, cellular localization via in situ hybridization also confirmed robust reduction in the expression of SNORD115-32 and SNORD114-22 in CCM tissues as compared to the healthy controls. By applying high-throughput sequencing and cellular localization analyses, we report here for the first time the genome-wide expression profile of snoRNAs in CCM tissues and a robust downregulation of candidate snoRNAs in CCM conditions. Future studies should warrant the screening in large CCM patient cohorts and will be helpful in the development of potential biomarkers and improved clinical diagnosis.     

Genome-Wide Analysis Reveals Novel Regulators of Growth in Drosophila melanogaster

Organismal size depends on the interplay between genetic and environmental factors. Genome-wide association (GWA) analyses in humans have implied many genes in the control of height but suffer from the inability to control the environment. Genetic analyses in Drosophila have identified conserved signaling pathways controlling size; however, how these pathways control phenotypic diversity is unclear. We performed GWA of size traits using the Drosophila Genetic Reference Panel of inbred, sequenced lines. We find that the top associated variants differ between traits and sexes; do not map to canonical growth pathway genes, but can be linked to these by epistasis analysis; and are enriched for genes and putative enhancers. Performing GWA on well-studied developmental traits under controlled conditions expands our understanding of developmental processes underlying phenotypic diversity.     

Genome-Wide Association Study of Metabolic Traits Reveals Novel Gene-Metabolite-Disease Links

Metabolic traits are molecular phenotypes that can drive clinical phenotypes and may predict disease progression. Here, we report results from a metabolome- and genome-wide association study on ¹H-NMR urine metabolic profiles. The study was conducted within an untargeted approach, employing a novel method for compound identification. From our discovery cohort of 835 Caucasian individuals who participated in the CoLaus study, we identified 139 suggestively significant (P<5×10⁻⁸) and independent associations between single nucleotide polymorphisms (SNP) and metabolome features. Fifty-six of these associations replicated in the TasteSensomics cohort, comprising 601 individuals from São Paulo of vastly diverse ethnic background. They correspond to eleven gene-metabolite associations, six of which had been previously identified in the urine metabolome and three in the serum metabolome. Our key novel findings are the associations of two SNPs with NMR spectral signatures pointing to fucose (rs492602, P = 6.9×10⁻⁴⁴) and lysine (rs8101881, P = 1.2×10⁻³³), respectively. Fine-mapping of the first locus pinpointed the FUT2 gene, which encodes a fucosyltransferase enzyme and has previously been associated with Crohn''s disease. This implicates fucose as a potential prognostic disease marker, for which there is already published evidence from a mouse model. The second SNP lies within the SLC7A9 gene, rare mutations of which have been linked to severe kidney damage. The replication of previous associations and our new discoveries demonstrate the potential of untargeted metabolomics GWAS to robustly identify molecular disease markers.     

Genome-Wide Double-Stranded RNA Sequencing Reveals the Functional Significance of Base-Paired RNAs in Arabidopsis

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA–producing substrates of RNA–DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA–DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.     

Genome-Wide Mutagenesis Reveals That ORF7 Is a Novel VZV Skin-Tropic Factor

The Varicella Zoster Virus (VZV) is a ubiquitous human alpha-herpesvirus that is the causative agent of chicken pox and shingles. Although an attenuated VZV vaccine (v-Oka) has been widely used in children in the United States, chicken pox outbreaks are still seen, and the shingles vaccine only reduces the risk of shingles by 50%. Therefore, VZV still remains an important public health concern. Knowledge of VZV replication and pathogenesis remains limited due to its highly cell-associated nature in cultured cells, the difficulty of generating recombinant viruses, and VZV''s almost exclusive tropism for human cells and tissues. In order to circumvent these hurdles, we cloned the entire VZV (p-Oka) genome into a bacterial artificial chromosome that included a dual-reporter system (GFP and luciferase reporter genes). We used PCR-based mutagenesis and the homologous recombination system in the E. coli to individually delete each of the genome''s 70 unique ORFs. The collection of viral mutants obtained was systematically examined both in MeWo cells and in cultured human fetal skin organ samples. We use our genome-wide deletion library to provide novel functional annotations to 51% of the VZV proteome. We found 44 out of 70 VZV ORFs to be essential for viral replication. Among the 26 non-essential ORF deletion mutants, eight have discernable growth defects in MeWo. Interestingly, four ORFs were found to be required for viral replication in skin organ cultures, but not in MeWo cells, suggesting their potential roles as skin tropism factors. One of the genes (ORF7) has never been described as a skin tropic factor. The global profiling of the VZV genome gives further insights into the replication and pathogenesis of this virus, which can lead to improved prevention and therapy of chicken pox and shingles.     

Joint Genome-Wide Profiling of miRNA and mRNA Expression in Alzheimer's Disease Cortex Reveals Altered miRNA Regulation

Although microRNAs are being extensively studied for their involvement in cancer and development, little is known about their roles in Alzheimer''s disease (AD). In this study, we used microarrays for the first joint profiling and analysis of miRNAs and mRNAs expression in brain cortex from AD and age-matched control subjects. These data provided the unique opportunity to study the relationship between miRNA and mRNA expression in normal and AD brains. Using a non-parametric analysis, we showed that the levels of many miRNAs can be either positively or negatively correlated with those of their target mRNAs. Comparative analysis with independent cancer datasets showed that such miRNA-mRNA expression correlations are not static, but rather context-dependent. Subsequently, we identified a large set of miRNA-mRNA associations that are changed in AD versus control, highlighting AD-specific changes in the miRNA regulatory system. Our results demonstrate a robust relationship between the levels of miRNAs and those of their targets in the brain. This has implications in the study of the molecular pathology of AD, as well as miRNA biology in general.     

Whole Genome Transcript Profiling of Drug Induced Steatosis in Rats Reveals a Gene Signature Predictive of Outcome

Drug induced steatosis (DIS) is characterised by excess triglyceride accumulation in the form of lipid droplets (LD) in liver cells. To explore mechanisms underlying DIS we interrogated the publically available microarray data from the Japanese Toxicogenomics Project (TGP) to study comprehensively whole genome gene expression changes in the liver of treated rats. For this purpose a total of 17 and 12 drugs which are diverse in molecular structure and mode of action were considered based on their ability to cause either steatosis or phospholipidosis, respectively, while 7 drugs served as negative controls. In our efforts we focused on 200 genes which are considered to be mechanistically relevant in the process of lipid droplet biogenesis in hepatocytes as recently published (Sahini and Borlak, 2014). Based on mechanistic considerations we identified 19 genes which displayed dose dependent responses while 10 genes showed time dependency. Importantly, the present study defined 9 genes (ANGPTL4, FABP7, FADS1, FGF21, GOT1, LDLR, GK, STAT3, and PKLR) as signature genes to predict DIS. Moreover, cross tabulation revealed 9 genes to be regulated ≥10 times amongst the various conditions and included genes linked to glucose metabolism, lipid transport and lipogenesis as well as signalling events. Additionally, a comparison between drugs causing phospholipidosis and/or steatosis revealed 26 genes to be regulated in common including 4 signature genes to predict DIS (PKLR, GK, FABP7 and FADS1). Furthermore, a comparison between in vivo single dose (3, 6, 9 and 24 h) and findings from rat hepatocyte studies (2 h, 8 h, 24 h) identified 10 genes which are regulated in common and contained 2 DIS signature genes (FABP7, FGF21). Altogether, our studies provide comprehensive information on mechanistically linked gene expression changes of a range of drugs causing steatosis and phospholipidosis and encourage the screening of DIS signature genes at the preclinical stage.     

mRNA Profiling Reveals Determinants of Trastuzumab Efficiency in HER2-Positive Breast Cancer

Intrinsic and acquired resistance to the monoclonal antibody drug trastuzumab is a major problem in the treatment of HER2-positive breast cancer. A deeper understanding of the underlying mechanisms could help to develop new agents. Our intention was to detect genes and single nucleotide polymorphisms (SNPs) affecting trastuzumab efficiency in cell culture. Three HER2-positive breast cancer cell lines with different resistance phenotypes were analyzed. We chose BT474 as model of trastuzumab sensitivity, HCC1954 as model of intrinsic resistance, and BTR50, derived from BT474, as model of acquired resistance. Based on RNA-Seq data, we performed differential expression analyses on these cell lines with and without trastuzumab treatment. Differentially expressed genes between the resistant cell lines and BT474 are expected to contribute to resistance. Differentially expressed genes between untreated and trastuzumab treated BT474 are expected to contribute to drug efficacy. To exclude false positives from the candidate gene set, we removed genes that were also differentially expressed between untreated and trastuzumab treated BTR50. We further searched for SNPs in the untreated cell lines which could contribute to trastuzumab resistance. The analysis resulted in 54 differentially expressed candidate genes that might be connected to trastuzumab efficiency. 90% of 40 selected candidates were validated by RT-qPCR. ALPP, CALCOCO1, CAV1, CYP1A2 and IGFBP3 were significantly higher expressed in the trastuzumab treated than in the untreated BT474 cell line. GDF15, IL8, LCN2, PTGS2 and 20 other genes were significantly higher expressed in HCC1954 than in BT474, while NCAM2, COLEC12, AFF3, TFF3, NRCAM, GREB1 and TFF1 were significantly lower expressed. Additionally, we inferred SNPs in HCC1954 for CAV1, PTGS2, IL8 and IGFBP3. The latter also had a variation in BTR50. 20% of the validated subset have already been mentioned in literature. For half of them we called and analyzed SNPs. These results contribute to a better understanding of trastuzumab action and resistance mechanisms.     

Genome-Wide Analysis Reveals Novel Genes Essential for Heme Homeostasis in Caenorhabditis elegans

Heme is a cofactor in proteins that function in almost all sub-cellular compartments and in many diverse biological processes. Heme is produced by a conserved biosynthetic pathway that is highly regulated to prevent the accumulation of heme—a cytotoxic, hydrophobic tetrapyrrole. Caenorhabditis elegans and related parasitic nematodes do not synthesize heme, but instead require environmental heme to grow and develop. Heme homeostasis in these auxotrophs is, therefore, regulated in accordance with available dietary heme. We have capitalized on this auxotrophy in C. elegans to study gene expression changes associated with precisely controlled dietary heme concentrations. RNA was isolated from cultures containing 4, 20, or 500 µM heme; derived cDNA probes were hybridized to Affymetrix C. elegans expression arrays. We identified 288 heme-responsive genes (hrgs) that were differentially expressed under these conditions. Of these genes, 42% had putative homologs in humans, while genomes of medically relevant heme auxotrophs revealed homologs for 12% in both Trypanosoma and Leishmania and 24% in parasitic nematodes. Depletion of each of the 288 hrgs by RNA–mediated interference (RNAi) in a transgenic heme-sensor worm strain identified six genes that regulated heme homeostasis. In addition, seven membrane-spanning transporters involved in heme uptake were identified by RNAi knockdown studies using a toxic heme analog. Comparison of genes that were positive in both of the RNAi screens resulted in the identification of three genes in common that were vital for organismal heme homeostasis in C. elegans. Collectively, our results provide a catalog of genes that are essential for metazoan heme homeostasis and demonstrate the power of C. elegans as a genetic animal model to dissect the regulatory circuits which mediate heme trafficking in both vertebrate hosts and their parasites, which depend on environmental heme for survival.