A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH).

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Identification of the sex chromosome system in a sand fly species, Lutzomyia longipalpis s.l.

In many animal species, sex determination is accomplished by heterogamety i.e., one of the sexes produces two types of gametes, which upon fertilization will direct the development toward males or females. Both male ("XY") and female ("ZW") heterogamety are known to occur and can be easily distinguished when the sex-chromosomes are morphologically different. However, this approach fails in cases of homomorphic sex chromosomes, such as the sand fly Lutzomyia longipalpis s.l. (Psychodidae, Diptera), which is the main vector of visceral leishmaniosis in Brazil. In order to identify the heterogametic sex in L. longipalpis s.l., we did a whole-genome sequencing of males and females separately and used the "Y chromosome Genome Scan" (YGS) method to find sex-specific sequences. Our results, which were confirmed by PCR, show that L. longipalpis s.l. has XY system. The YGS method can be especially useful in situations in which no morphological difference is observed in the sex-chromosomes or when fresh specimens are not readily available.

A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly.

BACKGROUND: The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics. RESULTS: Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS). CONCLUSIONS: Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.

Dinner and a show: The role of male copulatory courtship song and female blood-feeding in the reproductive success of Lutzomyia longipalpis from Lapinha, Brazil.

Lutzomyia longipalpis is the main vector of visceral Leishmaniasis in the Americas and is composed of a species complex. Males of this sand-fly produce acoustic signals during copulation and different patterns are observed among Brazilian populations. Such acoustic signals are commonly involved in species recognition. However, since the song is only produced during copulation it is not clear how it affects mating success or contributes to sexual isolation. Another aspect that may affect reproductive success is the presence of food. Since hematophagy is such an important aspect of L. longipalpis biology, we wanted to test if blood-feeding can influence the reproductive behaviour of this insect. We performed crossing experiments removing males' wings (silencing them) and playing back either the homo-specific or the hetero-specific song to either unfed or blood-fed females. Our results showed that both songs and blood-feeding affect insemination success, but not the frequency of copulation. In trials where females were not blood-fed song clearly affected insemination; males with wings, and males with homo-specific song playback had a higher insemination success than wingless males (no song) and trials with hetero-specific song. Blood-feeding females prior to the trials increased insemination in all groups including the control group which suggests that mating happens simultaneously with, or immediately after, the blood meal. Blood-fed females also seemed to discriminate less against the wrong song or the lack of song (wingless) one day after feeding, however trials with the correct song still had higher insemination rates. Altogether, our results show that both the male copulatory courtship songs and female blood-feeding are important for reproductive success and as such are important components of the sexual behaviour of L. longipalpis.

An Improved Genome Assembly for Drosophila navojoa, the Basal Species in the mojavensis Cluster.

Three North American cactophilic Drosophila species, D. mojavensis, D. arizonae, and D. navojoa, are of considerable evolutionary interest owing to the shift from breeding in Opuntia cacti to columnar species. The 3 species form the "mojavensis cluster" of Drosophila. The genome of D. mojavensis was sequenced in 2007 and the genomes of D. navojoa and D. arizonae were sequenced together in 2016 using the same technology (Illumina) and assembly software (AllPaths-LG). Yet, unfortunately, the D. navojoa genome was considerably more fragmented and incomplete than its sister species, rendering it less useful for evolutionary genetic studies. The D. navojoa read dataset does not fully meet the strict insert size required by the assembler used (AllPaths-LG) and this incompatibility might explain its assembly problems. Accordingly, when we re-assembled the genome of D. navojoa with the SPAdes assembler, which does not have the strict AllPaths-LG requirements, we obtained a substantial improvement in all quality indicators such as N50 (from 84 kb to 389 kb) and BUSCO coverage (from 77% to 97%). Here we share a new, improved reference assembly for D. navojoa genome, along with a RNAseq transcriptome. Given the basal relationship of the Opuntia breeding D. navojoa to the columnar breeding D. arizonae and D. mojavensis, the improved assembly and annotation will allow researchers to address a range of questions associated with the genomics of host shifts, chromosomal rearrangements and speciation in this group.

An investigation of Y chromosome incorporations in 400 species of Drosophila and related genera.

