Genome Size Changes by Duplication, Divergence, and Insertion in Caenorhabditis Worms

Genome size has been measurable since the 1940s but we still do not understand genome size variation. Caenorhabditis nematodes show strong conservation of chromosome number but vary in genome size between closely related species. Androdioecy, where populations are composed of males and self-fertile hermaphrodites, evolved from outcrossing, female-male dioecy, three times in this group. In Caenorhabditis, androdioecious genomes are 10-30% smaller than dioecious species, but in the nematode Pristionchus, androdioecy evolved six times and does not correlate with genome size. Previous hypotheses include genome size evolution through: 1) Deletions and "genome shrinkage" in androdioecious species; 2) Transposable element (TE) expansion and DNA loss through large deletions (the "accordion model"); and 3) Differing TE dynamics in androdioecious and dioecious species. We analyzed nematode genomes and found no evidence for these hypotheses. Instead, nematode genome sizes had strong phylogenetic inertia with increases in a few dioecious species, contradicting the "genome shrinkage" hypothesis. TEs did not explain genome size variation with the exception of the DNA transposon Mutator which was twice as abundant in dioecious genomes. Across short and long evolutionary distances Caenorhabditis genomes evolved through small structural mutations including gene-associated duplications and insertions. Seventy-one protein families had significant, parallel decreases across androdioecious Caenorhabditis including genes involved in the sensory system, regulatory proteins and membrane-associated immune responses. Our results suggest that within a dynamic landscape of frequent small rearrangements in Caenorhabditis, reproductive mode mediates genome evolution by altering the precise fates of individual genes, proteins, and the phenotypes they underlie.

Complete Genome Sequence of Rhodococcus opacus Strain MoAcy1 (DSM 44186), an Aerobic Acetylenotroph Isolated from Soil

We report the genome of Rhodococcus opacus strain MoAcy1 (DSM 44186), an aerobic soil isolate capable of using acetylene as its primary carbon and energy source (acetylenotrophy). The genome is composed of a single circular chromosome of ∼8 Mbp and two closed plasmids, with a G+C content of 67.3%.

Slow Recovery from Inbreeding Depression Generated by the Complex Genetic Architecture of Segregating Deleterious Mutations

The deleterious effects of inbreeding have been of extreme importance to evolutionary biology, but it has been difficult to characterize the complex interactions between genetic constraints and selection that lead to fitness loss and recovery after inbreeding. Haploid organisms and selfing organisms like the nematode Caenorhabditis elegans are capable of rapid recovery from the fixation of novel deleterious mutation; however, the potential for recovery and genomic consequences of inbreeding in diploid, outcrossing organisms are not well understood. We sought to answer two questions: 1) Can a diploid, outcrossing population recover from inbreeding via standing genetic variation and new mutation? and 2) How does allelic diversity change during recovery? We inbred C. remanei, an outcrossing relative of C. elegans, through brother-sister mating for 30 generations followed by recovery at large population size. Inbreeding reduced fitness but, surprisingly, recovery from inbreeding at large populations sizes generated only very moderate fitness recovery after 300 generations. We found that 65% of ancestral single nucleotide polymorphisms (SNPs) were fixed in the inbred population, far fewer than the theoretical expectation of ∼99%. Under recovery, 36 SNPs across 30 genes involved in alimentary, muscular, nervous, and reproductive systems changed reproducibly across replicates, indicating that strong selection for fitness recovery does exist. Our results indicate that recovery from inbreeding depression via standing genetic variation and mutation is likely to be constrained by the large number of segregating deleterious variants present in natural populations, limiting the capacity for recovery of small populations.

REVIEW OF THE DAUER HYPOTHESIS: WHAT NON-PARASITIC SPECIES CAN TELL US ABOUT THE EVOLUTION OF PARASITISM

Parasitic lineages have acquired suites of new traits compared to their nearest free-living relatives. When and why did these traits arise? We can envision lineages evolving through multiple stable intermediate steps such as a series of increasingly exploitative species interactions. This view allows us to use non-parasitic species that approximate those intermediate steps to uncover the timing and original function of parasitic traits, knowledge critical to understanding the evolution of parasitism. The dauer hypothesis proposes that free-living nematode lineages evolved into parasites through two intermediate steps, phoresy and necromeny. Here we delve into the proposed steps of the dauer hypothesis by collecting and organizing data from genetic, behavioral, and ecological studies in a range of nematode species. We argue that hypotheses on the evolution of parasites will be strengthened by complementing comparative genomic studies with ecological studies on non-parasites that approximate intermediate steps.

