The chromatin landscape of Drosophila: comparisons between species, sexes, and chromosomes

Chromatin Landscape Is Associated With Sex-Biased Expression and Drosophila-Like Dosage Compensation of the Z Chromosome in Artemia franciscana

The males and females of the brine shrimp Artemia franciscana are highly dimorphic, and this dimorphism is associated with substantial sex-biased gene expression in heads and gonads. How these sex-specific patterns of expression are regulated at the molecular level is unknown. A. franciscana also has differentiated ZW sex chromosomes, with complete dosage compensation, but the molecular mechanism through which compensation is achieved is unknown. Here, we conducted CUT&TAG assays targeting 7 post-translational histone modifications (H3K27me3, H3K9me2, H3K9me3, H3K36me3, H3K27ac, H3K4me3, and H4K16ac) in heads and gonads of A. franciscana, allowing us to divide the genome into 12 chromatin states. We further defined functional chromatin signatures for all genes, which were correlated with transcript level abundances. Differences in the occupancy of the profiled epigenetic marks between sexes were associated with differential gene expression between males and females. Finally, we found a significant enrichment of the permissive H4K16ac histone mark in the Z-specific region in both tissues of females but not males, supporting the role of this histone mark in mediating dosage compensation of the Z chromosome.

Interactions between U and V sex chromosomes during the life cycle of Ectocarpus

In many animals and flowering plants, sex determination occurs in the diploid phase of the life cycle with XX/XY or ZW/ZZ sex chromosomes. However, in early diverging plants and most macroalgae, sex is determined by female (U) or male (V) sex chromosomes in a haploid phase called the gametophyte. Once the U and V chromosomes unite at fertilization to produce a diploid sporophyte, sex determination no longer occurs, raising key questions about the fate of the U and V sex chromosomes in the sporophyte phase. Here, we investigate genetic and molecular interactions of the UV sex chromosomes in both the haploid and diploid phases of the brown alga Ectocarpus. We reveal extensive developmental regulation of sex chromosome genes across its life cycle and implicate the TALE-HD transcription factor OUROBOROS in suppressing sex determination in the diploid phase. Small RNAs may also play a role in the repression of a female sex-linked gene, and transition to the diploid sporophyte coincides with major reconfiguration of histone H3K79me2, suggesting a more intricate role for this histone mark in Ectocarpus development than previously appreciated.

Dynamic Evolution of Repetitive Elements and Chromatin States in Apis mellifera Subspecies

In this study, we elucidate the contribution of repetitive DNA sequences to the establishment of social structures in honeybees (Apis mellifera). Despite recent advancements in understanding the molecular mechanisms underlying the formation of honeybee castes, primarily associated with Notch signaling, the comprehensive identification of specific genomic cis-regulatory sequences remains elusive. Our objective is to characterize the repetitive landscape within the genomes of two honeybee subspecies, namely A. m. mellifera and A. m. ligustica. An observed recent burst of repeats in A. m. mellifera highlights a notable distinction between the two subspecies. After that, we transitioned to identifying differentially expressed DNA elements that may function as cis-regulatory elements. Nevertheless, the expression of these sequences showed minimal disparity in the transcriptome during caste differentiation, a pivotal process in honeybee eusocial organization. Despite this, chromatin segmentation, facilitated by ATAC-seq, ChIP-seq, and RNA-seq data, revealed a distinct chromatin state associated with repeats. Lastly, an analysis of sequence divergence among elements indicates successive changes in repeat states, correlating with their respective time of origin. Collectively, these findings propose a potential role of repeats in acquiring novel regulatory functions.

Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans

In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.

Sex-biased expression is associated with chromatin state in D. melanogaster and D. simulans

We propose a new model for the association of chromatin state and sex-bias in expression. We hypothesize enrichment of open chromatin in the sex where we see expression bias (OS) and closed chromatin in the opposite sex (CO). In this study of D. melanogaster and D. simulans head tissue, sex-bias in expression is associated with H3K4me3 (open mark) in males for male-biased genes and in females for female-biased genes in both species. Sex-bias in expression is also largely conserved in direction and magnitude between the two species on the X and autosomes. In male-biased orthologs, the sex-bias ratio is more divergent between species if both species have H3K27me2me3 marks in females compared to when either or neither species has H3K27me2me3 in females. H3K27me2me3 marks in females are associated with male-bias in expression on the autosomes in both species, but on the X only in D. melanogaster . In female-biased orthologs the relationship between the species for the sex-bias ratio is similar regardless of the H3K27me2me3 marks in males. Female-biased orthologs are more similar in the ratio of sex-bias than male-biased orthologs and there is an excess of male-bias in expression in orthologs that gain/lose sex-bias. There is an excess of male-bias in sex-limited expression in both species suggesting excess male-bias is due to rapid evolution between the species. The X chromosome has an enrichment in male-limited H3K4me3 in both species and an enrichment of sex-bias in expression compared to the autosomes.

