Now you see it: genome methylation makes a comeback in Drosophila

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates axon guidance genes in Drosophila

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases and catalyze the transition of 5mC to 5hmC in DNA. These enzymes have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila as Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by mapping one, Tet DNA-binding sites throughout the genome and two, the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC modifications are distributed along the entire transcript, while Tet DNA-binding sites are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are preferentially involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs primarily in cells of the nervous system.

High contiguity de novo genome assembly and DNA modification analyses for the fungus fly, Sciara coprophila, using single-molecule sequencing

BACKGROUND

The lower Dipteran fungus fly, Sciara coprophila, has many unique biological features that challenge the rule of genome DNA constancy. For example, Sciara undergoes paternal chromosome elimination and maternal X chromosome nondisjunction during spermatogenesis, paternal X elimination during embryogenesis, intrachromosomal DNA amplification of DNA puff loci during larval development, and germline-limited chromosome elimination from all somatic cells. Paternal chromosome elimination in Sciara was the first observation of imprinting, though the mechanism remains a mystery. Here, we present the first draft genome sequence for Sciara coprophila to take a large step forward in addressing these features.

RESULTS

We assembled the Sciara genome using PacBio, Nanopore, and Illumina sequencing. To find an optimal assembly using these datasets, we generated 44 short-read and 50 long-read assemblies. We ranked assemblies using 27 metrics assessing contiguity, gene content, and dataset concordance. The highest-ranking assemblies were scaffolded using BioNano optical maps. RNA-seq datasets from multiple life stages and both sexes facilitated genome annotation. A set of 66 metrics was used to select the first draft assembly for Sciara. Nearly half of the Sciara genome sequence was anchored into chromosomes, and all scaffolds were classified as X-linked or autosomal by coverage.

CONCLUSIONS

We determined that X-linked genes in Sciara males undergo dosage compensation. An entire bacterial genome from the Rickettsia genus, a group known to be endosymbionts in insects, was co-assembled with the Sciara genome, opening the possibility that Rickettsia may function in sex determination in Sciara. Finally, the signal level of the PacBio and Nanopore data support the presence of cytosine and adenine modifications in the Sciara genome, consistent with a possible role in imprinting.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

DNA methylation patterns respond to thermal stress in the viviparous cockroach Diploptera punctata

It is increasingly recognized that epigenetic mechanisms play a key role in acclimatization and adaptation to thermal stress in invertebrates. DNA methylation and its response to temperature variation has been poorly studied in insects. Here, we investigated DNA methylation and hydroxymethylation patterns in the viviparous cockroach Diploptera punctata at a global and gene specific level in response to variation in temperature. We specifically studied methylation percentage in the heat shock protein 70 (Hsp70), whose function is linked to thermal plasticity and resistance. We found high levels of DNA methylation in several tissues but only low levels of DNA hydroxymethylation in the brain. Hsp70 methylation patterns showed significant differences in response to temperature. We further found that global DNA methylation variation was considerably lower at 28°C compared to higher or lower temperatures, which may be indicative of the optimal temperature for this species. Our results demonstrate that DNA methylation could provide a mechanism for insects to dynamically respond to changing temperature conditions in their environment.

Can epigenetics translate environmental cues into phenotypes?

Living organisms are constantly exposed to wide ranges of environmental cues. They react to these cues by undergoing a battery of phenotypic responses, such as by altering their physiological and behavioral traits, in order to adapt and survive in the changed environments. The adaptive response of a species induced by environmental cues is typically thought to be associated with its genetic diversity such that higher genetic diversity provides increased adaptive potential. This originates from the general consensus that phenotypic traits have a genetic basis and are subject to Darwinian natural selection and Mendelian inheritance. There is no doubt about the validity of these principles, supported by the successful introgression of specific traits during (selective) breeding. However, a range of recent studies provided fascinating evidences suggesting that environmental effects experienced by an organism during its lifetime can have marked influences on its phenotype, and additionally the organism can pass on the acquired phenotypes to its subsequent generations through non-genetic mechanisms (also termed as epigenetic mechanism) - a notion that dates back to Lamarck and has been controversial ever since. In this review, we describe how the epigenetics has reshaped our long perception about the inheritance/development of phenotypes within organisms, contrasting with the classical gene-based view of inheritance. We particularly highlighted recent developments in our understanding of inheritance of parental environmental induced phenotypic traits in multicellular organisms under different environmental conditions, and discuss how modifications of the epigenome contribute to the determination of the adult phenotype of future generations.

