Intrapopulation genome size variation in D. melanogaster reflects life history variation and plasticity

Genetic variation in recalcitrant repetitive regions of the Drosophila melanogaster genome

Many essential functions of organisms are encoded in highly repetitive genomic regions, including histones involved in DNA packaging, centromeres that are core components of chromosome segregation, ribosomal RNA comprising the protein translation machinery, telomeres that ensure chromosome integrity, piRNA clusters encoding host defenses against selfish elements, and virtually the entire Y chromosome. These regions, formed by highly similar tandem arrays, pose significant challenges for experimental and informatic study, impeding sequence-level descriptions essential for understanding genetic variation. Here, we report the assembly and variation analysis of such repetitive regions in Drosophila melanogaster, offering significant improvements to the existing community reference assembly. Our work successfully recovers previously elusive segments, including complete reconstructions of the histone locus and the pericentric heterochromatin of the X chromosome, spanning the Stellate locus to the distal flank of the rDNA cluster. To infer structural changes in these regions where alignments are often not practicable, we introduce landmark anchors based on unique variants that are putatively orthologous. These regions display considerable structural variation between different D. melanogaster strains, exhibiting differences in copy number and organization of homologous repeat units between haplotypes. In the histone cluster, although we observe minimal genetic exchange indicative of crossing over, the variation patterns suggest mechanisms such as unequal sister chromatid exchange. We also examine the prevalence and scale of concerted evolution in the histone and Stellate clusters and discuss the mechanisms underlying these observed patterns.

No single rescue recipe: genome complexities modulate insect response to climate change

Declines in insect populations have gained formidable attention. Given their crucial role in the ecosystem, the causes of declining insect populations must be investigated. However, the insect clade has been associated with low extinction and high diversification rates. It is unlikely that insects underwent mass extinctions in the past. However, the pace of current climate change could make insect populations vulnerable to extinction. We propose genome size (GS) and transposable elements (TEs) to be rough estimates to assess extinction risk. Larger GS and/or proliferating TEs have been associated with adaptation in rapid climate change scenarios. We speculate that unstable, stressful environmental conditions are strongly associated with GS and TE expansion, which could be further correlated with adaptations. Alternately, stressful conditions trigger TE bursts that are not purged in smaller populations. GS and TE loads could be indicators of small effective populations in the wild, likely experiencing bottlenecks or drastic climatic perturbations, which calls for an urgent assessment of extinction risk.

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 size variation and evolution during invasive range expansion in an introduced plant

Plants demonstrate exceptional variation in genome size across species, and their genome sizes can also vary dramatically across individuals and populations within species. This aspect of genetic variation can have consequences for traits and fitness, but few studies attributed genome size differentiation to ecological and evolutionary processes. Biological invasions present particularly useful natural laboratories to infer selective agents that might drive genome size shifts across environments and population histories. Here, we test hypotheses for the evolutionary causes of genome size variation across 14 invading populations of yellow starthistle, Centaurea solstitialis, in California, United States. We use a survey of genome sizes and trait variation to ask: (1) Is variation in genome size associated with developmental trait variation? (2) Are genome sizes smaller toward the leading edge of the expansion, consistent with selection for "colonizer" traits? Or alternatively, does genome size increase toward the leading edge of the expansion, consistent with predicted consequences of founder effects and drift? (3) Finally, are genome sizes smaller at higher elevations, consistent with selection for shorter development times? We found that 2C DNA content varied 1.21-fold among all samples, and was associated with flowering time variation, such that plants with larger genomes reproduced later, with lower lifetime capitula production. Genome sizes increased toward the leading edge of the invasion, but tended to decrease at higher elevations, consistent with genetic drift during range expansion but potentially strong selection for smaller genomes and faster development time at higher elevations. These results demonstrate how genome size variation can contribute to traits directly tied to reproductive success, and how selection and drift can shape that variation. We highlight the influence of genome size on dynamics underlying a rapid range expansion in a highly problematic invasive plant.

Pleiotropic fitness effects across sexes and ages in the Drosophila genome and transcriptome

Selection varies between categories of individuals, with far-reaching ramifications: Sex-specific selection can impede or accelerate adaptation, and differences in selection between young and old individuals are ultimately responsible for senescence. Here, we measure early- and late-life fitness in adults of both sexes from the Drosophila genetic reference panel and perform quantitative genetic and transcriptomic analyses. Fitness was heritable, showed positive pleiotropy across sexes and age classes, and appeared to be influenced by very large numbers of loci with small effects plus a smaller number with moderate effects. Most loci affected male and female fitness in the same direction; relatively few candidate sexually antagonistic loci were found, though these were enriched on the X chromosome as predicted by theory. The expression level of many genes showed an opposite correlation with fitness in males and females, consistent with unresolved sexual conflict over transcription. The load of deleterious mutations correlated negatively with fitness across genotypes, and we found some evidence for the mutation accumulation (but not the antagonistic pleiotropy) theory of aging.

Identification of upper thermal thresholds during development in the endangered Nechako white sturgeon with management implications for a regulated river

Climate change-induced warming effects are already evident in river ecosystems, and projected increases in temperature will continue to amplify stress on fish communities. In addition, many rivers globally are impacted by dams, which have many negative effects on fishes by altering flow, blocking fish passage, and changing sediment composition. However, in some systems, dams present an opportunity to manage river temperature through regulated releases of cooler water. For example, there is a government mandate for Kenney dam operators in the Nechako river, British Columbia, Canada, to maintain river temperature <20°C in July and August to protect migrating sockeye salmon (Oncorhynchus nerka). However, there is another endangered fish species inhabiting the same river, Nechako white sturgeon (Acipenser transmontanus), and it is unclear if these current temperature regulations, or timing of the regulations, are suitable for spawning and developing sturgeon. In this study, we aimed to identify upper thermal thresholds in white sturgeon embryos and larvae to investigate if exposure to current river temperatures are playing a role in recruitment failure. We incubated embryos and yolk-sac larvae in three environmentally relevant temperatures (14, 18 and 21°C) throughout development to identify thermal thresholds across different levels of biological organization. Our results demonstrate upper thermal thresholds at 21°C across physiological measurements in embryo and yolk-sac larvae white sturgeon. Before hatch, both embryo survival and metabolic rate were reduced at 21°C. After hatch, sublethal consequences continued at 21°C because larval sturgeon had decreased thermal plasticity and a dampened transcriptional response during development. In recent years, the Nechako river has reached 21°C by the end of June, and at this temperature, a decrease in sturgeon performance is evident in most of the traits measured. As such, the thermal thresholds identified here suggest current temperature regulations may not be suitable for developing white sturgeon and future recruitment.