Y chromosomes are widely believed to evolve from a normal autosome through a process of massive gene loss (with preservation of some male genes), shaped by sex-antagonistic selection and complemented by occasional gains of male-related genes. The net result of these processes is a male-specialized chromosome. This might be expected to be an irreversible process, but it was found in 2005 that the Drosophila pseudoobscura Y chromosome was incorporated into an autosome. Y chromosome incorporations have important consequences: a formerly male-restricted chromosome reverts to autosomal inheritance, and the species may shift from an XY/XX to X0/XX sex-chromosome system. In order to assess the frequency and causes of this phenomenon we searched for Y chromosome incorporations in 400 species from Drosophila and related genera. We found one additional large scale event of Y chromosome incorporation, affecting the whole montium subgroup (40 species in our sample); overall 13% of the sampled species (52/400) have Y incorporations. While previous data indicated that after the Y incorporation the ancestral Y disappeared as a free chromosome, the much larger data set analyzed here indicates that a copy of the Y survived as a free chromosome both in montium and pseudoobscura species, and that the current Y of the pseudoobscura lineage results from a fusion between this free Y and the neoY. The 400 species sample also showed that the previously suggested causal connection between X-autosome fusions and Y incorporations is, at best, weak: the new case of Y incorporation (montium) does not have X-autosome fusion, whereas nine independent cases of X-autosome fusions were not followed by Y incorporations. Y incorporation is an underappreciated mechanism affecting Y chromosome evolution; our results show that at least in Drosophila it plays a relevant role and highlight the need of similar studies in other groups.

Two new species of Drosophila (Diptera, Drosophilidae) associated with inflorescences of Goeppertia monophylla (Marantaceae) in the city of São Paulo, state of São Paulo, Brazil

Abstract Two new Brazilian species of Drosophila (subgenus Drosophila) are described and illustrated: Drosophila asymmetrica sp. nov. and Drosophila peixotoi sp. nov. Both species were collected, and emerged, from inflorescences of Goeppertia monophylla (Marantaceae) in the urban Forest Reserve of the Instituto de Biociências da Universidade de São Paulo and their types will be deposited in the Museu de Zoologia da USP. The former species, which could not be assigned to any known group, has a conspicuously asymmetric aedeagus and a narrow oviscapt valve. The latter species belongs to the guarani group and is closely related to D. guaru, D. ornatifrons and D. subbadia, from which it can be distinguished by the presence of just one conspicuous large black spine at inner lower tip of cercus instead of two spines.

Cryptic diversity in an Atlantic Forest malaria vector from the mountains of South-East Brazil.

BACKGROUND: Anopheles (Kerteszia) cruzii is the primary vector of human and simian malarias in Brazilian regions covered by the Atlantic Rainforest. Previous studies found that An. cruzii presents high levels of behavioural, chromosomal and molecular polymorphisms, which led to the hypothesis that it may be a complex of cryptic species. Here, An. cruzii specimens were collected in five sites in South-East Brazil located at different altitudes on the inner and coastal slopes of two mountain ranges covered by Atlantic Rainforest, known as Serra do Mar and Serra da Mantiqueria. Partial sequences for two genes (Clock and cpr) were generated and compared with previously published sequences from Florianópolis (southern Brazil). Genetic diversity was analysed with estimates of population structure (F ST ) and haplotype phylogenetic trees in order to understand how many species of the complex may occur in this biome and how populations across the species distribution are related. RESULTS: The sequences from specimens collected at sites located on the lower coastal slopes of Serra do Mar (Guapimirim, Tinguá and Sana) clustered together in the phylogenetic analysis, while the major haplotypes from sites located on higher altitude and at the continental side of the same mountains (Bocaina) clustered with those from Serra da Mantiqueira (Itatiaia), an inner mountain range. These two An. cruzii lineages showed statistically significant genetic differentiation and fixed characters, and have high F ST values typical of between species comparisons. Finally, in Bocaina, where the two lineages occur in sympatry, we found deviations from Hardy-Weinberg equilibrium due to a deficit of heterozygotes, indicating partial reproductive isolation. These results strongly suggest that at least two distinct lineages of An. cruzii (provisorily named "Group 1" and "Group 2") occur in the mountains of South-East Brazil. CONCLUSIONS: At least two genetically distinct An. cruzii lineages occur in the Atlantic Forest covered mountains of South-East Brazil. The co-occurrence of distinct lineages of An. cruzii (possibly incipient species) in those mountains is an interesting biological phenomenon and may have important implications for malaria prevalence, Plasmodium transmission dynamics and control.

