Molecular phylogeny and divergence times of drosophilid species

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Cis element length variability does not confer differential transcription factor occupancy at the D. melanogaster histone locus

Histone genes require precise regulation to maintain histone homeostasis and ensure nucleosome stoichiometry. Animal histone genes often have unique clustered genomic organization. However, there is variability of histone gene number and organization as well as differential regulation of the histone genes across species. The Drosophila melanogaster histone locus has unique organizational characteristics as it exists as a series of ∼100 highly regular, tandemly repeated arrays of the 5 replication-dependent histone genes at a single locus. Yet D. melanogaster are viable with only 12 transgenic histone gene arrays. We hypothesized that the histone genes across the locus are differentially regulated. We discovered that the GA-repeat within the H3/H4 promoter is the only variable sequence across the histone gene arrays. The H3/H4 promoter GA-repeat is targeted by CLAMP to promote histone gene expression. We also show two additional GA-binding transcription factors, GAGA Factor and Pipsqueak, target the GA-repeat. When we further examined CLAMP and GAF targeting, we determined that neither CLAMP nor GAF show bias for any GA-repeat lengths. Furthermore, we found that the distribution of GA-repeats targeted by both CLAMP and GAF do not change throughout early development. Together our results suggest that the transcription factors targeting the H3/H4 GA-repeat do not impact differential regulation of the histone genes, but indicate that future studies should interrogate additional cis elements or factors that impact histone gene regulation.

Ancestral neural circuits potentiate the origin of a female sexual behavior

Courtship interactions are remarkably diverse in form and complexity among species. How neural circuits evolve to encode new behaviors that are functionally integrated into these dynamic social interactions is unknown. Here we report a recently originated female sexual behavior in the island endemic Drosophila species D. santomea, where females signal receptivity to male courtship songs by spreading their wings, which in turn promotes prolonged songs in courting males. Copulation success depends on this female signal and correlates with males' ability to adjust his singing in such a social feedback loop. Functional comparison of sexual circuitry across species suggests that a pair of descending neurons, which integrates male song stimuli and female internal state to control a conserved female abdominal behavior, drives wing spreading in D. santomea. This co-option occurred through the refinement of a pre-existing, plastic circuit that can be optogenetically activated in an outgroup species. Combined, our results show that the ancestral potential of a socially-tuned key circuit node to engage the wing motor program facilitates the expression of a new female behavior in appropriate sensory and motivational contexts. More broadly, our work provides insights into the evolution of social behaviors, particularly female behaviors, and the underlying neural mechanisms.

Sperm Cyst "Looping": A Developmental Novelty Enabling Extreme Male Ornament Evolution

Postcopulatory sexual selection is credited as a principal force behind the rapid evolution of reproductive characters, often generating a pattern of correlated evolution between interacting, sex-specific traits. Because the female reproductive tract is the selective environment for sperm, one taxonomically widespread example of this pattern is the co-diversification of sperm length and female sperm-storage organ dimension. In Drosophila, having testes that are longer than the sperm they manufacture was believed to be a universal physiological constraint. Further, the energetic and time costs of developing long testes have been credited with underlying the steep evolutionary allometry of sperm length and constraining sperm length evolution in Drosophila. Here, we report on the discovery of a novel spermatogenic mechanism—sperm cyst looping—that enables males to produce relatively long sperm in short testis. This phenomenon (restricted to members of the saltans and willistoni species groups) begins early during spermatogenesis and is potentially attributable to heterochronic evolution, resulting in growth asynchrony between spermatid tails and the surrounding spermatid and somatic cyst cell membranes. By removing the allometric constraint on sperm length, this evolutionary innovation appears to have enabled males to evolve extremely long sperm for their body mass while evading delays in reproductive maturation time. On the other hand, sperm cyst looping was found to exact a cost by requiring greater total energetic investment in testes and a pronounced reduction in male lifespan. We speculate on the ecological selection pressures underlying the evolutionary origin and maintenance of this unique adaptation.

Between virus correlations in the outcome of infection across host species: Evidence of virus by host species interactions

Virus host shifts are a major source of outbreaks and emerging infectious diseases, and predicting the outcome of novel host and virus interactions remains a key challenge for virus research. The evolutionary relationships between host species can explain variation in transmission rates, virulence, and virus community composition between hosts, but it is unclear if correlations exist between related viruses in infection traits across novel hosts. Here, we measure correlations in viral load of four Cripavirus isolates across experimental infections of 45 Drosophilidae host species. We find positive correlations between every pair of viruses tested, suggesting that some host clades show broad susceptibility and could act as reservoirs and donors for certain types of viruses. Additionally, we find evidence of virus by host species interactions, highlighting the importance of both host and virus traits in determining the outcome of virus host shifts. Of the four viruses tested here, those that were more closely related tended to be more strongly correlated, providing tentative evidence that virus evolutionary relatedness may be a useful proxy for determining the likelihood of novel virus emergence, which warrants further research.

Interspecies variation of larval locomotion kinematics in the genus Drosophila and its relation to habitat temperature

BACKGROUND

Speed and trajectory of locomotion are the characteristic traits of individual species. Locomotion kinematics may have been shaped during evolution towards increased survival in the habitats of each species. Although kinematics of locomotion is thought to be influenced by habitats, the quantitative relation between the kinematics and environmental factors has not been fully revealed. Here, we performed comparative analyses of larval locomotion in 11 Drosophila species.

RESULTS

We found that larval locomotion kinematics are divergent among the species. The diversity is not correlated to the body length but is correlated instead to the habitat temperature of the species. Phylogenetic analyses using Bayesian inference suggest that the evolutionary rate of the kinematics is diverse among phylogenetic tree branches.

CONCLUSIONS

The results of this study imply that the kinematics of larval locomotion has diverged in the evolutionary history of the genus Drosophila and evolved under the effects of the ambient temperature of habitats.

Evidence of Adaptive Evolution in Wolbachia-Regulated Gene DNMT2 and Its Role in the Dipteran Immune Response and Pathogen Blocking

Eukaryotic nucleic acid methyltransferase (MTase) proteins are essential mediators of epigenetic and epitranscriptomic regulation. DNMT2 belongs to a large, conserved family of DNA MTases found in many organisms, including holometabolous insects such as fruit flies and mosquitoes, where it is the lone MTase. Interestingly, despite its nomenclature, DNMT2 is not a DNA MTase, but instead targets and methylates RNA species. A growing body of literature suggests that DNMT2 mediates the host immune response against a wide range of pathogens, including RNA viruses. Curiously, although DNMT2 is antiviral in Drosophila, its expression promotes virus replication in mosquito species. We, therefore, sought to understand the divergent regulation, function, and evolution of these orthologs. We describe the role of the Drosophila-specific host protein IPOD in regulating the expression and function of fruit fly DNMT2. Heterologous expression of these orthologs suggests that DNMT2's role as an antiviral is host-dependent, indicating a requirement for additional host-specific factors. Finally, we identify and describe potential evidence of positive selection at different times throughout DNMT2 evolution within dipteran insects. We identify specific codons within each ortholog that are under positive selection and find that they are restricted to four distinct protein domains, which likely influence substrate binding, target recognition, and adaptation of unique intermolecular interactions. Collectively, our findings highlight the evolution of DNMT2 in Dipteran insects and point to structural, regulatory, and functional differences between mosquito and fruit fly homologs.

Comparative transcriptomics between Drosophila mojavensis and D. arizonae reveals transgressive gene expression and underexpression of spermatogenesis-related genes in hybrid testes

