Genome size diversity in the family Drosophilidae

The implications of satellite DNA instability on cellular function and evolution

Abundant tandemly repeated satellite DNA is present in most eukaryotic genomes. Previous limitations including a pervasive view that it was uninteresting junk DNA, combined with challenges in studying it, are starting to dissolve - and recent studies have found important functions for satellite DNAs. The observed rapid evolution and implied instability of satellite DNA now has important significance for their functions and maintenance within the genome. In this review, we discuss the processes that lead to satellite DNA copy number instability, and the importance of mechanisms to manage the potential negative effects of instability. Satellite DNA is vulnerable to challenges during replication and repair, since it forms difficult-to-process secondary structures and its homology within tandem arrays can result in various types of recombination. Satellite DNA instability may be managed by DNA or chromatin-binding proteins ensuring proper nuclear localization and repair, or by proteins that process aberrant structures that satellite DNAs tend to form. We also discuss the pattern of satellite DNA mutations from recent mutation accumulation (MA) studies that have tracked changes in satellite DNA for up to 1000 generations with minimal selection. Finally, we highlight examples of satellite evolution from studies that have characterized satellites across millions of years of Drosophila fruit fly evolution, and discuss possible ways that selection might act on the satellite DNA composition.

Genome sizes and repeatome evolution in zoantharians (Cnidaria: Hexacorallia: Zoantharia)

