Two male-killing Wolbachia from Drosophila birauraia that are closely related but distinct in genome structure.

Insects harbour diverse maternally inherited bacteria and viruses, some of which have evolved to kill the male progeny of their hosts (male killing: MK). The fly species Drosophila biauraria carries a maternally transmitted MK-inducing partiti-like virus, but it was unknown if it carries other MK-inducing endosymbionts. Here, we identified two male-killing Wolbachia strains (wBiau1 and wBiau2) from D. biauraria and compared their genomes to elucidate their evolutionary processes. The two strains were genetically closely related but had exceptionally different genome structures with considerable rearrangements compared with combinations of other Wolbachia strains. Despite substantial changes in the genome structure, the two Wolbachia strains did not experience gene losses that would disrupt the male-killing expression or persistence in the host population. The two Wolbachia-infected matrilines carried distinct mitochondrial haplotypes, suggesting that wBiau1 and wBiau2 have invaded D. biauraria independently and undergone considerable genome changes owing to unknown selective pressures in evolutionary history. This study demonstrated the presence of three male-killers from two distinct origins in one fly species and highlighted the diverse and rapid genome evolution of MK Wolbachia in the host.

Identification of key yeast species and microbe-microbe interactions impacting larval growth of Drosophila in the wild.

Microbiota consisting of various fungi and bacteria have a significant impact on the physiological functions of the host. However, it is unclear which species are essential to this impact and how they affect the host. This study analyzed and isolated microbes from natural food sources of Drosophila larvae, and investigated their functions. Hanseniaspora uvarum is the predominant yeast responsible for larval growth in the earlier stage of fermentation. As fermentation progresses, Acetobacter orientalis emerges as the key bacterium responsible for larval growth, although yeasts and lactic acid bacteria must coexist along with the bacterium to stabilize this host-bacterial association. By providing nutrients to the larvae in an accessible form, the microbiota contributes to the upregulation of various genes that function in larval cell growth and metabolism. Thus, this study elucidates the key microbial species that support animal growth under microbial transition.

A male-killing gene encoded by a symbiotic virus of Drosophila.

In most eukaryotes, biparentally inherited nuclear genomes and maternally inherited cytoplasmic genomes have different evolutionary interests. Strongly female-biased sex ratios that are repeatedly observed in various arthropods often result from the male-specific lethality (male-killing) induced by maternally inherited symbiotic bacteria such as Spiroplasma and Wolbachia. However, despite some plausible case reports wherein viruses are raised as male-killers, it is not well understood how viruses, having much smaller genomes than bacteria, are capable of inducing male-killing. Here we show that a maternally inherited double-stranded RNA (dsRNA) virus belonging to the family Partitiviridae (designated DbMKPV1) induces male-killing in Drosophila. DbMKPV1 localizes in the cytoplasm and possesses only four genes, i.e., one gene in each of the four genomic segments (dsRNA1-dsRNA4), in contrast to ca. 1000 or more genes possessed by Spiroplasma or Wolbachia. We also show that a protein (designated PVMKp1; 330 amino acids in size), encoded by a gene on the dsRNA4 segment, is necessary and sufficient for inducing male-killing. Our results imply that male-killing genes can be easily acquired by symbiotic viruses through reassortment and that symbiotic viruses are hidden players in arthropod evolution. We anticipate that host-manipulating genes possessed by symbiotic viruses can be utilized for controlling arthropods.

Highly contiguous assemblies of 101 drosophilid genomes.

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.

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.

Divergence of Drosophila species: Longevity and reproduction under different nutrient balances.

How nutrition impacts growth, reproduction and longevity is complex because relationships between these life events are difficult to disentangle. As a first step in sorting out these processes, we carried out a comparative analysis of related species of Drosophila with distinct feeding habits. In particular, we examined life spans and egg laying of two generalists and three specialists on diets with distinct protein-to-carbohydrate ratios. In contrast to the generalist D. melanogaster, adult males of two specialists, D. sechellia and D. elegans, lived longer on a protein-rich diet. These results and our previous studies collectively show that the diet to which larvae of each specialist species have adapted ensures a longer life span of adult males of that same species. We also found a species-specific sexual dimorphism of life span in the above two specialists regardless of the diets, which was in sharp contrast to D. melanogaster. In D. melanogaster, males lived longer than females, whereas females of D. sechellia and D. elegans were longer-lived than males, and those specialist females were exceedingly low in egg production, relative to the other species. We discuss our findings from perspectives of mechanisms, including a possible contribution of egg production to life span.

