Hawaiian Drosophila genomes: size variation and evolutionary expansions

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai, Drosophila takahashii, Drosophila bipectinata, and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma. However, the assembled portion of the Muller F element (orthologous to the fourth chromosome in Drosophila melanogaster) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to >18 Mb. Here, we present chromosome-level long-read genome assemblies for 4 Drosophila species with expanded F elements ranging in size from 2.3 to 20.5 Mb. Each Muller element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai , Drosophila takahashii , Drosophila bipectinata , and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma . However, the assembled portion of the Muller F Element (orthologous to the fourth chromosome in Drosophila melanogaster ) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to > 18 Mb. Here, we present chromosome-level long read genome assemblies for four Drosophila species with expanded F Elements ranging in size from 2.3 Mb to 20.5 Mb. Each Muller Element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

In Silico Identification and Characterization of Satellite DNAs in 23 Drosophila Species from the Montium Group

Satellite DNA (satDNA) is a class of tandemly repeated non-protein coding DNA sequences which can be found in abundance in eukaryotic genomes. They can be functional, impact the genomic architecture in many ways, and their rapid evolution has consequences for species diversification. We took advantage of the recent availability of sequenced genomes from 23 Drosophila species from the montium group to study their satDNA landscape. For this purpose, we used publicly available whole-genome sequencing Illumina reads and the TAREAN (tandem repeat analyzer) pipeline. We provide the characterization of 101 non-homologous satDNA families in this group, 93 of which are described here for the first time. Their repeat units vary in size from 4 bp to 1897 bp, but most satDNAs show repeat units < 100 bp long and, among them, repeats ≤ 10 bp are the most frequent ones. The genomic contribution of the satDNAs ranges from ~1.4% to 21.6%. There is no significant correlation between satDNA content and genome sizes in the 23 species. We also found that at least one satDNA originated from an expansion of the central tandem repeats (CTRs) present inside a Helitron transposon. Finally, some satDNAs may be useful as taxonomic markers for the identification of species or subgroups within the group.

Chromosome Comparisons of Australian Scaptodrosophila Species

The Scaptodrosophila represent a diverse group of Diptera closely related to Drosophila. Although they have radiated extensively in Australia, they have been the focus of few studies. Here, we characterized the karyotypes of 12 Scaptodrosophila species from several species groups and showed that they have undergone similar types of karyotypic change to those seen in Drosophila. This includes heterochromatin amplification involved in length changes of the sex and ‘dot’ chromosomes as well as the autosomes, particularly in the coracina group of species. Numerous weak points along the arms of the polytene chromosomes suggest the presence of internal repetitive sequence DNA, but these regions did not C-band in mitotic chromosomes, and their analysis will depend on DNA sequencing. The nucleolar organizing regions (NORs) are at the same chromosome positions in Scaptodrosophila as in Drosophila, and the various mechanisms responsible for changing arm configurations also appear to be the same. These chromosomal studies provide a complementary resource to other investigations of this group, with several species currently being sequenced.

Phylotranscriptomics reveals discordance in the phylogeny of Hawaiian Drosophila and Scaptomyza (Diptera: Drosophilidae)

Island radiations present natural laboratories for studying the evolutionary process. The Hawaiian Drosophilidae are one such radiation, with nearly 600 described species and substantial morphological and ecological diversification. These species are largely divided into a few major clades, but the relationship between clades remains uncertain. Here, we present new assembled transcriptomes from 12 species across these clades, and use these transcriptomes to resolve the base of the evolutionary radiation. We recover a new hypothesis for the relationship between clades, and demonstrate its support over previously published hypotheses. We then use the evolutionary radiation to explore dynamics of concordance in phylogenetic support, by analyzing the gene and site concordance factors for every possible topological combination of major groups. We show that high bootstrap values mask low evolutionary concordance, and we demonstrate that the most likely topology is distinct from the topology with the highest support across gene trees and from the topology with highest support across sites. We then combine all previously published genetic data for the group to estimate a time-calibrated tree for over 300 species of drosophilids. Finally, we digitize dozens of published Hawaiian Drosophilidae descriptions, and use this to pinpoint probable evolutionary shifts in reproductive ecology as well as body, wing, and egg size. We show that by examining the entire landscape of tree and trait space, we can gain a more complete understanding of how evolutionary dynamics play out across an island radiation.

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.

