Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Phylogeny of Drosophila saltans group (Diptera: Drosophilidae) based on morphological and molecular evidence

Drosophila saltans group belongs to the subgenus Sophophora (family Drosophilidae), and it is subdivided into five subgroups, with 23 species. The species in this group are widely distributed in the Americas, primarily in the Neotropics. In the literature, the phylogenetic reconstruction of this group has been performed with various markers, but many inconsistencies remain. Here, we present a phylogenetic reconstruction of the saltans group with a greater number of species, 16 species, which is the most complete to date for the saltans group and includes all subgroups, in a combined analysis with morphological and molecular markers. We incorporated 48 morphological characters of male terminalia, the highest number used to date, and molecular markers based on mitochondrial genes COI and COII. Based on the results, which have recovered the five subgroups as distinct lineages, we propose a new hypothesis regarding the phylogenetic relationships among the subgroups of the saltans group. The relationships of the species within the sturtevanti and elliptica subgroups were well supported. The saltans subgroup showed several polytomies, but the relationship between the sibling species D. austrosaltans and D. saltans and their close relation with D. nigrosaltans were well supported in the molecular and total evidence analyses. The morphological analysis additionally supported the formation of the clade D. nigrosaltans—D. pseudosaltans. The observed polytomies may represent synchronous radiations or have resulted from speciation rates that have been too fast relative to the pace of substitution accumulation.

Evolution of genome structure in the Drosophila simulans species complex

The rapid evolution of repetitive DNA sequences, including satellite DNA, tandem duplications, and transposable elements, underlies phenotypic evolution and contributes to hybrid incompatibilities between species. However, repetitive genomic regions are fragmented and misassembled in most contemporary genome assemblies. We generated highly contiguous de novo reference genomes for the Drosophila simulans species complex (D. simulans, D. mauritiana, and D. sechellia), which speciated ∼250,000 yr ago. Our assemblies are comparable in contiguity and accuracy to the current D. melanogaster genome, allowing us to directly compare repetitive sequences between these four species. We find that at least 15% of the D. simulans complex species genomes fail to align uniquely to D. melanogaster owing to structural divergence-twice the number of single-nucleotide substitutions. We also find rapid turnover of satellite DNA and extensive structural divergence in heterochromatic regions, whereas the euchromatic gene content is mostly conserved. Despite the overall preservation of gene synteny, euchromatin in each species has been shaped by clade- and species-specific inversions, transposable elements, expansions and contractions of satellite and tRNA tandem arrays, and gene duplications. We also find rapid divergence among Y-linked genes, including copy number variation and recent gene duplications from autosomes. Our assemblies provide a valuable resource for studying genome evolution and its consequences for phenotypic evolution in these genetic model species.

Heterochromatin-Enriched Assemblies Reveal the Sequence and Organization of the Drosophila melanogaster Y Chromosome

Heterochromatic regions of the genome are repeat-rich and poor in protein coding genes, and are therefore underrepresented in even the best genome assemblies. One of the most difficult regions of the genome to assemble are sex-limited chromosomes. The Drosophila melanogaster Y chromosome is entirely heterochromatic, yet has wide-ranging effects on male fertility, fitness, and genome-wide gene expression. The genetic basis of this phenotypic variation is difficult to study, in part because we do not know the detailed organization of the Y chromosome. To study Y chromosome organization in D. melanogaster, we develop an assembly strategy involving the in silico enrichment of heterochromatic long single-molecule reads and use these reads to create targeted de novo assemblies of heterochromatic sequences. We assigned contigs to the Y chromosome using Illumina reads to identify male-specific sequences. Our pipeline extends the D. melanogaster reference genome by 11.9 Mb, closes 43.8% of the gaps, and improves overall contiguity. The addition of 10.6 MB of Y-linked sequence permitted us to study the organization of repeats and genes along the Y chromosome. We detected a high rate of duplication to the pericentric regions of the Y chromosome from other regions in the genome. Most of these duplicated genes exist in multiple copies. We detail the evolutionary history of one sex-linked gene family, crystal-Stellate While the Y chromosome does not undergo crossing over, we observed high gene conversion rates within and between members of the crystal-Stellate gene family, Su(Ste), and PCKR, compared to genome-wide estimates. Our results suggest that gene conversion and gene duplication play an important role in the evolution of Y-linked genes.

The Drosophila bipectinata species complex: phylogenetic relationship among different members based on chromosomal variations