Mode and Tempo of 3D Genome Evolution in Drosophila

Topologically associating domains (TADs) are thought to play an important role in preventing gene misexpression by spatially constraining enhancer-promoter contacts. The deleterious nature of gene misexpression implies that TADs should, therefore, be conserved among related species. Several early studies comparing chromosome conformation between species reported high levels of TAD conservation; however, more recent studies have questioned these results. Furthermore, recent work suggests that TAD reorganization is not associated with extensive changes in gene expression. Here, we investigate the evolutionary conservation of TADs among 11 species of Drosophila. We use Hi-C data to identify TADs in each species and employ a comparative phylogenetic approach to derive empirical estimates of the rate of TAD evolution. Surprisingly, we find that TADs evolve rapidly. However, we also find that the rate of evolution depends on the chromatin state of the TAD, with TADs enriched for developmentally regulated chromatin evolving significantly slower than TADs enriched for broadly expressed, active chromatin. We also find that, after controlling for differences in chromatin state, highly conserved TADs do not exhibit higher levels of gene expression constraint. These results suggest that, in general, most TADs evolve rapidly and their divergence is not associated with widespread changes in gene expression. However, higher levels of evolutionary conservation and gene expression constraints in TADs enriched for developmentally regulated chromatin suggest that these TAD subtypes may be more important for regulating gene expression, likely due to the larger number of long-distance enhancer-promoter contacts associated with developmental genes.

Transposable element accumulation drives size differences among polymorphic Y Chromosomes in Drosophila

Y Chromosomes of many species are gene poor and show low levels of nucleotide variation, yet often display high amounts of structural diversity. Dobzhansky cataloged several morphologically distinct Y Chromosomes in Drosophila pseudoobscura that differ in size and shape, but the molecular causes of their dramatic size differences are unclear. Here we use cytogenetics and long-read sequencing to study the sequence content of polymorphic Y Chromosomes in D. pseudoobscura We show that Y Chromosomes differ almost 2-fold in size, ranging from 30 to 60 Mb. Most of this size difference is caused by a handful of active transposable elements (TEs) that have recently expanded on the largest Y Chromosome, with different elements being responsible for Y expansion on differently sized D. pseudoobscura Y's. We show that Y Chromosomes differ in their heterochromatin enrichment, expression of Y-enriched TEs, and also influence expression of dozens of autosomal and X-linked genes. The same helitron element that showed the most drastic amplification on the largest Y in D. pseudoobscura independently amplified on a polymorphic large Y Chromosome in D. affinis, suggesting that some TEs are inherently more prone to become deregulated on Y Chromosomes.

Chromatin landscape associated with sexual differentiation in a UV sex determination system

In many eukaryotes, such as dioicous mosses and many algae, sex is determined by UV sex chromosomes and is expressed during the haploid phase of the life cycle. In these species, the male and female developmental programs are initiated by the presence of the U- or V-specific regions of the sex chromosomes but, as in XY and ZW systems, sexual differentiation is largely driven by autosomal sex-biased gene expression. The mechanisms underlying the regulation of sex-biased expression of genes during sexual differentiation remain elusive. Here, we investigated the extent and nature of epigenomic changes associated with UV sexual differentiation in the brown alga Ectocarpus, a model UV system. Six histone modifications were quantified in near-isogenic lines, leading to the identification of 16 chromatin signatures across the genome. Chromatin signatures correlated with levels of gene expression and histone PTMs changes in males versus females occurred preferentially at genes involved in sex-specific pathways. Despite the absence of chromosome scale dosage compensation and the fact that UV sex chromosomes recombine across most of their length, the chromatin landscape of these chromosomes was remarkably different to that of autosomes. Hotspots of evolutionary young genes in the pseudoautosomal regions appear to drive the exceptional chromatin features of UV sex chromosomes.

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

A transposon expression burst accompanies the activation of Y-chromosome fertility genes during Drosophila spermatogenesis

Transposable elements (TEs) must replicate in germline cells to pass novel insertions to offspring. In Drosophila melanogaster ovaries, TEs can exploit specific developmental windows of opportunity to evade host silencing and increase their copy numbers. However, TE activity and host silencing in the distinct cell types of Drosophila testis are not well understood. Here, we reanalyze publicly available single-cell RNA-seq datasets to quantify TE expression in the distinct cell types of the Drosophila testis. We develop a method for identification of TE and host gene expression modules and find that a distinct population of early spermatocytes expresses a large number of TEs at much higher levels than other germline and somatic components of the testes. This burst of TE expression coincides with the activation of Y chromosome fertility factors and spermatocyte-specific transcriptional regulators, as well as downregulation of many components of the piRNA pathway. The TEs expressed by this cell population are specifically enriched on the Y chromosome and depleted on the X chromosome, relative to other active TEs. These data suggest that some TEs may achieve high insertional activity in males by exploiting a window of opportunity for mobilization created by the activation of spermatocyte-specific and Y chromosome-specific transcriptional programs.