Larger genomes show improved buffering of adult fitness against environmental stress in seed beetles

Our general understanding of the evolution of genome size (GS) is incomplete, and it has long been clear that GS does not reflect organismal complexity. Here, we assess the hypothesis that larger genomes may allow organisms to better cope with environmental variation. It is, for example, possible that genome expansion due to proliferation of transposable elements or gene duplications may affect the ability to regulate and fine-tune transcriptional profiles. We used 18 populations of the seed beetle Callosobruchus maculatus, which differ in GS by up to 4.5%, and exposed adults and juveniles to environmental stress in a series of experiments where stage-specific fitness was assayed. We found that populations with larger genomes were indeed better buffered against environmental stress for adult, but not for juvenile, fitness. The genetic correlation across populations between GS and canalization of adult fitness is consistent with a role for natural selection in the evolution of GS.

Transposable element accumulation drives size differences among polymorphic Y Chromosomes in Drosophila

Y Chromosomes of many species are gene poor and show low levels of nucleotide variation, yet often display high amounts of structural diversity. Dobzhansky cataloged several morphologically distinct Y Chromosomes in Drosophila pseudoobscura that differ in size and shape, but the molecular causes of their dramatic size differences are unclear. Here we use cytogenetics and long-read sequencing to study the sequence content of polymorphic Y Chromosomes in D. pseudoobscura We show that Y Chromosomes differ almost 2-fold in size, ranging from 30 to 60 Mb. Most of this size difference is caused by a handful of active transposable elements (TEs) that have recently expanded on the largest Y Chromosome, with different elements being responsible for Y expansion on differently sized D. pseudoobscura Y's. We show that Y Chromosomes differ in their heterochromatin enrichment, expression of Y-enriched TEs, and also influence expression of dozens of autosomal and X-linked genes. The same helitron element that showed the most drastic amplification on the largest Y in D. pseudoobscura independently amplified on a polymorphic large Y Chromosome in D. affinis, suggesting that some TEs are inherently more prone to become deregulated on Y Chromosomes.

Linking genome size variation to population phenotypic variation within the rotifer, Brachionus asplanchnoidis

Eukaryotic organisms usually contain much more genomic DNA than expected from their biological complexity. In explaining this pattern, selection-based hypotheses suggest that genome size evolves through selection acting on correlated life history traits, implicitly assuming the existence of phenotypic effects of (extra) genomic DNA that are independent of its information content. Here, we present conclusive evidence of such phenotypic effects within a well-mixed natural population that shows heritable variation in genome size. We found that genome size is positively correlated with body size, egg size, and embryonic development time in a population of the monogonont rotifer Brachionus asplanchnoidis. The effect on embryonic development time was mediated partly by an indirect effect (via egg size), and a direct effect, the latter indicating an increased replication cost of the larger amounts of DNA during mitosis. Our results suggest that selection-based change of genome size can operate in this population, provided it is strong enough to overcome drift or mutational change of genome size.

The Drosophila Y Chromosome Affects Heterochromatin Integrity Genome-Wide

The Drosophila Y chromosome is gene poor and mainly consists of silenced, repetitive DNA. Nonetheless, the Y influences expression of hundreds of genes genome-wide, possibly by sequestering key components of the heterochromatin machinery away from other positions in the genome. To test the influence of the Y chromosome on the genome-wide chromatin landscape, we assayed the genomic distribution of histone modifications associated with gene activation (H3K4me3) or heterochromatin (H3K9me2 and H3K9me3) in fruit flies with varying sex chromosome complements (X0, XY, and XYY males; XX and XXY females). Consistent with the general deficiency of active chromatin modifications on the Y, we find that Y gene dose has little influence on the genomic distribution of H3K4me3. In contrast, both the presence and the number of Y chromosomes strongly influence genome-wide enrichment patterns of repressive chromatin modifications. Highly repetitive regions such as the pericentromeres, the dot, and the Y chromosome (if present) are enriched for heterochromatic modifications in wildtype males and females, and even more strongly in X0 flies. In contrast, the additional Y chromosome in XYY males and XXY females diminishes the heterochromatic signal in these normally silenced, repeat-rich regions, which is accompanied by an increase in expression of Y-linked repeats. We find hundreds of genes that are expressed differentially between individuals with aberrant sex chromosome karyotypes, many of which also show sex-biased expression in wildtype Drosophila. Thus, Y chromosomes influence heterochromatin integrity genome-wide, and differences in the chromatin landscape of males and females may also contribute to sex-biased gene expression and sexual dimorphisms.

Excess histone H3 is a competitive Chk1 inhibitor that controls cell-cycle remodeling in the early Drosophila embryo

The DNA damage checkpoint is crucial to protect genome integrity.^1,2 However, the early embryos of many metazoans sacrifice this safeguard to allow for rapid cleavage divisions that are required for speedy development. At the mid-blastula transition (MBT), embryos switch from rapid cleavage divisions to slower, patterned divisions with the addition of gap phases and acquisition of DNA damage checkpoints. The timing of the MBT is dependent on the nuclear-to-cytoplasmic (N/C ratio)^3-7 and the activation of the checkpoint kinase, Chk1.^8-17 How Chk1 activity is coupled to the N/C ratio has remained poorly understood. Here, we show that dynamic changes in histone H3 availability in response to the increasing N/C ratio control Chk1 activity and thus time the MBT in the Drosophila embryo. We show that excess H3 in the early cycles interferes with cell-cycle slowing independent of chromatin incorporation. We find that the N-terminal tail of H3 acts as a competitive inhibitor of Chk1 in vitro and reduces Chk1 activity in vivo. Using a H3-tail mutant that has reduced Chk1 inhibitor activity, we show that the amount of available Chk1 sites in the H3 pool controls the dynamics of cell-cycle progression. Mathematical modeling quantitatively supports a mechanism where titration of H3 during early cleavage cycles regulates Chk1-dependent cell-cycle slowing. This study defines Chk1 regulation by H3 as a key mechanism that coordinates cell-cycle remodeling with developmental progression.

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.

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).