Making interspecific hybridizations, where possible remains an unparalleled option for studying the intricacies of speciation. In the Drosophila bipectinata species complex comprising of four species, namely D. bipectinata, D. parabipectinata, D. malerkotliana and D. pseudoananassae, interspecific hybrids can be obtained in the laboratory, thus bequeathing an ideal opportunity for studying speciation and phylogeny. With the view of investigating the degree of divergence between each species pair, we planned to study the polytene chromosomes of the F₁ hybrids, as it would mirror the level of compatibility between the genomes of the parental species. Two sets of crosses were made, one involving homozygous strains of all four species from India and the other including homozygous strains from different places across the globe. Polytene chromosomes of F₁ larvae from both sets of crosses had similar configurations. In F₁ larvae from crosses involving D. bipectinata, D. parabipectinata and D. malerkotliana, complex configurations (depicting overlapping inversions) could be detected in different arms. However, they were fairly synapsed, indicating that the differences are only at the level of gene arrangements. The polytene chromosomes of larvae obtained by crossing D. pseudoananassae with the other three species were very thin with gross asynapsis in all the arms, demonstrating that the genome of D. pseudoananassae is widely diverged from rest of the species. The overlapping inversions (reflected in complex configuration), are inferred in the light of earlier chromosomal studies performed in this complex.

Interspecific hybridization does not affect the level of fluctuating asymmetry (FA) in the Drosophila bipectinata species complex

The Drosophila bipectinata species complex comprises of four very closely related species namely D. bipectinata, D. parabipectinata, D. malerkotliana and D. pseudoananassae. It was found that irrespective of the evolutionary divergence among the species, FA which is reflective of the developmental precision remains nearly same in four species. During the present study, the level of FA in different morphological traits was studied in interspecific hybrids and compared with that of parental species with the view that it would throw light on the degree of divergence between the parental species. If they have not diverged much, the interspecific hybrids may have a similar FA level, incompatibilities between their genomes being negligible. On the other hand, if there is substantial divergence, the level of FA may be higher due to incompatibility between the genomes of the parental species. The morphological traits taken were sternopleural bristle number and wing length in both males and females and ovariole number and sex-comb tooth number in females and males respectively. However, except in a few cases, we could not detect any significant differences in the level of FA in hybrids as compared to pure species. On the other hand, a number of abnormalities like poor viability, dystrophied ovaries, asymmetrical eyes etc., could be detected in hybrids from crosses involving D. pseudoananassae as one of the parents. Therefore, we conclude that specific developmental pathways are more susceptible to developmental disturbances due to genomic incompatibilities than the large complex system bringing about developmental stability.

Genome diversity and divergence in Drosophila mauritiana: multiple signatures of faster X evolution

Drosophila mauritiana is an Indian Ocean island endemic species that diverged from its two sister species, Drosophila simulans and Drosophila sechellia, approximately 240,000 years ago. Multiple forms of incomplete reproductive isolation have evolved among these species, including sexual, gametic, ecological, and intrinsic postzygotic barriers, with crosses among all three species conforming to Haldane’s rule: F₁ hybrid males are sterile and F₁ hybrid females are fertile. Extensive genetic resources and the fertility of hybrid females have made D. mauritiana, in particular, an important model for speciation genetics. Analyses between D. mauritiana and both of its siblings have shown that the X chromosome makes a disproportionate contribution to hybrid male sterility. But why the X plays a special role in the evolution of hybrid sterility in these, and other, species remains an unsolved problem. To complement functional genetic analyses, we have investigated the population genomics of D. mauritiana, giving special attention to differences between the X and the autosomes. We present a de novo genome assembly of D. mauritiana annotated with RNAseq data and a whole-genome analysis of polymorphism and divergence from ten individuals. Our analyses show that, relative to the autosomes, the X chromosome has reduced nucleotide diversity but elevated nucleotide divergence; an excess of recurrent adaptive evolution at its protein-coding genes; an excess of recent, strong selective sweeps; and a large excess of satellite DNA. Interestingly, one of two centimorgan-scale selective sweeps on the D. mauritiana X chromosome spans a region containing two sex-ratio meiotic drive elements and a high concentration of satellite DNA. Furthermore, genes with roles in reproduction and chromosome biology are enriched among genes that have histories of recurrent adaptive protein evolution. Together, these genome-wide analyses suggest that genetic conflict and frequent positive natural selection on the X chromosome have shaped the molecular evolutionary history of D. mauritiana, refining our understanding of the possible causes of the large X-effect in speciation.

Meta-analysis reveals that genes regulated by the Y chromosome in Drosophila melanogaster are preferentially localized to repressive chromatin