Analysis of cell-type-specific chromatin modifications and gene expression in Drosophila neurons that direct reproductive behavior

Examining the role of chromatin modifications and gene expression in neurons is critical for understanding how the potential for behaviors are established and maintained. We investigate this question by examining Drosophila melanogaster fru P1 neurons that underlie reproductive behaviors in both sexes. We developed a method to purify cell-type-specific chromatin (Chromatag), using a tagged histone H2B variant that is expressed using the versatile Gal4/UAS gene expression system. Here, we use Chromatag to evaluate five chromatin modifications, at three life stages in both sexes. We find substantial changes in chromatin modification profiles across development and fewer differences between males and females. Additionally, we find chromatin modifications that persist in different sets of genes from pupal to adult stages, which may point to genes important for cell fate determination in fru P1 neurons. We generated cell-type-specific RNA-seq data sets, using translating ribosome affinity purification (TRAP). We identify actively translated genes in fru P1 neurons, revealing novel stage- and sex-differences in gene expression. We also find chromatin modification enrichment patterns that are associated with gene expression. Next, we use the chromatin modification data to identify cell-type-specific super-enhancer-containing genes. We show that genes with super-enhancers in fru P1 neurons differ across development and between the sexes. We validated that a set of genes are expressed in fru P1 neurons, which were chosen based on having a super-enhancer and TRAP-enriched expression in fru P1 neurons.

Epigenetics drive the evolution of sex chromosomes in animals and plants

We review how epigenetics affect sex chromosome evolution in animals and plants. In a few species, sex is determined epigenetically through the action of Y-encoded small RNAs. Epigenetics is also responsible for changing the sex of individuals through time, even in species that carry sex chromosomes, and could favour species adaptation through breeding system plasticity. The Y chromosome accumulates repeats that become epigenetically silenced which leads to an epigenetic conflict with the expression of Y genes and could accelerate Y degeneration. Y heterochromatin can be lost through ageing, which activates transposable elements and lowers male longevity. Y chromosome degeneration has led to the evolution of meiotic sex chromosome inactivation in eutherians (placentals) and marsupials, and dosage compensation mechanisms in animals and plants. X-inactivation convergently evolved in eutherians and marsupials via two independently evolved non-coding RNAs. In Drosophila, male X upregulation by the male specific lethal (MSL) complex can spread to neo-X chromosomes through the transposition of transposable elements that carry an MSL-binding motif. We discuss similarities and possible differences between plants and animals and suggest future directions for this dynamic field of research. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

The Drosophila Y Chromosome Affects Heterochromatin Integrity Genome-Wide

The Drosophila Y chromosome is gene poor and mainly consists of silenced, repetitive DNA. Nonetheless, the Y influences expression of hundreds of genes genome-wide, possibly by sequestering key components of the heterochromatin machinery away from other positions in the genome. To test the influence of the Y chromosome on the genome-wide chromatin landscape, we assayed the genomic distribution of histone modifications associated with gene activation (H3K4me3) or heterochromatin (H3K9me2 and H3K9me3) in fruit flies with varying sex chromosome complements (X0, XY, and XYY males; XX and XXY females). Consistent with the general deficiency of active chromatin modifications on the Y, we find that Y gene dose has little influence on the genomic distribution of H3K4me3. In contrast, both the presence and the number of Y chromosomes strongly influence genome-wide enrichment patterns of repressive chromatin modifications. Highly repetitive regions such as the pericentromeres, the dot, and the Y chromosome (if present) are enriched for heterochromatic modifications in wildtype males and females, and even more strongly in X0 flies. In contrast, the additional Y chromosome in XYY males and XXY females diminishes the heterochromatic signal in these normally silenced, repeat-rich regions, which is accompanied by an increase in expression of Y-linked repeats. We find hundreds of genes that are expressed differentially between individuals with aberrant sex chromosome karyotypes, many of which also show sex-biased expression in wildtype Drosophila. Thus, Y chromosomes influence heterochromatin integrity genome-wide, and differences in the chromatin landscape of males and females may also contribute to sex-biased gene expression and sexual dimorphisms.

Genome-wide analysis of H3K4me3 and H3K27me3 modifications due to Lr28 for leaf rust resistance in bread wheat (Triticum aestivum)

