Natural variation in infection specificity of Caenorhabditis briggsae isolates by two RNA viruses.

Antagonistic relationships such as host-virus interactions potentially lead to rapid evolution and specificity in interactions. The Orsay virus is so far the only horizontal virus naturally infecting the nematode C. elegans. In contrast, several related RNA viruses infect its congener C. briggsae, including Santeuil (SANTV) and Le Blanc (LEBV) viruses. Here we focus on the host's intraspecific variation in sensitivity to these two intestinal viruses. Many temperate-origin C. briggsae strains, including JU1264 and JU1498, are sensitive to both, while many tropical strains, such as AF16, are resistant to both. Interestingly, some C. briggsae strains exhibit a specific resistance, such as the HK104 strain, specifically resistant to LEBV. The viral sensitivity pattern matches the strains' geographic and genomic relationships. The heavily infected strains mount a seemingly normal small RNA response that is insufficient to suppress viral infection, while the resistant strains show no small RNA response, suggesting an early block in viral entry or replication. We use a genetic approach from the host side to map genomic regions participating in viral resistance polymorphisms. Using Advanced Intercrossed Recombinant Inbred Lines (RILs) between virus-resistant AF16 and SANTV-sensitive HK104, we detect Quantitative Trait Loci (QTLs) on chromosomes IV and III. Building RILs between virus-sensitive JU1498 and LEBV-resistant HK104 followed by bulk segregant analysis, we identify a chromosome II QTL. In both cases, further introgressions of the regions confirmed the QTLs. This diversity provides an avenue for studying virus entry, replication, and exit mechanisms, as well as host-virus specificity and the host response to a specific virus infection.

Genome-wide association and environmental suppression of the mortal germline phenotype of wild C. elegans.

The animal germline lineage needs to be maintained along generations. However, some Caenorhabditis elegans wild isolates display a mortal germline phenotype, leading to sterility after several generations at 25°C. Using a genome-wide association approach, we detect a significant peak on chromosome III around 5 Mb, confirmed by introgressions. Thus, a seemingly deleterious genotype is maintained at intermediate frequency in the species. Environmental rescue is a likely explanation, and indeed associated bacteria and microsporidia suppress the phenotype of wild isolates as well as mutants in small RNA inheritance (nrde-2) and histone modifications (set-2). Escherichia coli strains of the K-12 lineage suppress the phenotype compared to B strains. By shifting a wild strain from E. coli K-12 to E. coli B, we find that memory of the suppressing condition is maintained over several generations. Thus, the mortal germline phenotype of wild C. elegans is in part revealed by laboratory conditions and may represent variation in epigenetic inheritance and environmental interactions. This study also points to the importance of non-genetic memory in the face of environmental variation.

Vertical transmission in Caenorhabditis nematodes of RNA molecules encoding a viral RNA-dependent RNA polymerase.

Here, we report on the discovery in Caenorhabditis nematodes of multiple vertically transmitted RNAs coding for putative RNA-dependent RNA polymerases. Their sequences share similarity to distinct RNA viruses, including bunyaviruses, narnaviruses, and sobemoviruses. The sequences are present exclusively as RNA and are not found in DNA form. The RNAs persist in progeny after bleach treatment of adult animals, indicating vertical transmission of the RNAs. We tested one of the infected strains for transmission to an uninfected strain and found that mating of infected animals with uninfected animals resulted in infected progeny. By in situ hybridization, we detected several of these RNAs in the cytoplasm of the male and female germline of the nematode host. The Caenorhabditis hosts were found defective in degrading exogenous double-stranded RNAs, which may explain retention of viral-like RNAs. Strikingly, one strain, QG551, harbored three distinct virus-like RNA elements. Specific patterns of small RNAs complementary to the different viral-like RNAs were observed, suggesting that the different RNAs are differentially recognized by the RNA interference (RNAi) machinery. While vertical transmission of viruses in the family Narnaviridae, which are known as capsidless viruses, has been described in fungi, these observations provide evidence that multicellular animal cells harbor similar viruses.

