Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8–38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group—a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.     

The chromosomes of two Drosophila races: D. nasuta nasuta and D. nasuta albomicana

The microchromosomes of the totally cross fertile Drosophila races, D. nasuta nasuta and D. nasuta albomicana have been studied in nietaphase and polytene nuclei. In metaphase the microchromosome of D. n. albomicana is nearly five times longer than the homologous chromosome in D. n. nasuta. As shown by C-banding these length differences are mainly due to a massive addition of heterochromatin to the D. n. albomicana chromosome. In polytene nuclei these striking heterochromatin differences between the microchromosomes of the two Drosophila races cannot be observed. Analysis of the polytene banding pattern shows that the microchromosomes of both races differ by an inversion and by a duplication, present only in D. n, albomicana. The location and orientation of the duplicated regions in D. n. albomicana leads to a specific loop like chromosome configuration. On the basis of these differences within the Drosophila races studied it is assumed that the karyotype of D. n. albomicana is a more recent evolutionary product.     

Sex‐biased transcriptome divergence along a latitudinal gradient

Sex‐dependent gene expression is likely an important genomic mechanism that allows sex‐specific adaptation to environmental changes. Among Drosophila species, sex‐biased genes display remarkably consistent evolutionary patterns; male‐biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex‐biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex‐specific selection and the evolution of sex‐biased gene expression. We assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11 631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male‐biased genes, there was no overrepresentation of X‐linked genes in males. By contrast, X‐linked divergence was elevated in females, especially for female‐biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. Our results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro‐ and micro‐ecological spatial scales.     

Taxonomic identification of Drosophila nasuta subgroup strains by glue protein analysis

Protein fractions of salivary glands were analyzed from 30 wildtype strains of eight species belonging to the Drosophila nasuta subgroup by SDS-polyacrylamide gel electrophoresis. The electrophoretic patterns indicated several prominent bands which could be shown to represent the major glue protein fractions. The glue protein fractions are species-specific as well as wildtype strain-specific. Wildtype strain specificities are characterized by variations of the species-specific patterns. The patterns of the different wildtypes, species, and hybrids were used for taxonomic identification within the nasuta subgroup, in which the females are morphologically indistinguishable and the males differ only by the markings of their frons. The hybrids provide evidence for gonosomal as well as autosomal linkage of individual genes coding for the major glue protein fractions.     

Male accessory gland secretory protein polymorphism in natural populations of <Emphasis Type="Italic">Drosophila nasuta nasuta</Emphasis> and <Emphasis Type="Italic">Drosophila sulfurigaster neonasuta</Emphasis>

Male accessory gland secretory protein polymorphism was analysed in natural populations of Drosophila nasuta nasuta and D. sulfurigaster neonasuta for the first time, using SDS-PAGE to score polymorphism of these proteins in 2788 individuals of D. n. nasuta and 2232 individuals of D. s. neonasuta from 12 different populations from southern India. A total of 25 and 18 variant protein phenotypes were identified in D. n. nasuta and D. s. neonasuta, respectively. Protein fractions of group III were more polymorphic than those from groups I and II. The results show that accessory gland secretory proteins show high levels of polymorphism, irrespective of species or habitat. Moreover, we have used the variation in the accessory gland proteins to assess the extent of divergence between the species and to infer their population structure. The study suggests that though both D. n. nasuta and D. s. neonasuta belong to the same subgroup, they differ in population structure, as far as accessory gland protein polymorphism is concerned.     

Adaptive Evolution and the Birth of CTCF Binding Sites in the Drosophila Genome

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ∼2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes.     

In situ study of chorion gene amplification in ovarian follicle cells ofDrosophila nasuta

