The evolutionary novelty of insect defensins: from bacterial killing to toxin neutralization.

Insect host defense comprises two complementary dimensions, microbial killing-mediated resistance and microbial toxin neutralization-mediated resilience, both jointly providing protection against pathogen infections. Insect defensins are a class of effectors of innate immunity primarily responsible for resistance to Gram-positive bacteria. Here, we report a newly originated gene from an ancestral defensin via genetic deletion following gene duplication in Drosophila virilis, which confers an enhanced resilience to Gram-positive bacterial infection. This gene encodes an 18-mer arginine-rich peptide (termed DvirARP) with differences from its parent gene in its pattern of expression, structure and function. DvirARP specifically expresses in D. virilis female adults with a constitutive manner. It adopts a novel fold with a 310 helix and a two CXC motif-containing loop stabilized by two disulfide bridges. DvirARP exhibits no activity on the majority of microorganisms tested and only a weak activity against two Gram-positive bacteria. DvirARP knockout flies are viable and have no obvious defect in reproductivity but they are more susceptible to the DvirARP-resistant Staphylococcus aureus infection than the wild type files, which can be attributable to its ability in neutralization of the S. aureus secreted toxins. Phylogenetic distribution analysis reveals that DvirARP is restrictedly present in the Drosophila subgenus, but independent deletion variations also occur in defensins from the Sophophora subgenus, in support of the evolvability of this class of immune effectors. Our work illustrates for the first time how a duplicate resistance-mediated gene evolves an ability to increase the resilience of a subset of Drosophila species against bacterial infection.

Changes in the cellular makeup of motor patterning circuits drive courtship song evolution in Drosophila.

How evolutionary changes in genes and neurons encode species variation in complex motor behaviors is largely unknown. Here, we develop genetic tools that permit a neural circuit comparison between the model species Drosophila melanogaster and the closely related species D. yakuba, which has undergone a lineage-specific loss of sine song, one of the two major types of male courtship song in Drosophila. Neuroanatomical comparison of song-patterning neurons called TN1 across the phylogeny demonstrates a link between the loss of sine song and a reduction both in the number of TN1 neurons and the neurites supporting the sine circuit connectivity. Optogenetic activation confirms that TN1 neurons in D. yakuba have lost the ability to drive sine song, although they have maintained the ability to drive the singing wing posture. Single-cell transcriptomic comparison shows that D. yakuba specifically lacks a cell type corresponding to TN1A neurons, the TN1 subtype that is essential for sine song. Genetic and developmental manipulation reveals a functional divergence of the sex determination gene doublesex in D. yakuba to reduce TN1 number by promoting apoptosis. Our work illustrates the contribution of motor patterning circuits and cell type changes in behavioral evolution and uncovers the evolutionary lability of sex determination genes to reconfigure the cellular makeup of neural circuits.

First Records of the Invading Species Drosophila Nasuta (Diptera: Drosophilidae) in the Amazon.

Invading species pose a growing threat to biodiversity, ecosystemic systems, regional economies, and public health. In recent decades, South America has received five exotic drosophilids species, some of which have invaded natural ecosystems and caused harm to agriculture. The most recent case is the Asian fly Drosophila nasuta Lamb. In the present study, we record D. nasuta in the Amazon, which is the largest and most biodiverse tropical rainforest in the world. Sampling of drosophilids was carried out between 2012 and 2017 in the Brazilian state of Pará. Drosophila nasuta was first detected on 1st July 2017, with 145 individuals of this species sampled among the 11,496 drosophilids caught. Although at low abundance, D. nasuta was recorded in forest fragments, anthropized fields, and urban environment. The records of the species occurred in the six municipalities of the state of Pará investigated at locations separated by approximately 700 km. In less than 10 years, D. nasuta has occupied approximately 2.5 million km2 in South America. The present findings assist in understanding the susceptibility of tropical forests to biological invasions.

A framework for studying behavioral evolution by reconstructing ancestral repertoires.