High-resolution single-cell 3D-models of chromatin ensembles during Drosophila embryogenesis

Single-cell chromatin studies provide insights into how chromatin structure relates to functions of individual cells. However, balancing high-resolution and genome wide-coverage remains challenging. We describe a computational method for the reconstruction of large 3D-ensembles of single-cell (sc) chromatin conformations from population Hi-C that we apply to study embryogenesis in Drosophila. With minimal assumptions of physical properties and without adjustable parameters, our method generates large ensembles of chromatin conformations via deep-sampling. Our method identifies specific interactions, which constitute 5-6% of Hi-C frequencies, but surprisingly are sufficient to drive chromatin folding, giving rise to the observed Hi-C patterns. Modeled sc-chromatins quantify chromatin heterogeneity, revealing significant changes during embryogenesis. Furthermore, >50% of modeled sc-chromatin maintain topologically associating domains (TADs) in early embryos, when no population TADs are perceptible. Domain boundaries become fixated during development, with strong preference at binding-sites of insulator-complexes upon the midblastula transition. Overall, high-resolution 3D-ensembles of sc-chromatin conformations enable further in-depth interpretation of population Hi-C, improving understanding of the structure-function relationship of genome organization.

Genome Size of 17 Species From Caelifera (Orthoptera) and Determination of Internal Standards With Very Large Genome Size in Insecta

Comparative studies of insect genome size show that Orthoptera is a unique group of Insecta with a significantly enlarged genome. To determine a suitable internal standard for a large genome and to compare the effects of different internal standards on estimates of genome size, we used four internal standards to estimate nuclear DNA content in nine insect species with large genomes. The results showed that the combination of two internal standards, Locusta migratoria (♂1C = 6.20 pg, ♀1C = 6.60 pg) and Periplaneta americana♂ (1C = 3.41 pg), was suitable for estimating large genome of Caelifera by flow cytometry. Using these two internal standards, we estimated the genome sizes of 17 species of Caelifera (12 genera in Acrididae, 2 genera in Pamphagidae, 1 genus in Pyrgomorphidae) using flow cytometry. Genomes ranged from 6.57 pg (Shirakiacris shirakii) to 18.64 pg (Bryodemella holdereri), the largest described in insects to date. These species showed significant genomic dimorphism based on sex: females had a 0.56 pg larger genome than males on average, which might be due to the sex chromosome determinism mechanism of X0(♂)/XX(♀). To test the results obtained by flow cytometry, we used k-mers of Illumina sequencing data to gauge the C-value of Calliptamus abbreviatus and Haplotropis brunneriana. The results of the two methods are slightly different. Genomes were estimated to be about 0.28 and 0.26 pg smaller, respectively, than the flow cytometry values. Furthermore, we also reconstructed the evolutionary relationships of these taxa and discuss the genome size evolution in a phylogenetic framework.

Effects of acclimation temperature on the thermal physiology in two geographically distinct populations of lake sturgeon (Acipenser fulvescens)

Temperature is one of the most important abiotic factors regulating development and biological processes in ectotherms. By 2050, climate change may result in temperature increases of 2.1-3.4°C in Manitoba, Canada. Lake sturgeon, Acipenser fulvescens, from both northern and southern populations in Manitoba were acclimated to 16, 20 and 24°C for 30 days, after which critical thermal maximum (CT_max) trials were conducted to investigate their thermal plasticity. We also examined the effects of temperature on morphological and physiological indices. Acclimation temperature significantly influenced the CT_max, body mass, hepatosomatic index, metabolic rate and the mRNA expression of transcripts involved in the cellular response to heat shock and hypoxia (HSP70, HSP90a, HSP90b, HIF-1α) in the gill of lake sturgeon. Population significantly affected the above phenotypes, as well as the mRNA expression of Na⁺/K⁺ ATPase-α1 and the hepatic glutathione peroxidase enzyme activity. The southern population had an average CT_max that was 0.71 and 0.45°C higher than the northern population at 20 and 24°C, respectively. Immediately following CT_max trials, mRNA expression of HSP90a and HIF-1α was positively correlated with individual CT_max of lake sturgeon across acclimation treatments and populations (r = 0.7, r = 0.62, respectively; P < 0.0001). Lake sturgeon acclimated to 20 and 24°C had decreased hepatosomatic indices (93 and 244% reduction, respectively; P < 0.0001) and metabolic suppression (27.7 and 42.1% reduction, respectively; P < 0.05) when compared to sturgeon acclimated to 16°C, regardless of population. Glutathione peroxidase activity and mRNA expression Na⁺/K⁺ ATPase-α1 were elevated in the northern relative to the southern population. Acclimation to 24°C also induced mortality in both populations when compared to sturgeon acclimated to 16 and 20°C. Thus, increased temperatures have wide-ranging population-specific physiological consequences for lake sturgeon across biological levels of organization.

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.

The genetic adaptations of fall armyworm Spodoptera frugiperda facilitated its rapid global dispersal and invasion