Present study revealed a complex relationship among histone H3 methylation (examined using H3K4/K27me3 marks), cytosine DNA methylation and differential gene expression during Lr28 mediated leaf rust resistance in wheat. During the present study, genome-wide histone modifications were examined in a pair of near isogenic lines (NILs) (with and without Lr28 in the background of cv. HD2329). The two histone marks used included H3K4me3 (an activation mark) and H3K27me3 (a repression mark). The results were compared with levels of expression (using RNA-seq) and DNA methylation (MeDIP) data obtained using the same pair of NILs. Some of the salient features of the present study include the following: (i) large scale differential binding sites (DBS) were available for only H3K4me3 in the susceptible cultivar, but for both H3K4me3 and H3K27me3 in its resistant NIL; (ii) DBSs for H3K27me3 mark were more abundant (> 80%) in intergenic regions, whereas DBSs for H3K4me3 were distributed in all genomic regions including exons, introns, intergenic, TTS (transcription termination sites) and promoters; (iii) fourteen (14) genes associated with DBSs showed co-localization for both the marks; (iv) only a small fraction (7% for H3K4me3 and 12% for H3K27me3) of genes associated with DBSs matched with the levels of gene expression inferred from RNA-seq data; (v) validation studies using qRT-PCR were conducted on 26 selected representative genes; results for only 11 genes could be validated. The proteins encoded by important genes involved in promoting infection included domains generally carried by R gene proteins such as Mlo like protein, protein kinases and purple acid phosphatase. Similarly, proteins encoded by genes involved in resistance included those carrying domains for lectin kinase, R gene, aspartyl protease, etc. Overall, the results suggest a very complex network of downstream genes that are expressed during compatible and incompatible interactions; some of the genes identified during the present study may be used in future validation studies involving RNAi/overexpression approaches.

Sex and the TEs: transposable elements in sexual development and function in animals

Transposable elements are endogenous DNA sequences able to integrate into and multiply within genomes. They constitute a major source of genetic innovations, as they can not only rearrange genomes but also spread ready-to-use regulatory sequences able to modify host gene expression, and even can give birth to new host genes. As their evolutionary success depends on their vertical transmission, transposable elements are intrinsically linked to reproduction. In organisms with sexual reproduction, this implies that transposable elements have to manifest their transpositional activity in germ cells or their progenitors. The control of sexual development and function can be very versatile, and several studies have demonstrated the implication of transposable elements in the evolution of sex. In this review, we report the functional and evolutionary relationships between transposable elements and sexual reproduction in animals. In particular, we highlight how transposable elements can influence expression of sexual development genes, and how, reciprocally, they are tightly controlled in gonads. We also review how transposable elements contribute to the organization, expression and evolution of sexual development genes and sex chromosomes. This underscores the intricate co-evolution between host functions and transposable elements, which regularly shift from a parasitic to a domesticated status useful to the host.

Small Segmental Duplications in Drosophila-High Rate of Emergence and Elimination

Segmental duplications are an important class of mutations. Because a large proportion of segmental duplications may often be strongly deleterious, high frequency or fixed segmental duplications may represent only a tiny fraction of the mutational input. To understand the emergence and elimination of segmental duplications, we survey polymorphic duplications, including tandem and interspersed duplications, in natural populations of Drosophila by haploid embryo genomes. As haploid embryos are not expected to be heterozygous, the genome, sites of heterozygosity (referred to as pseudoheterozygous sites [PHS]), may likely represent recent duplications that have acquired new mutations. Among the 29 genomes of Drosophila melanogaster, we identify 2,282 polymorphic PHS duplications (linked PHS regions) in total or 154 PHS duplications per genome. Most PHS duplications are small (83.4% < 500 bp), Drosophila melanogaster lineage specific, and strain specific (72.6% singletons). The excess of the observed singleton PHS duplications deviates significantly from the neutral expectation, suggesting that most PHS duplications are strongly deleterious. In addition, these small segmental duplications are not evenly distributed in genomic regions and less common in noncoding functional element regions. The underrepresentation in RNA polymerase II binding sites and regions with active histone modifications is correlated with ages of duplications. In conclusion, small segmental duplications occur frequently in Drosophila but rapidly eliminated by natural selection.

Systematic Analysis of SIN3 Histone Modifying Complex Components During Development

Establishment and maintenance of histone acetylation levels are critical for metazoan development and viability. Disruption of the balance between acetylation and deacetylation by treatment with chemical histone deacetylase (HDAC) inhibitors results in loss of cell proliferation, differentiation and/or apoptosis. Histone deacetylation by the SIN3 complex is essential in Drosophila and mice, as loss of the scaffolding factor SIN3 or the associated HDAC results in lethality. The objective of this study is to elucidate contributions of SIN3 complex components to these essential processes. We used the Drosophila model organism to carry out a systematic functional analysis of the SIN3 complex. We find that SIN3 associated proteins are essential for viability and cell proliferation during development. Additionally, tissue specific reduction of SIN3 complex components results in abnormal wing development. Interestingly, while knockdown of each factor resulted in similar phenotypes, their individual effects on recruitment of SIN3 to polytene chromosomes are distinct. Reduction of some factors leads to large changes in the morphology of the chromosome and/or greatly reduced SIN3 binding. These findings suggest that while individual SIN3 complex components work through distinct molecular mechanisms, they each make a substantial contribution to the overall function of this highly conserved histone deacetylase complex.

Heterochromatin-Enriched Assemblies Reveal the Sequence and Organization of the Drosophila melanogaster Y Chromosome

