Accumulation of endosymbiont genomes in an insect autosome followed by endosymbiont replacement

Eukaryotic genomes can acquire bacterial DNA via lateral gene transfer (LGT).¹ A prominent source of LGT is Wolbachia,² a widespread endosymbiont of arthropods and nematodes that is transmitted maternally through female germline cells.^3,4 The DNA transfer from the Wolbachia endosymbiont wAna to Drosophila ananassae is extensive^5-7 and has been localized to chromosome 4, contributing to chromosome expansion in this lineage.⁶ As has happened frequently with claims of bacteria-to-eukaryote LGT, the contribution of wAna transfers to the expanded size of D. ananassae chromosome 4 has been specifically contested⁸ owing to an assembly where Wolbachia sequences were classified as contaminants and removed.⁹ Here, long-read sequencing with DNA from a Wolbachia-cured line enabled assembly of 4.9 Mbp of nuclear Wolbachia transfers (nuwts) in D. ananassae and a 24-kbp nuclear mitochondrial transfer. The nuwts are <8,000 years old in at least two locations in chromosome 4 with at least one whole-genome integration followed by rapid extensive duplication of most of the genome with regions that have up to 10 copies. The genes in nuwts are accumulating small indels and mobile element insertions. Among the highly duplicated genes are cifA and cifB, two genes associated with Wolbachia-mediated Drosophila cytoplasmic incompatibility. The wAna strain that was the source of nuwts was subsequently replaced by a different wAna endosymbiont. Direct RNA Nanopore sequencing of Wolbachia-cured lines identified nuwt transcripts, including spliced transcripts, but functionality, if any, remains elusive.

Three Quantitative Trait Loci Explain More than 60% of Variation for Chill Coma Recovery Time in a Natural Population of Drosophila ananassae

Ectothermic species such as insects are particularly vulnerable to climatic fluctuations. Nevertheless, many insects that evolved and diversified in the tropics have successfully colonized temperate regions all over the globe. To shed light on the genetic basis of cold tolerance in such species, we conducted a quantitative trait locus (QTL) mapping experiment for chill coma recovery time (CCRT) in Drosophila ananassae, a cosmopolitan species that has expanded its range from tropical to temperate regions. We created a mapping population of recombinant inbred advanced intercross lines (RIAILs) from two founder strains with diverging CCRT phenotypes. The RIAILs were phenotyped for their CCRT and, together with the founder strains, genotyped for polymorphic markers with double-digest restriction site-associated DNA (ddRAD) sequencing. Using a hierarchical mapping approach that combined standard interval mapping and a multiple-QTL model, we mapped three QTL which altogether explained 64% of the phenotypic variance. For two of the identified QTL, we found evidence of epistasis. To narrow down the list of cold tolerance candidate genes, we cross-referenced the QTL intervals with genes that we previously identified as differentially expressed in response to cold in D. ananassae, and with thermotolerance candidate genes of D. melanogaster. Among the 58 differentially expressed genes that were contained within the QTL, GF15058 showed a significant interaction of the CCRT phenotype and gene expression. Further, we identified the orthologs of four D. melanogaster thermotolerance candidate genes, MtnA, klarsicht, CG5246 (D.ana/GF17132) and CG10383 (D.ana/GF14829) as candidates for cold tolerance in D. ananassae.

Transcriptome Analysis Reveals Candidate Genes for Cold Tolerance in Drosophila ananassae

Coping with daily and seasonal temperature fluctuations is a key adaptive process for species to colonize temperate regions all over the globe. Over the past 18,000 years, the tropical species Drosophila ananassae expanded its home range from tropical regions in Southeast Asia to more temperate regions. Phenotypic assays of chill coma recovery time (CCRT) together with previously published population genetic data suggest that only a small number of genes underlie improved cold hardiness in the cold-adapted populations. We used high-throughput RNA sequencing to analyze differential gene expression before and after exposure to a cold shock in cold-tolerant lines (those with fast chill coma recovery, CCR) and cold-sensitive lines (slow CCR) from a population originating from Bangkok, Thailand (the ancestral species range). We identified two candidate genes with a significant interaction between cold tolerance and cold shock treatment: GF14647 and GF15058. Further, our data suggest that selection for increased cold tolerance did not operate through the increased activity of heat shock proteins, but more likely through the stabilization of the actin cytoskeleton and a delayed onset of apoptosis.

Tropical Drosophila ananassae of wet-dry seasons show cross resistance to heat, drought and starvation

Plastic responses to multiple environmental stressors in wet or dry seasonal populations of tropical Drosophila species have received less attention. We tested plastic effects of heat hardening, acclimation to drought or starvation, and changes in trehalose, proline and body lipids in Drosophila ananassae flies reared under wet or dry season-specific conditions. Wet season flies revealed significant increase in heat knockdown, starvation resistance and body lipids after heat hardening. However, accumulation of proline was observed only after desiccation acclimation of dry season flies while wet season flies elicited no proline but trehalose only. Therefore, drought-induced proline can be a marker metabolite for dry-season flies. Further, partial utilization of proline and trehalose under heat hardening reflects their possible thermoprotective effects. Heat hardening elicited cross-protection to starvation stress. Stressor-specific accumulation or utilization as well as rates of metabolic change for each energy metabolite were significantly higher in wet-season flies than dry-season flies. Energy metabolite changes due to inter-related stressors (heat versus desiccation or starvation) resulted in possible maintenance of energetic homeostasis in wet- or dry-season flies. Thus, low or high humidity-induced plastic changes in energy metabolites can provide cross-protection to seasonally varying climatic stressors.

Summary: In the tropical Drosophila ananassae, low or high humidity-induced plastic changes in energy metabolites provide cross-protection to seasonally varying climatic stressors.