The fall armyworm (Spodoptera frugiperda) is a lepidopteran insect pest that causes huge economic losses. This notorious insect pest has rapidly spread over the world in the past few years. However, the mechanisms of rapid dispersal are not well understood. Here, we report a chromosome-level assembled genome of the fall armyworm, named the ZJ-version, using PacBio and Hi-C technology. The sequenced individual was a female collected from the Zhejiang province of China and had high heterozygosity. The assembled genome size of ZJ-version was 486 Mb, containing 361 contigs with an N50 of 1.13 Mb. Hi-C scaffolding further assembled the genome into 31 chromosomes and a portion of W chromosome, representing 97.4% of all contigs and resulted in a chromosome-level genome with scaffold N50 of 16.3 Mb. The sex chromosomes were identified by genome resequencing of a single male pupa and a single female pupa. About 28% of the genome was annotated as repeat sequences, and 22,623 protein-coding genes were identified. Comparative genomics revealed the expansion of the detoxification-associated gene families, chemoreception-associated gene families, nutrition metabolism and transport system gene families in the fall armyworm. Transcriptomic and phylogenetic analyses focused on these gene families revealed the potential roles of the genes in polyphagia and invasion of fall armyworm. The high-quality of the fall armyworm genome provides an important genomic resource for further explorations of the mechanisms of polyphagia and insecticide resistance, as well as for pest management of fall armyworm.

The Y chromosome may contribute to sex-specific ageing in Drosophila

Heterochromatin suppresses repetitive DNA, and a loss of heterochromatin has been observed in aged cells of several species, including humans and Drosophila. Males often contain substantially more heterochromatic DNA than females, due to the presence of a large, repeat-rich Y chromosome, and male flies generally have shorter average life spans than females. Here we show that repetitive DNA becomes de-repressed more rapidly in old male flies relative to females, and repeats on the Y chromosome are disproportionally mis-expressed during aging. This is associated with a loss of heterochromatin at repetitive elements during aging in male flies, and a general loss of repressive chromatin in aged males away from pericentromeric regions and the Y. By generating flies with different sex chromosome karyotypes (XXY females; X0 and XYY males), we show that repeat de-repression and average lifespan is correlated with the number of Y chromosomes. This suggests that sex-specific chromatin differences may contribute to sex-specific aging in flies.

Triploidy in zebrafish larvae: Effects on gene expression, cell size and cell number, growth, development and swimming performance

There is renewed interest in the regulation and consequences of cell size adaptations in studies on understanding the ecophysiology of ectotherms. Here we test if induction of triploidy, which increases cell size in zebrafish (Danio rerio), makes for a good model system to study consequences of cell size. Ideally, diploid and triploid zebrafish should differ in cell size, but should otherwise be comparable in order to be suitable as a model. We induced triploidy by cold shock and compared diploid and triploid zebrafish larvae under standard rearing conditions for differences in genome size, cell size and cell number, development, growth and swimming performance and expression of housekeeping genes and hsp70.1. Triploid zebrafish have larger but fewer cells, and the increase in cell size matched the increase in genome size (+ 50%). Under standard conditions, patterns in gene expression, ontogenetic development and larval growth were near identical between triploids and diploids. However, under demanding conditions (i.e. the maximum swimming velocity during an escape response), triploid larvae performed poorer than their diploid counterparts, especially after repeated stimuli to induce swimming. This result is consistent with the idea that larger cells have less capacity to generate energy, which becomes manifest during repeated physical exertion resulting in increased fatigue. Triploidy induction in zebrafish appears a valid method to increase specifically cell size and this provides a model system to test for consequences of cell size adaptation for the energy budget and swimming performance of this ectothermic vertebrate.

The Genetic Basis of Natural Variation in Drosophila melanogaster Immune Defense against Enterococcus faecalis

Dissecting the genetic basis of natural variation in disease response in hosts provides insights into the coevolutionary dynamics of host-pathogen interactions. Here, a genome-wide association study of Drosophila melanogaster survival after infection with the Gram-positive entomopathogenic bacterium Enterococcus faecalis is reported. There was considerable variation in defense against E. faecalis infection among inbred lines of the Drosophila Genetics Reference Panel. We identified single nucleotide polymorphisms associated with six genes with a significant (p < 10^-08, corresponding to a false discovery rate of 2.4%) association with survival, none of which were canonical immune genes. To validate the role of these genes in immune defense, their expression was knocked-down using RNAi and survival of infected hosts was followed, which confirmed a role for the genes krishah and S6k in immune defense. We further identified a putative role for the Bomanin gene BomBc1 (also known as IM23), in E. faecalis infection response. This study adds to the growing set of association studies for infection in Drosophila melanogaster and suggests that the genetic causes of variation in immune defense differ for different pathogens.

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.

Quantification of Proteins and Histone Marks in Drosophila Embryos Reveals Stoichiometric Relationships Impacting Chromatin Regulation

Gene transcription in eukaryotes is regulated through dynamic interactions of a variety of different proteins with DNA in the context of chromatin. Here, we used mass spectrometry for absolute quantification of the nuclear proteome and methyl marks on selected lysine residues in histone H3 during two stages of Drosophila embryogenesis. These analyses provide comprehensive information about the absolute copy number of several thousand proteins and reveal unexpected relationships between the abundance of histone-modifying and -binding proteins and the chromatin landscape that they generate and interact with. For some histone modifications, the levels in Drosophila embryos are substantially different from those previously reported in tissue culture cells. Genome-wide profiling of H3K27 methylation during developmental progression and in animals with reduced PRC2 levels illustrates how mass spectrometry can be used for quantitatively describing and comparing chromatin states. Together, these data provide a foundation toward a quantitative understanding of gene regulation in Drosophila.

Small, but surprisingly repetitive genomes: transposon expansion and not polyploidy has driven a doubling in genome size in a metazoan species complex

BACKGROUND

The causes and consequences of genome size variation across Eukaryotes, which spans five orders of magnitude, have been hotly debated since before the advent of genome sequencing. Previous studies have mostly examined variation among larger taxonomic units (e.g., orders, or genera), while comparisons among closely related species are rare. Rotifers of the Brachionus plicatilis species complex exhibit a seven-fold variation in genome size and thus represent a unique opportunity to study such changes on a relatively short evolutionary timescale. Here, we sequenced and analysed the genomes of four species of this complex with nuclear DNA contents spanning 110-422 Mbp. To establish the likely mechanisms of genome size change, we analysed both sequencing read libraries and assemblies for signatures of polyploidy and repetitive element content. We also compared these genomes to that of B. calyciflorus, the closest relative with a sequenced genome (293 Mbp nuclear DNA content).