Heterochromatic regions of the genome are repeat-rich and poor in protein coding genes, and are therefore underrepresented in even the best genome assemblies. One of the most difficult regions of the genome to assemble are sex-limited chromosomes. The Drosophila melanogaster Y chromosome is entirely heterochromatic, yet has wide-ranging effects on male fertility, fitness, and genome-wide gene expression. The genetic basis of this phenotypic variation is difficult to study, in part because we do not know the detailed organization of the Y chromosome. To study Y chromosome organization in D. melanogaster, we develop an assembly strategy involving the in silico enrichment of heterochromatic long single-molecule reads and use these reads to create targeted de novo assemblies of heterochromatic sequences. We assigned contigs to the Y chromosome using Illumina reads to identify male-specific sequences. Our pipeline extends the D. melanogaster reference genome by 11.9 Mb, closes 43.8% of the gaps, and improves overall contiguity. The addition of 10.6 MB of Y-linked sequence permitted us to study the organization of repeats and genes along the Y chromosome. We detected a high rate of duplication to the pericentric regions of the Y chromosome from other regions in the genome. Most of these duplicated genes exist in multiple copies. We detail the evolutionary history of one sex-linked gene family, crystal-Stellate While the Y chromosome does not undergo crossing over, we observed high gene conversion rates within and between members of the crystal-Stellate gene family, Su(Ste), and PCKR, compared to genome-wide estimates. Our results suggest that gene conversion and gene duplication play an important role in the evolution of Y-linked genes.

Sex gap in aging and longevity: can sex chromosomes play a role?

It is well known that women live longer than men. This gap is observed in most human populations and can even reach 10-15 years. In addition, most of the known super centenarians (i.e., humans who lived for > 110 years) are women. The differences in life expectancy between men and women are often attributed to cultural differences in common thinking. However, sex hormones seem to influence differences in the prevalence of diseases, in the magnitude of aging, and in the longevity between men and women. Moreover, far from being human specific, the sex gap in longevity is extremely common in non-human animals, especially in mammals. Biological factors clearly contribute to such a sex gap in aging and longevity. Different hypotheses have been proposed to explain why males and females age and die differently. The cost of sexual selection and sexual dimorphism has long been considered the best explanation for the observed sex gap in aging/longevity. However, the way mitochondria are transmitted (i.e., through females in most species) could have an effect, called the mother's curse. Recent data suggest that sex chromosomes may also contribute to the sex gap in aging/longevity through several potential mechanisms, including the unguarded X/Z, the toxic Y/W and the loss of Y/W. We discuss future research directions to test these ideas.

Faster-X evolution: Theory and evidence from Drosophila

A faster rate of adaptive evolution of X-linked genes compared with autosomal genes can be caused by the fixation of recessive or partially recessive advantageous mutations, due to the full expression of X-linked mutations in hemizygous males. Other processes, including recombination rate and mutation rate differences between X chromosomes and autosomes, may also cause faster evolution of X-linked genes. We review population genetics theory concerning the expected relative values of variability and rates of evolution of X-linked and autosomal DNA sequences. The theoretical predictions are compared with data from population genomic studies of several species of Drosophila. We conclude that there is evidence for adaptive faster-X evolution of several classes of functionally significant nucleotides. We also find evidence for potential differences in mutation rates between X-linked and autosomal genes, due to differences in mutational bias towards GC to AT mutations. Many aspects of the data are consistent with the male hemizygosity model, although not all possible confounding factors can be excluded.

Chromatin-state discovery and genome annotation with ChromHMM

Noncoding DNA regions have central roles in human biology, evolution, and disease. ChromHMM helps to annotate the noncoding genome using epigenomic information across one or multiple cell types. It combines multiple genome-wide epigenomic maps, and uses combinatorial and spatial mark patterns to infer a complete annotation for each cell type. ChromHMM learns chromatin-state signatures using a multivariate hidden Markov model (HMM) that explicitly models the combinatorial presence or absence of each mark. ChromHMM uses these signatures to generate a genome-wide annotation for each cell type by calculating the most probable state for each genomic segment. ChromHMM provides an automated enrichment analysis of the resulting annotations to facilitate the functional interpretations of each chromatin state. ChromHMM is distinguished by its modeling emphasis on combinations of marks, its tight integration with downstream functional enrichment analyses, its speed, and its ease of use. Chromatin states are learned, annotations are produced, and enrichments are computed within 1 d.

Neo-sex Chromosomes in the Monarch Butterfly, Danaus plexippus

We report the discovery of a neo-sex chromosome in the monarch butterfly, Danaus plexippus, and several of its close relatives. Z-linked scaffolds in the D. plexippus genome assembly were identified via sex-specific differences in Illumina sequencing coverage. Additionally, a majority of the D. plexippus genome assembly was assigned to chromosomes based on counts of one-to-one orthologs relative to the butterfly Melitaea cinxia (with replication using two other lepidopteran species), in which genome scaffolds have been mapped to linkage groups. Sequencing coverage-based assessments of Z linkage combined with homology-based chromosomal assignments provided strong evidence for a Z-autosome fusion in the Danaus lineage, involving the autosome homologous to chromosome 21 in M. cinxia. Coverage analysis also identified three notable assembly errors resulting in chimeric Z-autosome scaffolds. Cytogenetic analysis further revealed a large W chromosome that is partially euchromatic, consistent with being a neo-W chromosome. The discovery of a neo-Z and the provisional assignment of chromosome linkage for >90% of D. plexippus genes lays the foundation for novel insights concerning sex chromosome evolution in this female-heterogametic model species for functional and evolutionary genomics.