Optimizing experimental design for genome sequencing and assembly with Oxford Nanopore Technologies

High quality reference genome sequences are the core of modern genomics. Oxford Nanopore Technologies (ONT) produces inexpensive DNA sequences, but has high error rates, which make sequence assembly and analysis difficult as genome size and complexity increases. Robust experimental design is necessary for ONT genome sequencing and assembly, but few studies have addressed eukaryotic organisms. Here, we present novel results using simulated and empirical ONT and DNA libraries to identify best practices for sequencing and assembly for several model species. We find that the unique error structure of ONT libraries causes errors to accumulate and assembly statistics plateau as sequence depth increases. High-quality assembled eukaryotic sequences require high-molecular-weight DNA extractions that increase sequence read length, and computational protocols that reduce error through pre-assembly correction and read selection. Our quantitative results will be helpful for researchers seeking guidance for de novo assembly projects.

Complex pleiotropic genetic architecture of evolved heat stress and oxidative stress resistance in the nematode Caenorhabditis remanei

The adaptation of complex organisms to changing environments has been a central question in evolutionary quantitative genetics since its inception. The structure of the genotype-phenotype maps is critical because pleiotropic effects can generate widespread correlated responses to selection and potentially restrict the extent of evolutionary change. In this study, we use experimental evolution to dissect the genetic architecture of natural variation for acute heat stress and oxidative stress response in the nematode Caenorhabiditis remanei. Previous work in the classic model nematode Caenorhabiditis elegans has found that abiotic stress response is controlled by a handful of genes of major effect and that mutations in any one of these genes can have widespread pleiotropic effects on multiple stress response traits. Here, we find that acute heat stress response and acute oxidative response in C. remanei are polygenic, complex traits, with hundreds of genomic regions responding to selection. In contrast to expectation from mutation studies, we find that evolved acute heat stress and acute oxidative stress response for the most part display independent genetic bases. This lack of correlation is reflected at the levels of phenotype, gene expression, and in the genomic response to selection. Thus, while these findings support the general view that rapid adaptation can be generated by changes at hundreds to thousands of sites in the genome, the architecture of segregating variation is likely to be determined by the pleiotropic structure of the underlying genetic networks.

Age-dependent impairment of disease tolerance is associated with a robust transcriptional response following RNA virus infection in Drosophila

Advanced age in humans is associated with greater susceptibility to and higher mortality rates from infections, including infections with some RNA viruses. The underlying innate immune mechanisms, which represent the first line of defense against pathogens, remain incompletely understood. Drosophila melanogaster is able to mount potent and evolutionarily conserved innate immune defenses against a variety of microorganisms including viruses and serves as an excellent model organism for studying host–pathogen interactions. With its relatively short lifespan, Drosophila also is an organism of choice for aging studies. Despite numerous advantages that this model offers, Drosophila has not been used to its full potential to investigate the response of the aged host to viral infection. Here, we show that, in comparison to younger flies, aged Drosophila succumb more rapidly to infection with the RNA-containing Flock House virus due to an age-dependent defect in disease tolerance. Relative to younger individuals, we find that older Drosophila mount transcriptional responses characterized by differential regulation of more genes and genes regulated to a greater extent. We show that loss of disease tolerance to Flock House virus with age associates with a stronger regulation of genes involved in apoptosis, some genes of the Drosophila immune deficiency NF-kB pathway, and genes whose products function in mitochondria and mitochondrial respiration. Our work shows that Drosophila can serve as a model to investigate host–virus interactions during aging and furthermore sets the stage for future analysis of the age-dependent mechanisms that govern survival and control of virus infections at older age.

Syntrophotalea acetylenivorans sp. nov., a diazotrophic, acetylenotrophic anaerobe isolated from intertidal sediments

A Gram-stain-negative, strictly anaerobic, non-motile, rod-shaped bacterium, designated SFB93^T, was isolated from the intertidal sediments of South San Francisco Bay, located near Palo Alto, CA, USA. SFB93^T was capable of acetylenotrophic and diazotrophic growth, grew at 22-37 °C, pH 6.3-8.5 and in the presence of 10-45 g l^-1 NaCl. Phylogenetic analyses based on 16S rRNA gene sequencing showed that SFB93^T represented a member of the genus Syntrophotalea with highest 16S rRNA gene sequence similarities to Syntrophotalea acetylenica DSM 3246^T (96.6 %), Syntrophotalea carbinolica DSM 2380^T (96.5 %), and Syntrophotalea venetiana DSM 2394^T (96.7 %). Genome sequencing revealed a genome size of 3.22 Mbp and a DNA G+C content of 53.4 %. SFB93^T had low genome-wide average nucleotide identity (81-87.5 %) and <70 % digital DNA-DNA hybridization value with other members of the genus Syntrophotalea. The phylogenetic position of SFB93^T within the family Syntrophotaleaceae and as a novel member of the genus Syntrophotalea was confirmed via phylogenetic reconstruction based on concatenated alignments of 92 bacterial core genes. On the basis of the results of phenotypic, genotypic and phylogenetic analyses, a novel species, Syntrophotalea acetylenivorans sp. nov., is proposed, with SFB93^T (=DSM 106009^T=JCM 33327^T=ATCC TSD-118^T) as the type strain.