Interspecific hybridization is a stressful condition that can lead to sterility and/or inviability through improper gene regulation in Drosophila species with a high divergence time. However, the extent of these abnormalities in hybrids of recently diverging species is not well known. Some studies have shown that in Drosophila, the mechanisms of postzygotic isolation may evolve more rapidly in males than in females and that the degree of viability and sterility is associated with the genetic distance between species. Here, we used transcriptomic comparisons between two Drosophila mojavensis subspecies and D. arizonae (repleta group, Drosophila) and identified greater differential gene expression in testes than in ovaries. We tested the hypothesis that the severity of the interspecies hybrid phenotype is associated with the degree of gene misregulation. We showed limited gene misregulation in fertile females and an increase in the amount of misregulation in males with more severe sterile phenotypes (motile vs. amotile sperm). In addition, for these hybrids, we identified candidate genes that were mostly associated with spermatogenesis dysfunction.

Stellate Genes and the piRNA Pathway in Speciation and Reproductive Isolation of Drosophila melanogaster

One of the main conditions of the species splitting from a common precursor lineage is the prevention of a gene flow between diverging populations. The study of Drosophila interspecific hybrids allows to reconstruct the speciation mechanisms and to identify hybrid incompatibility factors that maintain post-zygotic reproductive isolation between closely related species. The regulation, evolution, and maintenance of the testis-specific Ste-Su(Ste) genetic system in Drosophila melanogaster is the subject of investigation worldwide. X-linked tandem testis-specific Stellate genes encode proteins homologous to the regulatory β-subunit of protein kinase CK2, but they are permanently repressed in wild-type flies by the piRNA pathway via piRNAs originating from the homologous Y-linked Su(Ste) locus. Derepression of Stellate genes caused by Su(Ste) piRNA biogenesis disruption leads to the accumulation of crystalline aggregates in spermatocytes, meiotic defects and male sterility. In this review we summarize current data about the origin, organization, evolution of the Ste-Su(Ste) system, and piRNA-dependent regulation of Stellate expression. The Ste-Su(Ste) system is fixed only in the D. melanogaster genome. According to our hypothesis, the acquisition of the Ste-Su(Ste) system by a part of the ancient fly population appears to be the causative factor of hybrid sterility in crosses of female flies with males that do not carry Y-linked Su(Ste) repeats. To support this scenario, we have directly demonstrated Stellate derepression and the corresponding meiotic disorders in the testes of interspecies hybrids between D. melanogaster and D. mauritiana. This finding embraces our hypothesis about the contribution of the Ste-Su(Ste) system and the piRNA pathway to the emergence of reproductive isolation of D. melanogaster lineage from initial species.

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.

Evolution of sexual size dimorphism in the wing musculature of Drosophila

Male courtship songs in Drosophila are exceedingly diverse across species. While much of this variation is understood to have evolved from changes in the central nervous system, evolutionary transitions in the wing muscles that control the song may have also contributed to song diversity. Here, focusing on a group of four wing muscles that are known to influence courtship song in Drosophila melanogaster, we investigate the evolutionary history of wing muscle anatomy of males and females from 19 Drosophila species. We find that three of the wing muscles have evolved sexual dimorphisms in size multiple independent times, whereas one has remained monomorphic in the phylogeny. These data suggest that evolutionary changes in wing muscle anatomy may have contributed to species variation in sexually dimorphic wing-based behaviors, such as courtship song. Moreover, wing muscles appear to differ in their propensity to evolve size dimorphisms, which may reflect variation in the functional constraints acting upon different wing muscles.

Evolution of larval segment position across 12 Drosophila species

Many developmental traits that are critical to the survival of the organism are also robust. These robust traits are resistant to phenotypic change in the face of variation. This presents a challenge to evolution. In this article, we asked whether and how a well‐established robust trait, Drosophila segment patterning, changed over the evolutionary history of the genus. We compared segment position scaled to body length at the first‐instar larval stage among 12 Drosophila species. We found that relative segment position has changed many times across the phylogeny. Changes were frequent, but primarily small in magnitude. Phylogenetic analysis demonstrated that rates of change in segment position are variable along the Drosophila phylogenetic tree, and that these changes can occur in short evolutionary timescales. Correlation between position shifts of segments decreased as the distance between two segments increased, suggesting local control of segment position. The posterior‐most abdominal segment showed the highest magnitude of change on average, had the highest rate of evolution between species, and appeared to be evolving more independently as compared to the rest of the segments. This segment was exceptionally elongated in the cactophilic species in our dataset, raising questions as to whether this change may be adaptive.

Molecular Logic and Evolution of Bitter Taste in Drosophila

Taste systems detect a vast diversity of toxins, which are perceived as bitter. When a species adapts to a new environment, its taste system must adapt to detect new death threats. We deleted each of six commonly expressed bitter gustatory receptors (Grs) from Drosophila melanogaster. Systematic analysis revealed that requirements for these Grs differed for the same tastant in different neurons and for different tastants in the same neuron. Responses to some tastants in some neurons required four Grs, including Gr39a. Deletions also produced increased or novel responses, supporting a model of Gr-Gr inhibitory interactions. Coexpression of four Grs conferred several bitter responses to a sugar neuron. We then examined bitter coding in three other Drosophila species. We found major evolutionary shifts. One shift depended on the concerted activity of seven Grs. This work shows how the complex logic of bitter coding provides the capacity to detect innumerable hazards and the flexibility to adapt to new ones.

piRNA silencing contributes to interspecies hybrid sterility and reproductive isolation in Drosophila melanogaster

The piRNA pathway is an adaptive mechanism that maintains genome stability by repression of selfish genomic elements. In the male germline of Drosophila melanogaster repression of Stellate genes by piRNAs generated from Supressor of Stellate (Su(Ste)) locus is required for male fertility, but both Su(Ste) piRNAs and their targets are absent in other Drosophila species. We found that D. melanogaster genome contains multiple X-linked non-coding genomic repeats that have sequence similarity to the protein-coding host gene vasa. In the male germline, these vasa-related AT-chX repeats produce abundant piRNAs that are antisense to vasa; however, vasa mRNA escapes silencing due to imperfect complementarity to AT-chX piRNAs. Unexpectedly, we discovered AT-chX piRNAs target vasa of Drosophila mauritiana in the testes of interspecies hybrids. In the majority of hybrid flies, the testes were strongly reduced in size and germline content. A minority of hybrids maintained wild-type array of premeiotic germ cells in the testes, but in them harmful Stellate genes were derepressed due to the absence of Su(Ste) piRNAs, and meiotic failures were observed. Thus, the piRNA pathway contributes to reproductive isolation between D. melanogaster and closely related species, causing hybrid male sterility via misregulation of two different host protein factors.

High frequency of horizontal transfer in Jockey families (LINE order) of drosophilids

Background

The use of large-scale genomic analyses has resulted in an improvement of transposable element sampling and a significant increase in the number of reported HTT (horizontal transfer of transposable elements) events by expanding the sampling of transposable element sequences in general and of specific families of these elements in particular, which were previously poorly sampled. In this study, we investigated the occurrence of HTT events in a group of elements that, until recently, were uncommon among the HTT records in Drosophila – the Jockey elements, members of the LINE (long interspersed nuclear element) order of non-LTR (long terminal repeat) retrotransposons. The sequences of 111 Jockey families deposited in Repbase that met the criteria of the analysis were used to identify Jockey sequences in 48 genomes of Drosophilidae (genus Drosophila, subgenus Sophophora: melanogaster, obscura and willistoni groups; subgenus Drosophila: immigrans, melanica, repleta, robusta, virilis and grimshawi groups; subgenus Dorsilopha: busckii group; genus/subgenus Zaprionus and genus Scaptodrosophila).