Underrepresentation of active histone modification marks in evolutionarily young genes

It is known that evolutionarily new genes can rapidly evolve essential roles in fundamental biological processes. Nevertheless, the underlying molecular mechanism of how they acquire their novel transcriptional pattern is less characterized except for the role of cis-regulatory evolution. Epigenetic modification offers an alternative possibility. Here, we examined how histone modifications have changed among different gene age groups in Drosophila melanogaster by integrative analyses of an updated new gene dataset and published epigenomic data. We found a robust pattern across various datasets where both the coverage and intensity of active histone modifications, histone 3 lysine 4 trimethylation and lysine 36 trimethylation, increased with evolutionary age. Such a temporal correlation is negative and much weaker for the repressive histone mark, lysine 9 trimethylation, which is expected given its major association with heterochromatin. By further comparison with neighboring old genes, the depletion of active marks of new genes could be only partially explained by the local epigenetic context. All these data are consistent with the observation that older genes bear relatively higher expression levels and suggest that the evolution of histone modifications could be implicated in transcriptional evolution after gene birth.

Models of life: epigenetics, diversity and cycles

This review emphasizes aspects of biology that can be understood through repeated applications of simple causal rules. The selected topics include perspectives on gene regulation, phage lambda development, epigenetics, microbial ecology, as well as model approaches to diversity and to punctuated equilibrium in evolution. Two outstanding features are repeatedly described. One is the minimal number of rules to sustain specific states of complex systems for a long time. The other is the collapse of such states and the subsequent dynamical cycle of situations that restitute the system to a potentially new metastable state.

Helitrons shaping the genomic architecture of Drosophila: enrichment of DINE-TR1 in α- and β-heterochromatin, satellite DNA emergence, and piRNA expression

Drosophila INterspersed Elements (DINEs) constitute an abundant but poorly understood group of Helitrons present in several Drosophila species. The general structure of DINEs includes two conserved blocks that may or not contain a region with tandem repeats in between. These central tandem repeats (CTRs) are similar within species but highly divergent between species. It has been assumed that CTRs have independent origins. Herein, we identify a subset of DINEs, termed DINE-TR1, which contain homologous CTRs of approximately 150 bp. We found DINE-TR1 in the sequenced genomes of several Drosophila species and in Bactrocera tryoni (Acalyptratae, Diptera). However, interspecific high sequence identity (∼ 88 %) is limited to the first ∼ 30 bp of each tandem repeat, implying that evolutionary constraints operate differently over the monomer length. DINE-TR1 is unevenly distributed across the Drosophila phylogeny. Nevertheless, sequence analysis suggests vertical transmission. We found that CTRs within DINE-TR1 have independently expanded into satellite DNA-like arrays at least twice within Drosophila. By analyzing the genome of Drosophila virilis and Drosophila americana, we show that DINE-TR1 is highly abundant in pericentromeric heterochromatin boundaries, some telomeric regions and in the Y chromosome. It is also present in the centromeric region of one autosome from D. virilis and dispersed throughout several euchromatic sites in both species. We further found that DINE-TR1 is abundant at piRNA clusters, and small DINE-TR1-derived RNA transcripts (∼25 nt) are predominantly expressed in the testes and the ovaries, suggesting active targeting by the piRNA machinery. These features suggest potential piRNA-mediated regulatory roles for DINEs at local and genome-wide scales in Drosophila.

Expansion of GA Dinucleotide Repeats Increases the Density of CLAMP Binding Sites on the X-Chromosome to Promote Drosophila Dosage Compensation

Dosage compensation is an essential process that equalizes transcript levels of X-linked genes between sexes by forming a domain of coordinated gene expression. Throughout the evolution of Diptera, many different X-chromosomes acquired the ability to be dosage compensated. Once each newly evolved X-chromosome is targeted for dosage compensation in XY males, its active genes are upregulated two-fold to equalize gene expression with XX females. In Drosophila melanogaster, the CLAMP zinc finger protein links the dosage compensation complex to the X-chromosome. However, the mechanism for X-chromosome identification has remained unknown. Here, we combine biochemical, genomic and evolutionary approaches to reveal that expansion of GA-dinucleotide repeats likely accumulated on the X-chromosome over evolutionary time to increase the density of CLAMP binding sites, thereby driving the evolution of dosage compensation. Overall, we present new insight into how subtle changes in genomic architecture, such as expansions of a simple sequence repeat, promote the evolution of coordinated gene expression.