Noda-Like RNA Viruses Infecting Caenorhabditis Nematodes: Sympatry, Diversity, and Reassortment.

Three RNA viruses related to nodaviruses were previously described to naturally infect the nematode Caenorhabditis elegans and its relative, Caenorhabditis briggsae Here, we report on a collection of more than 50 viral variants from wild-caught Caenorhabditis. We describe the discovery of a new related virus, the Melník virus, infecting C. briggsae, which similarly infects intestinal cells. In France, a frequent pattern of coinfection of C. briggsae by the Santeuil virus and Le Blanc virus was observed at the level of an individual nematode and even a single cell. We do not find evidence of reassortment between the RNA1 and RNA2 molecules of Santeuil and Le Blanc viruses. However, by studying patterns of evolution of each virus, reassortments of RNA1 and RNA2 among variants of each virus were identified. We develop assays to test the relative infectivity and competitive ability of the viral variants and detect an interaction between host genotype and Santeuil virus genotype, such that the result depends on the host strain.IMPORTANCE The roundworm Caenorhabditis elegans is a laboratory model organism in biology. We study natural populations of this small animal and its relative, C. briggsae, and the viruses that infect them. We previously discovered three RNA viruses related to nodaviruses and here describe a fourth one, called the Melník virus. These viruses have a genome composed of two RNA molecules. We find that two viruses may infect the same animal and the same cell. The two RNA molecules may be exchanged between variants of a given viral species. We study the diversity of each viral species and devise an assay of their infectivity and competitive ability. Using this assay, we show that the outcome of the competition also depends on the host.

Chromosomal rearrangements maintain a polymorphic supergene controlling butterfly mimicry.

Supergenes are tight clusters of loci that facilitate the co-segregation of adaptive variation, providing integrated control of complex adaptive phenotypes. Polymorphic supergenes, in which specific combinations of traits are maintained within a single population, were first described for 'pin' and 'thrum' floral types in Primula and Fagopyrum, but classic examples are also found in insect mimicry and snail morphology. Understanding the evolutionary mechanisms that generate these co-adapted gene sets, as well as the mode of limiting the production of unfit recombinant forms, remains a substantial challenge. Here we show that individual wing-pattern morphs in the polymorphic mimetic butterfly Heliconius numata are associated with different genomic rearrangements at the supergene locus P. These rearrangements tighten the genetic linkage between at least two colour-pattern loci that are known to recombine in closely related species, with complete suppression of recombination being observed in experimental crosses across a 400-kilobase interval containing at least 18 genes. In natural populations, notable patterns of linkage disequilibrium (LD) are observed across the entire P region. The resulting divergent haplotype clades and inversion breakpoints are found in complete association with wing-pattern morphs. Our results indicate that allelic combinations at known wing-patterning loci have become locked together in a polymorphic rearrangement at the P locus, forming a supergene that acts as a simple switch between complex adaptive phenotypes found in sympatry. These findings highlight how genomic rearrangements can have a central role in the coexistence of adaptive phenotypes involving several genes acting in concert, by locally limiting recombination and gene flow.

Genome-wide mosaicism within Mycobacterium abscessus: evolutionary and epidemiological implications.

BACKGROUND: In mycobacteria, conjugation differs from the canonical Hfr model, but is still poorly understood. Here, we quantified this evolutionary processe in a natural mycobacterial population, taking advantage of a large clinical strain collection of the emerging pathogen Mycobacterium abscessus (MAB). RESULTS: Multilocus sequence typing confirmed the existence of three M. abscessus subspecies, and unravelled extensive allelic exchange between them. Furthermore, an asymmetrical gene flow occurring between these main lineages was detected, resulting in highly admixed strains. Intriguingly, these mosaic strains were significantly associated with cystic fibrosis patients with lung infections or chronic colonization. Genome sequencing of those hybrid strains confirmed that half of their genomic content was remodelled in large genomic blocks, leading to original tri-modal 'patchwork' architecture. One of these hybrid strains acquired a locus conferring inducible macrolide resistance, and a large genomic insertion from a slowly growing pathogenic mycobacteria, suggesting an adaptive gene transfer. This atypical genomic architecture of the highly recombinogenic strains is consistent with the distributive conjugal transfer (DCT) observed in M. smegmatis. Intriguingly, no known DCT function was found in M. abscessus chromosome, however, a p-RAW-like genetic element was detected in one of the highly admixed strains. CONCLUSION: Taken together, our results strongly suggest that MAB evolution is sporadically punctuated by dramatic genome wide remodelling events. These findings might have far reaching epidemiological consequences for emerging mycobacterial pathogens survey in the context of increasing numbers of rapidly growing mycobacteria and M. tuberculosis co-infections.