Identification and caste-dependent expression patterns of DNA methylation associated genes in Bombus terrestris

DNA methylation has been proposed to play critical roles in caste fate and behavioral plasticity in bumblebees, however, there is little information on its regulatory mechanisms. Here, we identified six important genes mediating the modification of DNA methylation and determined their expression patterns in the bumblebee Bombus terrestris. There is a complete functional DNA methylation system, including four DNA methyltransferases (DNMT1a, DNMT1b, DNMT2, and DNMT3), a DNA demethylase (Ten-eleven translocation), and a methyl-CpG-binding domain protein in B. terrestris. Most of these genes were highly expressed in fat bodies and gonads but lowly expressed in antennae and brains of bumblebee adults. Besides, these genes exhibited caste-specific expression patterns in bumblebees, with higher transcription levels in queens than workers and drones. Whereas their expression levels showed no remarkable difference in queenright and queenless workers. These results suggested potential roles of DNA methylation-related genes in caste differentiation in bumblebees.

Genome-wide DNA methylomes from discrete developmental stages reveal the predominance of non-CpG methylation in Tribolium castaneum

Cytosine DNA methylation is a vital epigenetic regulator of eukaryotic development. Whether this epigenetic modification occurs in Tribolium castaneum has been controversial, its distribution pattern and functions have not been established. Here, using bisulphite sequencing (BS-Seq), we confirmed the existence of DNA methylation and described the methylation profiles of the four life stages of T. castaneum. In the T. castaneum genome, both symmetrical CpG and non-CpG methylcytosines were observed. Symmetrical CpG methylation, which was catalysed by DNMT1 and occupied a small part in T. castaneum methylome, was primarily enriched in gene bodies and was positively correlated with gene expression levels. Asymmetrical non-CpG methylation, which was predominant in the methylome, was strongly concentrated in intergenic regions and introns but absent from exons. Gene body methylation was negatively correlated with gene expression levels. The distribution pattern and functions of this type of methylation were similar only to the methylome of Drosophila melanogaster, which further supports the existence of a novel methyltransferase in the two species responsible for this type of methylation. This first life-cycle methylome of T. castaneum reveals a novel and unique methylation pattern, which will contribute to the further understanding of the variety and functions of DNA methylation in eukaryotes.

Methyl-CpG binding domain proteins inhibit interspecies courtship and promote aggression in Drosophila

Reproductive isolation and speciation are driven by the convergence of environmental and genetic variation. The integration of these variation sources is thought to occur through epigenetic marks including DNA methylation. Proteins containing a methyl-CpG-binding domain (MBD) bind methylated DNA and interpret epigenetic marks, providing a dynamic yet evolutionarily adapted cellular output. Here, we report the Drosophila MBD-containing proteins, dMBD-R2 and dMBD2/3, contribute to reproductive isolation and survival behavioral strategies. Drosophila melanogaster males with a reduction in dMBD-R2 specifically in octopamine (OA) neurons exhibit courtship toward divergent interspecies D. virilis and D. yakuba females and a decrease in conspecific mating success. Conspecific male-male courtship is increased between dMBD-R2-deficient males while aggression is reduced. These changes in adaptive behavior are separable as males with a hypermethylated OA neuronal genome exhibited a decrease in aggression without altering male-male courtship. These results suggest Drosophila MBD-containing proteins are required within the OA neural circuitry to inhibit interspecies and conspecific male-male courtship and indicate that the genetically hard-wired neural mechanisms enforcing behavioral reproductive isolation include the interpretation of the epigenome.