Genome Size Covaries More Positively with Propagule Size than Adult Size: New Insights into an Old Problem

The body size and (or) complexity of organisms is not uniformly related to the amount of genetic material (DNA) contained in each of their cell nuclei ('genome size'). This surprising mismatch between the physical structure of organisms and their underlying genetic information appears to relate to variable accumulation of repetitive DNA sequences, but why this variation has evolved is little understood. Here, I show that genome size correlates more positively with egg size than adult size in crustaceans. I explain this and comparable patterns observed in other kinds of animals and plants as resulting from genome size relating strongly to cell size in most organisms, which should also apply to single-celled eggs and other reproductive propagules with relatively few cells that are pivotal first steps in their lives. However, since body size results from growth in cell size or number or both, it relates to genome size in diverse ways. Relationships between genome size and body size should be especially weak in large organisms whose size relates more to cell multiplication than to cell enlargement, as is generally observed. The ubiquitous single-cell 'bottleneck' of life cycles may affect both genome size and composition, and via both informational (genotypic) and non-informational (nucleotypic) effects, many other properties of multicellular organisms (e.g., rates of growth and metabolism) that have both theoretical and practical significance.

Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex

Heterochromatin is identified as a potential factor driving diversification of species. To understand the magnitude of heterochromatin variation within the Anopheles gambiae complex of malaria mosquitoes, we analyzed metaphase chromosomes in An. arabiensis, An. coluzzii, An. gambiae, An. merus, and An. quadriannulatus. Using fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA), a highly repetitive fraction of DNA, and heterochromatic Bacterial Artificial Chromosome (BAC) clones, we established the correspondence of pericentric heterochromatin between the metaphase and polytene X chromosomes of An. gambiae. We then developed chromosome idiograms and demonstrated that the X chromosomes exhibit qualitative differences in their pattern of heterochromatic bands and position of satellite DNA (satDNA) repeats among the sibling species with postzygotic isolation, An. arabiensis, An. merus, An. quadriannulatus, and An. coluzzii or An. gambiae. The identified differences in the size and structure of the X chromosome heterochromatin point to a possible role of repetitive DNA in speciation of mosquitoes. We found that An. coluzzii and An. gambiae, incipient species with prezygotic isolation, share variations in the relative positions of the satDNA repeats and the proximal heterochromatin band on the X chromosomes. This previously unknown genetic polymorphism in malaria mosquitoes may be caused by a differential amplification of DNA repeats or an inversion in the sex chromosome heterochromatin.

De novo identification of satellite DNAs in the sequenced genomes of Drosophila virilis and D. americana using the RepeatExplorer and TAREAN pipelines

Satellite DNAs are among the most abundant repetitive DNAs found in eukaryote genomes, where they participate in a variety of biological roles, from being components of important chromosome structures to gene regulation. Experimental methodologies used before the genomic era were insufficient, too laborious and time-consuming to recover the collection of all satDNAs from a genome. Today, the availability of whole sequenced genomes combined with the development of specific bioinformatic tools are expected to foster the identification of virtually all the "satellitome" of a particular species. While whole genome assemblies are important to obtain a global view of genome organization, most of them are incomplete and lack repetitive regions. We applied short-read sequencing and similarity clustering in order to perform a de novo identification of the most abundant satellite families in two Drosophila species from the virilis group: Drosophila virilis and D. americana, using the Tandem Repeat Analyzer (TAREAN) and RepeatExplorer pipelines. These species were chosen because they have been used as models to understand satDNA biology since the early 70's. We combined the computational approach with data from the literature and chromosome mapping to obtain an overview of the major tandem repeat sequences of these species. The fact that all of the abundant tandem repeats (TRs) we detected were previously identified in the literature allowed us to evaluate the efficiency of TAREAN in correctly identifying true satDNAs. Our results indicate that raw sequencing reads can be efficiently used to detect satDNAs, but that abundant tandem repeats present in dispersed arrays or associated with transposable elements are frequent false positives. We demonstrate that TAREAN with its parent method RepeatExplorer may be used as resources to detect tandem repeats associated with transposable elements and also to reveal families of dispersed tandem repeats.

The Drosophila Dot Chromosome: Where Genes Flourish Amidst Repeats