"Epidemic clones" of Listeria monocytogenes are widespread and ancient clonal groups.

The food-borne pathogen Listeria monocytogenes is genetically heterogeneous. Although some clonal groups have been implicated in multiple outbreaks, there is currently no consensus on how "epidemic clones" should be defined. The objectives of this work were to compare the patterns of sequence diversity on two sets of genes that have been widely used to define L. monocytogenes clonal groups: multilocus sequence typing (MLST) and multi-virulence-locus sequence typing (MvLST). Further, we evaluated the diversity within clonal groups by pulsed-field gel electrophoresis (PFGE). Based on 125 isolates of diverse temporal, geographical, and source origins, MLST and MvLST genes (i) had similar patterns of sequence polymorphisms, recombination, and selection, (ii) provided concordant phylogenetic clustering, and (iii) had similar discriminatory power, which was not improved when we combined both data sets. Inclusion of representative strains of previous outbreaks demonstrated the correspondence of epidemic clones with previously recognized MLST clonal complexes. PFGE analysis demonstrated heterogeneity within major clones, most of which were isolated decades before their involvement in outbreaks. We conclude that the "epidemic clone" denominations represent a redundant but largely incomplete nomenclature system for MLST-defined clones, which must be regarded as successful genetic groups that are widely distributed across time and space.

Rapid emergence of extensively drug-resistant Shigella sonnei in France.

Shigella sonnei, the main cause of bacillary dysentery in high-income countries, has become increasingly resistant to antibiotics. We monitored the antimicrobial susceptibility of 7121 S. sonnei isolates collected in France between 2005 and 2021. We detected a dramatic increase in the proportion of isolates simultaneously resistant to ciprofloxacin (CIP), third-generation cephalosporins (3GCs) and azithromycin (AZM) from 2015. Our genomic analysis of 164 such extensively drug-resistant (XDR) isolates identified 13 different clusters within CIP-resistant sublineage 3.6.1, which was selected in South Asia ∼15 years ago. AZM resistance was subsequently acquired, principally through IncFII (pKSR100-like) plasmids. The last step in the development of the XDR phenotype involved various extended-spectrum beta-lactamase genes (blaCTX-M-3, blaCTX-M-15, blaCTX-M-27, blaCTX-M-55, and blaCTX-M-134) carried by different plasmids (IncFII, IncI1, IncB/O/K/Z) or even integrated into the chromosome, and encoding resistance to 3GCs. This rapid emergence of XDR S. sonnei, including an international epidemic strain, is alarming, and good laboratory-based surveillance of shigellosis will be crucial for informed decision-making and appropriate public health action.

Local dispersal of Puccinia triticina and wheat canopy structure.

The development of dynamic models jointly to simulate host growth and disease spread necessitates a precise description of pathogen dispersal in relation to canopy structure. In this study, we measured disease spread from a single infected leaf positioned at different heights in wheat canopies. The resulting lesion distribution was described along crop rows and over three leaf layers. The spore sources, although limited to a single leaf, nearly saturated the host surface accessible to the spores. Most of the lesions were found within 30 to 40 cm of the source. The vertical position of the source influenced the lesion distribution and the steepness of the disease gradients. The leaf layer and the wheat row that contained the spore source were the most infected. Close to the source, a few heavily infected leaves produced steep disease gradients, whereas spore diffusion resulted in shallower gradients along the adjacent rows and on the other leaf layers. Depending on the precision needed, the lesion distribution can be described either at the level of leaf layers or by dispersal gradients for each row and leaf layer.

