The evolutionary history of Drosophila buzzatii XI. A new method for cytogenetic localization based on asynapsis of polytene chromosomes in interspecific hybrids of Drosophila

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

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

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

Genome-wide dissection of hybrid sterility in Drosophila confirms a polygenic threshold architecture

To date, different studies about the genetic basis of hybrid male sterility (HMS), a postzygotic reproductive barrier thoroughly investigated using Drosophila species, have demonstrated that no single major gene can produce hybrid sterility without the cooperation of several genetic factors. Early work using hybrids between Drosophila koepferae (Dk) and Drosophila buzzatii (Db) was consistent with the idea that HMS requires the cooperation of several genetic factors, supporting a polygenic threshold (PT) model. Here we present a genome-wide mapping strategy to test the PT model, analyzing serially backcrossed fertile and sterile males in which the Dk genome was introgressed into the Db background. We identified 32 Dk-specific markers significantly associated with hybrid sterility. Our results demonstrate 1) a strong correlation between the number of segregated sterility markers and males' degree of sterility, 2) the exchangeability among markers, 3) their tendency to cluster into low-recombining chromosomal regions, and 4) the requirement for a minimum number (threshold) of markers to elicit sterility. Although our findings do not contradict a role for occasional major hybrid-sterility genes, they conform more to the view that HMS primarily evolves by the cumulative action of many interacting genes of minor effect in a complex PT architecture.

Evaluation of the courtship and of the hybrid male sterility among Drosophila buzzatii cluster species (Diptera, Drosophilidae)

In the Drosophila repleta group the establishment of subgroups and complexes made on the basis of morphological and cytological evidences is supported by tests of reproductive isolation. Among species in the repleta group, the buzzatii cluster, due to its polymorphism and polytipism, is an excellent material for ecological and speciation studies. Some interspecific crosses involving Drosophila seriema, Drosophila sp. B, D. koepferae and D. buzzatii strains were completely sterile while others involving strains from these species produced F1 hybrids that did not yield F2. In the present work, data on courtship duration and copula occurrence obtained in the analysis of flies from parental sterile crosses and on spermatozoon mobility observed in F1 hybrids that did not yield F2 are presented. Copula did not occur during one hour of observation and the spermatozoon also did not show mobility at any of the analyzed stages (3, 7, 9 and 10 days old). There was a high variation in courtship average duration and in the percentage of males that courted the females. The reproductive isolation mechanisms indicated by these observations were pre and post-zygotic, as supported by the absence of copula and male sterility. Data obtained also showed the occurrence of different degrees of reproductive compatibility among the strains classified as the same species but from distinct geographic localities.

Chromosomal elements evolve at different rates in the Drosophila genome

Recent results indicate that the rate of chromosomal rearrangement in the genus Drosophila is the highest found so far in any eukaryote. This conclusion is based chiefly on the comparative mapping analysis of a single chromosomal element (Muller's element E) in two species, D. melanogaster and D. repleta, representing the two farthest lineages within the genus (the Sophophora and Drosophila subgenera, respectively). We have extended the analysis to two other chromosomal elements (Muller's elements A and D) and tested for differences in rate of evolution among chromosomes. With this purpose, detailed physical maps of chromosomes X and 4 of D. repleta were constructed by in situ hybridization of 145 DNA probes (gene clones, cosmids, and P1 phages) and their gene arrangements compared with those of the homologous chromosomes X and 3L of D. melanogaster. Both chromosomal elements have been extensively reshuffled over their entire length. The number of paracentric inversions fixed has been estimated as 118 +/- 17 for element A and 56 +/- 8 for element D. Comparison with previous data for elements E and B shows that there are fourfold differences in evolution rate among chromosomal elements, with chromosome X exhibiting the highest rate of rearrangement. Combining all results, we estimated that 393 paracentric inversions have been fixed in the whole genome since the divergence between D. repleta and D. melanogaster. This amounts to an average rate of 0.053 disruptions/Mb/myr, corroborating the high rate of rearrangement in the genus Drosophila.