The F element of the Drosophila karyotype (the fourth chromosome in Drosophila melanogaster) is often referred to as the "dot chromosome" because of its appearance in a metaphase chromosome spread. This chromosome is distinct from other Drosophila autosomes in possessing both a high level of repetitious sequences (in particular, remnants of transposable elements) and a gene density similar to that found in the other chromosome arms, ∼80 genes distributed throughout its 1.3-Mb "long arm." The dot chromosome is notorious for its lack of recombination and is often neglected as a consequence. This and other features suggest that the F element is packaged as heterochromatin throughout. F element genes have distinct characteristics (e.g, low codon bias, and larger size due both to larger introns and an increased number of exons), but exhibit expression levels comparable to genes found in euchromatin. Mapping experiments show the presence of appropriate chromatin modifications for the formation of DNaseI hypersensitive sites and transcript initiation at the 5' ends of active genes, but, in most cases, high levels of heterochromatin proteins are observed over the body of these genes. These various features raise many interesting questions about the relationships of chromatin structures with gene and chromosome function. The apparent evolution of the F element as an autosome from an ancestral sex chromosome also raises intriguing questions. The findings argue that the F element is a unique chromosome that occupies its own space in the nucleus. Further study of the F element should provide new insights into chromosome structure and function.

Hawaiian Drosophila as an Evolutionary Model Clade: Days of Future Past

The Hawaiian Drosophila have been a model system for evolutionary, ecological, and ethological studies since the inception of the Hawaiian Drosophila Project in the 1960s. Here we review the past and present research on this incredible lineage and provide a prospectus for future directions on genomics and microbial interactions. While the number of publications on this group has waxed and waned over the years, we assert that recent systematic, biogeographic, and ecological studies have reinvigorated Hawaiian Drosophila as an evolutionary model system. The characteristics that distinguish good model clades from good model organisms (e.g., Drosophila melanogaster) are somewhat different so we first define what constitutes a good evolutionary model. We argue that the Hawaiian Drosophila possess many desired aspects of a good evolutionary model, describe how this group of geographically isolated flies have been used in the past, and propose some exciting avenues for future evolutionary research on this diverse, dynamic clade of Drosophila.

Genomes of Diptera

Diptera (true flies) are among the most diverse holometabolan insect orders and were the first eukaryotic order to have a representative genome fully sequenced. 110 fly species have publically available genome assemblies and many hundreds of population-level genomes have been generated in the model organisms Drosophila melanogaster and the malaria mosquito Anopheles gambiae. Comparative genomics carried out in a phylogenetic context is illuminating many aspects of fly biology, providing unprecedented insight into variability in genome structure, gene content, genetic mechanisms, and rates and patterns of evolution in genes, populations, and species. Despite the rich availability of genomic resources in flies, there remain many fly lineages to which new genome sequencing efforts should be directed. Such efforts would be most valuable in fly families or clades that exhibit multiple origins of key fly behaviors such as blood feeding, phytophagy, parasitism, pollination, and mycophagy.

Retrotransposons Are the Major Contributors to the Expansion of the Drosophila ananassae Muller F Element

The discordance between genome size and the complexity of eukaryotes can partly be attributed to differences in repeat density. The Muller F element (∼5.2 Mb) is the smallest chromosome in Drosophila melanogaster, but it is substantially larger (>18.7 Mb) in D. ananassae To identify the major contributors to the expansion of the F element and to assess their impact, we improved the genome sequence and annotated the genes in a 1.4-Mb region of the D. ananassae F element, and a 1.7-Mb region from the D element for comparison. We find that transposons (particularly LTR and LINE retrotransposons) are major contributors to this expansion (78.6%), while Wolbachia sequences integrated into the D. ananassae genome are minor contributors (0.02%). Both D. melanogaster and D. ananassae F-element genes exhibit distinct characteristics compared to D-element genes (e.g., larger coding spans, larger introns, more coding exons, and lower codon bias), but these differences are exaggerated in D. ananassae Compared to D. melanogaster, the codon bias observed in D. ananassae F-element genes can primarily be attributed to mutational biases instead of selection. The 5' ends of F-element genes in both species are enriched in dimethylation of lysine 4 on histone 3 (H3K4me2), while the coding spans are enriched in H3K9me2. Despite differences in repeat density and gene characteristics, D. ananassae F-element genes show a similar range of expression levels compared to genes in euchromatic domains. This study improves our understanding of how transposons can affect genome size and how genes can function within highly repetitive domains.

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.