RESULTS

Despite the very large differences in genome size, we saw no evidence of ploidy level changes across the B. plicatilis complex. However, repetitive element content explained a large portion of genome size variation (at least 54%). The species with the largest genome, B. asplanchnoidis, has a strikingly high 44% repetitive element content, while the smaller B. plicatilis genomes contain between 14 and 25% repetitive elements. According to our analyses, the B. calyciflorus genome contains 39% repetitive elements, which is substantially higher than previously reported (21%), and suggests that high repetitive element load could be widespread in monogonont rotifers.

CONCLUSIONS

Even though the genome sizes of these species are at the low end of the metazoan spectrum, their genomes contain substantial amounts of repetitive elements. Polyploidy does not appear to play a role in genome size variations in these species, and these variations can be mostly explained by changes in repetitive element content. This contradicts the naïve expectation that small genomes are streamlined, or less complex, and that large variations in nuclear DNA content between closely related species are due to polyploidy.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

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 Estimation and Quantitative Cytogenetics in Insects

With care, it is possible using flow cytometry to create a precise and accurate estimate of the genome size of an insect that is useful for genomics, genetics, molecular/cell biology, or systematics. Genome size estimation is a useful first step in a complete genome sequencing project. The number of sequencing reads required to produce a given level of coverage depends directly upon the 1C amount of DNA per cell, while an even more critical need is an accurate 1C genome size estimate to compare against the final assembly. Here we present a detailed protocol to estimate genome size using flow cytometry. Published genome size estimates should be submitted to genomesize.com so that they are available to all.

A New Standard for Crustacean Genomes: The Highly Contiguous, Annotated Genome Assembly of the Clam Shrimp Eulimnadia texana Reveals HOX Gene Order and Identifies the Sex Chromosome

Vernal pool clam shrimp (Eulimnadia texana) are a promising model system due to their ease of lab culture, short generation time, modest sized genome, a somewhat rare stable androdioecious sex determination system, and a requirement to reproduce via desiccated diapaused eggs. We generated a highly contiguous genome assembly using 46× of PacBio long read data and 216× of Illumina short reads, and annotated using Illumina RNAseq obtained from adult males or hermaphrodites. Of the 120 Mb genome 85% is contained in the largest eight contigs, the smallest of which is 4.6 Mb. The assembly contains 98% of transcripts predicted via RNAseq. This assembly is qualitatively different from scaffolded Illumina assemblies: It is produced from long reads that contain sequence data along their entire length, and is thus gap free. The contiguity of the assembly allows us to order the HOX genes within the genome, identifying two loci that contain HOX gene orthologs, and which approximately maintain the order observed in other arthropods. We identified a partial duplication of the Antennapedia complex adjacent to the few genes homologous to the Bithorax locus. Because the sex chromosome of an androdioecious species is of special interest, we used existing allozyme and microsatellite markers to identify the E. texana sex chromosome, and find that it comprises nearly half of the genome of this species. Linkage patterns indicate that recombination is extremely rare and perhaps absent in hermaphrodites, and as a result the location of the sex determining locus will be difficult to refine using recombination mapping.

The Sleep Inbred Panel, a Collection of Inbred Drosophila melanogaster with Extreme Long and Short Sleep Duration

Understanding how genomic variation causes differences in observable phenotypes remains a major challenge in biology. It is difficult to trace the sequence of events originating from genomic variants to changes in transcriptional responses or protein modifications. Ideally, one would conduct experiments with individuals that are at either extreme of the trait of interest, but such resources are often not available. Further, advances in genome editing will enable testing of candidate polymorphisms individually and in combination. Here we have created a resource for the study of sleep with 39 inbred lines of Drosophila-the Sleep Inbred Panel (SIP). SIP lines have stable long- and short-sleeping phenotypes developed from naturally occurring polymorphisms. These lines are fully sequenced, enabling more accurate targeting for genome editing and transgenic constructs. This panel facilitates the study of intermediate transcriptional and proteomic correlates of sleep, and supports genome editing studies to verify polymorphisms associated with sleep duration.

Genotype Influences Day-to-Day Variability in Sleep in Drosophila melanogaster

Patterns of sleep often vary among individuals. But sleep and activity may also vary within an individual, fluctuating in pattern across time. One possibility is that these daily fluctuations in sleep are caused by the underlying genotype of the individual. However, differences attributable to genetic causes are difficult to distinguish from environmental factors in outbred populations such as humans. We therefore employed Drosophila as a model of intra-individual variability in sleep using previously collected sleep and activity data from the Drosophila Genetic Reference Panel, a collection of wild-derived inbred lines. Individual flies had significant daily fluctuations in their sleep patterns, and these fluctuations were heritable. Using the standard deviation of sleep parameters as a metric, we conducted a genome-wide association study. We found 663 polymorphisms in 104 genes associated with daily fluctuations in sleep. We confirmed the effects of 12 candidate genes on the standard deviation of sleep parameters. Our results suggest that daily fluctuations in sleep patterns are due in part to gene activity.

Dynamic Control of dNTP Synthesis in Early Embryos

Exponential increase of cell numbers in early embryos requires large amounts of DNA precursors (deoxyribonucleoside triphosphates (dNTPs)). Little is understood about how embryos satisfy this demand. We examined dNTP metabolism in the early Drosophila embryo, in which gastrulation is preceded by 13 sequential nuclear cleavages within only 2 hr of fertilization. Surprisingly, despite the breakneck speed at which Drosophila embryos synthesize DNA, maternally deposited dNTPs can generate less than half of the genomes needed to reach gastrulation. The rest of the dNTPs are synthesized "on the go." The rate-limiting enzyme of dNTP synthesis, ribonucleotide reductase, is inhibited by endogenous levels of deoxyATP (dATP) present at fertilization and is activated as dATP is depleted via DNA polymerization. This feedback inhibition renders the concentration of dNTPs at gastrulation robust, with respect to large variations in maternal supplies, and is essential for normal progression of embryogenesis.

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.

