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.

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.

Acetylenotrophic and Diazotrophic Bradyrhizobium sp. Strain I71 from TCE-Contaminated Soils.

Acetylene (C2H2) is a molecule rarely found in nature, with very few known natural sources, but acetylenotrophic microorganisms can use acetylene as their primary carbon and energy source. As of 2018 there were 15 known strains of aerobic and anaerobic acetylenotrophs; however, we hypothesize there may yet be unrecognized diversity of acetylenotrophs in nature. This study expands the known diversity of acetylenotrophs by isolating the aerobic acetylenotroph, Bradyrhizobium sp. strain I71, from trichloroethylene (TCE)-contaminated soils. Strain I71 is a member of the class Alphaproteobacteria and exhibits acetylenotrophic and diazotrophic activities, the only two enzymatic reactions known to transform acetylene. This unique capability in the isolated strain may increase the genus' economic impact beyond agriculture as acetylenotrophy is closely linked to bioremediation of chlorinated contaminants. Computational analyses indicate that the Bradyrhizobium sp. strain I71 genome contains 522 unique genes compared to close relatives. Moreover, applying a novel hidden Markov model of known acetylene hydratase (AH) enzymes identified a putative AH enzyme. Protein annotation with I-TASSER software predicted the AH from the microbe Syntrophotalea acetylenica as the closest structural and functional analog. Furthermore, the putative AH was flanked by horizontal gene transfer (HGT) elements, like that of AH in anaerobic acetylenotrophs, suggesting an unknown source of acetylene or acetylenic substrate in the environment that is selecting for the presence of AH. IMPORTANCE The isolation of Bradyrhizobium strain I71 expands the distribution of acetylene-consuming microbes to include a group of economically important microorganisms. Members of Bradyrhizobium are well studied for their abilities to improve plant health and increase crop yields by providing bioavailable nitrogen. Additionally, acetylene-consuming microbes have been shown to work in tandem with other microbes to degrade soil contaminants. Based on genome, cultivation, and protein prediction analysis, the ability to consume acetylene is likely not widespread within the genus Bradyrhizobium. These findings suggest that the suite of phenotypic capabilities of strain I71 may be unique and make it a good candidate for further study in several research avenues.

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

A Gram-stain-negative, strictly anaerobic, non-motile, rod-shaped bacterium, designated SFB93T, was isolated from the intertidal sediments of South San Francisco Bay, located near Palo Alto, CA, USA. SFB93T 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 SFB93T represented a member of the genus Syntrophotalea with highest 16S rRNA gene sequence similarities to Syntrophotalea acetylenica DSM 3246T (96.6 %), Syntrophotalea carbinolica DSM 2380T (96.5 %), and Syntrophotalea venetiana DSM 2394T (96.7 %). Genome sequencing revealed a genome size of 3.22 Mbp and a DNA G+C content of 53.4 %. SFB93T 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 SFB93T 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 SFB93T (=DSM 106009T=JCM 33327T=ATCC TSD-118T) as the type strain.

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.

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.

Detection of Diazotrophy in the Acetylene-Fermenting Anaerobe Pelobacter sp. Strain SFB93.