Across eukaryotes, large variations of genome sizes have been observed even between closely related species. Transposable elements as part of the repeated DNA have been proposed and confirmed as one of the most important contributors to genome size variation. However, the evolutionary implications of genome size variation and transposable element dynamics are not well understood. Together with phenotypic traits, they are commonly referred to as the "C-value enigma". The order Zoantharia are benthic cnidarians found from intertidal zones to the deep sea, and some species are particularly abundant in coral reefs. Despite their high ecological relevance, zoantharians have yet to be largely studied from the genomic point of view. This study aims at investigating the role of the repeatome (total content of repeated elements) in genome size variations across the order Zoantharia. To this end, whole-genomes of 32 zoantharian species representing five families were sequenced. Genome sizes were estimated and the abundances of different repeat classes were assessed. In addition, the repeat overlap between species was assessed by a sequence clustering method. The genome sizes in the dataset varied up to 2.4 fold magnitude. Significant correlations between genome size, repeated DNA content and transposable elements, respectively (Pearson's correlation test R2 = 0.47, p = 0.0016; R2 = 0.22, p = 0.05) were found, suggesting their involvement in the dynamics of genome expansion and reduction. In all species, long interspersed nuclear elements and DNA transposons were the most abundant identified elements. These transposable elements also appeared to have had a recent expansion event. This was in contrast to the comparative clustering analysis which revealed species-specific patterns of satellite elements' amplification. In summary, the genome sizes of zoantharians likely result from the complex dynamics of repeated elements. Finally, the majority of repeated elements (up to 70%) could not be annotated to a known repeat class, highlighting the need to further investigate non-model cnidarian genomes. More research is needed to understand how repeated DNA dynamics relate to zoantharian evolution and their biology.

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.

Genome size variation and polyploidy prevalence in the genus Eragrostis are associated with the global dispersal in arid area

BACKGROUND

Biologists have long debated the drivers of the genome size evolution and variation ever since Darwin. Assumptions for the adaptive or maladaptive consequences of the associations between genome sizes and environmental factors have been proposed, but the significance of these hypotheses remains controversial. Eragrostis is a large genus in the grass family and is often used as crop or forage during the dry seasons. The wide range and complex ploidy levels make Eragrostis an excellent model for investigating how the genome size variation and evolution is associated with environmental factors and how these changes can ben interpreted.

METHODS

We reconstructed the Eragrostis phylogeny and estimated genome sizes through flow cytometric analyses. Phylogenetic comparative analyses were performed to explore how genome size variation and evolution is related to their climatic niches and geographical ranges. The genome size evolution and environmental factors were examined using different models to study the phylogenetic signal, mode and tempo throughout evolutionary history.

RESULTS

Our results support the monophyly of Eragrostis. The genome sizes in Eragrostis ranged from ~0.66 pg to ~3.80 pg. We found that a moderate phylogenetic conservatism existed in terms of the genome sizes but was absent from environmental factors. In addition, phylogeny-based associations revealed close correlations between genome sizes and precipitation-related variables, indicating that the genome size variation mainly caused by polyploidization may have evolved as an adaptation to various environments in the genus Eragrostis.

CONCLUSION

This is the first study to take a global perspective on the genome size variation and evolution in the genus Eragrostis. Our results suggest that the adaptation and conservatism are manifested in the genome size variation, allowing the arid species of Eragrostis to spread the xeric area throughout the world.

Genome Survey Sequencing of the Mole Cricket Gryllotalpa orientalis

The mole cricket Gryllotalpa orientalis is an evolutionarily, medicinal, and agriculturally significant insect that inhabits underground environments and is distributed globally. This study measured genome size by flow cytometry and k-mer based on low-coverage sequencing, and nuclear repetitive elements were also identified. The haploid genome size estimate is 3.14 Gb by flow cytometry, 3.17 Gb, and 3.77 Gb-based two k-mer methods, respectively, which is well within the range previously reported for other species of the suborder Ensifera. 56% of repetitive elements were found in G. orientalis, similar to 56.83% in Locusta migratoria. However, the great size of repetitive sequences could not be annotated to specific repeat element families. For the repetitive elements that were annotated, Class I-LINE retrotransposon elements were the most common families and more abundant than satellite and Class I-LTR. These results based on the newly developed genome survey could be used in the taxonomic study and whole genome sequencing to improve the understanding of the biology of G. orientalis.

New estimates of genome size in Orthoptera and their evolutionary implications

Animal genomes vary widely in size, and much of their architecture and content remains poorly understood. Even among related groups, such as orders of insects, genomes may vary in size by orders of magnitude-for reasons unknown. The largest known insect genomes were repeatedly found in Orthoptera, e.g., Podisma pedestris (1C = 16.93 pg), Stethophyma grossum (1C = 18.48 pg) and Bryodemella holdereri (1C = 18.64 pg). While all these species belong to the suborder of Caelifera, the ensiferan Deracantha onos (1C = 19.60 pg) was recently found to have the largest genome. Here, we present new genome size estimates of 50 further species of Ensifera (superfamilies Gryllidea, Tettigoniidea) and Caelifera (Acrididae, Tetrigidae) based on flow cytometric measurements. We found that Bryodemella tuberculata (Caelifera: Acrididae) has the so far largest measured genome of all insects with 1C = 21.96 pg (21.48 gBp). Species of Orthoptera with 2n = 16 and 2n = 22 chromosomes have significantly larger genomes than species with other chromosome counts. Gryllidea genomes vary between 1C = 0.95 and 2.88 pg, and Tetrigidae between 1C = 2.18 and 2.41, while the genomes of all other studied Orthoptera range in size from 1C = 1.37 to 21.96 pg. Reconstructing ancestral genome sizes based on a phylogenetic tree of mitochondrial genomic data, we found genome size values of >15.84 pg only for the nodes of Bryodemella holdereri / B. tuberculata and Chrysochraon dispar / Euthystira brachyptera. The predicted values of ancestral genome sizes are 6.19 pg for Orthoptera, 5.37 pg for Ensifera, and 7.28 pg for Caelifera. The reasons for the large genomes in Orthoptera remain largely unknown, but a duplication or polyploidization seems unlikely as chromosome numbers do not differ much. Sequence-based genomic studies may shed light on the underlying evolutionary mechanisms.

The genome of the freshwater monogonont rotifer Brachionus rubens: Identification of phase I, II, and III detoxification genes

Monogonont rotifers are common species in aquatic environments and make model species for ecotoxicology studies. Whole genomes of several species of the genus Brachionus have been assembled, but no information on the freshwater rotifer Brachionus rubens has been reported. In this study, the whole-genome sequence of B. rubens was successfully assembled using NextDenovo. The total length of the genome was 132.7 Mb (N50 = 2.51 Mb), including 122 contigs. The GC contents accounted for 29.96% of the genome. Aquatic organisms are always exposed to various external stresses, and a comprehensive genomic analysis is needed to better understand the adverse effects on organisms. This paper focuses on the ecotoxicological aspect and conducted genome analysis of representative gene families involved in detoxification mechanisms against environmental stressors. Specifically, we identified cytochrome P450 genes (CYPs) of phase I, glutathione S-transferase genes (GSTs) of phase II, and ATP-binding cassette transporter genes (ABCs) of phase III in the genome of B. rubens. Gene duplications were found in CYP, GST, and ABC genes, as is the case for other Brachionus rotifers. Our results suggest that these detoxification-related gene families have evolved in a species-specific and/or lineage-specific manner. This paper improves our understanding of how the freshwater Brachionus rotifers respond to environmental stressors in a molecular ecotoxicology context.

High quality genome assembly of the anhydrobiotic midge provides insights on a single chromosome-based emergence of extreme desiccation tolerance

Non-biting midges (Chironomidae) are known to inhabit a wide range of environments, and certain species can tolerate extreme conditions, where the rest of insects cannot survive. In particular, the sleeping chironomid Polypedilum vanderplanki is known for the remarkable ability of its larvae to withstand almost complete desiccation by entering a state called anhydrobiosis. Chromosome numbers in chironomids are higher than in other dipterans and this extra genomic resource might facilitate rapid adaptation to novel environments. We used improved sequencing strategies to assemble a chromosome-level genome sequence for P. vanderplanki for deep comparative analysis of genomic location of genes associated with desiccation tolerance. Using whole genome-based cross-species and intra-species analysis, we provide evidence for the unique functional specialization of Chromosome 4 through extensive acquisition of novel genes. In contrast to other insect genomes, in the sleeping chironomid a uniquely high degree of subfunctionalization in paralogous anhydrobiosis genes occurs in this chromosome, as well as pseudogenization in a highly duplicated gene family. Our findings suggest that the Chromosome 4 in Polypedilum is a site of high genetic turnover, allowing it to act as a 'sandbox' for evolutionary experiments, thus facilitating the rapid adaptation of midges to harsh environments.

Biology and ecology of the Oriental flower-breeding Drosophila elegans and related species

Animals adapt to their environments in the course of evolution. One effective approach to elucidate mechanisms of adaptive evolution is to compare closely related species with model organisms in which knowledge of the molecular and physiological bases of various traits has been accumulated. Drosophila elegans and its close relatives, belonging to the same species group as the model organism D. melanogaster, exhibit various unique characteristics such as flower-breeding habit, courtship display, territoriality, sexual dimorphism, and colour polymorphism. Their ease of culturing and availability of genomic information makes them a useful model for understanding mechanisms of adaptive evolution. Here, we review the morphology, distribution, and phylogenetic relationships of D. elegans and related species, as well as their characteristic flower-dependent biology, food habits, and life-history traits. We also describe their unique mating and territorial behaviours and note their distinctive karyotype and the genetic mechanisms of morphological diversity that have recently been revealed.

A haploid diamondback moth (Plutella xylostella L.) genome assembly resolves 31 chromosomes and identifies a diamide resistance mutation

The diamondback moth, Plutella xylostella (L.), is a highly mobile brassica crop pest with worldwide distribution and can rapidly evolve resistance to insecticides, including group 28 diamides. Reference genomes assembled using Illumina sequencing technology have provided valuable resources to advance our knowledge regarding the biology, origin and movement of diamondback moth, and more recently with its sister species, Plutella australiana. Here we apply a trio binning approach to sequence and annotate a chromosome level reference genome of P. xylostella using PacBio Sequel and Dovetail Hi-C sequencing technology and identify a point mutation that causes resistance to commercial diamides. A P. xylostella population collected from brassica crops in the Lockyer Valley, Australia (LV-R), was reselected for chlorantraniliprole resistance then a single male was crossed to a P. australiana female and a hybrid pupa sequenced. A chromosome level 328 Mb P. xylostella genome was assembled with 98.1% assigned to 30 autosomes and the Z chromosome. The genome was highly complete with 98.4% of BUSCO Insecta genes identified and RNAseq informed protein prediction annotated 19,002 coding genes. The LV-R strain survived recommended field application doses of chlorantraniliprole, flubendiamide and cyclaniliprole. Some hybrids also survived these doses, indicating significant departure from recessivity, which has not been previously documented for diamides. Diamide chemicals modulate insect Ryanodine Receptors (RyR), disrupting calcium homeostasis, and we identified an amino acid substitution (I4790K) recently reported to cause diamide resistance in a strain from Japan. This chromosome level assembly provides a new resource for insect comparative genomics and highlights the emergence of diamide resistance in Australia. Resistance management plans need to account for the fact that resistance is not completely recessive.

Correlated evolution of larval development, egg size and genome size across two genera of snapping shrimp

Across plants and animals, genome size is often correlated with life-history traits: large genomes are correlated with larger seeds, slower development, larger body size and slower cell division. Among decapod crustaceans, caridean shrimps are among the most variable both in terms of genome size variation and life-history characteristics such as larval development mode and egg size, but the extent to which these traits are associated in a phylogenetic context is largely unknown. In this study, we examine correlations among egg size, larval development and genome size in two different genera of snapping shrimp, Alpheus and Synalpheus, using phylogenetically informed analyses. In both Alpheus and Synalpheus, egg size is strongly linked to larval development mode: species with abbreviated development had significantly larger eggs than species with extended larval development. We produced the first comprehensive dataset of genome size in Alpheus (n = 37 species) and demonstrated that genome size was strongly and positively correlated with egg size in both Alpheus and Synalpheus. Correlated trait evolution analyses showed that in Alpheus, changes in genome size were clearly dependent on egg size. In Synalpheus, evolutionary path analyses suggest that changes in development mode (from extended to abbreviated) drove increases in egg volume; larger eggs, in turn, resulted in larger genomes. These data suggest that variation in reproductive traits may underpin the high degree of variation in genome size seen in a wide variety of caridean shrimp groups more generally.

The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights

BACKGROUND AND AIMS

Genome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species.

METHODS

We generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples.

KEY RESULTS

We found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small.

CONCLUSIONS

We show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.

Genome Size Dynamics in Marine Ribbon Worms (Nemertea, Spiralia)

Nemertea is a phylum consisting of 1300 mostly marine species. Nemertea is distinguished by an eversible muscular proboscis, and most of the species are venomous. Genomic resources for this phylum are scarce despite their value in understanding biodiversity. Here, we present genome size estimates of Nemertea based on flow cytometry and their relationship to different morphological and developmental traits. Ancestral genome size estimations were done across the nemertean phylogeny. The results increase the available genome size estimates for Nemertea three-fold. Our analyses show that Nemertea has a narrow genome size range (0.43-3.89 pg) compared to other phyla in Lophotrochozoa. A relationship between genome size and evolutionary rate, developmental modes, and habitat was found. Trait analyses show that the highest evolutionary rate of genome size is found in upper intertidal, viviparous species with direct development. Despite previous findings, body size in nemerteans was not correlated with genome size. A relatively small genome (1.18 pg) is assumed for the most recent common ancestor of all extant nemerteans. The results provide an important basis for future studies in nemertean genomics, which will be instrumental to understanding the evolution of this enigmatic and often neglected phylum.

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 Estimation and Full-Length Transcriptome of Sphingonotus tsinlingensis: Genetic Background of a Drought-Adapted Grasshopper

Sphingonotus Fieber, 1852 (Orthoptera: Acrididae), is a grasshopper genus comprising approximately 170 species, all of which prefer dry environments such as deserts, steppes, and stony benchlands. In this study, we aimed to examine the adaptation of grasshopper species to arid environments. The genome size of Sphingonotus tsinlingensis was estimated using flow cytometry, and the first high-quality full-length transcriptome of this species was produced. The genome size of S. tsinlingensis is approximately 12.8 Gb. Based on 146.98 Gb of PacBio sequencing data, 221.47 Mb full-length transcripts were assembled. Among these, 88,693 non-redundant isoforms were identified with an N50 value of 2,726 bp, which was markedly longer than previous grasshopper transcriptome assemblies. In total, 48,502 protein-coding sequences were identified, and 37,569 were annotated using public gene function databases. Moreover, 36,488 simple tandem repeats, 12,765 long non-coding RNAs, and 414 transcription factors were identified. According to gene functions, 61 cytochrome P450 (CYP450) and 66 heat shock protein (HSP) genes, which may be associated with drought adaptation of S. tsinlingensis, were identified. We compared the transcriptomes of S. tsinlingensis and two other grasshopper species which were less tolerant to drought, namely Mongolotettix japonicus and Gomphocerus licenti. We observed the expression of CYP450 and HSP genes in S. tsinlingensis were higher. We produced the first full-length transcriptome of a Sphingonotus species that has an ultra-large genome. The assembly characteristics were better than those of all known grasshopper transcriptomes. This full-length transcriptome may thus be used to understand the genetic background and evolution of grasshoppers.

The genome of the euryhaline rotifer Brachionus paranguensis: Potential use in molecular ecotoxicology

Brachionus spp. rotifers have been proposed as model organisms for ecotoxicological studies. We analyzed the whole-genome sequence of B. paranguensis through NextDenovo, resulting in a total length of 106.2 Mb and 71 contigs. The N50 and the GC content were 4.13 Mb and 28%, respectively. A total of 18,501 genes were predicted within the genome of B. paranguensis. Prominent detoxification-related gene families of phase I and II detoxifications have been investigated. In parallel with other Brachionus rotifers, high gene expansion was observed in CYP clan 3 and GST sigma class in B. paranguensis. Moreover, species-specific expansion of sulfotransferase (SULTs) and gain of UDP-glucuronosyltransferases (UGTs) through horizontal gene transfer has been specifically found within B. plicatilis complex. This whole-genome analysis of B. paranguensis provides a basis for molecular ecotoxicological studies and provides useful information for comparative studies of the evolution of detoxification mechanisms in Brachionus spp.

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.

Evolutionary Dynamics of Abundant 7-bp Satellites in the Genome of Drosophila virilis

The factors that drive the rapid changes in abundance of tandem arrays of highly repetitive sequences, known as satellite DNA, are not well understood. Drosophila virilis has one of the highest relative amounts of simple satellites of any organism that has been studied, with an estimated >40% of its genome composed of a few related 7-bp satellites. Here, we use D. virilis as a model to understand technical biases affecting satellite sequencing and the evolutionary processes that drive satellite composition. By analyzing sequencing data from Illumina, PacBio, and Nanopore platforms, we identify platform-specific biases and suggest best practices for accurate characterization of satellites by sequencing. We use comparative genomics and cytogenetics to demonstrate that the highly abundant AAACTAC satellite family arose from a related satellite in the branch leading to the virilis phylad 4.5-11 Ma before exploding in abundance in some species of the clade. The most abundant satellite is conserved in sequence and location in the pericentromeric region but has diverged widely in abundance among species, whereas the satellites nearest the centromere are rapidly turning over in sequence composition. By analyzing multiple strains of D. virilis, we saw that the abundances of two centromere-proximal satellites are anticorrelated along a geographical gradient, which we suggest could be caused by ongoing conflicts at the centromere. In conclusion, we illuminate several key attributes of satellite evolutionary dynamics that we hypothesize to be driven by processes including selection, meiotic drive, and constraints on satellite sequence and abundance.

Transposable Element Mobilization in Interspecific Yeast Hybrids

Barbara McClintock first hypothesized that interspecific hybridization could provide a “genomic shock” that leads to the mobilization of transposable elements (TEs). This hypothesis is based on the idea that regulation of TE movement is potentially disrupted in hybrids. However, the handful of studies testing this hypothesis have yielded mixed results. Here, we set out to identify if hybridization can increase transposition rate and facilitate colonization of TEs in Saccharomyces cerevisiae × Saccharomyces uvarum interspecific yeast hybrids. Saccharomyces cerevisiae have a small number of active long terminal repeat retrotransposons (Ty elements), whereas their distant relative S. uvarum have lost the Ty elements active in S. cerevisiae. Although the regulation system of Ty elements is known in S. cerevisiae, it is unclear how Ty elements are regulated in other Saccharomyces species, and what mechanisms contributed to the loss of most classes of Ty elements in S. uvarum. Therefore, we first assessed whether TEs could insert in the S. uvarum sub-genome of a S. cerevisiae × S. uvarum hybrid. We induced transposition to occur in these hybrids and developed a sequencing technique to show that Ty elements insert readily and nonrandomly in the S. uvarum genome. We then used an in vivo reporter construct to directly measure transposition rate in hybrids, demonstrating that hybridization itself does not alter rate of mobilization. However, we surprisingly show that species-specific mitochondrial inheritance can change transposition rate by an order of magnitude. Overall, our results provide evidence that hybridization can potentially facilitate the introduction of TEs across species boundaries and alter transposition via mitochondrial transmission, but that this does not lead to unrestrained proliferation of TEs suggested by the genomic shock theory.

Shrinking body sizes in response to warming: explanations for the temperature-size rule with special emphasis on the role of oxygen

Body size is central to ecology at levels ranging from organismal fecundity to the functioning of communities and ecosystems. Understanding temperature-induced variations in body size is therefore of fundamental and applied interest, yet thermal responses of body size remain poorly understood. Temperature-size (T-S) responses tend to be negative (e.g. smaller body size at maturity when reared under warmer conditions), which has been termed the temperature-size rule (TSR). Explanations emphasize either physiological mechanisms (e.g. limitation of oxygen or other resources and temperature-dependent resource allocation) or the adaptive value of either a large body size (e.g. to increase fecundity) or a short development time (e.g. in response to increased mortality in warm conditions). Oxygen limitation could act as a proximate factor, but we suggest it more likely constitutes a selective pressure to reduce body size in the warm: risks of oxygen limitation will be reduced as a consequence of evolution eliminating genotypes more prone to oxygen limitation. Thus, T-S responses can be explained by the 'Ghost of Oxygen-limitation Past', whereby the resulting (evolved) T-S responses safeguard sufficient oxygen provisioning under warmer conditions, reflecting the balance between oxygen supply and demands experienced by ancestors. T-S responses vary considerably across species, but some of this variation is predictable. Body-size reductions with warming are stronger in aquatic taxa than in terrestrial taxa. We discuss whether larger aquatic taxa may especially face greater risks of oxygen limitation as they grow, which may be manifested at the cellular level, the level of the gills and the whole-organism level. In contrast to aquatic species, terrestrial ectotherms may be less prone to oxygen limitation and prioritize early maturity over large size, likely because overwintering is more challenging, with concomitant stronger end-of season time constraints. Mechanisms related to time constraints and oxygen limitation are not mutually exclusive explanations for the TSR. Rather, these and other mechanisms may operate in tandem. But their relative importance may vary depending on the ecology and physiology of the species in question, explaining not only the general tendency of negative T-S responses but also variation in T-S responses among animals differing in mode of respiration (e.g. water breathers versus air breathers), genome size, voltinism and thermally associated behaviour (e.g. heliotherms).

The genome of the marine monogonont rotifer Brachionus rotundiformis and insight into species-specific detoxification components in Brachionus spp

The monogonont rotifer Brachionus spp. have been widely used for ecotoxicological studies because of their advantages as one of the most suitable laboratory experimental species. In the present study, we obtained and assembled the whole genome sequence of the rotifer Brachionus rotundiformis, consisting of 13,612 annotated genes with 213 scaffolds and 58 Mb in total length. Focusing on ecotoxicological aspects, we conducted a comparative genome analysis on the gene families involved in detoxification, including four to six sulfotransferase gene families, seven uridine 5'-diphospho-glucuronosyltransferase gene families, and 58, 61, or 70 ATP-binding cassette genes in the genus Brachionus including Brachionus koreanus and Brachionus plicatilis. Our results suggest that these gene families have undergone a species- and/or lineage-specific evolution in response to the surrounding environmental pressure. Our genome resource for B. rotundiformis would be highly useful for future ecotoxicological studies and also provides a better understanding on the view of evolutionary mechanism of detoxification in the genus Brachionus spp.

Intraspecific variation in the karyotype length and genome size of fungus-farming ants (genus Mycetophylax), with remarks on procedures for the estimation of genome size in the Formicidae by flow cytometry

Ants (Formicidae) present considerable diversity in chromosome numbers, which vary from n = 1 to n = 60, although this variation is not proportional to that in genome size, for which estimates range from 0.18 pg to 0.77 pg. Intraspecific variation in the chromosome number and karyotype structure has been reported among species, although the variation among populations of the same species has received much less attention, and there are few data on genome size. Here, we studied the karyotype length and genome size of different populations of the fungus-farming ants Mycetophylax conformis (Mayr, 1884) and Mycetophylax morschi (Emery, 1888). We also provide remarks on procedure for the estimation of ant genome size by Flow Cytometry (FCM) analysis. Chromosome number and morphology did not vary among the populations of M. conformis or the cytotypes of M. morschi, but karyotype length and genome size were significantly distinct among the populations of these ants. Our results on the variation in karyotype length and genome size among M. morschi and M. conformis populations reveal considerable diversity that would be largely overlooked by more traditional descriptions of karyotypes, which were also supported by the estimates of genome size obtained using flow cytometry. Changes in the amount of DNA reflect variation in the fine structure of the chromosomes, which may represent the first steps of karyotype evolution and may occur previously to any changes in the chromosome number.

Near-chromosome level genome assembly of the fruit pest Drosophila suzukii using long-read sequencing

Over the past decade, the spotted wing Drosophila, Drosophila suzukii, has invaded Europe and America and has become a major agricultural pest in these areas, thereby prompting intense research activities to better understand its biology. Two draft genome assemblies already exist for this species but contain pervasive assembly errors and are highly fragmented, which limits their values. Our purpose here was to improve the assembly of the D. suzukii genome and to annotate it in a way that facilitates comparisons with D. melanogaster. For this, we generated PacBio long-read sequencing data and assembled a novel, high-quality D. suzukii genome assembly. It is one of the largest Drosophila genomes, notably because of the expansion of its repeatome. We found that despite 16 rounds of full-sib crossings the D. suzukii strain that we sequenced has maintained high levels of polymorphism in some regions of its genome. As a consequence, the quality of the assembly of these regions was reduced. We explored possible origins of this high residual diversity, including the presence of structural variants and a possible heterogeneous admixture pattern of North American and Asian ancestry. Overall, our assembly and annotation constitute a high-quality genomic resource that can be used for both high-throughput sequencing approaches, as well as manipulative genetic technologies to study D. suzukii.

The genome of the harpacticoid copepod Tigriopus japonicus: Potential for its use in marine molecular ecotoxicology

The copepod Tigriopus japonicus has been widely used as an experimental species in the field of ecotoxicology. We have sequenced and assembled the whole genome of T. japonicus with comparative analysis of gene families that represent detoxification phases in two additional public genomes of Tigriopus spp., namely, T. californicus and T. kingsejongensis. The total length of the T. japonicus assembled genome was 196.6 Mb with an N50 value of 10.65 Mb and consisted of 339 scaffolds and 25,143 annotated genes. The detoxification gene families encoding cytochrome P450s (CYP450s), glutathione S-transferases (GSTs), and ATP-binding cassette (ABC) proteins in Tigriopus spp. have shown species-dependent diversity in several gene sets, suggesting that these genes have undergone a species-specific expansion to increase their fitness to different marine habitats and environmental pressures. Our study will provide a better understanding of the detoxification system in Tigriopus spp. and will contribute to various areas of research, including ecotoxicology.

Effect of Phenotype Selection on Genome Size Variation in Two Species of Diptera

Genome size varies widely across organisms yet has not been found to be related to organismal complexity in eukaryotes. While there is no evidence for a relationship with complexity, there is evidence to suggest that other phenotypic characteristics, such as nucleus size and cell-cycle time, are associated with genome size, body size, and development rate. However, what is unknown is how the selection for divergent phenotypic traits may indirectly affect genome size. Drosophila melanogaster were selected for small and large body size for up to 220 generations, while Cochliomyia macellaria were selected for 32 generations for fast and slow development. Size in D. melanogaster significantly changed in terms of both cell-count and genome size in isolines, but only the cell-count changed in lines which were maintained at larger effective population sizes. Larger genome sizes only occurred in a subset of D. melanogaster isolines originated from flies selected for their large body size. Selection for development time did not change average genome size yet decreased the within-population variation in genome size with increasing generations of selection. This decrease in variation and convergence on a similar mean genome size was not in correspondence with phenotypic variation and suggests stabilizing selection on genome size in laboratory conditions.

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.

Absence of W Chromosome in Psychidae Moths and Implications for the Theory of Sex Chromosome Evolution in Lepidoptera

Moths and butterflies (Lepidoptera) are the largest group with heterogametic females. Although the ancestral sex chromosome system is probably Z0/ZZ, most lepidopteran species have the W chromosome. When and how the W chromosome arose remains elusive. Existing hypotheses place the W origin either at the common ancestor of Ditrysia and Tischeriidae, or prefer independent origins of W chromosomes in these two groups. Due to their phylogenetic position at the base of Ditrysia, bagworms (Psychidae) play an important role in investigating the W chromosome origin. Therefore, we examined the W chromosome status in three Psychidae species, namely Proutia betulina, Taleporia tubulosa, and Diplodoma laichartingella, using both classical and molecular cytogenetic methods such as sex chromatin assay, comparative genomic hybridization (CGH), and male vs. female genome size comparison by flow cytometry. In females of all three species, no sex chromatin was found, no female-specific chromosome regions were revealed by CGH, and a Z-chromosome univalent was observed in pachytene oocytes. In addition, the genome size of females was significantly smaller than males. Overall, our study provides strong evidence for the absence of the W chromosome in Psychidae, thus supporting the hypothesis of two independent W chromosome origins in Tischeriidae and in advanced Ditrysia.

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.

Long-read based assembly and synteny analysis of a reference Drosophila subobscura genome reveals signatures of structural evolution driven by inversions recombination-suppression effects

BACKGROUND

Drosophila subobscura has long been a central model in evolutionary genetics. Presently, its use is hindered by the lack of a reference genome. To bridge this gap, here we used PacBio long-read technology, together with the available wealth of genetic marker information, to assemble and annotate a high-quality nuclear and complete mitochondrial genome for the species. With the obtained assembly, we performed the first synteny analysis of genome structure evolution in the subobscura subgroup.

RESULTS

We generated a highly-contiguous ~ 129 Mb-long nuclear genome, consisting of six pseudochromosomes corresponding to the six chromosomes of a female haploid set, and a complete 15,764 bp-long mitogenome, and provide an account of their numbers and distributions of codifying and repetitive content. All 12 identified paracentric inversion differences in the subobscura subgroup would have originated by chromosomal breakage and repair, with some associated duplications, but no evidence of direct gene disruptions by the breakpoints. Between lineages, inversion fixation rates were 10 times higher in continental D. subobscura than in the two small oceanic-island endemics D. guanche and D. madeirensis. Within D. subobscura, we found contrasting ratios of chromosomal divergence to polymorphism between the A sex chromosome and the autosomes.

CONCLUSIONS

We present the first high-quality, long-read sequencing of a D. subobscura genome. Our findings generally support genome structure evolution in this species being driven indirectly, through the inversions' recombination-suppression effects in maintaining sets of adaptive alleles together in the face of gene flow. The resources developed will serve to further establish the subobscura subgroup as model for comparative genomics and evolutionary indicator of global change.

Genome Size Evolution within and between the Sexes

Genome sizes are known to vary between closely related species, but the patterns behind this variation have yet to be fully understood. Although this variation has been evaluated between species and within sexes, unknown is the extent to which this variation is driven by differentiation in sex chromosomes. To address this longstanding question, we examine the mode and tempo of genome size evolution for a total of 87 species of Drosophilidae, estimating and updating male genome size values for 44 of these species. We compare the evolution of genome size within each sex to the evolution of the differences between the sexes. Utilizing comparative phylogenetic methods, we find that male and female genome size evolution is largely a neutral process, reflective of phylogenetic relatedness between species, which supports the newly proposed accordion model for genome size change. When similarly analyzed, the difference between the sexes due to heteromorphic sex chromosomes is a dynamic process; the male-female genome size difference increases with time with or without known neo-Y events or complete loss of the Y. Observed instances of rapid change match theoretical expectations and known neo-Y and Y loss events in individual species.

Genome Size Reversely Correlates With Host Plant Range in Helicoverpa Species

In organisms with very low percentages of transposable elements (TEs), genome size may positively or negatively correlate with host range, depending on whether host adaptation or host modification is the main route to host generalism. To test if this holds true for insect herbivores with greater percentages of TEs, we conducted flow cytometry to measure the endopolyploidy levels and C-values of the host modification (salivary gland and mandibular gland in head), host adaptation (midgut), and host use-independent tissues (male gonad, hemolymph, and Malpighian tubules) of the generalist Helicoverpa armigera and the head of its older specialist sister H. assulta. Larval salivary gland displayed a consecutive chain of endopolyploidy particles from 8Cx to higher than 32Cx and larval head and midgut had endopolyploidy nuclei clusters of 16Cx and 32Cx, whereas larval male gonad, hemolymph, and Malpighian tubules possessed no endopolyploidy nuclei of higher than 8Cx. The estimated genome size of the Solanaceae plant specialist H. assulta is 430 Mb, significantly larger than that of its older generalist sister Heliothis virescens (408 Mb) and those of its two generalist descendants H. armigera (394 Mb) and H. zea (363 Mb). These data not only reveal a negative correlation between host plant range and genome size in this terminal lineage, but also imply that Helicoverpa species appear to depend more on host modification than on host adaptation to achieve polyphagy.

Ten things you should know about transposable elements

Transposable elements (TEs) are major components of eukaryotic genomes. However, the extent of their impact on genome evolution, function, and disease remain a matter of intense interrogation. The rise of genomics and large-scale functional assays has shed new light on the multi-faceted activities of TEs and implies that they should no longer be marginalized. Here, we introduce the fundamental properties of TEs and their complex interactions with their cellular environment, which are crucial to understanding their impact and manifold consequences for organismal biology. While we draw examples primarily from mammalian systems, the core concepts outlined here are relevant to a broad range of organisms.

Highly Contiguous Genome Assemblies of 15 Drosophila Species Generated Using Nanopore Sequencing

The Drosophila genus is a unique group containing a wide range of species that occupy diverse ecosystems. In addition to the most widely studied species, Drosophila melanogaster, many other members in this genus also possess a well-developed set of genetic tools. Indeed, high-quality genomes exist for several species within the genus, facilitating studies of the function and evolution of cis-regulatory regions and proteins by allowing comparisons across at least 50 million years of evolution. Yet, the available genomes still fail to capture much of the substantial genetic diversity within the Drosophila genus. We have therefore tested protocols to rapidly and inexpensively sequence and assemble the genome from any Drosophila species using single-molecule sequencing technology from Oxford Nanopore. Here, we use this technology to present highly contiguous genome assemblies of 15 Drosophila species: 10 of the 12 originally sequenced Drosophila species (ananassae, erecta, mojavensis, persimilis, pseudoobscura, sechellia, simulans, virilis, willistoni, and yakuba), four additional species that had previously reported assemblies (biarmipes, bipectinata, eugracilis, and mauritiana), and one novel assembly (triauraria). Genomes were generated from an average of 29x depth-of-coverage data that after assembly resulted in an average contig N50 of 4.4 Mb. Subsequent alignment of contigs from the published reference genomes demonstrates that our assemblies could be used to close over 60% of the gaps present in the currently published reference genomes. Importantly, the materials and reagents cost for each genome was approximately $1,000 (USD). This study demonstrates the power and cost-effectiveness of long-read sequencing for genome assembly in Drosophila and provides a framework for the affordable sequencing and assembly of additional Drosophila genomes.

Loss of males from mixed-sex societies in termites

Background

Sexual reproduction is the norm in almost all animal species, and in many advanced animal societies, both males and females participate in social activities. To date, the complete loss of males from advanced social animal lineages has been reported only in ants and honey bees (Hymenoptera), whose workers are always female and whose males display no helping behaviors even in normal sexual species. Asexuality has not previously been observed in colonies of another major group of social insects, the termites, where the ubiquitous presence of both male and female workers and soldiers indicate that males play a critical role beyond that of reproduction.

Results

Here, we report asexual societies in a lineage of the termite Glyptotermes nakajimai. We investigated the composition of mature colonies from ten distinct populations in Japan, finding six asexual populations characterized by a lack of any males in the reproductive, soldier, and worker castes of their colonies, an absence of sperm in the spermathecae of their queens, and the development of unfertilized eggs at a level comparable to that for the development of fertilized eggs in sexual populations of this species. Phylogenetic analyses indicated a single evolutionary origin of the asexual populations, with divergence from sampled sexual populations occurring about 14 million years ago. Asexual colonies differ from sexual colonies in having a more uniform head size in their all-female soldier caste, and fewer soldiers in proportion to other individuals, suggesting increased defensive efficiencies arising from uniform soldier morphology. Such efficiencies may have contributed to the persistence and spread of the asexual lineage. Cooperative colony foundation by multiple queens, the single-site nesting life history common to both the asexual and sexual lineages, and the occasional development of eggs without fertilization even in the sexual lineage are traits likely to have been present in the ancestors of the asexual lineage that may have facilitated the transition to asexuality.

Conclusions

Our findings demonstrate that completely asexual social lineages can evolve from mixed-sex termite societies, providing evidence that males are dispensable for the maintenance of advanced animal societies in which they previously played an active social role.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0563-y) contains supplementary material, which is available to authorized users.