Results

Phylogenetic analyses revealed 72 Jockey families in 41 genomes. Combined analyses revealed 15 potential HTT events between species belonging to different genera and species groups of Drosophilidae, providing evidence for the flow of genetic material favoured by the spatio-temporal sharing of these species present in the Palaeartic or Afrotropical region.

Conclusions

Our results provide phylogenetic, biogeographic and temporal evidence of horizontal transfers of the Jockey elements, increase the number of rare records of HTT in specific families of LINE elements, increase the number of known occurrences of these events, and enable a broad understanding of the evolutionary dynamics of these elements and the host species.

Thermosensitive alternative splicing senses and mediates temperature adaptation in Drosophila

Circadian rhythms are generated by the cyclic transcription, translation, and degradation of clock gene products, including timeless (tim), but how the circadian clock senses and adapts to temperature changes is not completely understood. Here, we show that temperature dramatically changes the splicing pattern of tim in Drosophila. We found that at 18°C, TIM levels are low because of the induction of two cold-specific isoforms: tim-cold and tim-short and cold. At 29°C, another isoform, tim-medium, is upregulated. Isoform switching regulates the levels and activity of TIM as each isoform has a specific function. We found that tim-short and cold encodes a protein that rescues the behavioral defects of tim⁰¹ mutants, and that flies in which tim-short and cold is abrogated have abnormal locomotor activity. In addition, miRNA-mediated control limits the expression of some of these isoforms. Finally, data that we obtained using minigenes suggest that tim alternative splicing might act as a thermometer for the circadian clock.

Hi-C guided assemblies reveal conserved regulatory topologies on X and autosomes despite extensive genome shuffling

In this study, Renschler et al. set out to analyze the impact of genomic rearrangements on genome topology using the Drosophila genus and X chromosome dosage compensation as a model. The authors developed a scaffolding algorithm and generated chromosome-length assemblies from Hi-C data for studying genome topology in three distantly related Drosophila species. Their data provides unique insights into genome topology evolution. RA

Genome rearrangements that occur during evolution impose major challenges on regulatory mechanisms that rely on three-dimensional genome architecture. Here, we developed a scaffolding algorithm and generated chromosome-length assemblies from Hi-C data for studying genome topology in three distantly related Drosophila species. We observe extensive genome shuffling between these species with one synteny breakpoint after approximately every six genes. A/B compartments, a set of large gene-dense topologically associating domains (TADs), and spatial contacts between high-affinity sites (HAS) located on the X chromosome are maintained over 40 million years, indicating architectural conservation at various hierarchies. Evolutionary conserved genes cluster in the vicinity of HAS, while HAS locations appear evolutionarily flexible, thus uncoupling functional requirement of dosage compensation from individual positions on the linear X chromosome. Therefore, 3D architecture is preserved even in scenarios of thousands of rearrangements highlighting its relevance for essential processes such as dosage compensation of the X chromosome.

Genome Size Evolution Differs Between Drosophila Subgenera with Striking Differences in Male and Female Genome Size in Sophophora

Genome size varies across the tree of life, with no clear correlation to organismal complexity or coding sequence, but with differences in non-coding regions. Phylogenetic methods have recently been incorporated to further disentangle this enigma, yet most of these studies have focused on widely diverged species. Few have compared patterns of genome size change in closely related species with known structural differences in the genome. As a consequence, the relationship between genome size and differences in chromosome number or inter-sexual differences attributed to XY systems are largely unstudied. We hypothesize that structural differences associated with chromosome number and X-Y chromosome differentiation, should result in differing rates and patterns of genome size change. In this study, we utilize the subgenera within the Drosophila to ask if patterns and rates of genome size change differ between closely related species with differences in chromosome numbers and states of the XY system. Genome sizes for males and females of 152 species are used to answer these questions (with 92 newly added or updated estimates). While we find no relationship between chromosome number and genome size or chromosome number and inter-sexual differences in genome size, we find evidence for differing patterns of genome size change between the subgenera, and increasing rates of change throughout time. Estimated shifts in rates of change in sex differences in genome size occur more often in Sophophora and correspond to known neo-sex events.

Quantitative Comparison of the Anterior-Posterior Patterning System in the Embryos of Five Drosophila Species

Complex spatiotemporal gene expression patterns direct the development of the fertilized egg into an adult animal. Comparisons across species show that, in spite of changes in the underlying regulatory DNA sequence, developmental programs can be maintained across millions of years of evolution. Reciprocally, changes in gene expression can be used to generate morphological novelty. Distinguishing between changes in regulatory DNA that lead to changes in gene expression and those that do not is therefore a central goal of evolutionary developmental biology. Quantitative, spatially-resolved measurements of developmental gene expression patterns play a crucial role in this goal, enabling the detection of subtle phenotypic differences between species and the development of computations models that link the sequence of regulatory DNA to expression patterns. Here we report the generation of two atlases of cellular resolution gene expression measurements for the primary anterior-posterior patterning genes in Drosophila simulans and Drosophila virilis By combining these data sets with existing atlases for three other Drosophila species, we detect subtle differences in the gene expression patterns and dynamics driving the highly conserved axis patterning system and delineate inter-species differences in the embryonic morphology. These data sets will be a resource for future modeling studies of the evolution of developmental gene regulatory networks.

Haploid selection drives new gene male germline expression

New genes are a major source of novelties, and a disproportionate amount of them are known to show testis expression in later phases of male gametogenesis in different groups such as mammals and plants. Here, we propose that this enhanced expression is a consequence of haploid selection during the latter stages of male gametogenesis. Because emerging adaptive mutations will be fixed faster if their phenotypes are expressed by haploid rather than diploid genotypes, new genes with advantageous functions arising during this unique stage of development have a better chance to become fixed. To test this hypothesis, expression levels of genes of differing evolutionary age were examined at various stages of Drosophila spermatogenesis. We found, consistent with a model based on haploid selection, that new Drosophila genes are both expressed in later haploid phases of spermatogenesis and harbor a significant enrichment of adaptive mutations. Additionally, the observed overexpression of new genes in the latter phases of spermatogenesis was limited to the autosomes. Because all male cells exhibit hemizygous expression for X-linked genes (and therefore effectively haploid), there is no expectation that selection acting on late spermatogenesis will have a different effect on X-linked genes in comparison to initial diploid phases. Together, our proposed hypothesis and the analyzed data suggest that natural selection in haploid cells elucidates several aspects of the origin of new genes by explaining the general prevalence of their testis expression, and a parsimonious solution for new alleles to avoid being lost by genetic drift or pseudogenization.

Hidden genomic evolution in a morphospecies-The landscape of rapidly evolving genes in Tetrahymena

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei—the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])—within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp–exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes.

Genomic comparison of ten morphologically very similar species of ciliate from the genus Tetrahymena reveals how parallel microevolutionary processes have shaped their genomes and created unique genes through retrotransposition.

Repeated evolution of asymmetric genitalia and right-sided mating behavior in the Drosophila nannoptera species group

BACKGROUND

Male genitals have repeatedly evolved left-right asymmetries, and the causes of such evolution remain unclear. The Drosophila nannoptera group contains four species, among which three exhibit left-right asymmetries of distinct genital organs. In the most studied species, Drosophila pachea, males display asymmetric genital lobes and they mate right-sided on top of the female. Copulation position of the other species is unknown.

RESULTS