A Natural Population Derived from Species Hybridizationin the Drosophila ananassae Species Complexon Penang Island, Malaysia

We surveyed natural population of the Drosophila ananassae species complex on Penang Island, Malaysia. Analyses of phenotypic traits, chromosome arrangements, molecular markers, and reproductive isolation suggest the existence of two species: D. ananassae and D. cf. parapallidosa. Molecular marker analysis indicates that D. cf. parapallidosa carries chromosome Y and 4 introgressions from D. ananassae. Thus, D. cf. parapallidosa seems to be a hybrid descendant that recently originated from a natural D. parapallidosa♀× D. ananassae♂ cross. Furthermore, D. cf. parapallidosa behaves differently from authentic D. parapallidosa with respect to its reproductive isolation from D. ananassae. Premating isolation is usually seen in only the D. ananassae♀× D. parapallidosa♂ cross, but we observed it in crosses of both directions between D. ananassae and D. cf. parapallidosa. In addition, hybrid males from the D. ananassae♀× D. parapallidosa♂ cross are usually sterile, but they were fertile when D. ananassae♀ were mated with D. cf. parapallidosa ♂. We attempted an artificial reconstruction of the hybrid species to simulate the evolutionary process(es) that produced D. cf. parapallidosa. This is a rare case of natural hybrid population in Drosophila and may be a useful system for elucidating speciation with gene flow.

Recent and Long-Term Selection Across Synonymous Sites in Drosophila ananassae

In Drosophila, many studies have examined the short- or long-term evolution occurring across synonymous sites. Few, however, have examined both the recent and long-term evolution to gain a complete view of this selection. Here we have analyzed Drosophila ananassae DNA polymorphism and divergence data using several different methods, and have identified evidence of positive selection favoring preferred codons in both recent and long-term evolutionary time scale. Further in D. ananassae, the strength of selection for preferred codons was stronger on the X chromosome compared to the autosomes. We show that this stronger selection is not due to higher gene expression of X-linked genes. Analysis of the selectively neutral introns indicated that the X chromosome also had a preference for GC over AT nucleotides, potentially from GC-biased gene conversions (gcBGCs) that can also affect the base composition of synonymous sites. Thus selection for preferred codons and gcBGC both seem to be partially responsible for shaping the D. ananassae synonymous site evolution.

The coevolutionary period of Wolbachia pipientis infecting Drosophila ananassae and its impact on the evolution of the host germline stem cell regulating genes

The endosymbiotic bacteria Wolbachia pipientis is known to infect a wide range of arthropod species yet less is known about the coevolutionary history it has with its hosts. Evidence of highly identical W. pipientis strains in evolutionary divergent hosts suggests horizontal transfer between hosts. For example, Drosophila ananassae is infected with a W. pipientis strain that is nearly identical in sequence to a strain that infects both D. simulans and D. suzukii, suggesting recent horizontal transfer among these three species. However, it is unknown whether the W. pipientis strain had recently invaded all three species or a more complex infectious dynamic underlies the horizontal transfers. Here, we have examined the coevolutionary history of D. ananassae and its resident W. pipientis to infer its period of infection. Phylogenetic analysis of D. ananassae mitochondrial DNA and W. pipientis DNA sequence diversity revealed the current W. pipientis infection is not recent. In addition, we examined the population genetics and molecular evolution of several germline stem cell (GSC) regulating genes of D. ananassae. These studies reveal significant evidence of recent and long-term positive selection at stonewall in D. ananassae, whereas pumillio showed patterns of variation consistent with only recent positive selection. Previous studies had found evidence for adaptive evolution of two key germline differentiation genes, bag of marbles (bam) and benign gonial cell neoplasm (bgcn), in D. melanogaster and D. simulans and proposed that the adaptive evolution at these two genes was driven by arms race between the host GSC and W. pipientis. However, we did not find any statistical departures from a neutral model of evolution for bam and bgcn in D. ananassae despite our new evidence that this species has been infected with W. pipientis for a period longer than the most recent infection in D. melanogaster. In the end, analyzing the GSC regulating genes individually showed two of the seven genes to have evidence of selection. However, combining the data set and fitting a specific population genetic model significant proportion of the nonsynonymous sites across the GSC regulating genes were driven to fixation by positive selection. Clearly the GSC system is under rapid evolution and potentially multiple drivers are causing the rapid evolution.

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.

Insights into differential activity patterns of drosophilids under semi-natural conditions

We showed recently that Drosophila ananassae, a closely related and sympatric species of the commonly studied fruitfly D. melanogaster, shows distinctly deviant patterns in circadian activity/rest rhythm from the latter under a variety of laboratory conditions. To examine whether such differences extend to more natural conditions where a variety of time cues and similar environmental pressures might force different species to adopt similar temporal patterns, we examined these two species under semi-natural conditions over a span of 1.5 years. Furthermore, we asked to what extent features of activity/rest rhythm of flies are conserved across species under changing environmental conditions encountered across seasons, and to do so, we studied two more drosophilid species. We found that while each species exhibits seasonality in activity patterns, this seasonality is marked by interesting inter-specific differences. Similar to laboratory studies, D. ananassae showed activity mostly during the day, while D. melanogaster and D. malerkotliana exhibited almost similar activity patterns across seasons, with predominantly two peaks of activity, one in the morning and another in the evening. Throughout the year, Zaprionus indianus displayed very low levels of activity compared with D. melanogaster, yet, compared with those seen in standard laboratory assays, this species exhibited more robust rhythms under semi-natural conditions. We hypothesise that different ecological factors may have influenced these species to adopt different temporal niches.