Caenorhabditis elegans dauers vary recovery in response to bacteria from natural habitat

Many species use dormant stages for habitat selection by tying recovery to informative external cues. Other species have an undiscerning strategy in which they recover randomly despite having advanced sensory systems. We investigated whether elements of a species' habitat structure and life history can bar it from developing a discerning recovery strategy. The nematode Caenorhabditis elegans has a dormant stage called the dauer larva that disperses between habitat patches. On one hand, C. elegans colonization success is profoundly influenced by the bacteria found in its habitat patches, so we might expect this to select for a discerning strategy. On the other hand, C. elegans' habitat structure and life history suggest that there is no fitness benefit to varying recovery, which might select for an undiscerning strategy. We exposed dauers of three genotypes to a range of bacteria acquired from the worms' natural habitat. We found that C. elegans dauers recover in all conditions but increase recovery on certain bacteria depending on the worm's genotype, suggesting a combination of undiscerning and discerning strategies. Additionally, the worms' responses did not match the bacteria's objective quality, suggesting that their decision is based on other characteristics.

Environmental stress maintains trioecy in nematode worms

Sex is determined by chromosomes in mammals but it can be influenced by the environment in many worms, crustaceans, and vertebrates. Despite this, there is little understanding of the relationship between ecology and the evolution of sexual systems. The nematode Auanema freiburgensis has a unique sex determination system in which individuals carrying one X chromosome develop into males while XX individuals develop into females in stress-free environments and self-fertile hermaphrodites in stressful environments. Theory predicts that trioecious populations with coexisting males, females, and hermaphrodites should be unstable intermediates in evolutionary transitions between mating systems. In this article, we study a mathematical model of reproductive evolution based on the unique life history and sex determination of A. freiburgensis. We develop the model in two scenarios, one where the relative production of hermaphrodites and females is entirely dependent on the environment and one based on empirical measurements of a population that displays incomplete, "leaky" environmental dependence. In the first scenario environmental conditions can push the population along an evolutionary continuum and result in the stable maintenance of multiple reproductive systems. The second "leaky" scenario results in the maintenance of three sexes for all environmental conditions. Theoretical investigations of reproductive system transitions have focused on the evolutionary costs and benefits of sex. Here, we show that the flexible sex determination system of A. freiburgensis may contribute to population-level resilience in the microscopic nematode's patchy, ephemeral natural habitat. Our results demonstrate that life history, ecology, and environment may play defining roles in the evolution of sexual systems.

Acetylenotrophy: a hidden but ubiquitous microbial metabolism?

Acetylene (IUPAC name: ethyne) is a colorless, gaseous hydrocarbon, composed of two triple bonded carbon atoms attached to hydrogens (C2H2). When microbiologists and biogeochemists think of acetylene, they immediately think of its use as an inhibitory compound of certain microbial processes and a tracer for nitrogen fixation. However, what is less widely known is that anaerobic and aerobic microorganisms can degrade acetylene, using it as a sole carbon and energy source and providing the basis of a microbial food web. Here, we review what is known about acetylene degrading organisms and introduce the term 'acetylenotrophs' to refer to the microorganisms that carry out this metabolic pathway. In addition, we review the known environmental sources of acetylene and postulate the presence of an hidden acetylene cycle. The abundance of bacteria capable of using acetylene and other alkynes as an energy and carbon source suggests that there are energy cycles present in the environment that are driven by acetylene and alkyne production and consumption that are isolated from atmospheric exchange. Acetylenotrophs may have developed to leverage the relatively high concentrations of acetylene in the pre-Cambrian atmosphere, evolving later to survive in specialized niches where acetylene and other alkynes were produced.

Decontaminating eukaryotic genome assemblies with machine learning

BACKGROUND

High-throughput sequencing has made it theoretically possible to obtain high-quality de novo assembled genome sequences but in practice DNA extracts are often contaminated with sequences from other organisms. Currently, there are few existing methods for rigorously decontaminating eukaryotic assemblies. Those that do exist filter sequences based on nucleotide similarity to contaminants and risk eliminating sequences from the target organism.