To assess whether the evolution of genital asymmetry could be linked to the evolution of one-sided mating, we examined phallus morphology and copulation position in D. pachea and closely related species. The phallus was found to be symmetric in all investigated species except D. pachea, which displays an asymmetric phallus with a right-sided gonopore, and D. acanthoptera, which harbors an asymmetrically bent phallus. In all examined species, males were found to position themselves symmetrically on top of the female, except in D. pachea and D. nannoptera, where males mated right-sided, in distinctive, species-specific positions. In addition, the copulation duration was found to be increased in the nannoptera group species compared to closely related outgroup species.

CONCLUSION

Our study shows that gains, and possibly losses, of asymmetry in genital morphology and mating position have evolved repeatedly in the nannoptera group. Current data does not allow us to conclude whether genital asymmetry has evolved in response to changes in mating position, or vice versa.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Phylogeny of the Genus Drosophila

Understanding phylogenetic relationships among taxa is key to designing and implementing comparative analyses. The genus Drosophila, which contains over 1600 species, is one of the most important model systems in the biological sciences. For over a century, one species in this group, Drosophila melanogaster, has been key to studies of animal development and genetics, genome organization and evolution, and human disease. As whole-genome sequencing becomes more cost-effective, there is increasing interest in other members of this morphologically, ecologically, and behaviorally diverse genus. Phylogenetic relationships within Drosophila are complicated, and the goal of this paper is to provide a review of the recent taxonomic changes and phylogenetic relationships in this genus to aid in further comparative studies.

Helena and BS: Two Travellers between the Genera Drosophila and Zaprionus

The frequency of horizontal transfers of transposable elements (HTTs) varies among the types of elements according to the transposition mode and the geographical and temporal overlap of the species involved in the transfer. The drosophilid species of the genus Zaprionus and those of the melanogaster, obscura, repleta, and virilis groups of the genus Drosophila investigated in this study shared space and time at some point in their evolutionary history. This is particularly true of the subgenus Zaprionus and the melanogaster subgroup, which overlapped both geographically and temporally in Tropical Africa during their period of origin and diversification. Here, we tested the hypothesis that this overlap may have facilitated the transfer of retrotransposons without long terminal repeats (non-LTRs) between these species. We estimated the HTT frequency of the non-LTRs BS and Helena at the genome-wide scale by using a phylogenetic framework and a vertical and horizontal inheritance consistence analysis (VHICA). An excessively low synonymous divergence among distantly related species and incongruities between the transposable element and species phylogenies allowed us to propose at least four relatively recent HTT events of Helena and BS involving ancestors of the subgroup melanogaster and ancestors of the subgenus Zaprionus during their concomitant diversification in Tropical Africa, along with older possible events between species of the subgenera Drosophila and Sophophora. This study provides the first evidence for HTT of non-LTRs retrotransposons between Drosophila and Zaprionus, including an in-depth reconstruction of the time frame and geography of these events.

Comparative transcriptomics across 14 Drosophila species reveals signatures of longevity

Lifespan varies dramatically among species, but the biological basis is not well understood. Previous studies in model organisms revealed the importance of nutrient sensing, mTOR, NAD/sirtuins, and insulin/IGF1 signaling in lifespan control. By studying life-history traits and transcriptomes of 14 Drosophila species differing more than sixfold in lifespan, we explored expression divergence and identified genes and processes that correlate with longevity. These longevity signatures suggested that longer-lived flies upregulate fatty acid metabolism, downregulate neuronal system development and activin signaling, and alter dynamics of RNA splicing. Interestingly, these gene expression patterns resembled those of flies under dietary restriction and several other lifespan-extending interventions, although on the individual gene level, there was no significant overlap with genes previously reported to have lifespan-extension effects. We experimentally tested the lifespan regulation potential of several candidate genes and found no consistent effects, suggesting that individual genes generally do not explain the observed longevity patterns. Instead, it appears that lifespan regulation across species is modulated by complex relationships at the system level represented by global gene expression.

Diverse Cis-Regulatory Mechanisms Contribute to Expression Evolution of Tandem Gene Duplicates

Pairs of duplicated genes generally display a combination of conserved expression patterns inherited from their unduplicated ancestor and newly acquired domains. However, how the cis-regulatory architecture of duplicated loci evolves to produce these expression patterns is poorly understood. We have directly examined the gene-regulatory evolution of two tandem duplicates, the Drosophila Ly6 genes CG9336 and CG9338, which arose at the base of the drosophilids between 40 and 60 Ma. Comparing the expression patterns of the two paralogs in four Drosophila species with that of the unduplicated ortholog in the tephritid Ceratitis capitata, we show that they diverged from each other as well as from the unduplicated ortholog. Moreover, the expression divergence appears to have occurred close to the duplication event and also more recently in a lineage-specific manner. The comparison of the tissue-specific cis-regulatory modules (CRMs) controlling the paralog expression in the four Drosophila species indicates that diverse cis-regulatory mechanisms, including the novel tissue-specific enhancers, differential inactivation, and enhancer sharing, contributed to the expression evolution. Our analysis also reveals a surprisingly variable cis-regulatory architecture, in which the CRMs driving conserved expression domains change in number, location, and specificity. Altogether, this study provides a detailed historical account that uncovers a highly dynamic picture of how the paralog expression patterns and their underlying cis-regulatory landscape evolve. We argue that our findings will encourage studying cis-regulatory evolution at the whole-locus level to understand how interactions between enhancers and other regulatory levels shape the evolution of gene expression.

Deducing the presence of proteins and proteoforms in quantitative proteomics

The human genome harbors just 20,000 genes suggesting that the variety of possible protein products per gene plays a significant role in generating functional diversity. In bottom-up proteomics peptides are mapped back to proteins and proteoforms to describe a proteome; however, accurate quantitation of proteoforms is challenging due to incomplete protein sequence coverage and mapping ambiguities. Here, we demonstrate that a new software tool called ProteinClusterQuant (PCQ) can be used to deduce the presence of proteoforms that would have otherwise been missed, as exemplified in a proteomic comparison of two fly species, Drosophilamelanogaster and D. virilis. PCQ was used to identify reduced levels of serine/threonine protein kinases PKN1 and PKN4 in CFBE41o⁻ cells compared to HBE41o⁻ cells and to elucidate that shorter proteoforms of full-length caspase-4 and ephrin B receptor are differentially expressed. Thus, PCQ extends current analyses in quantitative proteomics and facilitates finding differentially regulated proteins and proteoforms.

An accurate quantitation of different proteoforms remains challenging due to incomplete protein sequence coverage in mass spectrometry datasets. Here the authors describe a method to facilitate the discovery of proteoforms that may otherwise not be considered due to incomplete protein coverage or ambiguities in mapping peptides back to proteins.

Evolution of DNMT2 in drosophilids: Evidence for positive and purifying selection and insights into new protein (pathways) interactions

The DNA methyltransferase 2 (DNMT2) protein is the most conserved member of the DNA methyltransferase family. Nevertheless, its substrate specificity is still controversial and elusive. The genomic role and determinants of DNA methylation are poorly understood in invertebrates, and several mechanisms and associations are suggested. In Drosophila, the only known DNMT gene is Dnmt2. Here we present our findings from a wide search for Dnmt2 homologs in 68 species of Drosophilidae. We investigated its molecular evolution, and in our phylogenetic analyses the main clades of Drosophilidae species were recovered. We tested whether the Dnmt2 has evolved neutrally or under positive selection along the subgenera Drosophila and Sophophora and investigated positive selection in relation to several physicochemical properties. Despite of a major selective constraint on Dnmt2, we detected six sites under positive selection. Regarding the DNMT2 protein, 12 sites under positive-destabilizing selection were found, which suggests a selection that favors structural and functional shifts in the protein. The search for new potential protein partners with DNMT2 revealed 15 proteins with high evolutionary rate covariation (ERC), indicating a plurality of DNMT2 functions in different pathways. These events might represent signs of molecular adaptation, with molecular peculiarities arising from the diversity of evolutionary histories experienced by drosophilids.