A system of complementary genes in hybrids between Drosophila koepferae and D. buzzatii: a Markov chain model allows to make inferences about their number and relationships

In backcrosses between D. koepferae and D. buzzatii, the disruption of a system of species-specific complementary factors brings about hybrid male inviability. This system consists of a lethal factor, hmi-1, linked to the X chromosome of D. koepferae, and several conspecific autosomal suppressors. However, hmi-1 hybrid males can also be rescued by factors present in some strains of D. buzzatii. The present work aims to estimate the number of hmi-1 suppressors in one of these strains by means of Markov chains. The obtained results allow discarding models with one or more chromosomes having independent suppressor effect. On the other hand, models having n chromosomes that interact in groups of r, being 1 < r < or = n, to produce rescue effect, provide good approximations to the observed results. The best fit to the data is obtained with four or five chromosomes with suppressor effect, interacting epistatically in groups of three to rescue the viability of hmi-1 males.

Toward a physical map of Drosophila buzzatii. Use of randomly amplified polymorphic dna polymorphisms and sequence-tagged site landmarks

We present a physical map based on RAPD polymorphic fragments and sequence-tagged sites (STSs) for the repleta group species Drosophila buzzatii. One hundred forty-four RAPD markers have been used as probes for in situ hybridization to the polytene chromosomes, and positive results allowing the precise localization of 108 RAPDs were obtained. Of these, 73 behave as effectively unique markers for physical map construction, and in 9 additional cases the probes gave two hybridization signals, each on a different chromosome. Most markers (68%) are located on chromosomes 2 and 4, which partially agree with previous estimates on the distribution of genetic variation over chromosomes. One RAPD maps close to the proximal breakpoint of inversion 2z(3) but is not included within the inverted fragment. However, it was possible to conclude from this RAPD that the distal breakpoint of 2z(3) had previously been wrongly assigned. A total of 39 cytologically mapped RAPDs were converted to STSs and yielded an aggregate sequence of 28,431 bp. Thirty-six RAPDs (25%) did not produce any detectable hybridization signal, and we obtained the DNA sequence from three of them. Further prospects toward obtaining a more developed genetic map than the one currently available for D. buzzatii are discussed.

Characterization of Gandalf, a new inverted-repeat transposable element of Drosophila koepferae

The cloning and characterization of Gandalf, a new DNA-transposing mobile element obtained from the Drosophila koepferae (repleta group) genome is described. A fragment of Gandalf was found in a middle repetitive clone that shows variable chromosomal localization. Restriction, Southern blot, PCR and sequencing analyses have shown that most Gandalf copies are about 1 kb long, are flanked by 12 bp inverted terminal repeats and contain subterminal repetitive regions on both sides of the element. As with other elements of the DNA-transposing type (known as the 'Ac family'), the Gandalf element generates 8 bp direct duplications at the insertion point. Coding region analysis has shown that the longer open reading frame found in Gandalf copies could encode part of a protein. However, whether or not the 1 kb copies of the element are actually the active transposons remains to be elucidated. Gandalf shows a very low copy number in D. buzzatii, a sibling species of D. koepferae. An attempt to induce interspecific hybrid dysgenesis in hybrids of these two species has been unsuccessful.

High transposition rates of Osvaldo, a new Drosophila buzzatii retrotransposon

Transposition of a new Drosophila retrotransposon was investigated. Total genomic Southern analysis and polytene in situ hybridizations in D. buzzatii strains and other related species using a 6 kb D. buzzatii clone (cDb314) showed a dispersed, repetitive DNA pattern, suggesting that this clone contains a transposable element (TE). We have sequenced the cDb314 clone and demonstrated that it contains all the conserved protein sequences and motifs typical of retrovirus-related sequences. Although cDb314 does not include the complete TE, the protein sequence alignment demonstrates that it includes a defective copy of a new long terminal repeat (LTR) retrotransposon, related to the gypsy family, which we have named Osvaldo. Using a D. buzzatii inbred line in which all insertion sites are known, we have measured Osvaldo transposition rates in hybrids between this D. buzzatii line and its sibling species D. koepferae. The results show that Osvaldo transposes in bursts at high rate, both in the D. buzzatii inbred line and in species hybrids.