Interspecies Comparative Analyses Reveal Distinct Carbohydrate-Responsive Systems among Drosophila Species.

During evolution, organisms have acquired variable feeding habits. Some species are nutritional generalists that adapt to various food resources, while others are specialists, feeding on specific resources. However, much remains to be discovered about how generalists adapt to diversified diets. We find that larvae of the generalists Drosophila melanogaster and D. simulans develop on three diets with different nutrient balances, whereas specialists D. sechellia and D. elegans cannot develop on carbohydrate-rich diets. The generalist D. melanogaster downregulates the expression of diverse metabolic genes systemically by transforming growth factor ß (TGF-ß)/Activin signaling, maintains metabolic homeostasis, and successfully adapts to the diets. In contrast, the specialist D. sechellia expresses those metabolic genes at higher levels and accumulates various metabolites on the carbohydrate-rich diet, culminating in reduced adaptation. Phenotypic similarities and differences strongly suggest that the robust carbohydrate-responsive regulatory systems are evolutionarily retained through genome-environment interactions in the generalists and contribute to their nutritional adaptabilities.

Evolution and diversity of the courtship repertoire in the Drosophila montium species group (Diptera: Drosophilidae).

Changes in elements of courtship behaviour can influence sexual isolation between species. Large-scale analyses of changes, including loss and gain of courtship elements, across a relatively complete phylogenetic group are rare but needed to understand the significance of such changes, for example whether the gain and loss of courtship elements are essentially arbitrary or equally reversible. In most species of Drosophila, courtship, including singing, mainly occurs before mounting as premounting courtship. The Drosophila montium species group is unusual because loss of premounting courtship and gain of post-mounting one has been detected in this group. Here, we provide an extensive analysis on the courtship repertoire and songs of 42 species in this group. Synchronously captured video and audio recordings were analysed to describe courtship patterns and male courtship songs, and changes were analysed in a phylogenetic context. Ancestral state reconstruction suggests that a gain of post-mounting courtship singing at the ancestor of this species group has been accompanied by a concurrent decrease in the incidence of premounting courtship singing and has led to subsequent further decrease and eventually complete loss of premounting courtship song in several lineages. Alongside this evolutionary trend towards post-mounting courtship, sine song and a special type of "high pulse repetition song" have become more widely used for courtship during species diversification in the montium group. It is likely that the elaboration of post-mounting courtship behaviours is associated with changes in the relative importance of pre- and post-mounting components of mating systems, such as sperm competition or cryptic female choice.

Maternally transmitted non-bacterial male killer in Drosophila biauraria.

A maternally inherited, all-female trait is widely found among arthropods, which is caused by bacterial endosymbionts such as Wolbachia, Rickettsia, Spiroplasma and Cardinium We discovered a single female of Drosophila biauraria, collected from Tomakomai, Hokkaido, Japan, that produced all-female offspring. This all-female trait was maternally inherited in the iso-female line (SP12F) by backcrossing with males of a normal line (SP11-20) with a 1 : 1 sex ratio derived from the same population. The all-female trait was not affected by tetracycline treatment performed for two consecutive generations. However, the microinjection of filter-sterilized homogenate of SP12F females into SP11-20 females established all-female matrilines. Our data suggest the role of transmissible agents, most likely viruses, but not bacteria or protists, as the possible cause of the all-female phenotype, which is likely to be achieved by killing of male embryos because egg hatch rates of SP12F were nearly half those of SP11-20. This is the first report in Diptera to demonstrate a maternally inherited virus-like element as the cause of the male-killing phenotype in D. biauraria.

Genetic Changes to a Transcriptional Silencer Element Confers Phenotypic Diversity within and between Drosophila Species.

The modification of transcriptional regulation has become increasingly appreciated as a major contributor to morphological evolution. However, the role of negative-acting control elements (e.g. silencers) in generating morphological diversity has been generally overlooked relative to positive-acting "enhancer" elements. The highly variable body coloration patterns among Drosophilid insects represents a powerful model system in which the molecular alterations that underlie phenotypic diversity can be defined. In a survey of pigment phenotypes among geographically disparate Japanese populations of Drosophila auraria, we discovered a remarkable degree of variation in male-specific abdominal coloration. In testing the expression patterns of the major pigment-producing enzymes, we found that phenotypes uniquely correlated with differences in the expression of ebony, a gene required for yellow-colored cuticle. Assays of ebony's transcriptional control region indicated that a lightly pigmented strain harbored cis-regulatory mutations that caused correlated changes in its expression. Through a series of chimeric reporter constructs between light and dark strain alleles, we localized function-altering mutations to a conserved silencer that mediates a male-specific pattern of ebony repression. This suggests that the light allele was derived through the loss of this silencer's activity. Furthermore, examination of the ebony gene of D. serrata, a close relative of D. auraria which secondarily lost male-specific pigmentation revealed the parallel loss of this silencer element. These results demonstrate how loss-of-function mutations in a silencer element resulted in increased gene expression. We propose that the mutational inactivation of silencer elements may represent a favored path to evolve gene expression, impacting morphological traits.