Genome methylation patterns across castes and generations in a parasitoid wasp

Environmental influences shape phenotypes within and across generations, often through DNA methylations that modify gene expression. Methylations were proposed to mediate caste and task allocation in some eusocial insects, but how an insect's environment affects DNA methylation in its offspring is yet unknown. We characterized parental effects on methylation profiles in the polyembryonic parasitoid wasp Copidosoma koehleri, as well as methylation patterns associated with its simple caste system. We used methylation‐sensitive amplified fragment length polymorphism (MS‐AFLP) to compare methylation patterns, among (1) reproductive and soldier larvae; and (2) offspring (larvae, pupae, and adults) of wasps that were reared at either high or low larval density and mated in the four possible combinations. Methylation frequencies were similar across castes, but the profiles of methylated fragments differed significantly. Parental rearing density did not affect methylation frequencies in the offspring at any developmental stage. Principal coordinate analysis indicated no significant differences in methylation profiles among the four crossbreeding groups and the three developmental stages. Nevertheless, a clustering analysis, performed on a subset of the fragments, revealed similar methylation patterns in larvae, pupae, and adults in two of the four parental crosses. Nine fragments were methylated at two cytosine sites in all larvae, and five others were methylated at two sites in all adults. Thus, DNA methylations correlate with within‐generation phenotypic plasticity due to caste. However, their association with developmental stage and with transgenerational epigenetic effects is not clearly supported.

Functional conservation of MBD proteins: MeCP2 and Drosophila MBD proteins alter sleep

Proteins containing a methyl-CpG-binding domain (MBD) bind 5mC and convert the methylation pattern information into appropriate functional cellular states. The correct readout of epigenetic marks is of particular importance in the nervous system where abnormal expression or compromised MBD protein function, can lead to disease and developmental disorders. Recent evidence indicates that the genome of Drosophila melanogaster is methylated and two MBD proteins, dMBD2/3 and dMBD-R2, are present. Are Drosophila MBD proteins required for neuronal function, and as MBD-containing proteins have diverged and evolved, does the MBD domain retain the molecular properties required for conserved cellular function across species? To address these questions, we expressed the human MBD-containing protein, hMeCP2, in distinct amine neurons and quantified functional changes in sleep circuitry output using a high throughput assay in Drosophila. hMeCP2 expression resulted in phase-specific sleep loss and sleep fragmentation with the hMeCP2-mediated sleep deficits requiring an intact MBD domain. Reducing endogenous dMBD2/3 and dMBD-R2 levels also generated sleep fragmentation, with an increase in sleep occurring upon dMBD-R2 reduction. To examine if hMeCP2 and dMBD-R2 are targeting common neuronal functions, we reduced dMBD-R2 levels in combination with hMeCP2 expression and observed a complete rescue of sleep deficits. Furthermore, chromosomal binding experiments indicate MBD-R2 and MeCP2 associate on shared genomic loci. Our results provide the first demonstration that Drosophila MBD-containing family members are required for neuronal function and suggest that the MBD domain retains considerable functional conservation at the whole organism level across species.

Transcription factors, chromatin proteins and the diversification of Hemiptera

