The evolution and mutational robustness of chromatin accessibility in Drosophila

BACKGROUND

The evolution of genomic regulatory regions plays a critical role in shaping the diversity of life. While this process is primarily sequence-dependent, the enormous complexity of biological systems complicates the understanding of the factors underlying regulation and its evolution. Here, we apply deep neural networks as a tool to investigate the sequence determinants underlying chromatin accessibility in different species and tissues of Drosophila.

RESULTS

We train hybrid convolution-attention neural networks to accurately predict ATAC-seq peaks using only local DNA sequences as input. We show that our models generalize well across substantially evolutionarily diverged species of insects, implying that the sequence determinants of accessibility are highly conserved. Using our model to examine species-specific gains in accessibility, we find evidence suggesting that these regions may be ancestrally poised for evolution. Using in silico mutagenesis, we show that accessibility can be accurately predicted from short subsequences in each example. However, in silico knock-out of these sequences does not qualitatively impair classification, implying that accessibility is mutationally robust. Subsequently, we show that accessibility is predicted to be robust to large-scale random mutation even in the absence of selection. Conversely, simulations under strong selection demonstrate that accessibility can be extremely malleable despite its robustness. Finally, we identify motifs predictive of accessibility, recovering both novel and previously known motifs.

CONCLUSIONS

These results demonstrate the conservation of the sequence determinants of accessibility and the general robustness of chromatin accessibility, as well as the power of deep neural networks to explore fundamental questions in regulatory genomics and evolution.

The evolution and mutational robustness of chromatin accessibility in Drosophila

The evolution of regulatory regions in the genome plays a critical role in shaping the diversity of life. While this process is primarily sequence-dependent, the enormous complexity of biological systems has made it difficult to understand the factors underlying regulation and its evolution. Here, we apply deep neural networks as a tool to investigate the sequence determinants underlying chromatin accessibility in different tissues of Drosophila. We train hybrid convolution-attention neural networks to accurately predict ATAC-seq peaks using only local DNA sequences as input. We show that a model trained in one species has nearly identical performance when tested in another species, implying that the sequence determinants of accessibility are highly conserved. Indeed, model performance remains excellent even in distantly-related species. By using our model to examine species-specific gains in chromatin accessibility, we find that their orthologous inaccessible regions in other species have surprisingly similar model outputs, suggesting that these regions may be ancestrally poised for evolution. We then use in silico saturation mutagenesis to reveal evidence of selective constraint acting specifically on inaccessible chromatin regions. We further show that chromatin accessibility can be accurately predicted from short subsequences in each example. However, in silico knock-out of these sequences does not qualitatively impair classification, implying that chromatin accessibility is mutationally robust. Subsequently, we demonstrate that chromatin accessibility is predicted to be robust to large-scale random mutation even in the absence of selection. We also perform in silico evolution experiments under the regime of strong selection and weak mutation (SSWM) and show that chromatin accessibility can be extremely malleable despite its mutational robustness. However, selection acting in different directions in a tissue-specific manner can substantially slow adaptation. Finally, we identify motifs predictive of chromatin accessibility and recover motifs corresponding to known chromatin accessibility activators and repressors. These results demonstrate the conservation of the sequence determinants of accessibility and the general robustness of chromatin accessibility, as well as the power of deep neural networks as tools to answer fundamental questions in regulatory genomics and evolution.

Environmental and Genetic Contributions to Imperfect wMel-Like Wolbachia Transmission and Frequency Variation

Maternally transmitted Wolbachia bacteria infect about half of all insect species. They usually show imperfect maternal transmission and often produce cytoplasmic incompatibility (CI). Irrespective of CI, Wolbachia frequencies tend to increase when rare only if they benefit host fitness. Several Wolbachia, including wMel that infects Drosophila melanogaster, cause weak or no CI and persist at intermediate frequencies. On the island of São Tomé off West Africa, the frequencies of wMel-like Wolbachia infecting Drosophila yakuba (wYak) and Drosophila santomea (wSan) fluctuate, and the contributions of imperfect maternal transmission, fitness effects, and CI to these fluctuations are unknown. We demonstrate spatial variation in wYak frequency and transmission on São Tomé. Concurrent field estimates of imperfect maternal transmission do not predict spatial variation in wYak frequencies, which are highest at high altitudes where maternal transmission is the most imperfect. Genomic and genetic analyses provide little support for D. yakuba effects on wYak transmission. Instead, rearing at cool temperatures reduces wYak titer and increases imperfect transmission to levels observed on São Tomé. Using mathematical models of Wolbachia frequency dynamics and equilibria, we infer that temporally variable imperfect transmission or spatially variable effects on host fitness and reproduction are required to explain wYak frequencies. In contrast, spatially stable wSan frequencies are plausibly explained by imperfect transmission, modest fitness effects, and weak CI. Our results provide insight into causes of wMel-like frequency variation in divergent hosts. Understanding this variation is crucial to explain Wolbachia spread and to improve wMel biocontrol of human disease in transinfected mosquito systems.

Recurrent specialization on a toxic fruit in an island Drosophila population