The High-Quality Genome Sequence of the Oceanic Island Endemic Species Drosophila guanche Reveals Signals of Adaptive Evolution in Genes Related to Flight and Genome Stability

Drosophila guanche is a member of the obscura group that originated in the Canary Islands archipelago upon its colonization by D. subobscura. It evolved into a new species in the laurisilva, a laurel forest present in wet regions that in the islands have only minor long-term weather fluctuations. Oceanic island endemic species such as D. guanche can become model species to investigate not only the relative role of drift and adaptation in speciation processes but also how population size affects nucleotide variation. Moreover, the previous identification of two satellite DNAs in D. guanche makes this species attractive for studying how centromeric DNA evolves. As a prerequisite for its establishment as a model species suitable to address all these questions, we generated a high-quality D. guanche genome sequence composed of 42 cytologically mapped scaffolds, which are assembled into six super-scaffolds (one per chromosome). The comparative analysis of the D. guanche proteome with that of twelve other Drosophila species identified 151 genes that were subject to adaptive evolution in the D. guanche lineage, with a subset of them being involved in flight and genome stability. For example, the Centromere Identifier (CID) protein, directly interacting with centromeric satellite DNA, shows signals of adaptation in this species. Both genomic analyses and FISH of the two satellites would support an ongoing replacement of centromeric satellite DNA in D. guanche.