Availability of complete genomes provides a means to explore the evolution of enormous developmental, morphological, and behavioral diversity among insects. Hemipterans in particular show great diversity of both morphology and life history within a single order. To better understand the role of transcription regulators in the diversification of hemipterans, using sequence profile searches and hidden Markov models we computationally analyzed transcription factors (TFs) and chromatin proteins (CPs) in the recently available Rhodnius prolixus genome along with 13 other insect and 4 non-insect arthropod genomes. We generated a comprehensive collection of TFs and CPs across arthropods including 303 distinct types of domains in TFs and 139 in CPs. This, along with the availability of two hemipteran genomes, R. prolixus and Acyrthosiphon pisum, helped us identify possible determinants for their dramatic morphological and behavioral divergence. We identified five domain families (i.e. Pipsqueak, SAZ/MADF, THAP, FLYWCH and BED finger) as having undergone differential patterns of lineage-specific expansion in hemipterans or within hemipterans relative to other insects. These expansions appear to be at least in part driven by transposons, with the DNA-binding domains of transposases having provided the raw material for emergence of new TFs. Our analysis suggests that while R. prolixus probably retains a state closer to the ancestral hemipteran, A. pisum represents a highly derived state, with the emergence of asexual reproduction potentially favoring genome duplication and transposon expansion. Both hemipterans are predicted to possess active DNA methylation systems. However, in the course of their divergence, aphids seem to have expanded the ancestral hemipteran DNA methylation along with a distinctive linkage to the histone methylation system, as suggested by expansion of SET domain methylases, including those fused to methylated CpG recognition domains. Thus, differential use of DNA methylation and histone methylation might have played a role in emergence of polyphenism and cyclic parthenogenesis from the ancestral hemipteran.

What Can Plasticity Contribute to Insect Responses to Climate Change?

Plastic responses figure prominently in discussions on insect adaptation to climate change. Here we review the different types of plastic responses and whether they contribute much to adaptation. Under climate change, plastic responses involving diapause are often critical for population persistence, but key diapause responses under dry and hot conditions remain poorly understood. Climate variability can impose large fitness costs on insects showing diapause and other life cycle responses, threatening population persistence. In response to stressful climatic conditions, insects also undergo ontogenetic changes including hardening and acclimation. Environmental conditions experienced across developmental stages or by prior generations can influence hardening and acclimation, although evidence for the latter remains weak. Costs and constraints influence patterns of plasticity across insect clades, but they are poorly understood within field contexts. Plastic responses and their evolution should be considered when predicting vulnerability to climate change-but meaningful empirical data lag behind theory.

Basic mechanisms of longevity: A case study of Drosophila pro-longevity genes

Drosophila is one of the most convenient model organisms in the genetics of aging and longevity. Unlike the nematodes, which allow for the detection of new pro-aging genes by knockout and RNAi-mediated knock-down, Drosophila also provides an opportunity to find new pro-longevity genes by driver-induced overexpression. Similar studies on other models are extremely rare. In this review, we focused on genes whose overexpression prolongs the life of fruit flies. The majority of longevity-associated genes regulates metabolism and stress resistance, and belongs to the IGF-1R, PI3K, PKB, AMPK and TOR metabolic regulation cluster and the FOXO, HDAC, p53 stress response cluster.

N6-methyladenine functions as a potential epigenetic mark in eukaryotes

N(6)-methyladenine (6mA) is one of the most abundant types of DNA methylation, and plays an important role in bacteria; however, its roles in higher eukaryotes, such as plants, insects, and mammals, have been considered less important. Recent studies highlight that 6mA does indeed occur, and that it plays an important role in eukaryotes, such as worm, fly, and green algae, and thus the regulation of 6mA has emerged as a novel epigenetic mechanism in higher eukaryotes. Despite this intriguing development, a number of important issues regarding its biological roles are yet to be addressed. In this review, we focus on the 5mC and 6mA modifications in terms of their production, distribution, and the erasure of 6mA in higher eukaryotes including mammals. We perform an analysis of the potential functions of 6mA, hence widening understanding of this new epigenetic mark in higher eukaryotes, and suggesting future studies in this field.

Transgenerational inheritance of metabolic disease

Metabolic disease encompasses several disorders including obesity, type 2 diabetes, and dyslipidemia. Recently, the incidence of metabolic disease has drastically increased, driven primarily by a worldwide obesity epidemic. Transgenerational inheritance remains controversial, but has been proposed to contribute to human metabolic disease risk based on a growing number of proof-of-principle studies in model organisms ranging from Caenorhabditis elegans to Mus musculus to Sus scrofa. Collectively, these studies demonstrate that heritable risk is epigenetically transmitted from parent to offspring over multiple generations in the absence of a continued exposure to the triggering stimuli. A diverse assortment of initial triggers can induce transgenerational inheritance including high-fat or high-sugar diets, low-protein diets, various toxins, and ancestral genetic variants. Although the mechanistic basis underlying the transgenerational inheritance of disease risk remains largely unknown, putative molecules mediating transmission include small RNAs, histone modifications, and DNA methylation. Due to the considerable impact of metabolic disease on human health, it is critical to better understand the role of transgenerational inheritance of metabolic disease risk to open new avenues for therapeutic intervention and improve upon the current methods for clinical diagnoses and treatment.

