Intra- and interspecific divergence in the nuclear sequences of the clock gene period in species of the Drosophila buzzatii cluster

Structure, Organization, and Evolution of Satellite DNAs: Insights from the Drosophila repleta and D. virilis Species Groups

The fact that satellite DNAs (satDNAs) in eukaryotes are abundant genomic components, can perform functional roles, but can also change rapidly across species while being homogenous within a species, makes them an intriguing and fascinating genomic component to study. It is also becoming clear that satDNAs represent an important piece in genome architecture and that changes in their structure, organization, and abundance can affect the evolution of genomes and species in many ways. Since the discovery of satDNAs more than 50 years ago, species from the Drosophila genus have continuously been used as models to study several aspects of satDNA biology. These studies have been largely concentrated in D. melanogaster and closely related species from the Sophophora subgenus, even though the vast majority of all Drosophila species belong to the Drosophila subgenus. This chapter highlights some studies on the satDNA structure, organization, and evolution in two species groups from the Drosophila subgenus: the repleta and virilis groups. We also discuss and review the classification of other abundant tandem repeats found in these species in the light of the current information available.

Genetics and Distribution Modeling: The Demographic History of the Cactophilic Drosophila buzzatii Species Cluster in Open Areas of South America

Biodiversity is the result of historical and recurrent events acting on populations and species. The Drosophila buzzatii species cluster is distributed along a diagonal of open areas in South America. Combining genetic analyses with species distribution models we evaluated the influence of climatic changes in the demography history of this cluster. We performed a Bayesian Skyline analysis and reconstructed the ancestral areas based on mitochondrial cytochrome oxidase I (mtCOI) gene. We modeled the current distribution and projected it to past (mid-Holocene and Last Glacial Maximum) and future. Our results demonstrate that climate change plays a critical role in historical demography and in defining the current and future geographic ranges of these species. The inter-Andean dry valleys and the Chiquitano Seasonally Dry Tropical Forests (SDTF) in Bolivia are considered the ancestral area for the D. buzzatii cluster. From this area, the migration route was through a west-east corridor from central Andes, throughout Bolivia, Paraguay and Argentina toward eastern and northeastern Brazil, along the edges of rain forest. The responses of the species to the climatic changes differ from the dominant assumptions of expansion during dry/cold weather and contraction during wet/warm weather that characterized the glacial cycles. We suggest that the influence of ecology and the potential responsiveness of each taxon to the environmental dynamics should be considered as well. Predictions for the future suggest a decline of suitable areas for the cluster, threatening biodiversity of these habitats. This work showed the importance of an integrative analysis of genetics and geography information to improve the inferences about demographic history hypotheses for the cluster.

What does mitogenomics tell us about the evolutionary history of the Drosophila buzzatii cluster (repleta group)?

The Drosophila repleta group is an array of more than 100 species endemic to the “New World”, many of which are cactophilic. The ability to utilize decaying cactus tissues as breeding and feeding sites is a key aspect that allowed the successful diversification of the repleta group in American deserts and arid lands. Within this group, the Drosophila buzzatii cluster is a South American clade of seven closely related species in different stages of divergence, making them a valuable model system for evolutionary research. Substantial effort has been devoted to elucidating the phylogenetic relationships among members of the D. buzzatii cluster, including molecular phylogenetic studies that have generated ambiguous results where different tree topologies have resulted dependent on the kinds of molecular marker used. Even though mitochondrial DNA regions have become useful markers in evolutionary biology and population genetics, none of the more than twenty Drosophila mitogenomes assembled so far includes this cluster. Here, we report the assembly of six complete mitogenomes of five species: D. antonietae, D. borborema, D. buzzatii, two strains of D. koepferae and D. seriema, with the aim of revisiting phylogenetic relationships and divergence times by means of mitogenomic analyses. Our recovered topology using complete mitogenomes supports the hypothesis of monophyly of the D. buzzatii cluster and shows two main clades, one including D. buzzatii and D. koepferae (both strains), and the other containing the remaining species. These results are in agreement with previous reports based on a few mitochondrial and/or nuclear genes, but conflict with the results of a recent large-scale nuclear phylogeny, indicating that nuclear and mitochondrial genomes depict different evolutionary histories.

Transposable Element Misregulation Is Linked to the Divergence between Parental piRNA Pathways in Drosophila Hybrids