Extensive exchange of transposable elements in the Drosophila pseudoobscura group

BACKGROUND

As species diverge, so does their transposable element (TE) content. Within a genome, TE families may eventually become dormant due to host-silencing mechanisms, natural selection and the accumulation of inactive copies. The transmission of active copies from a TE families, both vertically and horizontally between species, can allow TEs to escape inactivation if it occurs often enough, as it may allow TEs to temporarily escape silencing in a new host. Thus, the contribution of horizontal exchange to TE persistence has been of increasing interest.

RESULTS

Here, we annotated TEs in five species with sequenced genomes from the D. pseudoobscura species group, and curated a set of TE families found in these species. We found that, compared to host genes, many TE families showed lower neutral divergence between species, consistent with recent transmission of TEs between species. Despite these transfers, there are differences in the TE content between species in the group.

CONCLUSIONS

The TE content is highly dynamic in the D. pseudoobscura species group, frequently transferring between species, keeping TEs active. This result highlights how frequently transposable elements are transmitted between sympatric species and, despite these transfers, how rapidly species TE content can diverge.

Strong phylogenetic inertia on genome size and transposable element content among 26 species of flies

While the evolutionary mechanisms driving eukaryote genome size evolution are still debated, repeated element content appears to be crucial. Here, we reconstructed the phylogeny and identified repeats in the genome of 26 Drosophila exhibiting a twofold variation in genome size. The content in transposable elements (TEs) is highly correlated to genome size evolution among these closely related species. We detected a strong phylogenetic signal on the evolution of both genome size and TE content, and a genome contraction in the Drosophila melanogaster subgroup.