The Drosophila Y chromosome is a degenerated, heterochromatic chromosome with few functional genes. Despite this, natural variation on the Y chromosome in D. melanogaster has substantial trans-acting effects on the regulation of X-linked and autosomal genes. It is not clear, however, whether these genes simply represent a random subset of the genome or whether specific functional properties are associated with susceptibility to regulation by Y-linked variation. Here, we present a meta-analysis of four previously published microarray studies of Y-linked regulatory variation (YRV) in D. melanogaster. We show that YRV genes are far from a random subset of the genome: They are more likely to be in repressive chromatin contexts, be expressed tissue specifically, and vary in expression within and between species than non-YRV genes. Furthermore, YRV genes are more likely to be associated with the nuclear lamina than non-YRV genes and are generally more likely to be close to each other in the nucleus (although not along chromosomes). Taken together, these results suggest that variation on the Y chromosome plays a role in modifying how the genome is distributed across chromatin compartments, either via changes in the distribution of DNA-binding proteins or via changes in the spatial arrangement of the genome in the nucleus.

Surprising differences in the variability of Y chromosomes in African and cosmopolitan populations of Drosophila melanogaster

The nonrecombining Drosophila melanogaster Y chromosome is heterochromatic and has few genes. Despite these limitations, there remains ample opportunity for natural selection to act on the genes that are vital for male fertility and on Y factors that modulate gene expression elsewhere in the genome. Y chromosomes of many organisms have low levels of nucleotide variability, but a formal survey of D. melanogaster Y chromosome variation had yet to be performed. Here we surveyed Y-linked variation in six populations of D. melanogaster spread across the globe. We find surprisingly low levels of variability in African relative to Cosmopolitan (i.e., non-African) populations. While the low levels of Cosmopolitan Y chromosome polymorphism can be explained by the demographic histories of these populations, the staggeringly low polymorphism of African Y chromosomes cannot be explained by demographic history. An explanation that is entirely consistent with the data is that the Y chromosomes of Zimbabwe and Uganda populations have experienced recent selective sweeps. Interestingly, the Zimbabwe and Uganda Y chromosomes differ: in Zimbabwe, a European Y chromosome appears to have swept through the population.

Interspecific sexual isolation and phylogeny among different members of the Drosophila bipectinata species complex

The degree and pattern of sexual isolation among closely related species can be utilized for predicting the direction of evolution and deriving a phylogeny. The Drosophila bipectinata species complex is a group of four morphologically very similar species, belonging to the ananassae subgroup of the melanogaster species group. Sexual isolation among the members of this complex was studied by employing the male-choice technique. To test the difference between the homogamic and heterogamic matings, chi-square values were calculated under the assumption of random mating and the degree of sexual isolation was measured by estimating isolation index. To test the difference in proportions of heterogamic matings between reciprocal crosses, normal deviate (z) was calculated. The direction of evolution among these species has been discussed in the light of models proposed by different workers. Our results suggest that while D. bipectinata, Drosophila parabipectinata and Drosophila malerkotliana are closely related to each other, D. pseudoananassae is distantly related to these three species. The species pairs D. bipectinata and D. parabipectinata and D. parabipectinata and D. malerkotliana show asymmetrical sexual isolation. Based on this, it is suggested that D. bipectinata and D. malerkotliana share a common ancestor where as D. parabipectinata has been derived from D. bipectinata. Their phylogenetic relationship has been discussed in the light of phylogenies suggested by earlier workers.

Interspecific Y chromosome introgressions disrupt testis-specific gene expression and male reproductive phenotypes in Drosophila

The Drosophila Y chromosome is a degenerated, heterochromatic chromosome with few functional genes. Nonetheless, natural variation on the Y chromosome in Drosophila melanogaster has substantial trans-acting effects on the regulation of X-linked and autosomal genes. However, the contribution of Y chromosome divergence to gene expression divergence between species is unknown. In this study, we constructed a series of Y chromosome introgression lines, in which Y chromosomes from either Drosophila sechellia or Drosophila simulans are introgressed into a common D. simulans genetic background. Using these lines, we compared genome-wide gene expression and male reproductive phenotypes between heterospecific and conspecific Y chromosomes. We find significant differences in expression for 122 genes, or 2.84% of all genes analyzed. Genes down-regulated in males with heterospecific Y chromosomes are significantly biased toward testis-specific expression patterns. These same lines show reduced fecundity and sperm competitive ability. Taken together, these results imply a significant role for Y/X and Y/autosome interactions in maintaining proper expression of male-specific genes, either directly or via indirect effects on male reproductive tissue development or function.

A large palindrome with interchromosomal gene duplications in the pericentromeric region of the D. melanogaster Y chromosome

The non-recombining Y chromosome is expected to degenerate over evolutionary time, however, gene gain is a common feature of Y chromosomes of mammals and Drosophila. Here, we report that a large palindrome containing interchromosomal segmental duplications is located in the vicinity of the first amplicon detected in the Y chromosome of D. melanogaster. The recent appearance of such amplicons suggests that duplications to the Y chromosome, followed by the amplification of the segmental duplications, are a mechanism for the continuing evolution of Drosophila Y chromosomes.