The genome of Diuraphis noxia, a global aphid pest of small grains

BACKGROUND

The Russian wheat aphid, Diuraphis noxia Kurdjumov, is one of the most important pests of small grains throughout the temperate regions of the world. This phytotoxic aphid causes severe systemic damage symptoms in wheat, barley, and other small grains as a direct result of the salivary proteins it injects into the plant while feeding.

RESULTS

We sequenced and de novo assembled the genome of D. noxia Biotype 2, the strain most virulent to resistance genes in wheat. The assembled genomic scaffolds span 393 MB, equivalent to 93% of its 421 MB genome, and contains 19,097 genes. D. noxia has the most AT-rich insect genome sequenced to date (70.9%), with a bimodal CpG(O/E) distribution and a complete set of methylation related genes. The D. noxia genome displays a widespread, extensive reduction in the number of genes per ortholog group, including defensive, detoxification, chemosensory, and sugar transporter groups in comparison to the Acyrthosiphon pisum genome, including a 65% reduction in chemoreceptor genes. Thirty of 34 known D. noxia salivary genes were found in this assembly. These genes exhibited less homology with those salivary genes commonly expressed in insect saliva, such as glucose dehydrogenase and trehalase, yet greater conservation among genes that are expressed in D. noxia saliva but not detected in the saliva of other insects. Genes involved in insecticide activity and endosymbiont-derived genes were also found, as well as genes involved in virus transmission, although D. noxia is not a viral vector.

CONCLUSIONS

This genome is the second sequenced aphid genome, and the first of a phytotoxic insect. D. noxia's reduced gene content of may reflect the influence of phytotoxic feeding in shaping the D. noxia genome, and in turn in broadening its host range. The presence of methylation-related genes, including cytosine methylation, is consistent with other parthenogenetic and polyphenic insects. The D. noxia genome will provide an important contrast to the A. pisum genome and advance functional and comparative genomics of insects and other organisms.

The relative importance of DNA methylation and Dnmt2-mediated epigenetic regulation on Wolbachia densities and cytoplasmic incompatibility

Wolbachia pipientis is a worldwide bacterial parasite of arthropods that infects germline cells and manipulates host reproduction to increase the ratio of infected females, the transmitting sex of the bacteria. The most common reproductive manipulation, cytoplasmic incompatibility (CI), is expressed as embryonic death in crosses between infected males and uninfected females. Specifically, Wolbachia modify developing sperm in the testes by unknown means to cause a post-fertilization disruption of the sperm chromatin that incapacitates the first mitosis of the embryo. As these Wolbachia-induced changes are stable, reversible, and affect the host cell cycle machinery including DNA replication and chromosome segregation, we hypothesized that the host methylation pathway is targeted for modulation during cytoplasmic incompatibility because it accounts for all of these traits. Here we show that infection of the testes is associated with a 55% increase of host DNA methylation in Drosophila melanogaster, but methylation of the paternal genome does not correlate with penetrance of CI. Overexpression and knock out of the Drosophila DNA methyltransferase Dnmt2 neither induces nor increases CI. Instead, overexpression decreases Wolbachia titers in host testes by approximately 17%, leading to a similar reduction in CI levels. Finally, strength of CI induced by several different strains of Wolbachia does not correlate with levels of DNA methylation in the host testes. We conclude that DNA methylation mediated by Drosophila’s only known methyltransferase is not required for the transgenerational sperm modification that causes CI.