Genomics of Natural Populations: How Differentially Expressed Genes Shape the Evolution of Chromosomal Inversions in Drosophila pseudoobscura

Chromosomal rearrangements can shape the structure of genetic variation in the genome directly through alteration of genes at breakpoints or indirectly by holding combinations of genetic variants together due to reduced recombination. The third chromosome of Drosophila pseudoobscura is a model system to test hypotheses about how rearrangements are established in populations because its third chromosome is polymorphic for >30 gene arrangements that were generated by a series of overlapping inversion mutations. Circumstantial evidence has suggested that these gene arrangements are selected. Despite the expected homogenizing effects of extensive gene flow, the frequencies of arrangements form gradients or clines in nature, which have been stable since the system was first described >80 years ago. Furthermore, multiple arrangements exist at appreciable frequencies across several ecological niches providing the opportunity for heterokaryotypes to form. In this study, we tested whether genes are differentially expressed among chromosome arrangements in first instar larvae, adult females and males. In addition, we asked whether transcriptional patterns in heterokaryotypes are dominant, semidominant, overdominant, or underdominant. We find evidence for a significant abundance of differentially expressed genes across the inverted regions of the third chromosome, including an enrichment of genes involved in sensory perception for males. We find the majority of loci show additivity in heterokaryotypes. Our results suggest that multiple genes have expression differences among arrangements that were either captured by the original inversion mutation or accumulated after it reached polymorphic frequencies, providing a potential source of genetic variation for selection to act upon. These data suggest that the inversions are favored because of their indirect effect of recombination suppression that has held different combinations of differentially expressed genes together in the various gene arrangement backgrounds.

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.

RNA-sequencing and whole-genome bisulfite sequencing in the hybrid offspring of two Nasonia parasitoid wasp species revealed strong cis-regulation of methylation and allele-specific expression. No gene was found to display genomic imprinting.

The relationship between methylation of genomic DNA and expression of the genes that it encodes—and how this relationship changes during evolution—has been widely studied in mammals, but remains less well understood for insects. Here we analyze the expressed mRNA transcripts and genomic DNA methylation of the hybrid offspring of a pair of Nasonia parasitoid wasp species, producing a wealth of information about the regulation of gene expression. We find that variation in DNA sequence impacts expression on the same strand (called “cis-regulation”), and that cytosine methylation state is also associated in cis with the regulatory consequences of this base alteration. We show that these wasp species lack differential expression dependent on parent-of-origin (called “genomic imprinting”), and that in the hybrids the alleles retain the methylation status of the parental species in a strong cis-regulated fashion. Transcript abundances were also largely driven in a cis-regulated manner, consistent with a correlation between methylation status and expression levels. Despite the many differences between Nasonia and mammals in the impact of genomic DNA methylation, in both groups the use of methylated cytosine has been co-opted in ways that help tune gene expression.

Dosage Compensation and the Distribution of Sex-Biased Gene Expression in Drosophila: Considerations and Genomic Constraints

Several studies in Drosophila have shown a paucity of male-biased genes (i.e., genes that express higher in males than in females) on the X chromosome. Dosage compensation (DC) is a regulatory mechanism of gene expression triggered in males that hypertranscribes the X-linked genes to the level of transcription in females. There are currently two different hypotheses about the effects of DC on the distribution of male-biased genes: (1) it might limit male-expression level, or (2) it might interfere with the male upregulation of gene expression. Here, we used previously published gene expression datasets to reevaluate both hypotheses and introduce a mutually exclusive prediction that helped us to reject the hypothesis that the paucity of male-biased genes in the X chromosome is due to a limit in the male-expression level. Our analysis also uncovers unanticipated details about how DC interferes with the genomic distribution of both, male-biased and female-biased genes. We suggest that DC actually interferes with female downregulation of gene expression and not male upregulation, as previously suggested.

The histone demethylase dKDM5/LID interacts with the SIN3 histone deacetylase complex and shares functional similarities with SIN3

Background

Regulation of gene expression by histone-modifying enzymes is essential to control cell fate decisions and developmental processes. Two histone-modifying enzymes, RPD3, a deacetylase, and dKDM5/LID, a demethylase, are present in a single complex, coordinated through the SIN3 scaffold protein. While the SIN3 complex has been demonstrated to have functional histone deacetylase activity, the role of the demethylase dKDM5/LID as part of the complex has not been investigated.

Results