Drosophila suzukii: the genetic footprint of a recent, worldwide invasion.

Native to Asia, the soft-skinned fruit pest Drosophila suzukii has recently invaded the United States and Europe. The eastern United States represents the most recent expansion of their range, and presents an opportunity to test alternative models of colonization history. Here, we investigate the genetic population structure of this invasive fruit fly, with a focus on the eastern United States. We sequenced six X-linked gene fragments from 246 individuals collected from a total of 12 populations. We examine patterns of genetic diversity within and between populations and explore alternative colonization scenarios using approximate Bayesian computation. Our results indicate high levels of nucleotide diversity in this species and suggest that the recent invasions of Europe and the continental United States are independent demographic events. More broadly speaking, our results highlight the importance of integrating population structure into demographic models, particularly when attempting to reconstruct invasion histories. Finally, our simulation results illustrate the general challenge in reconstructing invasion histories using genetic data and suggest that genome-level data are often required to distinguish among alternative demographic scenarios.

Courtship patterns in the Drosophila montium species subgroup: repeated loss of precopulatory courtship?

During precopulatory courtship, male Drosophila typically produce wing vibration to generate species-specific songs before mounting females. Three species in the lini clade of the montium species subgroup have been found to produce species-specific sine song only after mounting and during copulation. Here we investigated and analyzed the courtship behavior of 29 species in the montium subgroup from video and song recordings and measured the duration of wing vibration. We describe a great diversity of courtship behavior in the montium subgroup. The courtship patterns can be categorized into four types in the montium subgroup: 1) type P/C, species with both precopulatory and copulatory courtship, such as D. parvula and D. nikananu, 2) type P-/C, species with sporadic precopulatory and mainly copulatory courtship, such as D. auraria and D. triauraria. 3) type C, species with only copulatory courtship, such as D. tani and D. pectinifera, 4) type C-, species with only very brief copulatory courtship, such as D. rufa and D. asahinai. According to a phylogenetic tree based on sequences of mitochondrial COI and COII, and the nuclear Adh, both precopulatory courtship and copulatory courtship were present in the most basal species D. parvula. Each of two branches in the montium subgroup contains four types of courtship behavior. Type C is present in each sub-branch. These results suggest that the courtship behavior initially involved both precopulatory and copulatory courtship, but that subsequently precopulatory courtship has gradually been lost in the montium subgroup. We suggest reasons why precopulatory behavior might come to be lost in the montium subgroup.

Negative effects of low temperatures on the vertical transmission and infection density of a spiroplasma endosymbiont in Drosophila hydei.

Maternally transmitted endosymbionts of the genus Spiroplasma infecting several species of Drosophila are known to cause selective death of male offspring (male killing). The male-killing trait is considered to be advantageous for maternally transmitted endosymbionts. However, a non-male-killing spiroplasma is present in Japanese populations of Drosophila hydei at high frequencies (23-66%). This spiroplasma is phylogenetically closely related to the male-killing spiroplasma infecting other Drosophila species. It is unknown why this spiroplasma is maintained in its host populations despite its inability to cause male killing. We examined the susceptibilities of the spiroplasma in D. hydei to four different temperatures (28, 25, 18, and 15 degrees C). Diagnostic PCR revealed that vertical transmission of the spiroplasma was nearly perfect at 28 and 25 degrees C, partially suppressed at 18 degrees C, and completely blocked at 15 degrees C. Furthermore, quantitative PCR demonstrated that offspring treated at 18 degrees C exhibited dramatically lower densities of spiroplasma (i.e., approximately one-tenth) compared to offspring treated at 28 and 25 degrees C. Considering the low temperatures during winter in Japan, some unknown advantageous effects of the spiroplasma that compensate for the failure of vertical transmission are suggested to act in natural populations of D. hydei.