Although different animal species often exhibit extensive variation in many behaviors, typically scientists examine one or a small number of behaviors in any single study. Here, we propose a new framework to simultaneously study the evolution of many behaviors. We measured the behavioral repertoire of individuals from six species of fruit flies using unsupervised techniques and identified all stereotyped movements exhibited by each species. We then fit a Generalized Linear Mixed Model to estimate the intra- and inter-species behavioral covariances, and, by using the known phylogenetic relationships among species, we estimated the (unobserved) behaviors exhibited by ancestral species. We found that much of intra-specific behavioral variation has a similar covariance structure to previously described long-time scale variation in an individual's behavior, suggesting that much of the measured variation between individuals of a single species in our assay reflects differences in the status of neural networks, rather than genetic or developmental differences between individuals. We then propose a method to identify groups of behaviors that appear to have evolved in a correlated manner, illustrating how sets of behaviors, rather than individual behaviors, likely evolved. Our approach provides a new framework for identifying co-evolving behaviors and may provide new opportunities to study the mechanistic basis of behavioral evolution.

The contribution of mutation and selection to multivariate quantitative genetic variance in an outbred population of Drosophila serrata.

Genetic variance is not equal for all multivariate combinations of traits. This inequality, in which some combinations of traits have abundant genetic variation while others have very little, biases the rate and direction of multivariate phenotypic evolution. However, we still understand little about what causes genetic variance to differ among trait combinations. Here, we investigate the relative roles of mutation and selection in determining the genetic variance of multivariate phenotypes. We accumulated mutations in an outbred population of Drosophila serrata and analyzed wing shape and size traits for over 35,000 flies to simultaneously estimate the additive genetic and additive mutational (co)variances. This experimental design allowed us to gain insight into the phenotypic effects of mutation as they arise and come under selection in naturally outbred populations. Multivariate phenotypes associated with more (less) genetic variance were also associated with more (less) mutational variance, suggesting that differences in mutational input contribute to differences in genetic variance. However, mutational correlations between traits were stronger than genetic correlations, and most mutational variance was associated with only one multivariate trait combination, while genetic variance was relatively more equal across multivariate traits. Therefore, selection is implicated in breaking down trait covariance and resulting in a different pattern of genetic variance among multivariate combinations of traits than that predicted by mutation and drift. Overall, while low mutational input might slow evolution of some multivariate phenotypes, stabilizing selection appears to reduce the strength of evolutionary bias introduced by pleiotropic mutation.

Integrative taxonomy and a new species description in the sturtevanti subgroup of the Drosophila saltans group (Diptera: Drosophilidae).

Although the biological concept of species is well established in animals, sometimes the decision about the specific status of a new species is difficult and hence requires support of an integrative analysis of several character sets. To date, the species Drosophila sturtevanti, D. magalhaesi, D. milleri and D. dacunhai, belonging to the sturtevanti subgroup of the Neotropical saltans species group, are identified mainly by the aedeagus morphology, but also present some differences in spot coloration and patterning of the female sixth tergite and in the shape and size of the spermathecae, parallel to a pattern of reproductive isolation. In the present study, we describe a novel saltans group species from French Guiana belonging to the sturtevanti subgroup. Our species designation is based on an integrative approach covering (i) aedeagi and spermathecae morphology by scanning electron microscopy, (ii) analysis of female sixth-tergite color, (iii) morphometrical analysis of aedeagi and wings, (iv) analysis of partial sequence of the COI, COII and ND4 mitochondrial genes as well as (v) intercrosses for analysis of reproductive isolation. The comparative analysis of the results on these markers with those of D. sturtevanti, D. milleri and D. dacunhai supports that this line belongs to a new species of the sturtevanti subgroup that we name Drosophila lehrmanae sp. nov. in honor of Prof. Lee Ehrman´s 85th birthday.

Karyotype Description of Two Andean Species of the guarani Group of Drosophila (Díptera: Drosophilidae) and Cytological Notes.

