G-OnRamp: Generating genome browsers to facilitate undergraduate-driven collaborative genome annotation.

Scientists are sequencing new genomes at an increasing rate with the goal of associating genome contents with phenotypic traits. After a new genome is sequenced and assembled, structural gene annotation is often the first step in analysis. Despite advances in computational gene prediction algorithms, most eukaryotic genomes still benefit from manual gene annotation. This requires access to good genome browsers to enable annotators to visualize and evaluate multiple lines of evidence (e.g., sequence similarity, RNA sequencing [RNA-Seq] results, gene predictions, repeats) and necessitates many volunteers to participate in the work. To address the technical barriers to creating genome browsers, the Genomics Education Partnership (GEP; https://gep.wustl.edu/) has partnered with the Galaxy Project (https://galaxyproject.org) to develop G-OnRamp (http://g-onramp.org), a web-based platform for creating UCSC Genome Browser Assembly Hubs and JBrowse genome browsers. G-OnRamp also converts a JBrowse instance into an Apollo instance for collaborative genome annotations in research and educational settings. The genome browsers produced can be transferred to the CyVerse Data Store for long-term access. G-OnRamp enables researchers to easily visualize their experimental results, educators to create Course-based Undergraduate Research Experiences (CUREs) centered on genome annotation, and students to participate in genomics research. In the process, students learn about genes/genomes and about how to utilize large datasets. Development of G-OnRamp was guided by extensive user feedback. Sixty-five researchers/educators from >40 institutions participated through in-person workshops, which produced >20 genome browsers now available for research and education. Genome browsers generated for four parasitoid wasp species have been used in a CURE engaging students at 15 colleges and universities. Our assessment results in the classroom demonstrate that the genome browsers produced by G-OnRamp are effective tools for engaging undergraduates in research and in enabling their contributions to the scientific literature in genomics. Expansion of such genomics research/education partnerships will be beneficial to researchers, faculty, and students alike.

The GEP: Crowd-Sourcing Big Data Analysis with Undergraduates.

The era of 'big data' is also the era of abundant data, creating new opportunities for student-scientist research partnerships. By coordinating undergraduate efforts, the Genomics Education Partnership produces high-quality annotated data sets and analyses that could not be generated otherwise, leading to scientific publications while providing many students with research experience.

G-OnRamp: a Galaxy-based platform for collaborative annotation of eukaryotic genomes.

SUMMARY: G-OnRamp provides a user-friendly, web-based platform for collaborative, end-to-end annotation of eukaryotic genomes using UCSC Assembly Hubs and JBrowse/Apollo genome browsers with evidence tracks derived from sequence alignments, ab initio gene predictors, RNA-Seq data and repeat finders. G-OnRamp can be used to visualize large genomics datasets and to perform collaborative genome annotation projects in both research and educational settings. AVAILABILITY AND IMPLEMENTATION: The virtual machine images and tutorials are available on the G-OnRamp web site (http://g-onramp.org/deployments). The source code is available under an Academic Free License version 3.0 through the goeckslab GitHub repository (https://github.com/goeckslab). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

HP1a: a structural chromosomal protein regulating transcription.

Heterochromatin protein 1 (HP1a in Drosophila) is a conserved eukaryotic chromosomal protein that is prominently associated with pericentric heterochromatin and mediates the concomitant gene silencing. Mechanistic studies implicate HP1 family proteins as 'hub proteins,' able to interact with a variety of chromosomal proteins through the chromo-shadow domain (CSD), as well as to recognize key histone modification sites [primarily histone H3 di/trimethyl Lys9 (H3K9me2/3)] through the chromodomain (CD). Consequently, HP1 has many important roles in chromatin architecture and impacts both gene expression and gene silencing, utilizing a variety of mechanisms. Clearly, HP1 function is altered by context, and potentially by post-translational modifications (PTMs). Here, we report on recent ideas as to how this versatile protein accomplishes its diverse functions.

Comprehensive analysis of the chromatin landscape in Drosophila melanogaster.

Chromatin is composed of DNA and a variety of modified histones and non-histone proteins, which have an impact on cell differentiation, gene regulation and other key cellular processes. Here we present a genome-wide chromatin landscape for Drosophila melanogaster based on eighteen histone modifications, summarized by nine prevalent combinatorial patterns. Integrative analysis with other data (non-histone chromatin proteins, DNase I hypersensitivity, GRO-Seq reads produced by engaged polymerase, short/long RNA products) reveals discrete characteristics of chromosomes, genes, regulatory elements and other functional domains. We find that active genes display distinct chromatin signatures that are correlated with disparate gene lengths, exon patterns, regulatory functions and genomic contexts. We also demonstrate a diversity of signatures among Polycomb targets that include a subset with paused polymerase. This systematic profiling and integrative analysis of chromatin signatures provides insights into how genomic elements are regulated, and will serve as a resource for future experimental investigations of genome structure and function.

Stimulating conversations between HP1a and histone demethylase dKDM4A.

An interaction between Drosophila HP1a and H3K36 histone demethylase dKDM4A reported in this issue of Molecular Cell (Lin et al., 2008) offers an explanation for association of HP1a with active genes, reinforcing the notion of dynamic interplay among epigenetic marks on histone tails during transcription.

Maternal depletion of Piwi, a component of the RNAi system, impacts heterochromatin formation in Drosophila.

A persistent question in epigenetics is how heterochromatin is targeted for assembly at specific domains, and how that chromatin state is faithfully transmitted. Stable heterochromatin is necessary to silence transposable elements (TEs) and maintain genome integrity. Both the RNAi system and heterochromatin components HP1 (Swi6) and H3K9me2/3 are required for initial establishment of heterochromatin structures in S. pombe. Here we utilize both loss of function alleles and the newly developed Drosophila melanogaster transgenic shRNA lines to deplete proteins of interest at specific development stages to dissect their roles in heterochromatin assembly in early zygotes and in maintenance of the silencing chromatin state during development. Using reporters subject to Position Effect Variegation (PEV), we find that depletion of key proteins in the early embryo can lead to loss of silencing assayed at adult stages. The piRNA component Piwi is required in the early embryo for reporter silencing in non-gonadal somatic cells, but knock-down during larval stages has no impact. This implies that Piwi is involved in targeting HP1a when heterochromatin is established at the late blastoderm stage and possibly also during embryogenesis, but that the silent chromatin state created is transmitted through cell division independent of the piRNA system. In contrast, heterochromatin structural protein HP1a is required for both initial heterochromatin assembly and the following mitotic inheritance. HP1a profiles in piwi mutant animals confirm that Piwi depletion leads to decreased HP1a levels in pericentric heterochromatin, particularly in TEs. The results suggest that the major role of the piRNA system in assembly of heterochromatin in non-gonadal somatic cells occurs in the early embryo during heterochromatin formation, and further demonstrate that failure of heterochromatin formation in the early embryo impacts the phenotype of the adult.

Nature and function of insulator protein binding sites in the Drosophila genome.

Chromatin insulator elements and associated proteins have been proposed to partition eukaryotic genomes into sets of independently regulated domains. Here we test this hypothesis by quantitative genome-wide analysis of insulator protein binding to Drosophila chromatin. We find distinct combinatorial binding of insulator proteins to different classes of sites and uncover a novel type of insulator element that binds CP190 but not any other known insulator proteins. Functional characterization of different classes of binding sites indicates that only a small fraction act as robust insulators in standard enhancer-blocking assays. We show that insulators restrict the spreading of the H3K27me3 mark but only at a small number of Polycomb target regions and only to prevent repressive histone methylation within adjacent genes that are already transcriptionally inactive. RNAi knockdown of insulator proteins in cultured cells does not lead to major alterations in genome expression. Taken together, these observations argue against the concept of a genome partitioned by specialized boundary elements and suggest that insulators are reserved for specific regulation of selected genes.

Ectopic assembly of heterochromatin in Drosophila melanogaster triggered by transposable elements.

A persistent question in biology is how cis-acting sequence elements influence trans-acting factors and the local chromatin environment to modulate gene expression. We reported previously that the DNA transposon 1360 can enhance silencing of a reporter in a heterochromatic domain of Drosophila melanogaster. We have now generated a collection of variegating phiC31 landing-pad insertion lines containing 1360 and a heat-shock protein 70 (hsp70)-driven white reporter to explore the mechanism of 1360-sensitive silencing. Many 1360-sensitive sites were identified, some in apparently euchromatic domains, although all are close to heterochromatic masses. One such site (line 1198; insertion near the base of chromosome arm 2L) has been investigated in detail. ChIP analysis shows 1360-dependent Heterochromatin Protein 1a (HP1a) accumulation at this otherwise euchromatic site. The phiC31 landing pad system allows different 1360 constructs to be swapped with the full-length element at the same genomic site to identify the sequences that mediate 1360-sensitive silencing. Short deletions over sites with homology to PIWI-interacting RNAs (piRNAs) are sufficient to compromise 1360-sensitive silencing. Similar results were obtained on replacing 1360 with Invader4 (a retrotransposon), suggesting that this phenomenon likely applies to a broader set of transposable elements. Our results suggest a model in which piRNA sequence elements behave as cis-acting targets for heterochromatin assembly, likely in the early embryo, where piRNA pathway components are abundant, with the heterochromatic state subsequently propagated by chromatin modifiers present in somatic tissue.

Sequence-specific targeting of dosage compensation in Drosophila favors an active chromatin context.

The Drosophila MSL complex mediates dosage compensation by increasing transcription of the single X chromosome in males approximately two-fold. This is accomplished through recognition of the X chromosome and subsequent acetylation of histone H4K16 on X-linked genes. Initial binding to the X is thought to occur at "entry sites" that contain a consensus sequence motif ("MSL recognition element" or MRE). However, this motif is only ∼2 fold enriched on X, and only a fraction of the motifs on X are initially targeted. Here we ask whether chromatin context could distinguish between utilized and non-utilized copies of the motif, by comparing their relative enrichment for histone modifications and chromosomal proteins mapped in the modENCODE project. Through a comparative analysis of the chromatin features in male S2 cells (which contain MSL complex) and female Kc cells (which lack the complex), we find that the presence of active chromatin modifications, together with an elevated local GC content in the surrounding sequences, has strong predictive value for functional MSL entry sites, independent of MSL binding. We tested these sites for function in Kc cells by RNAi knockdown of Sxl, resulting in induction of MSL complex. We show that ectopic MSL expression in Kc cells leads to H4K16 acetylation around these sites and a relative increase in X chromosome transcription. Collectively, our results support a model in which a pre-existing active chromatin environment, coincident with H3K36me3, contributes to MSL entry site selection. The consequences of MSL targeting of the male X chromosome include increase in nucleosome lability, enrichment for H4K16 acetylation and JIL-1 kinase, and depletion of linker histone H1 on active X-linked genes. Our analysis can serve as a model for identifying chromatin and local sequence features that may contribute to selection of functional protein binding sites in the genome.

Enrichment of HP1a on Drosophila chromosome 4 genes creates an alternate chromatin structure critical for regulation in this heterochromatic domain.

Chromatin environments differ greatly within a eukaryotic genome, depending on expression state, chromosomal location, and nuclear position. In genomic regions characterized by high repeat content and high gene density, chromatin structure must silence transposable elements but permit expression of embedded genes. We have investigated one such region, chromosome 4 of Drosophila melanogaster. Using chromatin-immunoprecipitation followed by microarray (ChIP-chip) analysis, we examined enrichment patterns of 20 histone modifications and 25 chromosomal proteins in S2 and BG3 cells, as well as the changes in several marks resulting from mutations in key proteins. Active genes on chromosome 4 are distinct from those in euchromatin or pericentric heterochromatin: while there is a depletion of silencing marks at the transcription start sites (TSSs), HP1a and H3K9me3, but not H3K9me2, are enriched strongly over gene bodies. Intriguingly, genes on chromosome 4 are less frequently associated with paused polymerase. However, when the chromatin is altered by depleting HP1a or POF, the RNA pol II enrichment patterns of many chromosome 4 genes shift, showing a significant decrease over gene bodies but not at TSSs, accompanied by lower expression of those genes. Chromosome 4 genes have a low incidence of TRL/GAGA factor binding sites and a low T(m) downstream of the TSS, characteristics that could contribute to a low incidence of RNA polymerase pausing. Our data also indicate that EGG and POF jointly regulate H3K9 methylation and promote HP1a binding over gene bodies, while HP1a targeting and H3K9 methylation are maintained at the repeats by an independent mechanism. The HP1a-enriched, POF-associated chromatin structure over the gene bodies may represent one type of adaptation for genes embedded in repetitive DNA.

Heterochromatin Protein 1a (HP1a) partner specificity is determined by critical amino acids in the chromo shadow domain and C-terminal extension.

Drosophila melanogaster Heterochromatin Protein 1a (HP1a) is an essential protein critical for heterochromatin assembly and regulation. Its chromo shadow domain (CSD) homodimerizes, a requirement for binding protein partners that contain a PXVXL motif. How does HP1a select among its many different PXVXL-containing partners? HP1a binds tightly to Heterochromatin Protein 2 (HP2), but weakly to PIWI. We investigated differences in homodimerization and the impact of the C-terminal extension (CTE) by contrasting HP1a to its paralogue, HP1b. HP1a and HP1b differ in the dimerization interface, with HP1a having an Arg at position 188 rather than Glu. We find that while this substitution reduces the dimerization constant, it does not impact the binding surface as demonstrated by unchanged partner binding affinities. However, the CTE (only 4 residues in HP1a as compared with 87 residues in HP1b) is critical; the charged residues in HP1a are necessary for tight peptide binding. Examining a panel of amino acid substitutions in the HP1a CSD, we find that Leu-165 in HP1a interacts with HP2 but not PIWI, supporting the conclusion that different sites in the binding surface provide discrimination for partner selection. Partner sequence is also critical for affinity, as the remaining difference in binding between HP2 and PIWI polypeptides is eliminated by swapping the PXVXL motifs between the two. Taken together, these studies indicate that the binding surface of the HP1a CSD plus its short CTE provide the needed discrimination among HP1a's partners, and that the CTE is important for differentiating the interactions of the Drosophila HP1 paralogs.

Plasticity in patterns of histone modifications and chromosomal proteins in Drosophila heterochromatin.

Eukaryotic genomes are packaged in two basic forms, euchromatin and heterochromatin. We have examined the composition and organization of Drosophila melanogaster heterochromatin in different cell types using ChIP-array analysis of histone modifications and chromosomal proteins. As anticipated, the pericentric heterochromatin and chromosome 4 are on average enriched for the "silencing" marks H3K9me2, H3K9me3, HP1a, and SU(VAR)3-9, and are generally depleted for marks associated with active transcription. The locations of the euchromatin-heterochromatin borders identified by these marks are similar in animal tissues and most cell lines, although the amount of heterochromatin is variable in some cell lines. Combinatorial analysis of chromatin patterns reveals distinct profiles for euchromatin, pericentric heterochromatin, and the 4th chromosome. Both silent and active protein-coding genes in heterochromatin display complex patterns of chromosomal proteins and histone modifications; a majority of the active genes exhibit both "activation" marks (e.g., H3K4me3 and H3K36me3) and "silencing" marks (e.g., H3K9me2 and HP1a). The hallmark of active genes in heterochromatic domains appears to be a loss of H3K9 methylation at the transcription start site. We also observe complex epigenomic profiles of intergenic regions, repeated transposable element (TE) sequences, and genes in the heterochromatic extensions. An unexpectedly large fraction of sequences in the euchromatic chromosome arms exhibits a heterochromatic chromatin signature, which differs in size, position, and impact on gene expression among cell types. We conclude that patterns of heterochromatin/euchromatin packaging show greater complexity and plasticity than anticipated. This comprehensive analysis provides a foundation for future studies of gene activity and chromosomal functions that are influenced by or dependent upon heterochromatin.

Drosophila Piwi functions downstream of piRNA production mediating a chromatin-based transposon silencing mechanism in female germ line.

Transposon control is a critical process during reproduction. The PIWI family proteins can play a key role, using a piRNA-mediated slicing mechanism to suppress transposon activity posttranscriptionally. In Drosophila melanogaster, Piwi is predominantly localized in the nucleus and has been implicated in heterochromatin formation. Here, we use female germ-line-specific depletion to study Piwi function. This depletion of Piwi leads to infertility and to axis specification defects in the developing egg chambers; correspondingly, widespread loss of transposon silencing is observed. Germ-line Piwi does not appear to be required for piRNA production. Instead, Piwi requires Aubergine (and presumably secondary piRNA) for proper localization. A subset of transposons that show significant overexpression in germ-line Piwi-depleted ovaries exhibit a corresponding loss of HP1a and H3K9me2. Germ-line HP1a depletion also leads to a loss of transposon silencing, demonstrating the functional requirement for HP1a enrichment at these loci. Considering our results and those of others, we infer that germ-line Piwi functions downstream of piRNA production to promote silencing of some transposons via recruitment of HP1a. Thus, in addition to its better-known function in posttranscriptional silencing, piRNA also appears to function in a targeting mechanism for heterochromatin formation mediated by Piwi.

Transcription and RNA interference in the formation of heterochromatin.

Transcription in heterochromatin seems to be an oxymoron--surely the 'silenced' form of chromatin should not be transcribed. But there have been frequent reports of low-level transcription in heterochromatic regions, and several hundred genes are found in these regions in Drosophila. Most strikingly, recent investigations implicate RNA interference mechanisms in targeting and maintaining heterochromatin, and these mechanisms are inherently dependent on transcription. Silencing of chromatin might involve trans-acting sources of the crucial small RNAs that carry out RNA interference, but in some cases, transcription of the region to be silenced seems to be required--an apparent contradiction.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai, Drosophila takahashii, Drosophila bipectinata, and Drosophila ananassae.

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma. However, the assembled portion of the Muller F element (orthologous to the fourth chromosome in Drosophila melanogaster) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to >18 Mb. Here, we present chromosome-level long-read genome assemblies for 4 Drosophila species with expanded F elements ranging in size from 2.3 to 20.5 Mb. Each Muller element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai , Drosophila takahashii , Drosophila bipectinata , and Drosophila ananassae.

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma . However, the assembled portion of the Muller F Element (orthologous to the fourth chromosome in Drosophila melanogaster ) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to > 18 Mb. Here, we present chromosome-level long read genome assemblies for four Drosophila species with expanded F Elements ranging in size from 2.3 Mb to 20.5 Mb. Each Muller Element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

The HP1a disordered C terminus and chromo shadow domain cooperate to select target peptide partners.

Drosophila melanogaster heterochromatin protein 1a (HP1a) is essential for compacted heterochromatin structure and the associated gene silencing. Its chromo shadow domain (CSD) is well known for binding to peptides that contain a PXVXL motif. Heterochromatin protein 2 (HP2) is a non-histone chromosomal protein that associates with HP1a in the pericentric heterochromatin, telomeres, and the fourth chromosome. Using NMR spectroscopy, fluorescence polarization, and site-directed mutagenesis, we identified an LCVKI motif in HP2 that binds to the HP1a CSD. The binding affinity of the HP2 fragment is approximately two orders of magnitude higher than that of peptides from PIWI (with a PRVKV motif), AF10 (with a PLVVL motif), or CG15356 (with LYPLL and LSIVA motifs). To delineate differential interactions of the HP1a CSD, we characterized its structure, backbone dynamics, and dimerization constant. We found that the dimerization constant is bracketed by the affinities of HP2 and PIWI, which dock to the same HP1a homodimer surface. This suggests that HP2, but not PIWI, interaction can drive the homodimerization of HP1a. Interestingly, the integrity of the disordered C-terminal extension (CTE) of HP1a is essential for discriminatory binding, whereas swapping the PXVXL motifs does not confer specificity. Serine phosphorylation at the peptide binding surface of the CSD is thought to regulate heterochromatin assembly. Glutamic acid substitution at these sites destabilizes HP1a dimers, but improves the interaction with both binding partners. Our studies underscore the importance of CSD dimerization and cooperation with the CTE in forming distinct complexes of HP1a.

Small RNA-directed heterochromatin formation in the context of development: what flies might learn from fission yeast.

A link between the RNAi system and heterochromatin formation has been established in several model organisms including Schizosaccharomyces pombe and Arabidopsis thaliana. However, the data to support a role for small RNAs and the associated machinery in transcriptional gene silencing in animal systems is more tenuous. Using the S. pombe system as a model, we analyze the role of small RNA pathway components and associated small RNAs in regulating transposable elements and potentially directing heterochromatin formation at these elements in Drosophila melanogaster.

Multiple SET methyltransferases are required to maintain normal heterochromatin domains in the genome of Drosophila melanogaster.

Methylation of histone H3 lysine 9 (H3K9) is a key feature of silent chromatin and plays an important role in stabilizing the interaction of heterochromatin protein 1 (HP1) with chromatin. Genomes of metazoans such as the fruit fly Drosophila melanogaster generally encode three types of H3K9-specific SET domain methyltransferases that contribute to chromatin homeostasis during the life cycle of the organism. SU(VAR)3-9, dG9a, and dSETDB1 all function in the generation of wild-type H3K9 methylation levels in the Drosophila genome. Two of these enzymes, dSETDB1 and SU(VAR)3-9, govern heterochromatin formation in distinct but overlapping patterns across the genome. H3K9 methylation in the small, heterochromatic fourth chromosome of D. melanogaster is governed mainly by dSETDB1, whereas dSETDB1 and SU(VAR)3-9 function in concert to methylate H3K9 in the pericentric heterochromatin of all chromosomes, with dG9a having little impact in these domains, as shown by monitoring position effect variegation. To understand how these distinct heterochromatin compartments may be differentiated, we examined the developmental timing of dSETDB1 function using a knockdown strategy. dSETDB1 acts to maintain heterochromatin during metamorphosis, at a later stage in development than the reported action of SU(VAR)3-9. Surprisingly, depletion of both of these enzymes has less deleterious effect than depletion of one. These results imply that dSETDB1 acts as a heterochromatin maintenance factor that may be required for the persistence of earlier developmental events normally governed by SU(VAR)3-9. In addition, the genetic interactions between dSETDB1 and Su(var)3-9 mutations emphasize the importance of maintaining the activities of these histone methyltransferases in balance for normal genome function.