An Enhancer's Length and Composition Are Shaped by Its Regulatory Task

Enhancers drive the gene expression patterns required for virtually every process in metazoans. We propose that enhancer length and transcription factor (TF) binding site composition—the number and identity of TF binding sites—reflect the complexity of the enhancer's regulatory task. In development, we define regulatory task complexity as the number of fates specified in a set of cells at once. We hypothesize that enhancers with more complex regulatory tasks will be longer, with more, but less specific, TF binding sites. Larger numbers of binding sites can be arranged in more ways, allowing enhancers to drive many distinct expression patterns, and therefore cell fates, using a finite number of TF inputs. We compare ~100 enhancers patterning the more complex anterior-posterior (AP) axis and the simpler dorsal-ventral (DV) axis in Drosophila and find that the AP enhancers are longer with more, but less specific binding sites than the (DV) enhancers. Using a set of ~3,500 enhancers, we find enhancer length and TF binding site number again increase with increasing regulatory task complexity. Therefore, to be broadly applicable, computational tools to study enhancers must account for differences in regulatory task.

Pleiotropy, constraint, and modularity in the evolution of life histories: insights from genomic analyses

Multicellular organisms display an enormous range of life history (LH) strategies and present an evolutionary conundrum; despite strong natural selection, LH traits are characterized by high levels of genetic variation. To understand the evolution of life histories and maintenance of this variation, the specific phenotypic effects of segregating alleles and the genetic networks in which they act need to be elucidated. In particular, the extent to which LH evolution is constrained by the pleiotropy of alleles contributing to LH variation is generally unknown. Here, we review recent empirical results that shed light on this question, with an emphasis on studies employing genomic analyses. While genome-scale analyses are increasingly practical and affordable, they face limitations of genetic resolution and statistical power. We describe new research approaches that we believe can produce new insights and evaluate their promise and applicability to different kinds of organisms. Two approaches seem particularly promising: experiments that manipulate selection in multiple dimensions and measure phenotypic and genomic response and analytical approaches that take into account genome-wide associations between markers and phenotypes, rather than applying a traditional marker-by-marker approach.

Genome sequence of a rice pest, the white-backed planthopper (Sogatella furcifera)

Background

Sogatella furcifera is an important phloem sap-sucking and plant virus-transmitting migratory insect of rice. Because of its high reproductive potential, dispersal capability and transmission of plant viral diseases, S. furcifera causes considerable damage to rice grain production and has great economical and agricultural impacts. Comprehensive studies into ecological aspects and virus-host interactions of S. furcifera have been limited because of the lack of a well-assembled genome sequence.

Findings

A total of 241.3 Gb of raw reads from the whole genome of S. furcifera were generated by Illumina sequencing using different combinations of mate-pair and paired-end libraries from 17 insert libraries ranging between 180 bp and 40 kbp. The final genome assembly (0.72 Gb), with average N50 contig size of 70.7 kb and scaffold N50 of 1.18 Mb, covers 98.6 % of the estimated genome size of S. furcifera . Genome annotation, assisted by eight different developmental stages (embryos, 1 st -5 th instar nymphs, 5-day-old adults and 10-day-old adults), generated 21 254 protein-coding genes, which captured 99.59 % (247/248) of core CEGMA genes and 91.7 % (2453/2675) of BUSCO genes.

Conclusions

We report the first assembled and annotated whole genome sequence and transcriptome of S. furcifera . The assembled draft genome of S. furcifera will be a valuable resource for ecological and virus-host interaction studies of this pest.

Effects of larval crowding on quantitative variation for development time and viability in Drosophila melanogaster

Competition between individuals belonging to the same species is a universal feature of natural populations and is the process underpinning organismal adaptation. Despite its importance, still comparatively little is known about the genetic variation responsible for competitive traits. Here, we measured the phenotypic variation and quantitative genetics parameters for two fitness-related traits-egg-to-adult viability and development time-across a panel of Drosophila strains under varying larval densities. Both traits exhibited substantial genetic variation at all larval densities, as well as significant genotype-by-environment interactions (GEIs). GEI was attributable to changes in the rank order of reaction norms for both traits, and additionally to differences in the between-line variance for development time. The coefficient of genetic variation increased under stress conditions for development time, while it was higher at both high and low densities for viability. While development time also correlated negatively with fitness at high larval densities-meaning that fast developers have high fitness-there was no correlation with fitness at low density. This result suggests that GEI may be a common feature of fitness-related genetic variation and, further, that trait values under noncompetitive conditions could be poor indicators of individual fitness. The latter point could have significant implications for animal and plant breeding programs, as well as for conservation genetics.

A novel method for quantifying the rate of embryogenesis uncovers considerable genetic variation for the duration of embryonic development in Drosophila melanogaster

BACKGROUND

Embryogenesis is a highly conserved, canalized process, and variation in the duration of embryogenesis (DOE), i.e., time from egg lay to hatching, has a potentially profound effect on the outcome of within- and between-species competition. There is both intra- and inter-specific variation in this trait, which may provide important fuel for evolutionary processes, particularly adaptation. However, while genetic variation underlying simpler morphological traits, or with large phenotypic effects is well described in the literature, less is known about the underlying genetics of traits, such as DOE, partly due to a lack of tools with which to study them.

RESULTS

Here, we establish a novel microscope-based assay to survey genetic variation for the duration of embryogenesis (DOE). First, to establish the potential importance of DOE in competitive fitness, we performed a set of experiments where we experimentally manipulated the time until hatching, and show that short hatching times result in priority effect in the form of improved larval competitive ability. We then use our assay to measure DOE for 43 strains from the Drosophila Genetic Reference Panel (DGRP). Our assay greatly simplifies the measurement of DOE, making it possible to precisely quantify this trait for 59,295 individual embryos (mean ± S.D. of 1103 ± 293 per DGRP strain, and 1002 ± 203 per control). We find extensive genetic variation in DOE, with a 15 % difference in rate between the slowest and fastest strains measured, and 89 % of phenotypic variation due to DGRP strain. Using sequence information from the DGRP, we perform a genome-wide association study, which suggests that some well-known developmental genes affect the speed of embryonic development.