The evolutionary history of Drosophila buzzatii. XXI. Cumulative action of multiple sterility factors on spermatogenesis in hybrids of D. buzzatii and D. koepferae

The genetic basis of sterility in male hybrids of Drosophila buzzatii and D. koepferae has been investigated by assessment of the effects on spermatogenesis of substituting separate chromosome segments of the recipient species with the homologous material from the donor species, either in heterozygous (autosomes) or hemizygous (X chromosome) condition, after successive backcrossing of hybrid females to either parental species. Introgressed segments were identified by the characteristic asynapsis of the polytene chromosomes in their heterospecific regions. Except for one case, the introgression of chromosome segments either from autosome 3, 4, or 5 brings about sterility only when the introgressed segment exceeds a minimum size (threshold size). Segments of equal size frequently produce similar abnormalities, whose severity increases with the size of the introgressed segment. Apparently, throughout these autosomes of D. buzzatii and D. koepferae there are many non-allelic, minor sterility genes, whose individual segregation cannot be recognized phenotypically, and which act cumulatively on the same characteristics of spermatogenesis, each contributing a small effect to the phenotype. Accordingly, these genes should be considered as polygenes, and the type of sterility they bring about should be properly designated polygenic sterility.

The evolutionary history of D. buzzatii. XXII. Chromosomal and genic sterility in male hybrids of Drosophila buzzatii and Drosophila koepferae

The genetic basis of sterility in F1 male hybrids of Drosophila buzzatii and D. koepferae has been investigated in two steps. (1) By successive backcrossing of hybrid females to either parental species. (2) By assessment of the effects on male fertility of selected segments of polytene chromosomes from the donor species on a background entirely derived from the recipient species. The length of introgressed segments producing sterility was progressively reduced through repeated backcrosses. This procedure sometimes led to an approximate mapping of major genes of hybrid sterility (genic sterility) on the polytene chromosome map. At other times it was found that sterility was produced only when the introgressed segment exceeded a certain threshold size (chromosomal sterility). The contribution of the autosomes to hybrid sterility seems to be mainly of the chromosomal type. The evidence concerning the X chromosome is equivocal. No fertile males were found following introgression with any of the investigated segments of this chromosome. These results are compatible both with the presence of at least six major genes of hybrid sterility (genic sterility) and with the existence of a rather small threshold size for the chromosome segments producing sterility (chromosomal sterility). The role of the Y chromosome was not investigated in this study.

Differential response to environmental alcohol among second-chromosome arrangements in experimental populations of Drosophila buzzatii

Drosophila buzzatii feeds and breeds on the decaying cladodes and fruits of several species of Opuntia (prickly pear) which contain significant levels of ethanol and isopropanol. The potential influence of these two alcohols on the inversion polymorphism of the second and fourth chromosomes was investigated in ten experimental populations with different amounts of alcohol (either ethanol or isopropanol) added to the culture medium. All populations were started with the offspring of 29 wild females collected at Adeje (Tenerife, Canary Islands) and their genetic composition was monitored for about two years (more than 30 generations). Consistent changes in the frequency of most second- and fourth-chromosome arrangements occurred in all populations including those without alcohol (control). The comparison of inversion frequency through the various treatments revealed a significant influence of the alcohol on the frequency changes of the four second-chromosome arrangements. Moreover, this influence was of a different type for every one of them, sometimes with opposite effects between alcohols and/or concentrations. These results indicate genetic differentiation among second-chromosome arrangements with regard to alcohol and suggest that the alcohol heterogeneity found in the species' trophic niche may play an important role in the maintenance of this polymorphism and also in the recent historical changes in the frequency of some arrangements associated with colonization.