Molecular phylogeny of the Drosophila auraria species complex and allied species of Japan based on nuclear and mitochondrial DNA sequences.

We studied the phylogenetic relationships among the six species of the Drosophila auraria species complex and three species of the montium species subgroup endemic to Japan. To do this, we examined the DNA sequences of the nuclear ITS1 and the mitochondrial COII of these species. Intraspecific variation in ITS1 was observed for most of the species studied. However, interspecific variation was small and of the same magnitude as intraspecific variation among D. auraria, D. triauraria, and D. quadraria. In addition, the amount of intraspecific variation between COII sequences from D. auraria and D. triauraria was almost the same as that of the interspecific variation between these species. Phylogenetic trees based on ITS1 or COII data reveal that, with the exception of D. pectinifera, these species can be separated into two taxonomic groups. D. rufa formed a taxonomic lineage with two Japanese endemic species, D. asahinai and D. lacteicornis. Physiological and behavioral differences between D. auraria and D. triauraria have been reported, but our results differed from those previously reported. Such differences may be better explained by hybridization between D. auraria and D. triauraria than by similarity due to recent common ancestry. Both the ITS1- and COII-based phylogenetic trees placed D. pectinifera outside the auraria complex. However, D. pectinifera sequences were more similar to those of the auraria complex than those of D. punjabiensis, another species thought to be closely related to members of the auraria complex.

Prevalence of a non-male-killing spiroplasma in natural populations of Drosophila hydei.

Male-killing phenotypes are found in a variety of insects and are often associated with maternally inherited endosymbiotic bacteria. In several species of Drosophila, male-killing endosymbionts of the genus Spiroplasma have been found at low frequencies (0.1 to 3%). In this study, spiroplasma infection without causing male-killing was shown to be prevalent (23 to 66%) in Japanese populations of Drosophila hydei. Molecular phylogenetic analyses showed that D. hydei was infected with a single strain of spiroplasma, which was closely related to male-killing spiroplasmas from other Drosophila species. Artificial-transfer experiments suggested that the spiroplasma genotype rather than the host genotype was responsible for the absence of the male-killing phenotype. Infection densities of the spiroplasma in the natural host, D. hydei, and in the artificial host, Drosophila melanogaster, were significantly lower than those of the male-killing spiroplasma NSRO, which was in accordance with the hypothesis that a threshold infection density is needed for the spiroplasma-induced male-killing expression.

Phenotypic stability of the P-M system in wild populations of Drosophila melanogaster.

The P element appears to be one of the most recently invaded transposons of D. melanogaster. To study the dynamics and long-term fate of P elements in natural populations of D. melanogaster, 472 isofemale lines newly collected from 27 localities of Japan were examined for the P element-associated characteristics (abilities to induce and repress of P element transposition) and genomic P element composition (size classes and their numbers). There was variation in the P element-related phenotypes among local populations, but genomic P composition did not correlate strongly with the phenotype of each line: full-size P and KP elements predominated in their genomes (FP+ KP predominance). Comparison with previous results suggests a stability in the P-M system in local populations over about 15 years. In some populations, phenotypic stability for particularly long times was found: for 30 years or more Q strains predominated in Hikone and Tanushimaru, P or Q strains around Inakadate, and M' or Q strains around Tozukawa. There was no clear evidence of structural destruction underlying functional variation of P elements during this period. These results suggest that the current evolutionary status of P elements in the gene pool of D. melanogaster is not intermediary stage predicted by the original recent invasion hypothesis, and that several other factors such as the position effect play important roles.

Inversion polymorphisms in populations of Drosophila melanogaster in the South-West islands of Japan: comparisons with the mainland populations.

Natural populations of D. melanogaster newly collected from the South-West Islands of Japan were examined for the frequencies of inversions, and the results, together with previous data, were compared with those of the mainland populations. The four Common cosmopolitan inversions (2Lt, 2RNS, 3LP, 3RP) were detected from most of the populations with higher frequencies. The three Rare cosmopolitan inversions (3LM, 3RC, 3RMo), three Quasi cosmopolitan inversions (2LA, 3LY, 3RK) and two new Endemic inversions were also found but with lower frequencies. Some of the South-West Island populations still maintained higher levels of inversion polymorphisms compared to the mainland populations, while others lost the polymorphisms to various degrees. The South-West Island populations were characterized by the higher frequency of In(3L)P. The ratios of In(3L)P /In(3R)P in the South-West Islands are significantly higher than those of the mainland populations, suggesting that the In(3L)P is relatively adaptive in the South-West Islands.