Interspecific hybridization is a genomic stress condition that leads to the activation of transposable elements (TEs) in both animals and plants. In hybrids between Drosophila buzzatii and Drosophila koepferae, mobilization of at least 28 TEs has been described. However, the molecular mechanisms underlying this TE release remain poorly understood. To give insight on the causes of this TE activation, we performed a TE transcriptomic analysis in ovaries (notorious for playing a major role in TE silencing) of parental species and their F1 and backcrossed (BC) hybrids. We find that 15.2% and 10.6% of the expressed TEs are deregulated in F1 and BC1 ovaries, respectively, with a bias toward overexpression in both cases. Although differences between parental piRNA (Piwi-interacting RNA) populations explain only partially these results, we demonstrate that piRNA pathway proteins have divergent sequences and are differentially expressed between parental species. Thus, a functional divergence of the piRNA pathway between parental species, together with some differences between their piRNA pools, might be at the origin of hybrid instabilities and ultimately cause TE misregulation in ovaries. These analyses were complemented with the study of F1 testes, where TEs tend to be less expressed than in D. buzzatii. This can be explained by an increase in piRNA production, which probably acts as a defence mechanism against TE instability in the male germline. Hence, we describe a differential impact of interspecific hybridization in testes and ovaries, which reveals that TE expression and regulation are sex-biased.

Dissecting the Satellite DNA Landscape in Three Cactophilic Drosophila Sequenced Genomes

Eukaryote genomes are replete with repetitive DNAs. This class includes tandemly repeated satellite DNAs (satDNA) which are among the most abundant, fast evolving (yet poorly studied) genomic components. Here, we used high-throughput sequencing data from three cactophilic Drosophila species, D. buzzatii, D. seriema, and D. mojavensis, to access and study their whole satDNA landscape. In total, the RepeatExplorer software identified five satDNAs, three previously described (pBuM, DBC-150 and CDSTR198) and two novel ones (CDSTR138 and CDSTR130). Only pBuM is shared among all three species. The satDNA repeat length falls within only two classes, between 130 and 200 bp or between 340 and 390 bp. FISH on metaphase and polytene chromosomes revealed the presence of satDNA arrays in at least one of the following genomic compartments: centromeric, telomeric, subtelomeric, or dispersed along euchromatin. The chromosomal distribution ranges from a single chromosome to almost all chromosomes of the complement. Fiber-FISH and sequence analysis of contigs revealed interspersion between pBuM and CDSTR130 in the microchromosomes of D. mojavensis Phylogenetic analyses showed that the pBuM satDNA underwent concerted evolution at both interspecific and intraspecific levels. Based on RNA-seq data, we found transcription activity for pBuM (in D. mojavensis) and CDSTR198 (in D. buzzatii) in all five analyzed developmental stages, most notably in pupae and adult males. Our data revealed that cactophilic Drosophila present the lowest amount of satDNAs (1.9-2.9%) within the Drosophila genus reported so far. We discuss how our findings on the satDNA location, abundance, organization, and transcription activity may be related to functional aspects.

Expression of the Retrotransposon Helena Reveals a Complex Pattern of TE Deregulation in Drosophila Hybrids

Transposable elements (TEs), repeated mobile sequences, are ubiquitous in the eukaryotic kingdom. Their mobilizing capacity confers on them a high mutagenic potential, which must be strongly regulated to guarantee genome stability. In the Drosophila germline, a small RNA-mediated silencing system, the piRNA (Piwi-interacting RNA) pathway, is the main responsible TE regulating mechanism, but some stressful conditions can destabilize it. For instance, during interspecific hybridization, genomic stress caused by the shock of two different genomes can lead, in both animals and plants, to higher transposition rates. A recent study in D. buzatii-D. koepferae hybrids detected mobilization of 28 TEs, yet little is known about the molecular mechanisms explaining this transposition release. We have characterized one of the mobilized TEs, the retrotransposon Helena, and used quantitative expression to assess whether its high transposition rates in hybrids are preceded by increased expression. We have also localized Helena expression in the gonads to see if cellular expression patterns have changed in the hybrids. To give more insight into changes in TE regulation in hybrids, we analysed Helena-specific piRNA populations of hybrids and parental species. Helena expression is not globally altered in somatic tissues, but male and female gonads have different patterns of deregulation. In testes, Helena is repressed in F1, increasing then its expression up to parental values. This is linked with a mislocation of Helena transcripts along with an increase of their specific piRNA levels. Ovaries have additive levels of Helena expression, but the ping-pong cycle efficiency seems to be reduced in F1 hybrids. This could be at the origin of new Helena insertions in hybrids, which would be transmitted to F1 hybrid female progeny.

Variations on a theme: diversification of cuticular hydrocarbons in a clade of cactophilic Drosophila

Background