Interspecific studies of circadian genes period and timeless in Drosophila

The level of rescue of clock function in genetically arrhythmic Drosophila melanogaster hosts using interspecific clock gene transformation was used to study the putative intermolecular coevolution between interacting clock proteins. Among them PER and TIM are the two important negative regulators of the circadian clock feedback loop. We transformed either the D. pseudoobscura per or tim transgenes into the corresponding arrhythmic D. melanogaster mutant (per01 or tim01) and observed >50% rhythmicity but the period of activity rhythm was either longer (D. pseudoobscura-per) or shorter than 24 h (D. pseudoobscura-tim) compared to controls. By introducing both transgenes simultaneously into double mutants, we observed that the period of the activity rhythm was rescued by the pair of hemizygous transgenes (~24 h). These flies also showed a more optimal level of temperature compensation for the period. Under LD 12:12 these flies have a D. pseudoobscura like activity profile with the absence of morning anticipation as well as a very prominent earlier evening peak of activity rhythm. These observation are consistent with the view that TIM and PER form a heterospecific coevolved module at least for the circadian period of activity rhythms. However the strength of rhythmicity was reduced by having both transgenes present, so while evidence for a coevolution between PER and TIM is observed for some characters it is not for others.

Combining morphology and molecular data to improve Drosophila paulistorum (Diptera, Drosophilidae) taxonomic status

The willistoni species subgroup has been the subject of several studies since the latter half of the past century and is considered a Neotropical model for evolutionary studies, given the many levels of reproductive isolation and different evolutionary stages occurring within them. Here we present for the first time a phylogenetic reconstruction combining morphological characters and molecular data obtained from 8 gene fragments (COI, COII, Cytb, Adh, Ddc, Hb, kl-3 and per). Some relationships were incongruent when comparing morphological and molecular data. Also, morphological data presented some unresolved polytomies, which could reflect the very recent divergence of the subgroup. The total evidence phylogenetic reconstruction presented well-supported relationships and summarized the results of all analyses. The diversification of the willistoni subgroup began about 7.3 Ma with the split of D. insularis while D.paulistorum complex has a much more recent diversification history, which began about 2.1 Ma and apparently has not completed the speciation process, since the average time to sister species separation is one million years, and some entities of the D. paulistorum complex diverge between 0.3 and 1 Ma. Based on the obtained data, we propose the categorization of the former "semispecies" of D. paulistorum as a subspecies and describe the subspecies D. paulistorum amazonian, D. paulistorum andeanbrazilian, D. paulistorum centroamerican, D. paulistorum interior, D. paulistorum orinocan and D. paulistorum transitional.

Fine scale mapping of genomic introgressions within the Drosophila yakuba clade

The process of speciation involves populations diverging over time until they are genetically and reproductively isolated. Hybridization between nascent species was long thought to directly oppose speciation. However, the amount of interspecific genetic exchange (introgression) mediated by hybridization remains largely unknown, although recent progress in genome sequencing has made measuring introgression more tractable. A natural place to look for individuals with admixed ancestry (indicative of introgression) is in regions where species co-occur. In west Africa, D. santomea and D. yakuba hybridize on the island of São Tomé, while D. yakuba and D. teissieri hybridize on the nearby island of Bioko. In this report, we quantify the genomic extent of introgression between the three species of the Drosophila yakuba clade (D. yakuba, D. santomea), D. teissieri). We sequenced the genomes of 86 individuals from all three species. We also developed and applied a new statistical framework, using a hidden Markov approach, to identify introgression. We found that introgression has occurred between both species pairs but most introgressed segments are small (on the order of a few kilobases). After ruling out the retention of ancestral polymorphism as an explanation for these similar regions, we find that the sizes of introgressed haplotypes indicate that genetic exchange is not recent (>1,000 generations ago). We additionally show that in both cases, introgression was rarer on X chromosomes than on autosomes which is consistent with sex chromosomes playing a large role in reproductive isolation. Even though the two species pairs have stable contemporary hybrid zones, providing the opportunity for ongoing gene flow, our results indicate that genetic exchange between these species is currently rare.

Histone locus regulation by the Drosophila dosage compensation adaptor protein CLAMP

Rieder et al. report that conserved GA repeat cis elements within the bidirectional histone3–histone4 promoter direct histone locus body (HLB) formation in Drosophila. In addition, the CLAMP zinc finger protein binds these GA repeat motifs, increases chromatin accessibility, enhances histone gene transcription, and promotes HLB formation.

The conserved histone locus body (HLB) assembles prior to zygotic gene activation early during development and concentrates factors into a nuclear domain of coordinated histone gene regulation. Although HLBs form specifically at replication-dependent histone loci, the cis and trans factors that target HLB components to histone genes remained unknown. Here we report that conserved GA repeat cis elements within the bidirectional histone3–histone4 promoter direct HLB formation in Drosophila. In addition, the CLAMP (chromatin-linked adaptor for male-specific lethal [MSL] proteins) zinc finger protein binds these GA repeat motifs, increases chromatin accessibility, enhances histone gene transcription, and promotes HLB formation. We demonstrated previously that CLAMP also promotes the formation of another domain of coordinated gene regulation: the dosage-compensated male X chromosome. Therefore, CLAMP binding to GA repeat motifs promotes the formation of two distinct domains of coordinated gene activation located at different places in the genome.

Patterns of transcriptional parallelism and variation in the developing olfactory system of Drosophila species

Organisms have evolved strikingly parallel phenotypes in response to similar selection pressures suggesting that there may be shared constraints limiting the possible evolutionary trajectories. For example, the behavioral adaptation of specialist Drosophila species to specific host plants can exhibit parallel changes in their adult olfactory neuroanatomy. We investigated the genetic basis of these parallel changes by comparing gene expression during the development of the olfactory system of two specialist Drosophila species to that of four other generalist species. Our results suggest that the parallelism observed in the adult olfactory neuroanatomy of ecological specialists extends more broadly to their developmental antennal expression profiles, and to the transcription factor combinations specifying olfactory receptor neuron (ORN) fates. Additionally, comparing general patterns of variation for the antennal transcriptional profiles in the adult and developing olfactory system of the six species suggest the possibility that specific, non-random components of the developmental programs underlying the Drosophila olfactory system harbor a disproportionate amount of interspecies variation. Further examination of these developmental components may be able to inform a deeper understanding of how traits evolve.

Mutation predicts 40 million years of fly wing evolution

Mutation enables evolution, but the idea that adaptation is also shaped by mutational variation is controversial. Simple evolutionary hypotheses predict such a relationship if the supply of mutations constrains evolution, but it is not clear that constraints exist, and, even if they do, they may be overcome by long-term natural selection. Quantification of the relationship between mutation and phenotypic divergence among species will help to resolve these issues. Here we use precise data on over 50,000 Drosophilid fly wings to demonstrate unexpectedly strong positive relationships between variation produced by mutation, standing genetic variation, and the rate of evolution over the last 40 million years. Our results are inconsistent with simple constraint hypotheses because the rate of evolution is very low relative to what both mutational and standing variation could allow. In principle, the constraint hypothesis could be rescued if the vast majority of mutations are so deleterious that they cannot contribute to evolution, but this also requires the implausible assumption that deleterious mutations have the same pattern of effects as potentially advantageous ones. Our evidence for a strong relationship between mutation and divergence in a slowly evolving structure challenges the existing models of mutation in evolution.