CONCLUSIONS

We showed that the duration of embryogenesis (DOE) can be efficiently and precisely measured in Drosophila, and that the DGRP strains show remarkable variation in DOE. A genome-wide analysis suggests that some well-known developmental genes are potentially associated with DOE. Further functional assays, or transcriptomic analysis of embryos from the DGRP, can validate the role of our candidates in early developmental processes.

Genome Stability and Evolution: Attempting a Holistic View

The reason why the DNA content, chromosome number and shape, and gene content of eukaryotic genomes vary independently remains a matter of speculation. The same is true for the questions of whether there is a general tendency for increase or decrease of genome size and chromosome number and whether genome size and/or chromosome number have an adaptive value and, if so, what this value is. Here we assume that three strategies of genome evolution (shrinkage, expansion, and equilibrium) have developed to find the optimal balance between genomic stability and plasticity. We suggest various modes of DNA double-strand break (DSB) repair in combination with whole-genome duplication (WGD) and dysploid chromosome number alteration to explain the different strategies of genome size and karyotype evolution.

Genome size correlates with reproductive fitness in seed beetles

The ultimate cause of genome size (GS) evolution in eukaryotes remains a major and unresolved puzzle in evolutionary biology. Large-scale comparative studies have failed to find consistent correlations between GS and organismal properties, resulting in the 'C-value paradox'. Current hypotheses for the evolution of GS are based either on the balance between mutational events and drift or on natural selection acting upon standing genetic variation in GS. It is, however, currently very difficult to evaluate the role of selection because within-species studies that relate variation in life-history traits to variation in GS are very rare. Here, we report phylogenetic comparative analyses of GS evolution in seed beetles at two distinct taxonomic scales, which combines replicated estimation of GS with experimental assays of life-history traits and reproductive fitness. GS showed rapid and bidirectional evolution across species, but did not show correlated evolution with any of several indices of the relative importance of genetic drift. Within a single species, GS varied by 4-5% across populations and showed positive correlated evolution with independent estimates of male and female reproductive fitness. Collectively, the phylogenetic pattern of GS diversification across and within species in conjunction with the pattern of correlated evolution between GS and fitness provide novel support for the tenet that natural selection plays a key role in shaping GS evolution.

Patterns of genome size variation in snapping shrimp

Although crustaceans vary extensively in genome size, little is known about how genome size may affect the ecology and evolution of species in this diverse group, in part due to the lack of large genome size datasets. Here we investigate interspecific, intraspecific, and intracolony variation in genome size in 39 species of Synalpheus shrimps, representing one of the largest genome size datasets for a single genus within crustaceans. We find that genome size ranges approximately 4-fold across Synalpheus with little phylogenetic signal, and is not related to body size. In a subset of these species, genome size is related to chromosome size, but not to chromosome number, suggesting that despite large genomes, these species are not polyploid. Interestingly, there appears to be 35% intraspecific genome size variation in Synalpheus idios among geographic regions, and up to 30% variation in Synalpheus duffyi genome size within the same colony.

Mitonuclear Epistasis for Development Time and Its Modification by Diet in Drosophila

Mitochondrial (mtDNA) and nuclear genes have to operate in a coordinated manner to maintain organismal function, and the regulation of this homeostasis presents a substantial source of potential epistatic (G × G) interactions. How these interactions shape the fitness landscape is poorly understood. Here we developed a novel mitonuclear epistasis model, using selected strains of the Drosophila Genetic Reference Panel (DGRP) and mitochondrial genomes from within Drosophila melanogaster and D. simulans to test the hypothesis that mtDNA × nDNA interactions influence fitness. In total we built 72 genotypes (12 nuclear backgrounds × 6 mtDNA haplotypes, with 3 from each species) to dissect the relationship between genotype and phenotype. Each genotype was assayed on four food environments. We found considerable variation in several phenotypes, including development time and egg-to-adult viability, and this variation was partitioned into genetic (G), environmental (E), and higher-order (G × G, G × E, and G × G × E) components. Food type had a significant impact on development time and also modified mitonuclear epistases, evidencing a broad spectrum of G × G × E across these genotypes. Nuclear background effects were substantial, followed by mtDNA effects and their G × G interaction. The species of mtDNA haplotype had negligible effects on phenotypic variation and there was no evidence that mtDNA variation has different effects on male and female fitness traits. Our results demonstrate that mitonuclear epistases are context dependent, suggesting the selective pressure acting on mitonuclear genotypes may vary with food environment in a genotype-specific manner.

Genome-Wide Mapping Targets of the Metazoan Chromatin Remodeling Factor NURF Reveals Nucleosome Remodeling at Enhancers, Core Promoters and Gene Insulators

