Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times. The transition from outcrossing to selfing significantly decreases the effective population size, effective recombination rate, and heterozygosity within a species. These changes lead to a reduction in genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection. Within the nematode genus Caenorhabditis , selfing has evolved at least three times and all three species, including in the model organism Caenorhabditis elegans , show substantially reduced genetic diversity relative to outcrossing species. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to these environments. Here, we examine the whole-genome sequences from 609 wild C. elegans strains isolated worldwide and show that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response, and xenobiotic stress response. Population genomic evidence suggests that genetic diversity in these regions has been maintained by long-term balancing selection. Using long- read genome assemblies for 15 wild strains, we show that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to levels found between Caenorhabditis species that diverged millions of years ago. These results provide an example for how species can avoid the evolutionary “dead end” associated with selfing.

Flaviviruses such as Yellow fever, Dengue, West Nile, and Zika generate disease-linked viral noncoding RNAs called subgenomic flavivirus RNAs. Subgenomic flavivirus RNAs result when the 5′–3′ progression of cellular exoribonuclease Xrn1 is blocked by RNA elements called Xrn1-resistant RNAs located within the viral genome’s 3′-untranslated region that operate without protein co-factors. Here, we show that Xrn1-resistant RNAs can halt diverse exoribonucleases, revealing a mechanism in which they act as general mechanical blocks that ‘brace’ against an enzyme’s surface, presenting an unfolding problem that confounds further enzyme progression. Further, we directly demonstrate that Xrn1-resistant RNAs exist in a diverse set of flaviviruses, including some specific to insects or with no known arthropod vector. These Xrn1-resistant RNAs comprise two secondary structural classes that mirror previously reported phylogenic analysis. Our discoveries have implications for the evolution of exoribonuclease resistance, the use of Xrn1-resistant RNAs in synthetic biology, and the development of new therapies.

Background: The nematode Caenorhabditis elegans has been extensively used to explore the relationships between complex traits, genotypes, and environments. Complex traits can vary across different genotypes of a species, and the genetic regulators of trait variation can be mapped on the genome using quantitative trait locus (QTL) analysis of recombinant inbred lines (RILs) derived from genetically and phenotypically divergent parents. Most RILs have been derived from crossing two parents from globally distant locations. However, the genetic diversity between local C. elegans populations can be as diverse as between global populations and could thus provide means of identifying genetic variation associated with complex traits relevant on a broader scale.Results: To investigate the effect of local genetic variation on heritable traits, we developed a new RIL population derived from 4 parental wild isolates collected from 2 closely located sites in France: Orsay and Santeuil. We crossed these 4 genetically diverse parental isolates to generate a population of 200 multi-parental RILs and used RNA-seq to obtain sequence polymorphisms identifying almost 9000 SNPs variable between the 4 genotypes with an average spacing of 11 kb, doubling the mapping resolution relative to currently available RIL panels for many loci. The SNPs were used to construct a genetic map to facilitate QTL analysis. We measured life history traits such as lifespan, stress resistance, developmental speed, and population growth in different environments, and found substantial variation for most traits. We detected multiple QTLs for most traits, including novel QTLs not found in previous QTL analysis, including those for lifespan and pathogen responses. This shows that recombining genetic variation across C. elegans populations that are in geographical close proximity provides ample variation for QTL mapping. Conclusion: Taken together, we show that using more parents than the classical two parental genotypes to construct a RIL population facilitates the detection of QTLs and that the use of wild isolates facilitates the detection of QTLs. The use of multi-parent RIL populations can further enhance our understanding of local adaptation and life history trade-offs.

There is considerable insight into pathways and genes associated with heat-stress conditions. Most genes involved in stress response have been identified using mutant screens or gene knockdowns. Yet, there is limited understanding of the temporal dynamics of global gene expression in stressful environments. Here, we studied global gene expression profiles during 12 hours of heat stress in the nematode C. elegans. Using a high-resolution time series of increasing stress exposures, we found a distinct shift in gene expression patterns between 3–4 hours into the stress response, separating an initially highly dynamic phase from a later relatively stagnant phase. This turning point in expression dynamics coincided with a phenotypic turning point, as shown by a strong decrease in movement, survival and, progeny count in the days following the stress. Both detectable at transcriptional and phenotypic level, this study pin-points a relatively small time frame during heat stress at which enough damage is accumulated, making it impossible to recover the next few days.

No abstract

Summary paragraphThe mating system of a species profoundly influences its evolutionary trajectory1–3. Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times4. The transition from outcrossing to selfing significantly decreases the effective population size, effective recombination rate, and heterozygosity within a species5. These changes lead to a reduction in the genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection6,7. Selfing has evolved at least three times independently in the nematode genus Caenorhabditis8, including in the model organism Caenorhabditis elegans, and all three selfing species show substantially reduced genetic diversity relative to outcrossing species8,9. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to and inhabit these same diverse environments. Here, we discovered previously uncharacterized levels and patterns of genetic diversity by examining the whole-genome sequences from 609 wild C. elegans strains isolated worldwide. We found that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response, and xenobiotic stress. Population genomic evidence suggests that these regions have been maintained by balancing selection. Using long-read genome assemblies for 15 wild isolates, we found that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to that found between Caenorhabditis species that diverged millions of years ago. Taken together, these results suggest that ancient genetic diversity present in the outcrossing ancestor of C. elegans has been maintained by long-term balancing selection since the evolution of selfing. These results provide an example for how species can avoid the evolutionary “dead end” associated with selfing by maintaining ancestral genetic diversity.

Summary Susceptibility to the root‐knot nematode Meloidogyne incognita in plants is thought to be a complex trait based on multiple genes involved in cell differentiation, growth and defence. Previous genetic analyses of susceptibility to M. incognita have mainly focused on segregating dominant resistance genes in crops. It is not known if plants harbour significant genetic variation in susceptibility to M. incognita independent of dominant resistance.To study the genetic architecture of susceptibility to M. incognita, we analysed nematode reproduction on a highly diverse set of 340 natural inbred lines of Arabidopsis thaliana with genome‐wide association mapping. We observed a surprisingly large variation in nematode reproduction among these lines.Genome‐wide association mapping revealed four quantitative trait loci (QTLs) located on chromosomes 1 and 5 of A. thaliana significantly associated with reproductive success of M. incognita, none of which harbours typical resistance gene homologues. Mutant analysis of three genes located in two QTLs showed that the transcription factor BRASSINAZOLE RESISTANT1 and an F‐box family protein may function as (co‐)regulators of susceptibility to M. incognita in Arabidopsis.Our data suggest that breeding for loss‐of‐susceptibility, based on allelic variants critically involved in nematode feeding, could be used to make crops more resilient to root‐knot nematodes.

BackgroundCryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population.ResultsWe measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40-57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL.ConclusionsThese results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 hours) and that CGV is mainly determined by response related trans regulatory eQTL.