Profuse evolutionary diversification and speciation on volcanic islands: transposon instability and amplification bursts explain the genetic paradox

BACKGROUND

Species-rich adaptive radiations arising from rare plant and animal colonizers are common on remote volcanic archipelagoes. However, they present a paradox. The severe genetic bottleneck of founder events and effects of inbreeding depression, coupled with the inherently stressful volcanic environment, would seem to predict reduced evolutionary potential and increased risk of extinction, rather than rapid adaptive divergence and speciation. Significantly, eukaryotic genomes harbor many families of transposable elements (TEs) that are mobilized by genome shock; these elements may be the primary drivers of genetic reorganization and speciation on volcanic islands.

PRESENTATION OF THE HYPOTHESIS

Here I propose that a central factor in the spectacular radiation and diversification of the endemic Hawaiian Drosophila and other terrestrial lineages on the Hawaiian and other oceanic islands has been repeated bursts of transposition of multiple TEs induced by the unique ecological features of volcanic habitats. Founder individuals and populations on remote volcanic islands experience significant levels of physiological and genomic stress as a consequence of both biotic and abiotic factors. This results in disruption of the usual epigenetic suppression of TEs, unleashing them to proliferate and spread, which in turn gives rise to novel genetic variation and remodels genomic regulatory circuits, facilitating rapid morphological, ecological and behavioral change, and adaptive radiation.

TESTING THE HYPOTHESIS

To obtain empirical support for the hypothesis, test organisms should be exposed to prolonged heat stress, high levels of carbon dioxide and other volcanic gases, along with inbreeding. Data from subsequent whole genome sequencing and bioinformatics screening for TE numbers and locations would then be compared with initial pre-exposure TE information for the test strains, a labor-intensive project. Several predicted outcomes arising from the hypothesis are discussed. Currently available data are consistent with the proposed concept of stress-induced TE mobilization as a trigger of evolutionary diversification and speciation on volcanic islands.

IMPLICATIONS OF THE HYPOTHESIS

The main implication is that both TEs and species should proliferate at a much higher rate on volcanic islands than elsewhere. Second, the evolvability of a lineage may correlate with the abundance and distribution of TEs in the genome. Successful colonizers of volcanic habitats with high genomic proportions of TEs may be best poised to found a speciose lineage that gives rise to a dramatic adaptive radiation. Colonizers that are depauperate in TEs are likely to be evolutionarily constrained and diversify little, if at all.

REVIEWERS

This article was reviewed by Dr. James Shapiro and Dr. Wolfgang Miller (nominated by Editorial Board member Dr. I. King Jordan).

Cytogenetic Analysis of the South American Fruit Fly Anastrepha fraterculus (Diptera:Tephritidae) Species Complex: Construction of Detailed Photographic Polytene Chromosome Maps of the Argentinian Af. sp.1 Member

Genetic and cytogenetic studies constitute a significant basis for understanding the biology of insect pests and the design and the construction of genetic tools for biological control strategies. Anastrepha fraterculus is an important pest of the Tephritidae family. It is distributed from southern Texas through eastern Mexico, Central America and South America causing significant crop damage and economic losses. Currently it is considered as a species complex; until now seven members have been described based on multidisciplinary approaches. Here we report the cytogenetic analysis of an Argentinian population characterized as Af. sp.1 member of the Anastrepha fraterculus species complex. The mitotic karyotype and the first detailed photographic maps of the salivary gland polytene chromosomes are presented. The mitotic metaphase complement consists of six (6) pairs of chromosomes, including one pair of heteromorphic sex chromosomes, with the male being the heterogametic sex. The analysis of the salivary gland polytene complement shows a total number of five long chromosomes that correspond to the five autosomes of the mitotic karyotype and a heterochromatic network corresponding to the sex chromosomes. Comparison of the polytene chromosome maps between this species and Anastrepha ludens shows significant similarity. The polytene maps presented here are suitable for cytogenetic studies that could shed light on the species limits within this species complex and support the development of genetic tools for sterile insect technique (SIT) applications.