NURF is a conserved higher eukaryotic ISWI-containing chromatin remodeling complex that catalyzes ATP-dependent nucleosome sliding. By sliding nucleosomes, NURF is able to alter chromatin dynamics to control transcription and genome organization. Previous biochemical and genetic analysis of the specificity-subunit of Drosophila NURF (Nurf301/Enhancer of Bithorax (E(bx)) has defined NURF as a critical regulator of homeotic, heat-shock and steroid-responsive gene transcription. It has been speculated that NURF controls pathway specific transcription by co-operating with sequence-specific transcription factors to remodel chromatin at dedicated enhancers. However, conclusive in vivo demonstration of this is lacking and precise regulatory elements targeted by NURF are poorly defined. To address this, we have generated a comprehensive map of in vivo NURF activity, using MNase-sequencing to determine at base pair resolution NURF target nucleosomes, and ChIP-sequencing to define sites of NURF recruitment. Our data show that, besides anticipated roles at enhancers, NURF interacts physically and functionally with the TRF2/DREF basal transcription factor to organize nucleosomes downstream of active promoters. Moreover, we detect NURF remodeling and recruitment at distal insulator sites, where NURF functionally interacts with and co-localizes with DREF and insulator proteins including CP190 to establish nucleosome-depleted domains. This insulator function of NURF is most apparent at subclasses of insulators that mark the boundaries of chromatin domains, where multiple insulator proteins co-associate. By visualizing the complete repertoire of in vivo NURF chromatin targets, our data provide new insights into how chromatin remodeling can control genome organization and regulatory interactions.

The Adaptive Significance of Natural Genetic Variation in the DNA Damage Response of Drosophila melanogaster

Despite decades of work, our understanding of the distribution of fitness effects of segregating genetic variants in natural populations remains largely incomplete. One form of selection that can maintain genetic variation is spatially varying selection, such as that leading to latitudinal clines. While the introduction of population genomic approaches to understanding spatially varying selection has generated much excitement, little successful effort has been devoted to moving beyond genome scans for selection to experimental analysis of the relevant biology and the development of experimentally motivated hypotheses regarding the agents of selection; it remains an interesting question as to whether the vast majority of population genomic work will lead to satisfying biological insights. Here, motivated by population genomic results, we investigate how spatially varying selection in the genetic model system, Drosophila melanogaster, has led to genetic differences between populations in several components of the DNA damage response. UVB incidence, which is negatively correlated with latitude, is an important agent of DNA damage. We show that sensitivity of early embryos to UVB exposure is strongly correlated with latitude such that low latitude populations show much lower sensitivity to UVB. We then show that lines with lower embryo UVB sensitivity also exhibit increased capacity for repair of damaged sperm DNA by the oocyte. A comparison of the early embryo transcriptome in high and low latitude embryos provides evidence that one mechanism of adaptive DNA repair differences between populations is the greater abundance of DNA repair transcripts in the eggs of low latitude females. Finally, we use population genomic comparisons of high and low latitude samples to reveal evidence that multiple components of the DNA damage response and both coding and non-coding variation likely contribute to adaptive differences in DNA repair between populations.

Drosophila Females Undergo Genome Expansion after Interspecific Hybridization

Genome size (or C-value) can present a wide range of values among eukaryotes. This variation has been attributed to differences in the amplification and deletion of different noncoding repetitive sequences, particularly transposable elements (TEs). TEs can be activated under different stress conditions such as interspecific hybridization events, as described for several species of animals and plants. These massive transposition episodes can lead to considerable genome expansions that could ultimately be involved in hybrid speciation processes. Here, we describe the effects of hybridization and introgression on genome size of Drosophila hybrids. We measured the genome size of two close Drosophila species, Drosophila buzzatii and Drosophila koepferae, their F1 offspring and the offspring from three generations of backcrossed hybrids; where mobilization of up to 28 different TEs was previously detected. We show that hybrid females indeed present a genome expansion, especially in the first backcross, which could likely be explained by transposition events. Hybrid males, which exhibit more variable C-values among individuals of the same generation, do not present an increased genome size. Thus, we demonstrate that the impact of hybridization on genome size can be detected through flow cytometry and is sex-dependent.

Comparison of life history and genetic properties of cowpea bruchid strains and their response to hypoxia

The cowpea bruchid (Callosobruchus maculatus) is the most important storage pest of grain legumes and comprises geographically distinct strains. Storage under a modified atmosphere with decreased O2 content represents an alternative to chemical fumigants for pest control of stored grains. In this study, we compared reproduction, development and survival, as well as genome size of bruchid strains from South India (SI), Burkina Faso (BF), Niger (CmNnC) and the United States (OH), reared on mung bean (Vigna radiata). Fecundity and egg-to-adult duration varied significantly among these strains. Notably, strain BF had the highest fecundity, and strain SI displayed the fastest development whereas strain OH was the slowest. Differences in adult lifespan among strains were only detected in unmated but not in the mated group. Genome size of SI females was significantly larger than that of OH females, and for all four strains, the female genomes were larger than those of their corresponding males. Furthermore, we studied effects of exposure to 1% O2+99% N2 on strains SI and BF. Mortality caused by hypoxia was influenced by not only developmental stage but also by insect strain. Eggs were most sensitive, particularly at the early stage, whereas the 3rd and 4th instar larvae were most tolerant and could survive up to 15 days of low O2. Strain SI was slightly more resistant than BF in egg and larval stages. Proteolytic activity prior to, during and after hypoxia treatment revealed remarkable metabolic plasticity of cowpea bruchids in response to modified atmosphere.

Correlated variation and population differentiation in satellite DNA abundance among lines of Drosophila melanogaster

Tandemly repeating satellite DNA elements in heterochromatin occupy a substantial portion of many eukaryotic genomes. Although often characterized as genomic parasites deleterious to the host, they also can be crucial for essential processes such as chromosome segregation. Adding to their interest, satellite DNA elements evolve at high rates; among Drosophila, closely related species often differ drastically in both the types and abundances of satellite repeats. However, due to technical challenges, the evolutionary mechanisms driving this rapid turnover remain unclear. Here we characterize natural variation in simple-sequence repeats of 2-10 bp from inbred Drosophila melanogaster lines derived from multiple populations, using a method we developed called k-Seek that analyzes unassembled Illumina sequence reads. In addition to quantifying all previously described satellite repeats, we identified many novel repeats of low to medium abundance. Many of the repeats show population differentiation, including two that are present in only some populations. Interestingly, the population structure inferred from overall satellite quantities does not recapitulate the expected population relationships based on the demographic history of D. melanogaster. We also find that some satellites of similar sequence composition are correlated across lines, revealing concerted evolution. Moreover, correlated satellites tend to be interspersed with each other, further suggesting that concerted change is partially driven by higher order structure. Surprisingly, we identified negative correlations among some satellites, suggesting antagonistic interactions. Our study demonstrates that current genome assemblies vastly underestimate the complexity, abundance, and variation of highly repetitive satellite DNA and presents approaches to understand their rapid evolutionary divergence.