Recurrent specialization on similar host plants offers a unique opportunity to unravel the evolutionary and genetic mechanisms underlying dietary shifts. Recent studies have focused on ecological races belonging to the same species, but it is hard in many cases to untangle the role of adaptive introgression versus distinct mutations in facilitating recurrent evolution. We discovered on the island of Mayotte a population of the generalist fly Drosophila yakuba that is strictly associated with noni (Morinda citrifolia). This case strongly resembles Drosophila sechellia, a genetically isolated insular relative of D. yakuba whose intensely studied specialization on toxic noni fruits has always been considered a unique event in insect evolution. Experiments revealed that unlike mainland D. yakuba strains, Mayotte flies showed strong olfactory attraction and significant toxin tolerance to noni. Island females strongly discriminated against mainland males, suggesting that dietary adaptation has been accompanied by partial reproductive isolation. Population genomic analysis indicated a recent colonization (∼29 kya), at a time when year-round noni fruits may have presented a predictable resource on the small island, with ongoing migration after colonization. This relatively recent time scale allowed us to search for putatively adaptive loci based on genetic variation. Strong signals of genetic differentiation were found for several detoxification genes, including a major toxin tolerance locus in D. sechellia Our results suggest that recurrent evolution on a toxic resource can involve similar historical events and common genetic bases, and they establish an important genetic system for the study of early stages of ecological specialization and speciation.

Hydrocarbon Patterns and Mating Behaviour in Populations of Drosophila yakuba

Drosophila yakuba is widespread in Africa. Here we compare the cuticular hydrocarbon (CHC) profiles and mating behavior of mainland (Kounden, Cameroon) and island (Mayotte, Sao-Tome, Bioko) populations. The strains each had different CHC profiles: Bioko and Kounden were the most similar, while Mayotte and Sao-Tome contained significant amounts of 7-heptacosene. The CHC profile of the Sao-Tome population differed the most, with half the 7-tricosene of the other populations and more 7-heptacosene and 7-nonacosene. We also studied the characteristics of the mating behavior of the four strains: copulation duration was similar but latency times were higher in Mayotte and Sao-Tome populations. We found partial reproductive isolation between populations, especially in male-choice experiments with Sao-Tome females.

Landscape of standing variation for tandem duplications in Drosophila yakuba and Drosophila simulans

We have used whole genome paired-end Illumina sequence data to identify tandem duplications in 20 isofemale lines of Drosophila yakuba and 20 isofemale lines of D. simulans and performed genome wide validation with PacBio long molecule sequencing. We identify 1,415 tandem duplications that are segregating in D. yakuba as well as 975 duplications in D. simulans, indicating greater variation in D. yakuba. Additionally, we observe high rates of secondary deletions at duplicated sites, with 8% of duplicated sites in D. simulans and 17% of sites in D. yakuba modified with deletions. These secondary deletions are consistent with the action of the large loop mismatch repair system acting to remove polymorphic tandem duplication, resulting in rapid dynamics of gain and loss in duplicated alleles and a richer substrate of genetic novelty than has been previously reported. Most duplications are present in only single strains, suggesting that deleterious impacts are common. Drosophila simulans shows larger numbers of whole gene duplications in comparison to larger proportions of gene fragments in D. yakuba. Drosophila simulans displays an excess of high-frequency variants on the X chromosome, consistent with adaptive evolution through duplications on the D. simulans X or demographic forces driving duplicates to high frequency. We identify 78 chimeric genes in D. yakuba and 38 chimeric genes in D. simulans, as well as 143 cases of recruited noncoding sequence in D. yakuba and 96 in D. simulans, in agreement with rates of chimeric gene origination in D. melanogaster. Together, these results suggest that tandem duplications often result in complex variation beyond whole gene duplications that offers a rich substrate of standing variation that is likely to contribute both to detrimental phenotypes and disease, as well as to adaptive evolutionary change.

A novel dataset for identifying sex-biased genes in Drosophila

Phenotypic differences between males and females of sexually dimorphic species are caused in large part by differences in gene expression between the sexes, most of which occurs in the gonads. To accurately identify genes differentially expressed between males and females in Drosophila, we sequenced the testis and ovary transcriptomes of D. yakuba, D. pseudoobscura, and D. ananassae and used them to identify sex-biased genes in the latter two species. We highlight the increased sensitivity and improved power of sex-biased gene detection methods when using our testis/ovary data versus male and female whole body transcriptome data. We thus provide a resource specifically designed to accurately identify and characterize sex-biased genes across Drosophila. This dataset is available through NCBI GEO accession GSE52058.

Comparative gene expression analysis of Dtg, a novel target gene of Dpp signaling pathway in the early Drosophila melanogaster embryo

In the early Drosophila melanogaster embryo, Dpp, a secreted molecule that belongs to the TGF-β superfamily of growth factors, activates a set of downstream genes to subdivide the dorsal region into amnioserosa and dorsal epidermis. Here, we examined the expression pattern and transcriptional regulation of Dtg, a new target gene of Dpp signaling pathway that is required for proper amnioserosa differentiation. We showed that the expression of Dtg was controlled by Dpp and characterized a 524-bp enhancer that mediated expression in the dorsal midline, as well as, in the differentiated amnioserosa in transgenic reporter embryos. This enhancer contained a highly conserved region of 48-bp in which bioinformatic predictions and in vitro assays identified three Mad binding motifs. Mutational analysis revealed that these three motifs were necessary for proper expression of a reporter gene in transgenic embryos, suggesting that short and highly conserved genomic sequences may be indicative of functional regulatory regions in D. melanogaster genes. Dtg orthologs were not detected in basal lineages of Dipterans, which unlike D. melanogaster develop two extra-embryonic membranes, amnion and serosa, nevertheless Dtg orthologs were identified in the transcriptome of Musca domestica, in which dorsal ectoderm patterning leads to the formation of a single extra-embryonic membrane. These results suggest that Dtg was recruited as a new component of the network that controls dorsal ectoderm patterning in the lineage leading to higher Cyclorrhaphan flies, such as D. melanogaster and M. domestica.