Genome size of 14 species of fireflies (Insecta, Coleoptera, Lampyridae)

Eukaryotic genome size data are important both as the basis for comparative research into genome evolution and as estimators of the cost and difficulty of genome sequencing programs for non-model organisms. In this study, the genome size of 14 species of fireflies (Lampyridae) (two genera in Lampyrinae, three genera in Luciolinae, and one genus in subfamily incertae sedis) were estimated by propidium iodide (PI)-based flow cytometry. The haploid genome sizes of Lampyridae ranged from 0. 42 to 1. 31 pg, a 3. 1-fold span. Genome sizes of the fireflies varied within the tested subfamilies and genera. Lamprigera and Pyrocoelia species had large and small genome sizes, respectively. No correlation was found between genome size and morphological traits such as body length, body width, eye width, and antennal length. Our data provide additional information on genome size estimation of the firefly family Lampyridae. Furthermore, this study will help clarify the cost and difficulty of genome sequencing programs for non-model organisms and will help promote studies on firefly genome evolution.

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.

Genome size in arthropods; different roles of phylogeny, habitat and life history in insects and crustaceans

Despite the major role of genome size for physiology, ecology, and evolution, there is still mixed evidence with regard to proximate and ultimate drivers. The main causes of large genome size are proliferation of noncoding elements and/or duplication events. The relative role and interplay between these proximate causes and the evolutionary patterns shaped by phylogeny, life history traits or environment are largely unknown for the arthropods. Genome size shows a tremendous variability in this group, and it has a major impact on a range of fitness-related parameters such as growth, metabolism, life history traits, and for many species also body size. In this study, we compared genome size in two major arthropod groups, insects and crustaceans, and related this to phylogenetic patterns and parameters affecting ambient temperature (latitude, depth, or altitude), insect developmental mode, as well as crustacean body size and habitat, for species where data were available. For the insects, the genome size is clearly phylogeny-dependent, reflecting primarily their life history and mode of development, while for crustaceans there was a weaker association between genome size and phylogeny, suggesting life cycle strategies and habitat as more important determinants. Maximum observed latitude and depth, and their combined effect, showed positive, and possibly phylogenetic independent, correlations with genome size for crustaceans. This study illustrate the striking difference in genome sizes both between and within these two major groups of arthropods, and that while living in the cold with low developmental rates may promote large genomes in marine crustaceans, there is a multitude of proximate and ultimate drivers of genome size.