The guarani group of Drosophila genus (Diptera: Drosophilidae) is formed by 24 species however the relationship of these species is not clear. In the present study are described the karyotypes of Drosophila sachapuyu Peñafiel and Rafael, 2018 and Drosophila zamorana Peñafiel and Rafael, 2018, two Andean species members of the guarani group. Mitotic chromosomes from cerebral ganglia of third stand larval were obtained by thermal shock and cell suspension techniques. The karyotype of D. sachapuyu, presents 2n = 10 (4R, 1V; X = R, Y = R) while D. zamorana exhibits karyotype 2n = 12 (5R, 1V; X = V, Y = R).

Plastic Variation in the Phyletic Lineages of Cactophilic Drosophila meridionalis and Relation to Hosts as Potential for Diversification.

The insect/plant interaction is known to be a trigger for diversification and even speciation. Experimental analyses on fitness traits and phenotypic variation using alternative host sites have been performed to understand the process of diversification relative to insect/plant interactions. For cactophilic species of Drosophila, the speciation process is considered an adaptive radiation in response to the exploration of species of the Cactaceae as breeding and feeding sites. In this work, we analyzed life history and morphological traits in individuals from two phyletic lineages (Evolutionarily Significant Units ESU) of the cactophilic species Drosophila meridionalis (Wasserman 1962) (Diptera: Drosophilidae) raised from media prepare. The characters analyzed corresponded to viability, developmental time, and four morphological measurements. The experiments were performed in a semi-natural medium prepared with fermenting tissues of the natural hosts, Cereus hildmaniannus and Opuntia monacantha. Viability, development time, and three morphological measurements were influenced by lineage, suggesting differentiation between the lineages. However, in O. monacantha, the mean viability was greater (~15%) and development time was longer (~336 h) than in C. hildmaniannus (~11% and ~301 h, respectively). Only the developmental time was significantly affected by the host cactus. In general, ESU group A had better values than ESU group BC for the evaluated traits. This finding suggested differentiation between the two lineages and different plastic responsiveness to the contrasting environments of the hosts, and that C. hildmaniannus may be a relatively stressful environment for the larvae, as for other Drosophila species.

Allometry constrains the evolution of sexual dimorphism in Drosophila across 33 million years of divergence.

Sexual dimorphism is widely viewed as adaptive, reflecting the evolution of males and females toward divergent fitness optima. Its evolution, however, may often be constrained by the shared genetic architecture of the sexes, and by allometry. Here, we investigated the evolution of sexual size dimorphism, shape dimorphism, and their allometric relationship, in the wings of 82 taxa in the family Drosophilidae that have been diverging for at least 33 million years. Shape dimorphism among species was remarkably similar, with males characterized by longer, thinner wings than females. There was, however, quantitative variation among species in both size and shape dimorphism, with evidence that they have adapted to different evolutionary optima in different clades on timescales of about 10 million years. Within species, shape dimorphism was predicted by size, and among species, there was a strong relationship between size dimorphism and shape dimorphism. Allometry constrained the evolution of shape dimorphism for the two most variable traits we studied, but dimorphism was evolutionary labile in other traits. The keys for disentangling alternative explanations for dimorphism evolution are studies of natural and sexual selection, together with a deeper understanding of how microevolutionary parameters of evolvability relate to macroevolutionary patterns of divergence.

Topologically associating domains and their role in the evolution of genome structure and function in Drosophila.