Here, we analyzed the developmental and transcriptional activities of dKDM5/LID in relation to SIN3. Knockdown of either Sin3A or lid resulted in decreased cell proliferation in S2 cells and wing imaginal discs. Conditional knockdown of either Sin3A or lid resulted in flies that displayed wing developmental defects. Interestingly, overexpression of dKDM5/LID rescued the wing developmental defect due to reduced levels of SIN3 in female flies, indicating a major role for dKDM5/LID in cooperation with SIN3 during development. Together, these observed phenotypes strongly suggest that dKDM5/LID as part of the SIN3 complex can impact previously uncharacterized transcriptional networks. Transcriptome analysis revealed that SIN3 and dKDM5/LID regulate many common genes. While several genes implicated in cell cycle and wing developmental pathways were affected upon altering the level of these chromatin factors, a significant affect was also observed on genes required to mount an effective stress response. Further, under conditions of induced oxidative stress, reduction of SIN3 and/or dKDM5/LID altered the expression of a greater number of genes involved in cell cycle-related processes relative to normal conditions. This highlights an important role for SIN3 and dKDM5/LID proteins to maintain proper progression through the cell cycle in environments of cellular stress. Further, we find that target genes are bound by both SIN3 and dKDM5/LID, however, histone acetylation, not methylation, plays a predominant role in gene regulation by the SIN3 complex.

Conclusions

We have provided genetic evidence to demonstrate functional cooperation between the histone demethylase dKDM5/LID and SIN3. Biochemical and transcriptome data further support functional links between these proteins. Together, the data provide a solid framework for analyzing the gene regulatory pathways through which SIN3 and dKDM5/LID control diverse biological processes in the organism.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0053-9) contains supplementary material, which is available to authorized users.

Chromatin state analysis of the barley epigenome reveals a higher-order structure defined by H3K27me1 and H3K27me3 abundance

Combinations of histones carrying different covalent modifications are a major component of epigenetic variation. We have mapped nine modified histones in the barley seedling epigenome by chromatin immunoprecipitation next-generation sequencing (ChIP-seq). The chromosomal distributions of the modifications group them into four different classes, and members of a given class also tend to coincide at the local DNA level, suggesting that global distribution patterns reflect local epigenetic environments. We used this peak sharing to define 10 chromatin states representing local epigenetic environments in the barley genome. Five states map mainly to genes and five to intergenic regions. Two genic states involving H3K36me3 are preferentially associated with constitutive gene expression, while an H3K27me3-containing genic state is associated with differentially expressed genes. The 10 states display striking distribution patterns that divide barley chromosomes into three distinct global environments. First, telomere-proximal regions contain high densities of H3K27me3 covering both genes and intergenic DNA, together with very low levels of the repressive H3K27me1 modification. Flanking these are gene-rich interior regions that are rich in active chromatin states and have greatly decreased levels of H3K27me3 and increasing amounts of H3K27me1 and H3K9me2. Lastly, H3K27me3-depleted pericentromeric regions contain gene islands with active chromatin states separated by extensive retrotransposon-rich regions that are associated with abundant H3K27me1 and H3K9me2 modifications. We propose an epigenomic framework for barley whereby intergenic H3K27me3 specifies facultative heterochromatin in the telomere-proximal regions and H3K27me1 is diagnostic for constitutive heterochromatin elsewhere in the barley genome.

Genome architecture: from linear organisation of chromatin to the 3D assembly in the nucleus

The genetic information is stored in the eukaryotic nucleus in the form of chromatin. This is a macromolecular entity that includes genomic DNA and histone proteins that form nucleosomes, plus a large variety of chromatin-associated non-histone proteins. Chromatin is structurally and functionally organised at various levels. One reveals the linear topography of DNA, histones and their post-translational modifications and non-histone proteins along each chromosome. This level provides regulatory information about the association of genomic elements with particular signatures that have been used to define chromatin states. Importantly, these chromatin states correlate with structural and functional genomic features. Another regulatory layer is established at the level of the 3D organisation of chromatin within the nucleus, which has been revealed clearly as non-random. Instead, a variety of intra- and inter-chromosomal genomic domains with specific epigenetic and functional properties has been identified. In this review, we discuss how the recent advances in genomic approaches have contributed to our understanding of these two levels of genome architecture. We have emphasised our analysis with the aim of integrating information available for yeast, Arabidopsis, Drosophila, and mammalian cells. We consider that this comparative study helps define common and unique features in each system, providing a basis to better understand the complexity of genome organisation.

Insights into the epigenomic landscape of the human malaria vector Anopheles gambiae

The epigenome of the human malaria vector Anopheles gambiae was characterized in midgut cells by mapping the distribution and levels of two post-translational histone modifications, H3K27ac and H3K27me3. These histone profiles were then correlated with levels of gene expression obtained by RNA-seq. Analysis of the transcriptome of A. gambiae midguts and salivary glands led to the discovery of 13,898 new transcripts not present in the most recent genome assembly. A subset of these transcripts is differentially expressed between midgut and salivary glands. The enrichment profiles of H3K27ac and H3K27me3 are mutually exclusive and associate with high and low levels of transcription, respectively. This distribution agrees with previous findings in Drosophila showing association of these two histone modifications with either active or inactive transcriptional states, including Polycomb-associated domains in silenced genes. This study provides a mosquito epigenomics platform for future comparative studies in other mosquito species, opening future investigations into the role of epigenetic processes in vector-borne systems of medical and economic importance.