Genome Size in North American Fireflies: Substantial Variation Likely Driven by Neutral Processes

Eukaryotic genomes show tremendous size variation across taxa. Proximate explanations for genome size variation include differences in ploidy and amounts of noncoding DNA, especially repetitive DNA. Ultimate explanations include selection on physiological correlates of genome size such as cell size, which in turn influence body size, resulting in the often-observed correlation between body size and genome size. In this study, we examined body size and repetitive DNA elements in relationship to the evolution of genome size in North American representatives of a single beetle family, the Lampyridae (fireflies). The 23 species considered represent an excellent study system because of the greater than 5-fold range of genome sizes, documented here using flow cytometry, and the 3-fold range in body size, measured using pronotum width. We also identified common genomic repetitive elements using low-coverage sequencing. We found a positive relationship between genome size and repetitive DNA, particularly retrotransposons. Both genome size and these elements were evolving as expected given phylogenetic relatedness. We also tested whether genome size varied with body size and found no relationship. Together, our results suggest that genome size is evolving neutrally in fireflies.

Transposable elements in Drosophila

Transposable elements (TEs) are mobile genetic elements that can mobilize within host genomes. As TEs comprise more than 40% of the human genome and are linked to numerous diseases, understanding their mechanisms of mobilization and regulation is important. Drosophila melanogaster is an ideal model organism for the study of eukaryotic TEs as its genome contains a diverse array of active TEs. TEs universally impact host genome size via transposition and deletion events, but may also adopt unique functional roles in host organisms. There are 2 main classes of TEs: DNA transposons and retrotransposons. These classes are further divided into subgroups of TEs with unique structural and functional characteristics, demonstrating the significant variability among these elements. Despite this variability, D. melanogaster and other eukaryotic organisms utilize conserved mechanisms to regulate TEs. This review focuses on the transposition mechanisms and regulatory pathways of TEs, and their functional roles in D. melanogaster.