Topologically associating domains (TADs) were recently identified as fundamental units of three-dimensional eukaryotic genomic organization, although our knowledge of the influence of TADs on genome evolution remains preliminary. To study the molecular evolution of TADs in Drosophila species, we constructed a new reference-grade genome assembly and accompanying high-resolution TAD map for D. pseudoobscura Comparison of D. pseudoobscura and D. melanogaster, which are separated by ∼49 million years of divergence, showed that ∼30%-40% of their genomes retain conserved TADs. Comparative genomic analysis of 17 Drosophila species revealed that chromosomal rearrangement breakpoints are enriched at TAD boundaries but depleted within TADs. Additionally, genes within conserved TADs show lower expression divergence than those located in nonconserved TADs. Furthermore, we found that a substantial proportion of long genes (>50 kbp) in D. melanogaster (42%) and D. pseudoobscura (26%) constitute their own TADs, implying transcript structure may be one of the deterministic factors for TAD formation. By using structural variants (SVs) identified from 14 D. melanogaster strains, its three closest sibling species from the D. simulans species complex, and two obscura clade species, we uncovered evidence of selection acting on SVs at TAD boundaries, but with the nature of selection differing between SV types. Deletions are depleted at TAD boundaries in both divergent and polymorphic SVs, suggesting purifying selection, whereas divergent tandem duplications are enriched at TAD boundaries relative to polymorphism, suggesting they are adaptive. Our findings highlight how important TADs are in shaping the acquisition and retention of structural mutations that fundamentally alter genome organization.

Interspecific introgression reveals a role of male genital morphology during the evolution of reproductive isolation in Drosophila.

Rapid divergence in genital structures among nascent species has been posited to be an early-evolving cause of reproductive isolation, although evidence supporting this idea as a widespread phenomenon remains mixed. Using a collection of interspecific introgression lines between two Drosophila species that diverged approximately 240,000 years ago, we tested the hypothesis that even modest divergence in genital morphology can result in substantial fitness losses. We studied the reproductive consequences of variation in the male epandrial posterior lobes between Drosophila mauritiana and Drosophila sechellia and found that divergence in posterior lobe morphology has significant fitness costs on several prefertilization and postcopulatory reproductive measures. Males with divergent posterior lobe morphology also significantly reduced the life span of their mates. Interestingly, one of the consequences of genital divergence was decreased oviposition and fertilization, which suggests that a sensory bias for posterior lobe morphology could exist in females, and thus, posterior lobe morphology may be the target of cryptic female choice in these species. Our results provide evidence that divergence in genitalia can in fact give rise to substantial reproductive isolation early during species divergence, and they also reveal novel reproductive functions of the external male genitalia in Drosophila.

A phylogeny for the Drosophila montium species group: A model clade for comparative analyses.

The Drosophila montium species group is a clade of 94 named species, closely related to the model species D. melanogaster. The montium species group is distributed over a broad geographic range throughout Asia, Africa, and Australasia. Species of this group possess a wide range of morphologies, mating behaviors, and endosymbiont associations, making this clade useful for comparative analyses. We use genomic data from 42 available species to estimate the phylogeny and relative divergence times within the montium species group, and its relative divergence time from D. melanogaster. To assess the robustness of our phylogenetic inferences, we use 3 non-overlapping sets of 20 single-copy coding sequences and analyze all 60 genes with both Bayesian and maximum likelihood methods. Our analyses support monophyly of the group. Apart from the uncertain placement of a single species, D. baimaii, our analyses also support the monophyly of all seven subgroups proposed within the montium group. Our phylograms and relative chronograms provide a highly resolved species tree, with discordance restricted to estimates of relatively short branches deep in the tree. In contrast, age estimates for the montium crown group, relative to its divergence from D. melanogaster, depend critically on prior assumptions concerning variation in rates of molecular evolution across branches, and hence have not been reliably determined. We discuss methodological issues that limit phylogenetic resolution - even when complete genome sequences are available - as well as the utility of the current phylogeny for understanding the evolutionary and biogeographic history of this clade.

PseudoBase: a genomic visualization and exploration resource for the Drosophila pseudoobscura subgroup.