Improved assembly of noisy long reads by k-mer validation.

Genome assembly depends critically on read length. Two recent technologies, from Pacific Biosciences (PacBio) and Oxford Nanopore, produce read lengths >20 kb, which yield de novo genome assemblies with vastly greater contiguity than those based on Sanger, Illumina, or other technologies. However, the very high error rates of these two new technologies (∼15% per base) makes assembly imprecise at repeats longer than the read length and computationally expensive. Here we show that the contiguity and quality of the assembly of these noisy long reads can be significantly improved at a minimal cost, by leveraging on the low error rate and low cost of Illumina short reads. Namely, k-mers from the PacBio raw reads that are not present in Illumina reads (which account for ∼95% of the distinct k-mers) are deemed sequencing errors and ignored at the seed alignment step. By focusing on the ∼5% of k-mers that are error free, read overlap sensitivity is dramatically increased. Of equal importance, the validation procedure can be extended to exclude repetitive k-mers, which prevents read miscorrection at repeats and further improves the resulting assemblies. We tested the k-mer validation procedure using one long-read technology (PacBio) and one assembler (MHAP/Celera Assembler), but it is very likely to yield analogous improvements with alternative long-read technologies and assemblers, such as Oxford Nanopore and BLASR/DALIGNER/Falcon, respectively.

Birth of a new gene on the Y chromosome of Drosophila melanogaster.

Contrary to the pattern seen in mammalian sex chromosomes, where most Y-linked genes have X-linked homologs, the Drosophila X and Y chromosomes appear to be unrelated. Most of the Y-linked genes have autosomal paralogs, so autosome-to-Y transposition must be the main source of Drosophila Y-linked genes. Here we show how these genes were acquired. We found a previously unidentified gene (flagrante delicto Y, FDY) that originated from a recent duplication of the autosomal gene vig2 to the Y chromosome of Drosophila melanogaster. Four contiguous genes were duplicated along with vig2, but they became pseudogenes through the accumulation of deletions and transposable element insertions, whereas FDY remained functional, acquired testis-specific expression, and now accounts for ∼20% of the vig2-like mRNA in testis. FDY is absent in the closest relatives of D. melanogaster, and DNA sequence divergence indicates that the duplication to the Y chromosome occurred ∼2 million years ago. Thus, FDY provides a snapshot of the early stages of the establishment of a Y-linked gene and demonstrates how the Drosophila Y has been accumulating autosomal genes.

Positive and purifying selection on the Drosophila Y chromosome.

Y chromosomes, with their reduced effective population size, lack of recombination, and male-limited transmission, present a unique collection of constraints for the operation of natural selection. Male-limited transmission may greatly increase the efficacy of selection for male-beneficial mutations, but the reduced effective size also inflates the role of random genetic drift. Together, these defining features of the Y chromosome are expected to influence rates and patterns of molecular evolution on the Y as compared with X-linked or autosomal loci. Here, we use sequence data from 11 genes in 9 Drosophila species to gain insight into the efficacy of natural selection on the Drosophila Y relative to the rest of the genome. Drosophila is an ideal system for assessing the consequences of Y-linkage for molecular evolution in part because the gene content of Drosophila Y chromosomes is highly dynamic, with orthologous genes being Y-linked in some species whereas autosomal in others. Our results confirm the expectation that the efficacy of natural selection at weakly selected sites is reduced on the Y chromosome. In contrast, purifying selection on the Y chromosome for strongly deleterious mutations does not appear to be compromised. Finally, we find evidence of recurrent positive selection for 4 of the 11 genes studied here. Our results thus highlight the variable nature of the mode and impact of natural selection on the Drosophila Y chromosome.

Efficient identification of Y chromosome sequences in the human and Drosophila genomes.