RESULTS

We introduce a novel application of an established machine learning method, a decision tree, that can rigorously classify sequences. The major strength of the decision tree is that it can take any measured feature as input and does not require a priori identification of significant descriptors. We use the decision tree to classify de novo assembled sequences and compare the method to published protocols.

CONCLUSIONS

A decision tree performs better than existing methods when classifying sequences in eukaryotic de novo assemblies. It is efficient, readily implemented, and accurately identifies target and contaminant sequences. Importantly, a decision tree can be used to classify sequences according to measured descriptors and has potentially many uses in distilling biological datasets.

Using linkage maps to correct and scaffold de novo genome assemblies: methods, challenges, and computational tools

Modern high-throughput DNA sequencing has made it possible to inexpensively produce genome sequences, but in practice many of these draft genomes are fragmented and incomplete. Genetic linkage maps based on recombination rates between physical markers have been used in biology for over 100 years and a linkage map, when paired with a de novo sequencing project, can resolve mis-assemblies and anchor chromosome-scale sequences. Here, I summarize the methodology behind integrating de novo assemblies and genetic linkage maps, outline the current challenges, review the available software tools, and discuss new mapping technologies.

Reproductive Mode and the Evolution of Genome Size and Structure in Caenorhabditis Nematodes

The self-fertile nematode worms Caenorhabditis elegans, C. briggsae, and C. tropicalis evolved independently from outcrossing male-female ancestors and have genomes 20-40% smaller than closely related outcrossing relatives. This pattern of smaller genomes for selfing species and larger genomes for closely related outcrossing species is also seen in plants. We use comparative genomics, including the first high quality genome assembly for an outcrossing member of the genus (C. remanei) to test several hypotheses for the evolution of genome reduction under a change in mating system. Unlike plants, it does not appear that reductions in the number of repetitive elements, such as transposable elements, are an important contributor to the change in genome size. Instead, all functional genomic categories are lost in approximately equal proportions. Theory predicts that self-fertilization should equalize the effective population size, as well as the resulting effects of genetic drift, between the X chromosome and autosomes. Contrary to this, we find that the self-fertile C. briggsae and C. elegans have larger intergenic spaces and larger protein-coding genes on the X chromosome when compared to autosomes, while C. remanei actually has smaller introns on the X chromosome than either self-reproducing species. Rather than being driven by mutational biases and/or genetic drift caused by a reduction in effective population size under self reproduction, changes in genome size in this group of nematodes appear to be caused by genome-wide patterns of gene loss, most likely generated by genomic adaptation to self reproduction per se.

Variance in epistasis links gene regulation and evolutionary rate in the yeast genetic interaction network

Organisms face a constantly shifting landscape of environmental conditions and internal physiological states. How gene regulation and cellular functions are maintained across genetic and environmental variation is therefore a fundamental question in biology. Here, we analyze the Saccharomyces cerevisiae genetic interaction network to understand how the yeast cell maintains regulatory capacity across genetic backgrounds and environmental conditions. We used the recently characterized synthetic sick/lethal network in yeast, which measures the fitness effects of knocking out pairs of genes, to analyze interactions among 4,364 genes. Genes with large variance in epistatic effects on fitness are highly and ubiquitously expressed (with open chromatin conformations in their promoter regions) and evolve more slowly than genes with weak effects on fitness. Thus, rather than being the elements responsible for the regulation and responsiveness of the genetic network, genes with large epistatic effects tend to be more mundane “housekeeping” genes whose consistent expression is critical to fitness under all environments and that are thereby deeply embedded within the regulatory structure of the network. Our analysis shows that the yeast cell has evolved a system whereby a physical mechanism of regulation (nucleosome occupancy) buffers key genes from the variability experienced by the cell as a whole.

A history of phenotypic plasticity accelerates adaptation to a new environment

Can a history of phenotypic plasticity increase the rate of adaptation to a new environment? Theory suggests it can be through two different mechanisms. Phenotypically plastic organisms can adapt rapidly to new environments through genetic assimilation, or the fluctuating environments that result in phenotypic plasticity can produce evolvable genetic architectures. In this article, I studied a model of a gene regulatory network that determined a phenotypic character in one population selected for phenotypic plasticity and a second population in a constant environment. A history of phenotypic plasticity increased the rate of adaptation in a new environment, but the amount of this increase was dependent on the strength of selection in the original environment. Phenotypic variance in the original environment predicted the adaptive capacity of the trait within, but not between, plastic and nonplastic populations. These results have implications for invasive species and ecological studies of rapid adaptation.