Acetylene (C2H2) is a trace constituent of the present Earth's oxidizing atmosphere, reflecting a mixture of terrestrial and marine emissions from anthropogenic, biomass-burning, and unidentified biogenic sources. Fermentation of acetylene was serendipitously discovered during C2H2 block assays of N2O reductase, and Pelobacter acetylenicus was shown to grow on C2H2 via acetylene hydratase (AH). AH is a W-containing, catabolic, low-redox-potential enzyme that, unlike nitrogenase (N2ase), is specific for acetylene. Acetylene fermentation is a rare metabolic process that is well characterized only in P. acetylenicus DSM3246 and DSM3247 and Pelobacter sp. strain SFB93. To better understand the genetic controls for AH activity, we sequenced the genomes of the three acetylene-fermenting Pelobacter strains. Genome assembly and annotation produced three novel genomes containing gene sequences for AH, with two copies being present in SFB93. In addition, gene sequences for all five compulsory genes for iron-molybdenum N2ase were also present in the three genomes, indicating the cooccurrence of two acetylene transformation pathways. Nitrogen fixation growth assays showed that DSM3426 could ferment acetylene in the absence of ammonium, but no ethylene was produced. However, SFB93 degraded acetylene and, in the absence of ammonium, produced ethylene, indicating an active N2ase. Diazotrophic growth was observed under N2 but not in experimental controls incubated under argon. SFB93 exhibits acetylene fermentation and nitrogen fixation, the only known biochemical mechanisms for acetylene transformation. Our results indicate complex interactions between N2ase and AH and suggest novel evolutionary pathways for these relic enzymes from early Earth to modern days.IMPORTANCE Here we show that a single Pelobacter strain can grow via acetylene fermentation and carry out nitrogen fixation, using the only two enzymes known to transform acetylene. These findings provide new insights into acetylene transformations and adaptations for nutrient (C and N) and energy acquisition by microorganisms. Enhanced understanding of acetylene transformations (i.e., extent, occurrence, and rates) in modern environments is important for the use of acetylene as a potential biomarker for extraterrestrial life and for degradation of anthropogenic contaminants.

Modeling the evolution of complex genetic systems: the gene network family tree.

In 1994 and 1996, Andreas Wagner introduced a novel model in two papers addressing the evolution of genetic regulatory networks. This work, and a suite of papers that followed using similar models, helped integrate network thinking into biology and motivate research focused on the evolution of genetic networks. The Wagner network has its mathematical roots in the Ising model, a statistical physics model describing the activity of atoms on a lattice, and in neural networks. These models have given rise to two branches of applications, one in physics and biology and one in artificial intelligence and machine learning. Here, we review development along these branches, outline similarities and differences between biological models of genetic regulatory circuits and neural circuits models used in machine learning, and identify ways in which these models can provide novel insights into biological systems.

Female mating preferences determine system-level evolution in a gene network model.

Environmental patterns of directional, stabilizing and fluctuating selection can influence the evolution of system-level properties like evolvability and mutational robustness. Intersexual selection produces strong phenotypic selection and these dynamics may also affect the response to mutation and the potential for future adaptation. In order to to assess the influence of mating preferences on these evolutionary properties, I modeled a male trait and female preference determined by separate gene regulatory networks. I studied three sexual selection scenarios: sexual conflict, a Gaussian model of the Fisher process described in Lande (in Proc Natl Acad Sci 78(6):3721-3725, 1981) and a good genes model in which the male trait signalled his mutational condition. I measured the effects these mating preferences had on the potential for traits and preferences to evolve towards new states, and mutational robustness of both the phenotype and the individual's overall viability. All types of sexual selection increased male phenotypic robustness relative to a randomly mating population. The Fisher model also reduced male evolvability and mutational robustness for viability. Under good genes sexual selection, males evolved an increased mutational robustness for viability. Females choosing their mates is a scenario that is sufficient to create selective forces that impact genetic evolution and shape the evolutionary response to mutation and environmental selection. These dynamics will inevitably develop in any population where sexual selection is operating, and affect the potential for future adaptation.

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

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.

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.

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.

Genome Evolution: On the Road to Parasitism.

Roughly 1 billion people are infected with nematode parasites, but there is little understanding of the genomic changes that accompany the evolution of parasitism. A new study analyzes the genome of Caenorhabditis bovis, a nematode that may be evolving a parasitic lifestyle.

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.

De Novo Genome Assemblies for Three North American Bumble Bee Species: Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii.

Bumble bees are ecologically and economically important insect pollinators. Three abundant and widespread species in western North America, Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii, have been the focus of substantial research relating to diverse aspects of bumble bee ecology and evolutionary biology. We present de novo genome assemblies for each of the three species using hybrid assembly of Illumina and Oxford Nanopore Technologies sequences. All three assemblies are of high quality with large N50s (> 2.2 Mb), BUSCO scores indicating > 98% complete genes, and annotations producing 13,325 - 13,687 genes, comparing favorably with other bee genomes. Analysis of synteny against the most complete bumble bee genome, Bombus terrestris, reveals a high degree of collinearity. These genomes should provide a valuable resource for addressing questions relating to functional genomics and evolutionary biology in these species.