Selection and gene flow shape niche-associated variation in pheromone response.

From quorum sensing in bacteria to pheromone signalling in social insects, chemical communication mediates interactions among individuals in local populations. In Caenorhabditis elegans, ascaroside pheromones can dictate local population density; high levels of pheromones inhibit the reproductive maturation of individuals. Little is known about how natural genetic diversity affects the pheromone responses of individuals from diverse habitats. Here, we show that a niche-associated variation in pheromone receptor genes contributes to natural differences in pheromone responses. We identified putative loss-of-function deletions that impair duplicated pheromone receptor genes (srg-36 and srg-37), which were previously shown to be lost in population-dense laboratory cultures. A common natural deletion in srg-37 arose recently from a single ancestral population that spread throughout the world; this deletion underlies reduced pheromone sensitivity across the global C. elegans population. We found that many local populations harbour individuals with a wild-type or a deletion allele of srg-37, suggesting that balancing selection has maintained the recent variation in this pheromone receptor gene. The two srg-37 genotypes are associated with niche diversity underlying boom-and-bust population dynamics. We hypothesize that human activities likely contributed to the gene flow and balancing selection of srg-37 variation through facilitating the migration of species and providing a favourable niche for the recently arisen srg-37 deletion.

Comparative genomics of 10 new Caenorhabditis species.

The nematode Caenorhabditis elegans has been central to the understanding of metazoan biology. However, C. elegans is but one species among millions and the significance of this important model organism will only be fully revealed if it is placed in a rich evolutionary context. Global sampling efforts have led to the discovery of over 50 putative species from the genus Caenorhabditis, many of which await formal species description. Here, we present species descriptions for 10 new Caenorhabditis species. We also present draft genome sequences for nine of these new species, along with a transcriptome assembly for one. We exploit these whole-genome data to reconstruct the Caenorhabditis phylogeny and use this phylogenetic tree to dissect the evolution of morphology in the genus. We reveal extensive variation in genome size and investigate the molecular processes that underlie this variation. We show unexpected complexity in the evolutionary history of key developmental pathway genes. These new species and the associated genomic resources will be essential in our attempts to understand the evolutionary origins of the C. elegans model.

Four years of DNA barcoding: current advances and prospects.

Research using cytochrome c oxidase barcoding techniques on zoological specimens was initiated by Hebert et al. [Hebert, P.D.N., Ratnasingham, S., deWaard, J.R., 2003. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. Lond. B 270, S96-S99]. By March 2004, the Consortium for the Barcode of Life started to promote the use of a standardized DNA barcoding approach, consisting of identifying a specimen as belonging to a certain animal species based on a single universal marker: the DNA barcode sequence. Over the last 4 years, this approach has become increasingly popular and advances as well as limitations have clearly emerged as increasing amounts of organisms have been studied. Our purpose is to briefly expose DNA Barcode of Life principles, pros and cons, relevance and universality. The initially proposed Barcode of life framework has greatly evolved, giving rise to a flexible description of DNA barcoding and a larger range of applications.

Adaptation of a Fungal Pathogen to Host Quantitative Resistance.

Impact of host quantitative resistance on pathogen evolution is still poorly documented. In our study, we characterized the adaptation of the pathogenic fungus Colletotrichum gloeosporioides, to the quantitative resistance of its host, the water yam (Dioscorea alata). Genetic and pathogenic diversities of C. gloeosporioides populations were specified at the field scale. We used nuclear markers to describe fungal population structuring within and between six fields of three cultivars differently susceptible to the fungus. Strain aggressiveness was then quantified in the laboratory through cross-inoculation tests. The high level of genetic diversity and significant linkage disequilibrium revealed a significant influence of clonal reproduction in the C. gloeosporioides evolution. The recorded fungal migration between fields was weak (evidence for a dispersion mode via tubers rather than splashing dispersal), which provides the first molecular evidence for limited C. gloeosporioides migration via yam tuber exchanges. C. gloeosporioides's populations are adapted to their host resistance. The aggressiveness of the fungal clones seems to have evolved toward an accumulation of components specific to each host cultivar. Despite the remaining marks of adaptation to the former widely cultivated host, adaptation to current cultivars was clearly depicted.