A history of phenotypic plasticity accelerates adaptation to a new environment

Can a history of phenotypic plasticity increase the rate of adaptation to a new environment? Theory suggests it can be through two different mechanisms. Phenotypically plastic organisms can adapt rapidly to new environments through genetic assimilation, or the fluctuating environments that result in phenotypic plasticity can produce evolvable genetic architectures. In this article, I studied a model of a gene regulatory network that determined a phenotypic character in one population selected for phenotypic plasticity and a second population in a constant environment. A history of phenotypic plasticity increased the rate of adaptation in a new environment, but the amount of this increase was dependent on the strength of selection in the original environment. Phenotypic variance in the original environment predicted the adaptive capacity of the trait within, but not between, plastic and nonplastic populations. These results have implications for invasive species and ecological studies of rapid adaptation.

Spatial distribution and reproductive phenology of sexual and asexual Mastocarpus papillatus (Rhodophyta)

Species of the genus Mastocarpus exhibit two distinct life cycles, a sexual alternation of generations and an obligate, asexual direct life cycle that produces only female upright fronds. In the intertidal red alga, M. papillatus (Kützing) sexual fronds dominate southern populations and asexual fronds dominate northern populations along the northeast Pacific coast, a pattern of spatial separation called geographic parthenogenesis. Along the central coast of California, sexual and asexual variants occur in mixed populations, but it is not known whether they are spatially separated within the intertidal zone at a given site. We investigated reproductive phenologies and analyzed patterns of spatial distributions of sexual and asexual M. papillatus at three sites in this region. Sexual M. papillatus were aggregated lower on the shore at two sites and only reproduced during part of a year, while asexual M. papillatus occurred throughout the intertidal range at all sites and reproduced throughout the year. The distribution patterns of sexual and asexual M. papillatus are consistent with a hypothesis of shoreline topography influencing their dynamics of dispersal and colonization. Spatial and temporal partitioning may contribute to the long-term coexistence of sexual and asexual life histories in this, and other, species of Mastocarpus. The occurrence of geographic parthenogenesis at multiple spatial scales in M. papillatus provides an opportunity to gain insight into the phenomenon.

Genetic architecture and postzygotic reproductive isolation: evolution of Bateson-Dobzhansky-Muller incompatibilities in a polygenic model

The Bateson-Dobzhansky-Muller model predicts that postzygotic isolation evolves due to the accumulation of incompatible epistatic interactions, but few studies have quantified the relationship between genetic architecture and patterns of reproductive divergence. We examined how the direction and magnitude of epistatic interactions in a polygenic trait under stabilizing selection influenced the evolution of hybrid incompatibilities. We found that populations evolving independently under stabilizing selection experienced suites of compensatory allelic changes that resulted in genetic divergence between populations despite the maintenance of a stable, high-fitness phenotype. A small number of loci were then incompatible with multiple alleles in the genetic background of the hybrid and the identity of these incompatibility loci changed over the evolution of the populations. For F(1) hybrids, reduced fitness evolved in a window of intermediate strengths of epistatic interactions, but F(2) and backcross hybrids evolved reduced fitness across weak and moderate strengths of epistasis due to segregation variance. Strong epistatic interactions constrained the allelic divergence of parental populations and prevented the development of reproductive isolation. Because many traits with varying genetic architectures must be under stabilizing selection, our results indicate that polygenetic drift is a plausible hypothesis for the evolution of postzygotic reproductive isolation.

Sexual selection accelerates the elimination of a deleterious mutant in Drosophila melanogaster

Although theory indicates that indirect genetic benefits through mate choice should be widespread, empirical work has often either failed to detect the operation of such benefits or shown a net cost to the presence of sexual selection. We tested whether sexual selection can increase the speed with which a conditionally deleterious allele is removed from a laboratory population of Drosophila melanogaster. The alcohol dehydrogenase null allele (Adh-) confers slightly lower viability than wild-type alleles in the absence of ethanol but is lethal in homozygotes when ethanol comprises 6% of the medium. We tracked the frequency of this allele in artificially constructed populations reared at three different levels of ethanol (0%, 2%, and 4%) that either experienced sexual selection or did not. Loss of the deleterious Adh- allele was more rapid when sexual selection was allowed to act, especially in the presence of ethanol. We also quantified the strength of both nonsexual and sexual selection against the Adh- allele using maximum-likelihood estimation. In contrast to recent experiments employing monogamy/polygamy designs, our results demonstrate a fitness benefit to sexual selection. This is consistent with the operation of good-genes female choice.