Drosophila pseudoobscura is a classic model system for the study of evolutionary genetics and genomics. Given this long-standing interest, many genome sequences have accumulated for D. pseudoobscura and closely related species D. persimilis, D. miranda, and D. lowei. To facilitate the exploration of genetic variation within species and comparative genomics across species, we present PseudoBase, a database that couples extensive publicly available genomic data with simple visualization and query tools via an intuitive graphical interface, amenable for use in both research and educational settings. All genetic variation (SNPs and indels) within the database is derived from the same workflow, so variants are easily comparable across data sets. Features include an embedded JBrowse interface, ability to pull out alignments of individual genes/regions, and batch access for gene lists. Here, we introduce PseudoBase, and we demonstrate how this resource facilitates use of extensive genomic data from flies of the Drosophila pseudoobscura subgroup.

The effects of the sex chromosomes on the inheritance of species-specific traits of the copulatory organ shape in Drosophila virilis and Drosophila lummei.

The shape of the male genitalia in many taxa is the most rapidly evolving morphological structure, often driving reproductive isolation, and is therefore widely used in systematics as a key character to distinguish between sibling species. However, only a few studies have used the genital arch of the male copulatory organ as a model to study the genetic basis of species-specific differences in the Drosophila copulatory system. Moreover, almost nothing is known about the effects of the sex chromosomes on the shape of the male mating organ. In our study, we used a set of crosses between D. virilis and D. lummei and applied the methods of quantitative genetics to assess the variability of the shape of the male copulatory organ and the effects of the sex chromosomes and autosomes on its variance. Our results showed that the male genital shape depends on the species composition of the sex chromosomes and autosomes. Epistatic interactions of the sex chromosomes with autosomes and the species origin of the Y-chromosome in a male in interspecific crosses also influenced the expression of species-specific traits in the shape of the male copulatory system. Overall, the effects of sex chromosomes were comparable to the effects of autosomes despite the great differences in gene numbers between them. It may be reasonably considered that sexual selection for specific genes associated with the shape of the male mating organ prevents the demasculinization of the X chromosome.

The application of rapid evaporative ionization mass spectrometry in the analysis of Drosophila species-a potential new tool in entomology.

There is increasing emphasis on the use of new analytical approaches in subject analysis and classification, particularly in respect to minimal sample preparation. Here, we demonstrate that rapid evaporative ionization mass spectrometry (REIMS), a method that captures metabolite mass spectra after rapid combustive degradation of an intact biological specimen, generates informative mass spectra from several arthropods, and more specifically, is capable of discerning differences between species and sex of several adult Drosophila species. A model including five Drosophila species, built using pattern recognition, achieves high correct classification rates (over 90%) using test datasets and is able to resolve closely related species. The ease of discrimination of male and female specimens also demonstrates that sex-specific differences reside in the REIMS metabolite patterns, whether analysed across all five species or specifically for D. melanogaster. Further, the same approach can correctly discriminate and assign Drosophila species at the larval stage, where these are morphologically highly similar or identical. REIMS offers a novel approach to insect typing and analysis, requiring a few seconds of data acquisition per sample and has considerable potential as a new tool for the field biologist.

Reverse and forward engineering of Drosophila corneal nanocoatings.

Insect eyes have an anti-reflective coating, owing to nanostructures on the corneal surface creating a gradient of refractive index between that of air and that of the lens material1,2. These nanocoatings have also been shown to provide anti-adhesive functionality3. The morphology of corneal nanocoatings are very diverse in arthropods, with nipple-like structures that can be organized into arrays or fused into ridge-like structures4. This diversity can be attributed to a reaction-diffusion mechanism4 and patterning principles developed by Alan Turing5, which have applications in numerous biological settings6. The nanocoatings on insect corneas are one example of such Turing patterns, and the first known example of nanoscale Turing patterns4. Here we demonstrate a clear link between the morphology and function of the nanocoatings on Drosophila corneas. We find that nanocoatings that consist of individual protrusions have better anti-reflective properties, whereas partially merged structures have better anti-adhesion properties. We use biochemical analysis and genetic modification techniques to reverse engineer the protein Retinin and corneal waxes as the building blocks of the nanostructures. In the context of Turing patterns, these building blocks fulfil the roles of activator and inhibitor, respectively. We then establish low-cost production of Retinin, and mix this synthetic protein with waxes to forward engineer various artificial nanocoatings with insect-like morphology and anti-adhesive or anti-reflective function. Our combined reverse- and forward-engineering approach thus provides a way to economically produce functional nanostructured coatings from biodegradable materials.