Single-Molecule Sequencing of the Drosophila serrata Genome

Long-read sequencing technology promises to greatly enhance de novo assembly of genomes for nonmodel species. Although the error rates of long reads have been a stumbling block, sequencing at high coverage permits the self-correction of many errors. Here, we sequence and de novo assemble the genome of Drosophila serrata, a species from the montium subgroup that has been well-studied for latitudinal clines, sexual selection, and gene expression, but which lacks a reference genome. Using 11 PacBio single-molecule real-time (SMRT cells), we generated 12 Gbp of raw sequence data comprising ∼65 × whole-genome coverage. Read lengths averaged 8940 bp (NRead50 12,200) with the longest read at 53 kbp. We self-corrected reads using the PBDagCon algorithm and assembled the genome using the MHAP algorithm within the PBcR assembler. Total genome length was 198 Mbp with an N50 just under 1 Mbp. Contigs displayed a high degree of chromosome arm-level conservation with the D. melanogaster genome and many could be sensibly placed on the D. serrata physical map. We also provide an initial annotation for this genome using in silico gene predictions that were supported by RNA-seq data.

The mode and tempo of genome size evolution in the subgenus Sophophora

Genome size varies widely across organisms, with no apparent tie to organismal complexity. While genome size is inherited, there is no established evolutionary model for this trait. Hypotheses have been postulated for the observed variation in genome sizes across species, most notably the effective population size hypothesis, the mutational equilibrium hypothesis, and the adaptive hypothesis. While much data has been collected on genome size, the above hypotheses have largely ignored impacts from phylogenetic relationships. In order to test these competing hypotheses, genome sizes of 87 Sophophora species were analyzed in a comparative phylogenetic approach using Pagel’s parameters of evolution, Blomberg’s K, Abouheif’s C_mean and Moran’s I. In addition to testing the mode and rate of genome size evolution in Sophophora species, the effect of number of taxa on detection of phylogenetic signal was analyzed for each of these comparative phylogenetic methods. Sophophora genome size was found to be dependent on the phylogeny, indicating that evolutionary time was important for predicting the variation among species. Genome size was found to evolve gradually on branches of the tree, with a rapid burst of change early in the phylogeny. These results suggest that Sophophora genome size has experienced gradual changes, which support the largely theoretical mutational equilibrium hypothesis. While some methods (Abouheif’s C_mean and Moran’s I) were found to be affected by increasing taxa numbers, more commonly used methods (λ and Blomberg’s K) were found to have increasing reliability with increasing taxa number, with significantly more support with fifteen or more taxa. Our results suggest that these comparative phylogenetic methods, with adequate taxon sampling, can be a powerful way to uncover the enigma that is genome size variation through incorporation of phylogenetic relationships.