Natural Genetic Variation in a Multigenerational Phenotype in C. elegans.

Although heredity mostly relies on the transmission of DNA sequence, additional molecular and cellular features are heritable across several generations. In the nematode Caenorhabditis elegans, insights into such unconventional inheritance result from two lines of work. First, the mortal germline (Mrt) phenotype was defined as a multigenerational phenotype whereby a selfing lineage becomes sterile after several generations, implying multigenerational memory [1, 2]. Second, certain RNAi effects are heritable over several generations in the absence of the initial trigger [3-5]. Both lines of work converged when the subset of Mrt mutants that are heat sensitive were found to closely correspond to mutants defective in the RNAi-inheritance machinery, including histone modifiers [6-9]. Here, we report the surprising finding that several C. elegans wild isolates display a heat-sensitive mortal germline phenotype in laboratory conditions: upon chronic exposure to higher temperatures, such as 25°C, lines reproducibly become sterile after several generations. This phenomenon is reversible, as it can be suppressed by temperature alternations at each generation, suggesting a non-genetic basis for the sterility. We tested whether natural variation in the temperature-induced Mrt phenotype was of genetic nature by building recombinant inbred lines between the isolates MY10 (Mrt) and JU1395 (non-Mrt). Using bulk segregant analysis, we detected two quantitative trait loci. After further recombinant mapping and genome editing, we identified the major causal locus as a polymorphism in the set-24 gene, encoding a SET- and SPK-domain protein. We conclude that C. elegans natural populations may harbor natural genetic variation in epigenetic inheritance phenomena.

Supergene Evolution Triggered by the Introgression of a Chromosomal Inversion.

Supergenes are groups of tightly linked loci whose variation is inherited as a single Mendelian locus and are a common genetic architecture for complex traits under balancing selection [1-8]. Supergene alleles are long-range haplotypes with numerous mutations underlying distinct adaptive strategies, often maintained in linkage disequilibrium through the suppression of recombination by chromosomal rearrangements [1, 5, 7-9]. However, the mechanism governing the formation of supergenes is not well understood and poses the paradox of establishing divergent functional haplotypes in the face of recombination. Here, we show that the formation of the supergene alleles encoding mimicry polymorphism in the butterfly Heliconius numata is associated with the introgression of a divergent, inverted chromosomal segment. Haplotype divergence and linkage disequilibrium indicate that supergene alleles, each allowing precise wing-pattern resemblance to distinct butterfly models, originate from over a million years of independent chromosomal evolution in separate lineages. These "superalleles" have evolved from a chromosomal inversion captured by introgression and maintained in balanced polymorphism, triggering supergene inheritance. This mode of evolution involving the introgression of a chromosomal rearrangement is likely to be a common feature of complex structural polymorphisms associated with the coexistence of distinct adaptive syndromes. This shows that the reticulation of genealogies may have a powerful influence on the evolution of genetic architectures in nature.

Pan-arthropod analysis reveals somatic piRNAs as an ancestral defence against transposable elements.

In animals, small RNA molecules termed PIWI-interacting RNAs (piRNAs) silence transposable elements (TEs), protecting the germline from genomic instability and mutation. piRNAs have been detected in the soma in a few animals, but these are believed to be specific adaptations of individual species. Here, we report that somatic piRNAs were probably present in the ancestral arthropod more than 500 million years ago. Analysis of 20 species across the arthropod phylum suggests that somatic piRNAs targeting TEs and messenger RNAs are common among arthropods. The presence of an RNA-dependent RNA polymerase in chelicerates (horseshoe crabs, spiders and scorpions) suggests that arthropods originally used a plant-like RNA interference mechanism to silence TEs. Our results call into question the view that the ancestral role of the piRNA pathway was to protect the germline and demonstrate that small RNA silencing pathways have been repurposed for both somatic and germline functions throughout arthropod evolution.