Adaptation of gene loci to heterochromatin in the course of Drosophila evolution is associated with insulator proteins.

Pericentromeric heterochromatin is generally composed of repetitive DNA forming a transcriptionally repressive environment. Dozens of genes were embedded into pericentromeric heterochromatin during evolution of Drosophilidae lineage while retaining activity. However, factors that contribute to insusceptibility of gene loci to transcriptional silencing remain unknown. Here, we find that the promoter region of genes that can be embedded in both euchromatin and heterochromatin exhibits a conserved structure throughout the Drosophila phylogeny and carries motifs for binding of certain chromatin remodeling factors, including insulator proteins. Using ChIP-seq data, we demonstrate that evolutionary gene relocation between euchromatin and pericentric heterochromatin occurred with preservation of sites of insulation of BEAF-32 in evolutionarily distant species, i.e. D. melanogaster and D. virilis. Moreover, promoters of virtually all protein-coding genes located in heterochromatin in D. melanogaster are enriched with insulator proteins BEAF-32, GAF and dCTCF. Applying RNA-seq of a BEAF-32 mutant, we show that the impairment of BEAF-32 function has a complex effect on gene expression in D. melanogaster, affecting even those genes that lack BEAF-32 association in their promoters. We propose that conserved intrinsic properties of genes, such as sites of insulation near the promoter regions, may contribute to adaptation of genes to the heterochromatic environment and, hence, facilitate the evolutionary relocation of genes loci between euchromatin and heterochromatin.

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.

The evolution of sex-biased gene expression in the Drosophila brain.

Genes with sex-biased expression in Drosophila are thought to underlie sexually dimorphic phenotypes and have been shown to possess unique evolutionary properties. However, the forces and constraints governing the evolution of sex-biased genes in the somatic tissues of Drosophila are largely unknown. By using population-scale RNA sequencing data, we show that sex-biased genes in the Drosophila brain are highly enriched on the X Chromosome and that most are biased in a species-specific manner. We show that X-linked male-biased genes, and to a lesser extent female-biased genes, are enriched for signatures of directional selection at the gene expression level. By examining the evolutionary properties of gene-flanking regions on the X Chromosome, we find evidence that adaptive cis-regulatory changes are more likely to drive the expression evolution of X-linked male-biased genes than other X-linked genes. Finally, we examine whether constraint owing to broad expression across multiple tissues and genetic constraint owing to the largely shared male and female genomes could be responsible for the observed patterns of gene expression evolution. We find that expression breadth does not constrain the directional evolution of gene expression in the brain. Additionally, we find that the shared genome between males and females imposes a substantial constraint on the expression evolution of sex-biased genes. Overall, these results significantly advance our understanding of the patterns and forces shaping the evolution of sexual dimorphism in the Drosophila brain.

Post-transcriptional gene regulation regulates germline stem cell to oocyte transition during Drosophila oogenesis.

During oogenesis, several developmental processes must be traversed to ensure effective completion of gametogenesis including, stem cell maintenance and asymmetric division, differentiation, mitosis and meiosis, and production of maternally contributed mRNAs, making the germline a salient model for understanding how cell fate transitions are mediated. Due to silencing of the genome during meiotic divisions, there is little instructive transcription, barring a few examples, to mediate these critical transitions. In Drosophila, several layers of post-transcriptional regulation ensure that the mRNAs required for these processes are expressed in a timely manner and as needed during germline differentiation. These layers of regulation include alternative splicing, RNA modification, ribosome production, and translational repression. Many of the molecules and pathways involved in these regulatory activities are conserved from Drosophila to humans making the Drosophila germline an elegant model for studying the role of post-transcriptional regulation during stem cell differentiation and meiosis.