Genome Evolution in Three Species of Cactophilic Drosophila

We report genomes of two species of cactophilic Drosophila: Drosophila arizonae and D. navojoa These two are the closest relatives of D. mojavensis, forming the D. mojavensis cluster. D. mojavensis and D. arizonae diverged from D. navojoa ∼5.8 Mya, while the split between D. arizonae and D. mojavensis is more recent, at 1.5 Mya. Together the three genomes provide opportunities to examine genomic changes associated with speciation and host shifts in this ecologically defined group of flies. The three species are also separated by fixed inversion differences in three of their six chromosomes. While the levels of nucleotide divergence in the colinear chromosomes are significantly lower than in the inverted chromosomes, consistent with a past role of the inversions in preventing gene flow, the patterns differ among the inverted chromosomes when the locations of nucleotides inside or outside of the inversions are considered. For Muller element E, there is greater divergence external to the inversion breakpoints. For Muller A, the divergence is slightly higher inside the inversions, while for Muller B, the breakpoints and hence the difference in substitutions in relation to the inversions could not be determined. The differences among the inverted chromosomes, especially once the breakpoints are clearly established, could aid in dating the origins of the inversions.

Potential Direct Regulators of the Drosophila yellow Gene Identified by Yeast One-Hybrid and RNAi Screens

The regulation of gene expression controls development, and changes in this regulation often contribute to phenotypic evolution. Drosophila pigmentation is a model system for studying evolutionary changes in gene regulation, with differences in expression of pigmentation genes such as yellow that correlate with divergent pigment patterns among species shown to be caused by changes in cis- and trans-regulation. Currently, much more is known about the cis-regulatory component of divergent yellow expression than the trans-regulatory component, in part because very few trans-acting regulators of yellow expression have been identified. This study aims to improve our understanding of the trans-acting control of yellow expression by combining yeast-one-hybrid and RNAi screens for transcription factors binding to yellow cis-regulatory sequences and affecting abdominal pigmentation in adults, respectively. Of the 670 transcription factors included in the yeast-one-hybrid screen, 45 showed evidence of binding to one or more sequence fragments tested from the 5′ intergenic and intronic yellow sequences from D. melanogaster, D. pseudoobscura, and D. willistoni, suggesting that they might be direct regulators of yellow expression. Of the 670 transcription factors included in the yeast-one-hybrid screen, plus another TF previously shown to be genetically upstream of yellow, 125 were also tested using RNAi, and 32 showed altered abdominal pigmentation. Nine transcription factors were identified in both screens, including four nuclear receptors related to ecdysone signaling (Hr78, Hr38, Hr46, and Eip78C). This finding suggests that yellow expression might be directly controlled by nuclear receptors influenced by ecdysone during early pupal development when adult pigmentation is forming.

Repeated Duplication of Argonaute2 Is Associated with Strong Selection and Testis Specialization in Drosophila

Argonaute2 (Ago2) is a rapidly evolving nuclease in the Drosophila melanogaster RNA interference (RNAi) pathway that targets viruses and transposable elements in somatic tissues. Here we reconstruct the history of Ago2 duplications across the D. obscura group and use patterns of gene expression to infer new functional specialization. We show that some duplications are old, shared by the entire species group, and that losses may be common, including previously undetected losses in the lineage leading to D. pseudoobscura We find that while the original (syntenic) gene copy has generally retained the ancestral ubiquitous expression pattern, most of the novel Ago2 paralogs have independently specialized to testis-specific expression. Using population genetic analyses, we show that most testis-specific paralogs have significantly lower genetic diversity than the genome-wide average. This suggests recent positive selection in three different species, and model-based analyses provide strong evidence of recent hard selective sweeps in or near four of the six D. pseudoobscura Ago2 paralogs. We speculate that the repeated evolution of testis specificity in obscura group Ago2 genes, combined with their dynamic turnover and strong signatures of adaptive evolution, may be associated with highly derived roles in the suppression of transposable elements or meiotic drive. Our study highlights the lability of RNAi pathways, even within well-studied groups such as Drosophila, and suggests that strong selection may act quickly after duplication in RNAi pathways, potentially giving rise to new and unknown RNAi functions in nonmodel species.

Genomic Mining of Phylogenetically Informative Nuclear Markers in Bark and Ambrosia Beetles

Deep level insect relationships are generally difficult to resolve, especially within taxa of the most diverse and species rich holometabolous orders. In beetles, the major diversity occurs in the Phytophaga, including charismatic groups such as leaf beetles, longhorn beetles and weevils. Bark and ambrosia beetles are wood boring weevils that contribute 12 percent of the diversity encountered in Curculionidae, one of the largest families of beetles with more than 50000 described species. Phylogenetic resolution in groups of Cretaceous age has proven particularly difficult and requires large quantity of data. In this study, we investigated 100 nuclear genes in order to select a number of markers with low evolutionary rates and high phylogenetic signal. A PCR screening using degenerate primers was applied to 26 different weevil species. We obtained sequences from 57 of the 100 targeted genes. Sequences from each nuclear marker were aligned and examined for detecting multiple copies, pseudogenes and introns. Phylogenetic informativeness (PI) and the capacity for reconstruction of previously established phylogenetic relationships were used as proxies for selecting a subset of the 57 amplified genes. Finally, we selected 16 markers suitable for large-scale phylogenetics of Scolytinae and related weevil taxa.

Activation of Latent Courtship Circuitry in the Brain of Drosophila Females Induces Male-like Behaviors

Courtship in Drosophila melanogaster offers a powerful experimental paradigm for the study of innate sexually dimorphic behaviors [1, 2]. Fruit fly males exhibit an elaborate courtship display toward a potential mate [1, 2]. Females never actively court males, but their response to the male’s display determines whether mating will actually occur. Sex-specific behaviors are hardwired into the nervous system via the actions of the sex determination genes doublesex (dsx) and fruitless (fru) [1]. Activation of male-specific dsx/fru⁺ P1 neurons in the brain initiates the male’s courtship display [3, 4], suggesting that neurons unique to males trigger this sex-specific behavior. In females, dsx⁺ neurons play a pivotal role in sexual receptivity and post-mating behaviors [1, 2, 5, 6, 7, 8, 9]. Yet it is still unclear how dsx⁺ neurons and dimorphisms in these circuits give rise to the different behaviors displayed by males and females. Here, we manipulated the function of dsx⁺ neurons in the female brain to investigate higher-order neurons that drive female behaviors. Surprisingly, we found that activation of female dsx⁺ neurons in the brain induces females to behave like males by promoting male-typical courtship behaviors. Activated females display courtship toward conspecific males or females, as well other Drosophila species. We uncovered specific dsx⁺ neurons critical for driving male courtship and identified pheromones that trigger such behaviors in activated females. While male courtship behavior was thought to arise from male-specific central neurons, our study shows that the female brain is equipped with latent courtship circuitry capable of inducing this male-specific behavioral program.

•

Activation of brain dsx-pC1 neurons promote male-like courtship in females

•

Activated females court conspecific males and females and other Drosophila species

•

Methyl pheromones trigger male courtship behaviors in activated females

•

The female brain is equipped with latent circuitry underlying male-like behavior

Comparative Genomic Analyses Provide New Insights into the Evolutionary Dynamics of Heterochromatin in Drosophila