We characterized variation and chemical composition of epicuticular hydrocarbons (CHCs) in the seven species of the Drosophila buzzatii cluster with gas chromatography/mass spectrometry. Despite the critical role of CHCs in providing resistance to desiccation and involvement in communication, such as courtship behavior, mating, and aggregation, few studies have investigated how CHC profiles evolve within and between species in a phylogenetic context. We analyzed quantitative differences in CHC profiles in populations of the D. buzzatii species cluster in order to assess the concordance of CHC differentiation with species divergence.

Results

Thirty-six CHC components were scored in single fly extracts with carbon chain lengths ranging from C₂₉ to C₃₉, including methyl-branched alkanes, n-alkenes, and alkadienes. Multivariate analysis of variance revealed that CHC amounts were significantly different among all species and canonical discriminant function (CDF) analysis resolved all species into distinct, non-overlapping groups. Significant intraspecific variation was found in different populations of D. serido suggesting that this taxon is comprised of at least two species. We summarized CHC variation using CDF analysis and mapped the first five CHC canonical variates (CVs) onto an independently derived period (per) gene + chromosome inversion + mtDNA COI gene for each sex. We found that the COI sequences were not phylogenetically informative due to introgression between some species, so only per + inversion data were used. Positive phylogenetic signal was observed mainly for CV1 when parsimony methods and the test for serial independence (TFSI) were used. These results changed when no outgroup species were included in the analysis and phylogenetic signal was then observed for female CV3 and/or CV4 and male CV4 and CV5. Finally, removal of divergent populations of D. serido significantly increased the amount of phylogenetic signal as up to four out of five CVs then displayed positive phylogenetic signal.

Conclusions

CHCs were conserved among species while quantitative differences in CHC profiles between populations and species were statistically significant. Most CHCs were species-, population-, and sex-specific. Mapping CHCs onto an independently derived phylogeny revealed that a significant portion of CHC variation was explained by species' systematic affinities indicating phylogenetic conservatism in the evolution of these hydrocarbon arrays, presumptive waterproofing compounds and courtship signals as in many other drosophilid species.

The non-regular orbit: three satellite DNAs in Drosophila martensis (buzzatii complex, repleta group) followed three different evolutionary pathways

The genome of species from the buzzatii cluster (buzzatii complex, repleta group) is hosted by a number of satellite DNAs (satDNAs) showing contrasting structural characteristics, genomic organization and evolution, such as pBuM-alpha (~190 bp repeats), pBuM-alpha/beta (~370 bp repeats) and the DBC-150 (~150 bp repeats). In the present study, we aimed to investigate the evolution of these three satDNAs by looking for homologous sequences in the genome of the closest outgroup species: Drosophila martensis (buzzatii complex). After PCR, we isolated and sequenced 9 alpha, 8 alpha/beta and 11 DBC-150 sequences from this species. The results were compared to all pBuM and DBC-150 sequences available in literature. After D. martensis split from the buzzatii cluster some 6 Mya, the three satDNAs evolved differently in the genome of D. martensis by: (1) maintenance of a collection of major types of ancestral repeats in the genome (alpha); (2) fixation for a single major type of ancestral repeats (alpha/beta) or (3) fixation for new divergent species-specific repeat types (DBC-150). Curiously, D. seriema and D. martensis, although belonging to different and allopatric clusters, became independently fixed for the same major type of alpha/beta ancestral repeats, illustrating a rare case of parallelism in satDNA evolution. The contrasting pictures illustrate the diversity of evolutionary pathways a satDNA can follow, defining a "non-regular orbit" with outcomes difficult to predict.

Male genital and wing morphology in the cactophilic sibling species Drosophila gouveai and Drosophila antonietae and their hybrids reared in different host plants

Cactophilic Drosophila flies are excellent models to study adaptation to a relatively narrow spectrum of potential host plants and host-driven evolutionary diversification. Previous studies suggested a complex genetic architecture of wing and male genital morphology in phylogenetically basal species of the D. buzzatii cluster. In this work, we investigate the effect of experimental hybridization and host plant shifts on male genital and wing morphology in D. gouveai Tidon-Sklorz and Sene and D. antonietae Tidon-Sklorz and Sene, a pair of more recently derived species. We explicitly tested the hypotheses that wing and male genital morphology in interspecific hybrids depend on the host plant in which flies were grown. Our study shows that cactus hosts exert a strong effect on genital and wing morphology and that hybrids can be clearly differentiated on the basis of wing and genital morphology from both parental species. However, the extent of morphological differentiation between hybrids and pure species as well as plasticity patterns varied across organs, suggesting a complex genetic architecture for the studied traits.