Molecular evolution of a Y chromosome to autosome gene duplication in Drosophila

In contrast to the rest of the genome, the Y chromosome is restricted to males and lacks recombination. As a result, Y chromosomes are unable to respond efficiently to selection, and newly formed Y chromosomes degenerate until few genes remain. The rapid loss of genes from newly formed Y chromosomes has been well studied, but gene loss from highly degenerate Y chromosomes has only recently received attention. Here, we identify and characterize a Y to autosome duplication of the male fertility gene kl-5 that occurred during the evolution of the testacea group species of Drosophila. The duplication was likely DNA based, as other Y-linked genes remain on the Y chromosome, the locations of introns are conserved, and expression analyses suggest that regulatory elements remain linked. Genetic mapping reveals that the autosomal copy of kl-5 resides on the dot chromosome, a tiny autosome with strongly suppressed recombination. Molecular evolutionary analyses show that autosomal copies of kl-5 have reduced polymorphism and little recombination. Importantly, the rate of protein evolution of kl-5 has increased significantly in lineages where it is on the dot versus Y linked. Further analyses suggest this pattern is a consequence of relaxed purifying selection, rather than adaptive evolution. Thus, although the initial fixation of the kl-5 duplication may have been advantageous, slightly deleterious mutations have accumulated in the dot-linked copies of kl-5 faster than in the Y-linked copies. Because the dot chromosome contains seven times more genes than the Y and is exposed to selection in both males and females, these results suggest that the dot suffers the deleterious effects of genetic linkage to more selective targets compared with the Y chromosome. Thus, a highly degenerate Y chromosome may not be the worst environment in the genome, as is generally thought, but may in fact be protected from the accumulation of deleterious mutations relative to other nonrecombining regions that contain more genes.

Functional copies of the Mst77F gene on the Y chromosome of Drosophila melanogaster

The Y chromosome of Drosophila melanogaster has <20 protein-coding genes. These genes originated from the duplication of autosomal genes and have male-related functions. In 1993, Russell and Kaiser found three Y-linked pseudogenes of the Mst77F gene, which is a testis-expressed autosomal gene that is essential for male fertility. We did a thorough search using experimental and computational methods and found 18 Y-linked copies of this gene (named Mst77Y-1-Mst77Y-18). Ten Mst77Y genes encode defective proteins and the other eight are potentially functional. These eight genes produce approximately 20% of the functional Mst77F-like mRNA, and molecular evolutionary analysis shows that they evolved under purifying selection. Hence several Mst77Y genes have all the features of functional genes. Mst77Y genes are present only in D. melanogaster, and phylogenetic analysis confirmed that the duplication is a recent event. The identification of functional Mst77Y genes reinforces the previous finding that gene gains play a prominent role in the evolution of the Drosophila Y chromosome.

Phylogenetic relationships among cetaceans revealed by Y-chromosome sequences

The Y chromosome has recently come into the spotlight as a new and efficient genetic marker for tracing paternal lineages. We reconstructed cetacean phylogeny using a 1.7-kbp fragment of the non-recombining Y chromosome (NRY), including the SRY gene and a flanking non-coding region. The topology of the Y-chromosome tree is robust to various methods of analysis and exhibits high branch-support values, possibly due to the absence of recombination, small effective population size, and low homoplasy. The Y-chromosome tree indicates monophyly of each suborder, Mysticeti and Odontoceti, with high branch support values (BS> or =86%; PP> or =98%). In the Odontoceti clade, three superfamilies, Physeteroidea, Ziphioidea, and Delphinoidea, diverged soon after the split between Mysticeti and Odontoceti. Our analysis allows resolution of this rapid radiation and indicates that Physeteroidea is basal in the Odontoceti clade (BS, 99%; PP, 100%; MBS, 61%). The major split within the superfamily Delphinoidea is between the Delphinidae clade and the Monodontidae+ Phocoenidae clade. The phylogenetic relationships among delphinid species are ambiguous, probably because of the rapid radiation of this family. In the Mysticeti clade, the first major split is between Balaenidae and Balaenopteridae; within Balaenopteridae, a Balaenoptera acutorostrata+B. bonaerensis (minke whales) clade forms a sister clade with the other balaenopterid species. Megaptera novaeangliae is nested within Balaenoptera, making the latter paraphyletic. The low homoplasy exhibited by the Y-chromosome data presented here suggests that an extended data set incorporating longer sequences would provide better resolution of cetacean lower-level pylogeny.