The term heterochromatin has been long considered synonymous with gene silencing, but it is now clear that the presence of transcribed genes embedded in pericentromeric heterochromatin is a conserved feature in the evolution of eukaryotic genomes. Several studies have addressed the epigenetic changes that enable the expression of genes in pericentric heterochromatin, yet little is known about the evolutionary processes through which this has occurred. By combining genome annotation analysis and high-resolution cytology, we have identified and mapped 53 orthologs of D. melanogaster heterochromatic genes in the genomes of two evolutionarily distant species, D. pseudoobscura and D. virilis. Our results show that the orthologs of the D. melanogaster heterochromatic genes are clustered at three main genomic regions in D. virilis and D. pseudoobscura. In D. virilis, the clusters lie in the middle of euchromatin, while those in D. pseudoobscura are located in the proximal portion of the chromosome arms. Some orthologs map to the corresponding Muller C element in D. pseudoobscura and D. virilis, while others localize on the Muller B element, suggesting that chromosomal rearrangements that have been instrumental in the fusion of two separate elements involved the progenitors of genes currently located in D. melanogaster heterochromatin. These results demonstrate an evolutionary repositioning of gene clusters from ancestral locations in euchromatin to the pericentromeric heterochromatin of descendent D. melanogaster chromosomes. Remarkably, in both D. virilis and D. pseudoobscura the gene clusters show a conserved association with the HP1a protein, one of the most highly evolutionarily conserved epigenetic marks. In light of these results, we suggest a new scenario whereby ancestral HP1-like proteins (and possibly other epigenetic marks) may have contributed to the evolutionary repositioning of gene clusters into heterochromatin.

Functional interplay between MSL1 and CDK7 controls RNA polymerase II Ser5 phosphorylation

Proper gene expression requires coordinated interplay among transcriptional coactivators, transcription factors and the general transcription machinery. We report here that MSL1, a central component of the dosage compensation complex in Drosophila melanogaster and Drosophila virilis, displays evolutionarily conserved sex-independent binding to promoters. Genetic and biochemical analyses reveal a functional interaction of MSL1 with CDK7, a subunit of the Cdk-activating kinase (CAK) complex of the general transcription factor TFIIH. Importantly, MSL1 depletion leads to decreased phosphorylation of Ser5 of RNA polymerase II. In addition, we demonstrate that MSL1 is a phosphoprotein, and transgenic flies expressing MSL1 phosphomutants show mislocalization of the histone acetyltransferase MOF and histone H4 K16 acetylation, thus ultimately causing male lethality due to a failure of dosage compensation. We propose that, by virtue of its interaction with components of the general transcription machinery, MSL1 exists in different phosphorylation states, thereby modulating transcription in flies.

Identification of Bari Transposons in 23 Sequenced Drosophila Genomes Reveals Novel Structural Variants, MITEs and Horizontal Transfer

Bari elements are members of the Tc1-mariner superfamily of DNA transposons, originally discovered in Drosophila melanogaster, and subsequently identified in silico in 11 sequenced Drosophila genomes and as experimentally isolated in four non-sequenced Drosophila species. Bari-like elements have been also studied for their mobility both in vivo and in vitro. We analyzed 23 Drosophila genomes and carried out a detailed characterization of the Bari elements identified, including those from the heterochromatic Bari1 cluster in D. melanogaster. We have annotated 401 copies of Bari elements classified either as putatively autonomous or inactive according to the structure of the terminal sequences and the presence of a complete transposase-coding region. Analyses of the integration sites revealed that Bari transposase prefers AT-rich sequences in which the TA target is cleaved and duplicated. Furthermore evaluation of transposon’s co-occurrence near the integration sites of Bari elements showed a non-random distribution of other transposable elements. We also unveil the existence of a putatively autonomous Bari1 variant characterized by two identical long Terminal Inverted Repeats, in D. rhopaloa. In addition, we detected MITEs related to Bari transposons in 9 species. Phylogenetic analyses based on transposase gene and the terminal sequences confirmed that Bari-like elements are distributed into three subfamilies. A few inconsistencies in Bari phylogenetic tree with respect to the Drosophila species tree could be explained by the occurrence of horizontal transfer events as also suggested by the results of dS analyses. This study further clarifies the Bari transposon’s evolutionary dynamics and increases our understanding on the Tc1-mariner elements’ biology.

Rapid evolutionary turnover underlies conserved lncRNA-genome interactions

Many long noncoding RNAs (lncRNAs) can regulate chromatin states, but the evolutionary origin and dynamics driving lncRNA-genome interactions are unclear. We adapted an integrative strategy that identifies lncRNA orthologs in different species despite limited sequence similarity, which is applicable to mammalian and insect lncRNAs. Analysis of the roX lncRNAs, which are essential for dosage compensation of the single X chromosome in Drosophila males, revealed 47 new roX orthologs in diverse Drosophilid species across ∼40 million years of evolution. Genetic rescue by roX orthologs and engineered synthetic lncRNAs showed that altering the number of focal, repetitive RNA structures determines roX ortholog function. Genomic occupancy maps of roX RNAs in four species revealed conserved targeting of X chromosome neighborhoods but rapid turnover of individual binding sites. Many new roX-binding sites evolved from DNA encoding a pre-existing RNA splicing signal, effectively linking dosage compensation to transcribed genes. Thus, dynamic change in lncRNAs and their genomic targets underlies conserved and essential lncRNA-genome interactions.

Hawaiian Drosophila genomes: size variation and evolutionary expansions

This paper reports genome sizes of one Hawaiian Scaptomyza and 16 endemic Hawaiian Drosophila species that include five members of the antopocerus species group, one member of the modified mouthpart group, and ten members of the picture wing clade. Genome size expansions have occurred independently multiple times among Hawaiian Drosophila lineages, and have resulted in an over 2.3-fold range of genome sizes among species, with the largest observed in Drosophila cyrtoloma (1C = 0.41 pg). We find evidence that these repeated genome size expansions were likely driven by the addition of significant amounts of heterochromatin and satellite DNA. For example, our data reveal that the addition of seven heterochromatic chromosome arms to the ancestral haploid karyotype, and a remarkable proportion of ~70 % satellite DNA, account for the greatly expanded size of the D. cyrtoloma genome. Moreover, the genomes of 13/17 Hawaiian picture wing species are composed of substantial proportions (22-70 %) of detectable satellites (all but one of which are AT-rich). Our results suggest that in this tightly knit group of recently evolved species, genomes have expanded, in large part, via evolutionary amplifications of satellite DNA sequences in centric and pericentric domains (especially of the X and dot chromosomes), which have resulted in longer acrocentric chromosomes or metacentrics with an added heterochromatic chromosome arm. We discuss possible evolutionary mechanisms that may have shaped these patterns, including rapid fixation of novel expanded genomes during founder-effect speciation.

Sex Bias and Maternal Contribution to Gene Expression Divergence in Drosophila Blastoderm Embryos