C. elegans outside the Petri dish.

The roundworm Caenorhabditis elegans has risen to the status of a top model organism for biological research in the last fifty years. Among laboratory animals, this tiny nematode is one of the simplest and easiest organisms to handle. And its life outside the laboratory is beginning to be unveiled. Like other model organisms, C. elegans has a boom-and-bust lifestyle. It feasts on ephemeral bacterial blooms in decomposing fruits and stems. After resource depletion, its young larvae enter a migratory diapause stage, called the dauer. Organisms known to be associated with C. elegans include migration vectors (such as snails, slugs and isopods) and pathogens (such as microsporidia, fungi, bacteria and viruses). By deepening our understanding of the natural history of C. elegans, we establish a broader context and improved tools for studying its biology.

Orsay, Santeuil and Le Blanc viruses primarily infect intestinal cells in Caenorhabditis nematodes.

The discoveries of Orsay, Santeuil and Le Blanc viruses, three viruses infecting either Caenorhabditis elegans or its relative Caenorhabditis briggsae, enable the study of virus-host interactions using natural pathogens of these two well-established model organisms. We characterized the tissue tropism of infection in Caenorhabditis nematodes by these viruses. Using immunofluorescence assays targeting proteins from each of the viruses, and in situ hybridization, we demonstrate viral proteins and RNAs localize to intestinal cells in larval stage Caenorhabditis nematodes. Viral proteins were detected in one to six of the 20 intestinal cells present in Caenorhabditis nematodes. In Orsay virus-infected C. elegans, viral proteins were detected as early as 6h post-infection. The RNA-dependent RNA polymerase and capsid proteins of Orsay virus exhibited different subcellular localization patterns. Collectively, these observations provide the first experimental insights into viral protein expression in any nematode host, and broaden our understanding of viral infection in Caenorhabditis nematodes.

A deletion polymorphism in the Caenorhabditis elegans RIG-I homolog disables viral RNA dicing and antiviral immunity.

RNA interference defends against viral infection in plant and animal cells. The nematode Caenorhabditis elegans and its natural pathogen, the positive-strand RNA virus Orsay, have recently emerged as a new animal model of host-virus interaction. Using a genome-wide association study in C. elegans wild populations and quantitative trait locus mapping, we identify a 159 base-pair deletion in the conserved drh-1 gene (encoding a RIG-I-like helicase) as a major determinant of viral sensitivity. We show that DRH-1 is required for the initiation of an antiviral RNAi pathway and the generation of virus-derived siRNAs (viRNAs). In mammals, RIG-I-domain containing proteins trigger an interferon-based innate immunity pathway in response to RNA virus infection. Our work in C. elegans demonstrates that the RIG-I domain has an ancient role in viral recognition. We propose that RIG-I acts as modular viral recognition factor that couples viral recognition to different effector pathways including RNAi and interferon responses. DOI:http://dx.doi.org/10.7554/eLife.00994.001.

Autoinfection in wheat leaf rust epidemics.

Autoinfection (within-host inoculum transmission) allows plant pathogens locally to increase their density on an infected host. Estimating autoinfection is of particular importance in understanding epidemic development in host mixtures. More generally, autoinfection influences the rate of host colonization by the pathogen, as well as pathogen evolution. Despite its importance in epidemiological models, autoinfection has not yet been directly quantified. It was measured here on wheat (Triticum aestivum) leaves infected by a pathogenic fungus (Puccinia triticina). Autoinfection was measured either on inoculated leaves or by assessing the local progeny of spontaneous infections, and was described by a model of the form y = microx(alpha), where alpha accounts for host saturation and micro represents the pathogen multiplication rate resulting from autoinfection. It was shown that autoinfection resulted in typical patterns of disease aggregation at the leaf level and influenced lesion distribution in the crop during the first epidemic stages. The parameter micro was calculated by taking overdispersion of the data and density dependence into account. It was found that a single lesion produced between 50 and 200 offspring by autoinfection, within a pathogen generation. By taking into account environmental variability, it was possible to estimate autoinfection under optimal conditions for epidemic development.