Notwithstanding their biological importance, Y chromosomes remain poorly known in most species. A major obstacle to their study is the identification of Y chromosome sequences; due to its high content of repetitive DNA, in most genome projects, the Y chromosome sequence is fragmented into a large number of small, unmapped scaffolds. Identification of Y-linked genes among these fragments has yielded important insights about the origin and evolution of Y chromosomes, but the process is labor intensive, restricting studies to a small number of species. Apart from these fragmentary assemblies, in a few mammalian species, the euchromatic sequence of the Y is essentially complete, owing to painstaking BAC mapping and sequencing. Here we use female short-read sequencing and k-mer comparison to identify Y-linked sequences in two very different genomes, Drosophila virilis and human. Using this method, essentially all D. virilis scaffolds were unambiguously classified as Y-linked or not Y-linked. We found 800 new scaffolds (totaling 8.5 Mbp), and four new genes in the Y chromosome of D. virilis, including JYalpha, a gene involved in hybrid male sterility. Our results also strongly support the preponderance of gene gains over gene losses in the evolution of the Drosophila Y. In the intensively studied human genome, used here as a positive control, we recovered all previously known genes or gene families, plus a small amount (283 kb) of new, unfinished sequence. Hence, this method works in large and complex genomes and can be applied to any species with sex chromosomes.

Low conservation of gene content in the Drosophila Y chromosome.

Chromosomal organization is sufficiently evolutionarily stable that large syntenic blocks of genes can be recognized even between species as distantly related as mammals and puffer fish (450 million years (Myr) of divergence). In Diptera, the gene content of the X chromosome and the autosomes is well conserved: in Drosophila more than 95% of the genes have remained on the same chromosome arm in the 12 sequenced species (63 Myr of divergence, traversing 400 Myr of evolution), and the same linkage groups are clearly recognizable in mosquito genomes (260 Myr of divergence). Here we investigate the conservation of Y-linked gene content among the 12 sequenced Drosophila species. We found that only a quarter of the Drosophila melanogaster Y-linked genes (3 out of 12) are Y-linked in all sequenced species, and that most of them (7 out of 12) were acquired less than 63 Myr ago. Hence, whereas the organization of other Drosophila chromosomes traces back to the common ancestor with mosquitoes, the gene content of the D. melanogaster Y chromosome is much younger. Gene losses are known to have an important role in the evolution of Y chromosomes, and we indeed found two such cases. However, the rate of gene gain in the Drosophila Y chromosomes investigated is 10.9 times higher than the rate of gene loss (95% confidence interval: 2.3-52.5), indicating a clear tendency of the Y chromosomes to increase in gene content. In contrast with the mammalian Y chromosome, gene gains have a prominent role in the evolution of the Drosophila Y chromosome.

Y chromosome of D. pseudoobscura is not homologous to the ancestral Drosophila Y.

We report a genome-wide search of Y-linked genes in Drosophila pseudoobscura. All six identifiable orthologs of the D. melanogaster Y-linked genes have autosomal inheritance in D. pseudoobscura. Four orthologs were investigated in detail and proved to be Y-linked in D. guanche and D. bifasciata, which shows that less than 18 million years ago the ancestral Drosophila Y chromosome was translocated to an autosome in the D. pseudoobscura lineage. We found 15 genes and pseudogenes in the current Y of D. pseudoobscura, and none are shared with the D. melanogaster Y. Hence, the Y chromosome in the D. pseudoobscura lineage appears to have arisen de novo and is not homologous to the D. melanogaster Y.

Evolution of autosomal suppression of the sex-ratio trait in Drosophila.

The sex-ratio trait is the production of female-biased progenies due to X-linked meiotic drive in males of several Drosophila species. The driving X chromosome (called SR) is not fixed due to at least two stabilizing factors: natural selection (favoring ST, the nondriving standard X) and drive suppression by either Y-linked or autosomal genes. The evolution of autosomal suppression is explained by Fisher's principle, a mechanism of natural selection that leads to equal proportion of males and females in a sexually reproducing population. In fact, sex-ratio expression is partially suppressed by autosomal genes in at least three Drosophila species. The population genetics of this system is not completely understood. In this article we develop a mathematical model for the evolution of autosomal suppressors of SR (sup alleles) and show that: (i). an autosomal suppressor cannot invade when SR is very deleterious in males (c < (1)/(3), where c is the fitness of SR/Y males); (ii). "SR/ST, sup/+" polymorphisms occur when SR is partially deleterious ( approximately 0.3 < c < 1); while (iii). SR neutrality (c = 1) results in sup fixation and thus in total abolishment of drive. So, surprisingly, as long as there is any selection against SR/Y males, neutral autosomal suppressors will not be fixed. In that case, when a polymorphic equilibrium exists, the average female proportion in SR/Y males' progeny is given approximately by ac + 1 - a + a (2) c + 1 (2) + 1 - 4ac /4ac, where a is the fitness of SR/ST females.