The temporal and spatial pattern of replication of chorion gene clusters in follicle cells during oogenesis inDrosophila melanogaster andDrosophila nasuta was examined by [^3H thymidine autoradiography and byin situ hybridization with chorion gene probes. When pulse labelled with [³H] thymidine, the follicle cells from stage 10–12 ovarian follicles of bothDrosophila melanogaster and,Drosophila nasuta often showed intense labelling at only one or two sites per nucleus.In situ hybridization of chorion gene probes derived fromDrosophila melanogaster with follicle cell nuclei ofDrosophila melanogaster andDrosophila nasuta revealed these discrete [³H] thymidine labelled sites to correspond to the two amplifying chorion gene clusters. It appears, therefore, that in spite of evolutionary divergence, the organization and programme of selective amplification of chorion genes in ovarian follicle cells have remained generally similar in these two species. The endoreplicated and amplified copies of each chorion gene cluster remain closely associated but the two clusters occupy separate sites in follicle cell nucleus.     

The speciation history of the Drosophila nasuta complex

The Drosophila nasuta subgroup of the immigrans species group is widely distributed throughout the South-East Asian region, consisting of morphologically similar species with varying degrees of reproductive isolation. Here, I report nucleotide variability data for five X-linked and two mtDNA loci in eight taxa from the nasuta subgroup, with deeper sampling from D. albomicans and its sister species D. nasuta. Phylogenetic relationships among these species vary among different genomic regions, and levels of genetic differentiation suggest that this species group diversified only about one million years ago. D. albomicans and D. nasuta share nucleotide polymorphisms and are distinguished by relatively few fixed differences. Patterns of genetic differentiation between this species pair are compatible with a simple isolation model with no gene flow. Nucleotide variability levels of species in the nasuta group are comparable to those in members of the melanogaster and pseudoobscura species groups, indicating effective population sizes on the order of several million. Population genetic analyses reveal that summaries of the frequency distribution of neutral polymorphisms in both D. albomicans and D. nasuta generally fit the assumptions of the standard neutral model. D. albomicans is of particular interest for evolutionary studies because of its recently formed neo-sex chromosomes, and our phylogenetic and population genetic analyses suggest that it might be an ideal model to study the very early stages of Y chromosome evolution.     

Hybridization and introgression of the genomes of Drosophila nasuta and Drosophila albomicans: evolution of new karyotypes.

Drosophila nasuta (2n = 8) and Drosophila albomicans (2n = 6) are cross-fertile allopatric sibling chromosomal races of the nasuta subgroup of Drosophila. Hybrids of these races can be maintained for any number of generations. Some of the introgressed hybrid lineages of D. nasuta and D. albomicans, after passing through a transient phase of karyotypic polymorphism, ended up with a stable karyotype whose composition is different from those of the parental races. Such hybrid populations were called cytoraces, in which the chromosomes of D. nasuta and D. albomicans are represented in different combinations. The karyotypic composition of 16 such cytoraces have been presented and discussed with reference to evolutionary strategies such as balancing selection, directional selection, and sex-specific effect on different components of the evolving karyotypes.     

Extensive introgression of mitochondrial DNA relative to nuclear genes in the Drosophila yakuba species group

Studies of gene flow between recently diverged species can illuminate the role of natural selection in the formation of new species. Drosophila santomea and D. yakuba are recently diverged, partially reproductively isolated species that continue to hybridize in the wild, and appear to be reproductively isolated from the more distantly related species D. teissieri. We examine patterns of nucleotide polymorphism and divergence in these three species at multiple X-linked, Y-linked, and mitochondrial markers. All three species harbor drastically reduced variability on the Y chromosome relative to the X, as expected for a nonrecombining chromosome subject to variation-reducing selection. The three species are generally well differentiated at the nuclear markers, with little evidence for recent introgression for either the X- or Y-linked genes. Based on the nuclear genes, we estimate that D. santomea and D. yakuba diverged about one-half million years ago and split from D. teissieri about one million years ago. In contrast to the pattern at nuclear loci, all three species share a very similar mtDNA haplotype. We show that the mtDNA must have recently introgressed across species boundaries in the D. yakuba subgroup and that its fixation was driven by either selection on the mitochondria itself or other cytoplasmic factors. These results demonstrate that different regions of the genome can have distinct evolutionary dynamics in the context of species formation. Although natural selection is usually thought of as accentuating divergence between species, our results imply that it can also act as a homogenizing force.     

Taxonomic identification of Drosophila nasuta subgroup strains by glue protein analysis

Protein fractions of salivary glands were analyzed from 30 wildtype strains of eight species belonging to the Drosophila nasuta subgroup by SDS-polyacrylamide gel electrophoresis. The electrophoretic patterns indicated several prominent bands which could be shown to represent the major glue protein fractions. The glue protein fractions are species-specific as well as wildtype strain-specific. Wildtype strain specificities are characterized by variations of the species-specific patterns. The patterns of the different wildtypes, species, and hybrids were used for taxonomic identification within the nasuta subgroup, in which the females are morphologically indistinguishable and the males differ only by the markings of their frons. The hybrids provide evidence for gonosomal as well as autosomal linkage of individual genes coding for the major glue protein fractions.     

Metaphase chromosomes of four species of the Drosophila nasuta subgroup

Metaphase chromosomes of four species of the Drosophila nasuta subgroup, D. nasuta, D. kepulauana, D. kohkoa, and D. neveifrons, were analysed by the C banding method to clarify the chromosomal differentiation during speciation. Four species were different in heterochromatic regions of their chromosomes, especially microchromosomes of the fourth chromosomes and the Y chromosomes. Intraspecific variations of heterochromatic regions were also found among isofemale lines from various localities of D. nasuta, D. kepulauana, and D. kohkoa. Intraspecific variations of heterochromatic regions were not different from interspecific variations morphologically. From these results, the evolutionary process of heterochromatic regions was discussed.     

Comparison of the gap segmentation gene hunchback between Drosophila melanogaster and Drosophila virilis reveals novel modes of evolutionary change.

We have cloned and sequenced a large portion of the hunchback (hb) locus from Drosophila virilis. Comparison with the Drosophila melanogaster hb sequence shows multiple strong homologies in the upstream and downstream regions of the gene, including most of the known functional parts. The coding sequence is highly conserved within the presumptive DNA-binding finger regions, but more diverged outside of them. The regions of high divergence are correlated with regions which are rich in short direct repeats (regions of high 'cryptic simplicity'), suggesting a significant influence of slippage-like mechanisms in the evolutionary divergence of the two genes. Staining of early D.virilis embryos with an hb antibody reveals conserved and divergent features of the spatial expression pattern at blastoderm stage. It appears that the basic expression pattern, which serves as the gap gene function of hb, is conserved, while certain secondary expression patterns, which have separate functions for the segmentation process, are partly diverged. Thus, both slippage driven mutations in the coding region, which are likely to occur at higher rates than point mutations and the evolutionary divergence of secondary expression patterns may contribute to the evolution of regulatory genes.     

果蝇Drosophila nasuta亚群系统发育问题的探讨

果蝇Ｄｒｏｓｏｐｈｉｌａｎａｓｕｔａ亚群由在１４个处于不同物种分化阶段的种,亚种和分类元组成,这个亚群的物种有许多进化上的独特之处,使得它在物种分化研究方面倍受关注,然而,在形态学,生殖隔离,染色体和同工酶多态,线粒体ＤＮＡＲＦＬＰ,求偶歌特征惟及线粒体和核基因序列分析等方面的研究都未能清楚地阐明这一亚群的系统进化关系,本文综合分析了关于这一亚群的进化遗传学的研究结果,并提出了有待进一步的一些问题     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

Functional divergence caused by ancient positive selection of a Drosophila hybrid incompatibility locus

Interspecific hybrid lethality and sterility are a consequence of divergent evolution between species and serve to maintain the discrete identities of species. The evolution of hybrid incompatibilities has been described in widely accepted models by Dobzhansky and Muller where lineage-specific functional divergence is the essential characteristic of hybrid incompatibility genes. Experimentally tractable models are required to identify and test candidate hybrid incompatibility genes. Several Drosophila melanogaster genes involved in hybrid incompatibility have been identified but none has yet been shown to have functionally diverged in accordance with the Dobzhansky-Muller model. By introducing transgenic copies of the X-linked Hybrid male rescue (Hmr) gene into D. melanogaster from its sibling species D. simulans and D. mauritiana, we demonstrate that Hmr has functionally diverged to cause F1 hybrid incompatibility between these species. Consistent with the Dobzhansky-Muller model, we find that Hmr has diverged extensively in the D. melanogaster lineage, but we also find extensive divergence in the sibling-species lineage. Together, these findings implicate over 13% of the amino acids encoded by Hmr as candidates for causing hybrid incompatibility. The exceptional level of divergence at Hmr cannot be explained by neutral processes because we use phylogenetic methods and population genetic analyses to show that the elevated amino-acid divergence in both lineages is due to positive selection in the distant past—at least one million generations ago. Our findings suggest that multiple substitutions driven by natural selection may be a general phenomenon required to generate hybrid incompatibility alleles.     

The chromosomes of two Drosophila races: D. nasuta nasuta and D. nasuta albomicana

Interracial hybridization between Drosophila nasuta nasuta (2n=8) and D. n. albomicana (2n=6) has resulted in the evolution of two new karyotypic strains, called Cytoraces I and II. Males and females of Cytorace I have 2n=7 and 2n=6 respectively. The reconstituted karyotype is totally new in its composition, the chromosomes being drawn from both the parental races. The individuals of Cytorace II have 2n=6. Even though the chromosomes of the parental races are duly represented in the F₁, there is selective retention/elimination of certain chromosomes in the succeeding generations during which repatterning of the karyotype has taken place. Dynamics of each one of the parental chromosomes are presented and its implications re discussed.We dedicate this paper to the memory of the founder of our Department, the late Prof. M.R. Rajasekarasetty on the occasion of the Silver Jubilee of our Department     

Inversion polymorphism in populations of Drosophila sulphurigaster albostrigata and Drosophila nasuta albomicans from Phuket,Thailand

Two species of the Drosophila nasuta subgroup of the Drosophila immigrans group, D. sulphurigaster albostrigata and D. nasuta albomicans were investigated in this study. Collections of both species were made from Phuket, Thailand. Both species have similar salivary chromosomes, with four autosomal arms and one sex chromosome arm, and both are highly polymorphic for paracentric inversions. D. s. albostrigata accounted for the majority of the isolines collected and exhibited the greater number of inversions. One inversion, C₁, was common to both species, indicating common ancestry.A non-random distribution of inversions was observed on the proximal end of chromosome II in both D. s. albostrigata and D. n. albomicans. An inter-collection comparison revealed that both rigid and flexible chromosomal polymorphism were operating in the two species, with a seasonal variation noted for one inversion in D. s. albostrigata. A non-random association of two inversions was observed in D. n. albomicans.Based on a comparison of the indices of crossing over, both D. s. albostrigata and D. n. albomicans were found to be more heterozygous than in previous studies, with D. n. albomicans appearing to have evolved further than D. s. albostrigata.Based on a thesis submitted for the degree of Ph. D. at the University of Queensland.