Early embryogenesis is a unique developmental stage where genetic control of development is handed off from mother to zygote. Yet the contribution of this transition to the evolution of gene expression is poorly understood. Here we study two aspects of gene expression specific to early embryogenesis in Drosophila: sex-biased gene expression prior to the onset of canonical X chromosomal dosage compensation, and the contribution of maternally supplied mRNAs. We sequenced mRNAs from individual unfertilized eggs and precisely staged and sexed blastoderm embryos, and compared levels between D. melanogaster, D. yakuba, D. pseudoobscura and D. virilis. First, we find that mRNA content is highly conserved for a given stage and that studies relying on pooled embryos likely systematically overstate the degree of gene expression divergence. Unlike studies done on larvae and adults where most species show a larger proportion of genes with male-biased expression, we find that transcripts in Drosophila embryos are largely female-biased in all species, likely due to incomplete dosage compensation prior to the activation of the canonical dosage compensation mechanism. The divergence of sex-biased gene expression across species is observed to be often due to lineage-specific decrease of expression; the most drastic example of which is the overall reduction of male expression from the neo-X chromosome in D. pseudoobscura, leading to a pervasive female-bias on this chromosome. We see no evidence for a faster evolution of expression on the X chromosome in embryos (no “faster-X” effect), unlike in adults, and contrary to a previous study on pooled non-sexed embryos. Finally, we find that most genes are conserved in regard to their maternal or zygotic origin of transcription, and present evidence that differences in maternal contribution to the blastoderm transcript pool may be due to species-specific divergence of transcript degradation rates.

The early embryonic stage of development is unique as this stage begins under the genetic control of the mother and transitions to genetic control of the zygote. In order to investigate how these developmental specificities have shaped gene expression evolution, we studied the contribution of maternal gene expression to the pool of transcripts in the early embryo and how gene expression differs in early female and male embryos. We find extensive female-bias in gene expression, something not observed in adults. This is largely due to the delay in the establishment of a mechanism that balances gene expression from the differing number of X chromosomes in females and males. The maternal contribution to the transcript pool in the embryo is largely conserved over species, with some exceptions. However, there may also be some species-specific differences in the degradation rate of maternal transcripts, an often-overlooked contributor to variation in transcript level.

cis-regulatory analysis of the Drosophila pdm locus reveals a diversity of neural enhancers

Background

One of the major challenges in developmental biology is to understand the regulatory events that generate neuronal diversity. During Drosophila embryonic neural lineage development, cellular temporal identity is established in part by a transcription factor (TF) regulatory network that mediates a cascade of cellular identity decisions. Two of the regulators essential to this network are the POU-domain TFs Nubbin and Pdm-2, encoded by adjacent genes collectively known as pdm. The focus of this study is the discovery and characterization of cis-regulatory DNA that governs their expression.

Results

Phylogenetic footprinting analysis of a 125 kb genomic region that spans the pdm locus identified 116 conserved sequence clusters. To determine which of these regions function as cis-regulatory enhancers that regulate the dynamics of pdm gene expression, we tested each for in vivo enhancer activity during embryonic development and postembryonic neurogenesis. Our screen revealed 77 unique enhancers positioned throughout the noncoding region of the pdm locus. Many of these activated neural-specific gene expression during different developmental stages and many drove expression in overlapping patterns. Sequence comparisons of functionally related enhancers that activate overlapping expression patterns revealed that they share conserved elements that can be predictive of enhancer behavior. To facilitate data accessibility, the results of our analysis are catalogued in cisPatterns, an online database of the structure and function of these and other Drosophila enhancers.

Conclusions

These studies reveal a diversity of modular enhancers that most likely regulate pdm gene expression during embryonic and adult development, highlighting a high level of temporal and spatial expression specificity. In addition, we discovered clusters of functionally related enhancers throughout the pdm locus. A subset of these enhancers share conserved elements including sequences that correspond to known TF DNA binding sites. Although comparative analysis of the nubbin and pdm-2 encoding sequences indicate that these two genes most likely arose from a duplication event, we found only partial evidence of sequence duplication between their enhancers, suggesting that after the putative duplication their cis-regulatory DNA diverged at a higher rate than their coding sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1897-2) contains supplementary material, which is available to authorized users.

Gonadal transcriptomics elucidate patterns of adaptive evolution within marine rockfishes (Sebastes)

Background

The genetic mechanisms of speciation and adaptation in the marine environment are not well understood. The rockfish genus Sebastes provides a unique model system for studying adaptive evolution because of the extensive diversity found within this group, which includes morphology, ecology, and a broad range of life spans. Examples of adaptive radiations within marine ecosystems are considered an anomaly due to the absence of geographical barriers and the presence of gene flow. Using marine rockfishes, we identified signatures of natural selection from transcriptomes developed from gonadal tissue of two rockfish species (Sebastes goodei and S. saxicola). We predicted orthologous transcript pairs, and estimated their distributions of nonsynonymous (Ka) and synonymous (Ks) substitution rates.

Results

We identified 144 genes out of 1079 orthologous pairs under positive selection, of which 11 are functionally annotated to reproduction based on gene ontologies (GOs). One orthologous pair of the zona pellucida gene family, which is known for its role in the selection of sperm by oocytes, out of ten was identified to be evolving under positive selection. In addition to our results in the protein coding-regions of transcripts, we found substitution rates in 3’ and 5’ UTRs to be significantly lower than Ks substitution rates implying negative selection in these regions.

Conclusions

We were able to identify a series of candidate genes that are useful for the assessment of the critical genes that diverged and are responsible for the radiation within this genus. Genes associated with longevity hold potential for understanding the molecular mechanisms that have contributed to the radiation within this genus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1870-0) contains supplementary material, which is available to authorized users.

Pheromones mediating copulation and attraction in Drosophila

Intraspecific olfactory signals known as pheromones play important roles in insect mating systems. In the model Drosophila melanogaster, a key part of the pheromone-detecting system has remained enigmatic through many years of research in terms of both its behavioral significance and its activating ligands. Here we show that Or47b-and Or88a-expressing olfactory sensory neurons (OSNs) detect the fly-produced odorants methyl laurate (ML), methyl myristate, and methyl palmitate. Fruitless (fru(M))-positive Or47b-expressing OSNs detect ML exclusively, and Or47b- and Or47b-expressing OSNs are required for optimal male copulation behavior. In addition, activation of Or47b-expressing OSNs in the male is sufficient to provide a competitive mating advantage. We further find that the vigorous male courtship displayed toward oenocyte-less flies is attributed to an oenocyte-independent sustained production of the Or47b ligand, ML. In addition, we reveal that Or88a-expressing OSNs respond to all three compounds, and that these neurons are necessary and sufficient for attraction behavior in both males and females. Beyond the OSN level, information regarding the three fly odorants is transferred from the antennal lobe to higher brain centers in two dedicated neural lines. Finally, we find that both Or47b- and Or88a-based systems and their ligands are remarkably conserved over a number of drosophilid species. Taken together, our results close a significant gap in the understanding of the olfactory background to Drosophila mating and attraction behavior; while reproductive isolation barriers between species are created mainly by species-specific signals, the mating enhancing signal in several Drosophila species is conserved.

Diversity of epithelial morphogenesis during eggshell formation in drosophilids

The eggshells of drosophilid species provide a powerful model for studying the origins of morphological diversity. The dorsal appendages, or respiratory filaments, of these eggshells display a remarkable interspecies variation in number and shape, and the epithelial patterning underlying the formation of these structures is an area of active research. To extend the analysis of dorsal appendage formation to include morphogenesis, we developed an improved 3D image reconstruction approach. This approach revealed considerable interspecies variation in the cell shape changes and neighbor exchanges underlying appendage formation. Specifically, although the appendage floor in Drosophila melanogaster is formed through spatially ordered neighbor exchanges, the same structure in Scaptodrosophila pattersoni is formed through extreme changes in cell shape, whereas Drosophila funebris appears to display a combination of both cellular mechanisms. Furthermore, localization patterns of Par3/Bazooka suggest a self-organized, cell polarity-based origin for the variability of appendage number in S. pattersoni. Our results suggest that species deploy different combinations of apically and basally driven mechanisms to convert a two-dimensional primordium into a three-dimensional structure, and provide new directions for exploring the molecular origins of interspecies morphological variation.