Variation in genome size and karyotype among closely related aphid parasitoids (Hymenoptera, Aphelinidae)

Genome sizes were measured and determined for the karyotypes of nine species of aphid parasitoids in the genus Aphelinus Dalman,1820. Large differences in genome size and karyotype were found between Aphelinus species, which is surprising given the similarity in their morphology and life history. Genome sizes estimated from flow cytometry were larger for species in the Aphelinus mali (Haldeman, 1851) complex than those for the species in the Aphelinus daucicola Kurdjumov, 1913 and Aphelinus varipes (Förster,1841) complexes. Haploid karyotypes of the Aphelinus daucicola and Aphelinus mali complexes comprised five metacentric chromosomes of similar size, whereas those of the Aphelinus varipes complex had four chromosomes, including a larger and a smaller metacentric chromosome and two small acrocentric chromosomes or a large metacentric and three smaller acrocentric chromosomes. Total lengths of female haploid chromosome sets correlated with genome sizes estimated from flow cytometry. Phylogenetic analysis of karyotypic variation revealed a chromosomal fusion together with pericentric inversions in the common ancestor of the Aphelinus varipes complex and further pericentric inversions in the clade comprising Aphelinus kurdjumovi Mercet, 1930 and Aphelinus hordei Kurdjumov, 1913. Fluorescence in situ hybridization with a 28S ribosomal DNA probe revealed a single site on chromosomes of the haploid karyotype of Aphelinus coreae Hopper & Woolley, 2012. The differences in genome size and total chromosome length between species complexes matched the phylogenetic divergence between them.