A conserved family of immune effectors cleaves cellular ATP upon viral infection.

During viral infection, cells can deploy immune strategies that deprive viruses of molecules essential for their replication. Here, we report a family of immune effectors in bacteria that, upon phage infection, degrade cellular adenosine triphosphate (ATP) and deoxyadenosine triphosphate (dATP) by cleaving the N-glycosidic bond between the adenine and sugar moieties. These ATP nucleosidase effectors are widely distributed within multiple bacterial defense systems, including cyclic oligonucleotide-based antiviral signaling systems (CBASS), prokaryotic argonautes, and nucleotide-binding leucine-rich repeat (NLR)-like proteins, and we show that ATP and dATP degradation during infection halts phage propagation. By analyzing homologs of the immune ATP nucleosidase domain, we discover and characterize Detocs, a family of bacterial defense systems with a two-component phosphotransfer-signaling architecture. The immune ATP nucleosidase domain is also encoded within diverse eukaryotic proteins with immune-like architectures, and we show biochemically that eukaryotic homologs preserve the ATP nucleosidase activity. Our findings suggest that ATP and dATP degradation is a cell-autonomous innate immune strategy conserved across the tree of life.

Innate immunity: the bacterial connection.

Pathogens have fueled the diversification of intracellular defense strategies that collectively define cell-autonomous innate immunity. In bacteria, innate immunity is manifested by a broad arsenal of defense systems that provide protection against bacterial viruses, called phages. The complexity of the bacterial immune repertoire has only been realized recently and is now suggesting that innate immunity has commonalities across the tree of life: many components of eukaryotic innate immunity are found in bacteria where they protect against phages, including the cGAS-STING pathway, gasdermins, and viperins. Here, I summarize recent findings on the conservation of innate immune pathways between prokaryotes and eukaryotes and hypothesize that bacterial defense mechanisms can catalyze the discovery of novel molecular players of eukaryotic innate immunity.

GSpace: an exact coalescence simulator of recombining genomes under isolation by distance.

MOTIVATION: Simulation-based inference can bypass the limitations of statistical methods based on analytical approximations, but software allowing simulation of structured population genetic data without the classical n-coalescent approximations (such as those following from assuming large population size) are scarce or slow. RESULTS: We present GSpace, a simulator for genomic data, based on a generation-by-generation coalescence algorithm taking into account small population size, recombination and isolation by distance. AVAILABILITY AND IMPLEMENTATION: Freely available at site web INRAe (http://www1.montpellier.inra.fr/CBGP/software/gspace/download.html).

Sex-specific spatial variation in fitness in the highly dimorphic Leucadendron rubrum.

Sexual dimorphism in plants may emerge as a result of sex-specific selection on traits enhancing access to nutritive resources and/or to sexual partners. Here we investigated sex-specific differences in selection of sexually dimorphic traits and in the spatial distribution of effective fecundity (our fitness proxy) in a highly dimorphic dioecious wind-pollinated shrub, Leucadendron rubrum. In particular, we tested for the effect of density on male and female effective fecundity. We used spatial and genotypic data of parent and offspring cohorts to jointly estimate individual male and female effective fecundity on the one hand and pollen and seed dispersal kernels on the other hand. This methodology was adapted to the case of dioecious species. Explicitly modelling dispersal avoids the confounding effects of heterogeneous spatial distribution of mates and sampled seedlings on the estimation of effective fecundity. We also estimated selection gradients on plant traits while modelling sex-specific spatial autocorrelation in fecundity. Males exhibited spatial autocorrelation in effective fecundity at a smaller scale than females. A higher local density of plants was associated with lower effective fecundity in males but was not related to female effective fecundity. These results suggest sex-specific sensitivities to environmental heterogeneity in L. rubrum. Despite these sexual differences, we found directional selection for wider canopies and smaller leaves in both sexes, and no sexually antagonistic selection on strongly dimorphic traits in L. rubrum. Many empirical studies in animals similarly failed to detect sexually antagonistic selection in species expressing strong sexual dimorphism, and we discuss reasons explaining this common pattern.

Genome-wide CRISPR-dCas9 screens in E. coli identify essential genes and phage host factors.

High-throughput genetic screens are powerful methods to identify genes linked to a given phenotype. The catalytic null mutant of the Cas9 RNA-guided nuclease (dCas9) can be conveniently used to silence genes of interest in a method also known as CRISPRi. Here, we report a genome-wide CRISPR-dCas9 screen using a starting pool of ~ 92,000 sgRNAs which target random positions in the chromosome of E. coli. To benchmark our method, we first investigate its utility to predict gene essentiality in the genome of E. coli during growth in rich medium. We could identify 79% of the genes previously reported as essential and demonstrate the non-essentiality of some genes annotated as essential. In addition, we took advantage of the intermediate repression levels obtained when targeting the template strand of genes to show that cells are very sensitive to the expression level of a limited set of essential genes. Our data can be visualized on CRISPRbrowser, a custom web interface available at crispr.pasteur.fr. We then apply the screen to discover E. coli genes required by phages λ, T4 and 186 to kill their host, highlighting the involvement of diverse host pathways in the infection process of the three tested phages. We also identify colanic acid capsule synthesis as a shared resistance mechanism to all three phages. Finally, using a plasmid packaging system and a transduction assay, we identify genes required for the formation of functional λ capsids, thus covering the entire phage cycle. This study demonstrates the usefulness and convenience of pooled genome-wide CRISPR-dCas9 screens in bacteria and paves the way for their broader use as a powerful tool in bacterial genomics.

When Do Individuals Maximize Their Inclusive Fitness?

Adaptation is often described in behavioral ecology as individuals maximizing their inclusive fitness. Under what conditions does this hold, and how does this relate to the gene-centered perspective of adaptation? We unify and extend the literature on these questions to class-structured populations. We demonstrate that the maximization (in the best-response sense) of class-specific inclusive fitness obtains in uninvadable population states (meaning that all deviating mutants become extinct). This defines a genuine actor-centered perspective on adaptation. But this inclusive fitness is assigned to all bearers of a mutant allele in a given class and depends on distributions of demographic and genetic contexts. These distributions, in turn, usually depend on events in previous generations and are thus not under individual control. This prevents, in general, envisioning individuals themselves as autonomous fitness maximizers, each with its own inclusive fitness. For weak selection, however, the dependence on earlier events can be neglected. We then show that each individual in each class appears to maximize its own inclusive fitness when all other individuals exhibit inclusive fitness-maximizing behavior. This defines a genuine individual-centered perspective of adaptation and justifies formally, as a first-order approximation, the long-heralded view of individuals appearing to maximize their own inclusive fitness.

Farming plant cooperation in crops.

Selection of the fittest can promote individual competitiveness but often results in the erosion of group performance. Recently, several authors revisited this idea in crop production and proposed new practices based on selection for cooperative phenotypes, i.e. phenotypes that increase crop yield through decreased competitiveness. These recommendations, however, remain difficult to evaluate without a formal description of crop evolutionary dynamics under different selection strategies. Here, we develop a theoretical framework to investigate the evolution of cooperation-related traits in crops, using plant height as a case study. Our model is tailored to realistic agricultural practices and shows that combining high plant density, high relatedness and selection among groups favours the evolution of shorter plants that maximize grain yield. Our model allows us to revisit past and current breeding practices in light of kin selection theory, and yields practical recommendations to increase cooperation among crops and promote sustainable agriculture.

Pollen dispersal slows geographical range shift and accelerates ecological niche shift under climate change.

Species may survive climate change by migrating to track favorable climates and/or adapting to different climates. Several quantitative genetics models predict that species escaping extinction will change their geographical distribution while keeping the same ecological niche. We introduce pollen dispersal in these models, which affects gene flow but not directly colonization. We show that plant populations may escape extinction because of both spatial range and ecological niche shifts. Exact analytical formulas predict that increasing pollen dispersal distance slows the expected spatial range shift and accelerates the ecological niche shift. There is an optimal distance of pollen dispersal, which maximizes the sustainable rate of climate change. These conclusions hold in simulations relaxing several strong assumptions of our analytical model. Our results imply that, for plants with long distance of pollen dispersal, models assuming niche conservatism may not accurately predict their future distribution under climate change.

Regulated in development and DNA damage responses 1 (REDD1) links stress with IL-1ß-mediated familial Mediterranean fever attack through autophagy-driven neutrophil extracellular traps.

BACKGROUND: Familial Mediterranean fever (FMF) is an IL-1ß-dependent autoinflammatory disease caused by mutations of Mediterranean fever (MEFV) encoding pyrin and characterized by inflammatory attacks induced by physical or psychological stress. OBJECTIVE: We investigated the underlying mechanism that links stress-induced inflammatory attacks with neutrophil activation and release of IL-1ß-bearing neutrophil extracellular traps (NETs) in patients with FMF. METHODS: RNA sequencing was performed in peripheral neutrophils from 3 patients with FMF isolated both during attacks and remission, 8 patients in remission, and 8 healthy subjects. NET formation and proteins were analyzed by using confocal immunofluorescence microscopy, immunoblotting, myeloperoxidase-DNA complex ELISA, and flow cytometry. Samples from patients with Still's disease and bacterial infections were used also. RESULTS: The stress-related protein regulated in development and DNA damage responses 1 (REDD1) is significantly overexpressed during FMF attacks. Neutrophils from patients with FMF during remission are resistant to autophagy-mediated NET release, which can be overcome through REDD1 induction. Stress-related mediators (eg, epinephrine) decrease this threshold, leading to autophagy-driven NET release, whereas the synchronous inflammatory environment of FMF attack leads to intracellular production of IL-1ß and its release through NETs. REDD1 in autolysosomes colocalizes with pyrin and nucleotide-binding domain, leucine-rich repeat/pyrin domain-containing 3. Mutated pyrin prohibits this colocalization, leading to higher IL-1ß levels on NETs. CONCLUSIONS: This study provides a link between stress and initiation of inflammatory attacks in patients with FMF. REDD1 emerges as a regulator of neutrophil function upstream to pyrin, is involved in NET release and regulation of IL-1ß, and might constitute an important piece in the IL-1ß-mediated inflammation puzzle.

The Evolution of Mutual Mate Choice under Direct Benefits.

In nature, the intensity of mate choice (i.e., choosiness) is highly variable within and between sexes. Despite growing empirical evidence of male and/or mutual mate choice, theoretical investigations of the joint evolution of female and male choosiness are few. In addition, previous approaches have often assumed an absence of trade-off between the direct benefits per mating and the lower mating rate that results from being choosy. Here we model the joint evolution of female and male choosiness when it is solely ruled by this fundamental trade-off. We show that this trade-off can generate a diversity of stable combinations of choosiness. Mutual mate choice can evolve only if both females and males exhibit long latency after mating. Furthermore, we show that an increase in choosiness in one sex does not necessarily prevent the evolution of mutual mate choice; the outcome depends on details shaping the trade-off: the life history, the decision rule for mate choice, and how the fecundity of a pair is shaped by the quality of both individuals. Last, we discuss the power of the sensitivity of the relative searching time (i.e., of the proportion of a lifetime spent searching for mates) as a predictor of the joint evolution of choosiness.

Maximum-likelihood inference of population size contractions from microsatellite data.

Understanding the demographic history of populations and species is a central issue in evolutionary biology and molecular ecology. In this work, we develop a maximum-likelihood method for the inference of past changes in population size from microsatellite allelic data. Our method is based on importance sampling of gene genealogies, extended for new mutation models, notably the generalized stepwise mutation model (GSM). Using simulations, we test its performance to detect and characterize past reductions in population size. First, we test the estimation precision and confidence intervals coverage properties under ideal conditions, then we compare the accuracy of the estimation with another available method (MSVAR) and we finally test its robustness to misspecification of the mutational model and population structure. We show that our method is very competitive compared with alternative ones. Moreover, our implementation of a GSM allows more accurate analysis of microsatellite data, as we show that the violations of a single step mutation assumption induce very high bias toward false contraction detection rates. However, our simulation tests also showed some limits, which most importantly are large computation times for strong disequilibrium scenarios and a strong influence of some form of unaccounted population structure. This inference method is available in the latest implementation of the MIGRAINE software package.

Stable coexistence of incompatible Wolbachia along a narrow contact zone in mosquito field populations.

In arthropods, the intracellular bacteria Wolbachia often induce cytoplasmic incompatibility (CI) between sperm and egg, which causes conditional embryonic death and promotes the spatial spread of Wolbachia infections into host populations. The ability of Wolbachia to spread in natural populations through CI has attracted attention for using these bacteria in vector-borne disease control. The dynamics of incompatible Wolbachia infections have been deeply investigated theoretically, whereas in natural populations, there are only few examples described, especially among incompatible infected hosts. Here, we have surveyed the distribution of two molecular Wolbachia strains (wPip11 and wPip31) infecting the mosquito Culex pipiens in Tunisia. We delineated a clear spatial structure of both infections, with a sharp contact zone separating their distribution areas. Crossing experiments with isofemale lines from different localities showed three crossing types: wPip11-infected males always sterilize wPip31-infected females; however, while most wPip31-infected males were compatible with wPip11-infected females, a few completely sterilize them. The wPip11 strain was thus expected to spread, but temporal dynamics over 7 years of monitoring shows the stability of the contact zone. We examined which factors may contribute to the observed stability, both theoretically and empirically. Population cage experiments, field samples and modelling did not support significant impacts of local adaptation or assortative mating on the stability of wPip infection structure. By contrast, low dispersal probability and metapopulation dynamics in the host Cx. pipiens probably play major roles. This study highlights the need of understanding CI dynamics in natural populations to design effective and sustainable Wolbachia-based control strategies.

How choosy should I be? The relative searching time predicts evolution of choosiness under direct sexual selection.

Most theoretical research in sexual selection has focused on indirect selection. However, empirical studies have not strongly supported indirect selection. A well-established finding is that direct benefits and costs exert a strong influence on the evolution of mate choice. We present an analytical model in which unilateral mate choice evolves solely by direct sexual selection on choosiness. We show this is sufficient to generate the evolution of all possible levels of choosiness, because of the fundamental trade-off between mating rate and mating benefits. We further identify the relative searching time (RST, i.e. the proportion of lifetime devoted to searching for mates) as a predictor of the effect of any variable affecting the mating rate on the evolution of choosiness. We show that the RST: (i) allows one to make predictions about the evolution of choosiness across a wide variety of mating systems; (ii) encompasses all alternative variables proposed thus far to explain the evolution of choosiness by direct sexual selection; and (iii) can be empirically used to infer qualitative differences in choosiness.

Matrix inversions for chromosomal inversions: a method to construct summary statistics in complex coalescent models.

Chromosomal inversions allow genetic divergence of locally adapted populations by reducing recombination between chromosomes with different arrangements. While patterns of genetic variation within inverted regions are increasingly documented, inferential methods are largely missing to analyze such data. Previous work has provided expectations for coalescence patterns of neutral sites linked to an inversion polymorphism in two locally adapted populations. Here, we define a method to construct summary statistics in such complex population structure models. Under a scenario of selection on the inversion breakpoints, we first construct estimators of the migration rate between the two habitats, and of the recombination rate of a nucleotide site between the two inversion backgrounds. Next, we analyze the disequilibrium between two sites within an inversion and provide an estimator of the distinct recombination rate between these two sites in homokaryotypes and heterokaryotypes. These estimators should be suitable summary statistics for simulation-based methods that can handle the complex dependences in the data.

Quantifying maternal investment in mammals using allometry.

Maternal investment influences the survival and reproduction of both mothers and their progeny and plays a crucial role in understanding individuals' life-history and population ecology. To reveal the complex mechanisms associated with reproduction and investment, it is necessary to examine variations in maternal investment across species. Comparisons across species call for a standardised method to quantify maternal investment, which remained to be developed. This paper addresses this limitation by introducing the maternal investment metric - MI - for mammalian species, established through the allometric scaling of the litter mass at weaning age by the adult mass and investment duration (i.e. gestation + lactation duration) of a species. Using a database encompassing hundreds of mammalian species, we show that the metric is not highly sensitive to the regression method used to fit the allometric relationship or to the proxy used for adult body mass. The comparison of the maternal investment metric between mammalian subclasses and orders reveals strong differences across taxa. For example, our metric confirms that Eutheria have a higher maternal investment than Metatheria. We discuss how further research could use the maternal investment metric as a valuable tool to understand variation in reproductive strategies.

Cassette recombination dynamics within chromosomal integrons are regulated by toxin-antitoxin systems.

Integrons are adaptive bacterial devices that rearrange promoter-less gene cassettes into variable ordered arrays under stress conditions, thereby sampling combinatorial phenotypic diversity. Chromosomal integrons often carry hundreds of silent gene cassettes, with integrase-mediated recombination leading to rampant DNA excision and integration, posing a potential threat to genome integrity. How this activity is regulated and controlled, particularly through selective pressures, to maintain such large cassette arrays is unknown. Here, we show a key role of promoter-containing toxin-antitoxin (TA) cassettes as systems that kill the cell when the overall cassette excision rate is too high. These results highlight the importance of TA cassettes regulating the cassette recombination dynamics and provide insight into the evolution and success of integrons in bacterial genomes.

Likelihood-based inferences under isolation by distance: two-dimensional habitats and confidence intervals.

Likelihood-based methods of inference of population parameters from genetic data in structured populations have been implemented but still little tested in large networks of populations. In this work, a previous software implementation of inference in linear habitats is extended to two-dimensional habitats, and the coverage properties of confidence intervals are analyzed in both cases. Both standard likelihood and an efficient approximation are considered. The effects of misspecification of mutation model and dispersal distribution, and of spatial binning of samples, are considered. In the absence of model misspecification, the estimators have low bias, low mean square error, and the coverage properties of confidence intervals are consistent with theoretical expectations. Inferences of dispersal parameters and of the mutation rate are sensitive to misspecification or to approximations inherent to the coalescent algorithms used. In particular, coalescent approximations are not appropriate to infer the shape of the dispersal distribution. However, inferences of the neighborhood parameter (or of the product of population density and mean square dispersal rate) are generally robust with respect to complicating factors, such as misspecification of the mutation process and of the shape of the dispersal distribution, and with respect to spatial binning of samples. Likelihood inferences appear feasible in moderately sized networks of populations (up to 400 populations in this work), and they are more efficient than previous moment-based spatial regression method in realistic conditions.

The evolutionary success of regulated cell death in bacterial immunity.

Bacteria employ a complex arsenal of immune mechanisms to defend themselves against phages. Recent studies demonstrate that these immune mechanisms frequently involve regulated cell death in response to phage infection. By sacrificing infected cells, this strategy prevents the spread of phages within the surrounding population. In this review, we discuss the principles of regulated cell death in bacterial defense, and show that over 70% of sequenced prokaryotes employ this strategy as part of their defensive arsenals. We highlight the modularity of defense systems involving regulated cell death, explaining how shuffling between phage-sensing and cell-killing protein domains dominates their evolution. Some of these defense systems are the evolutionary ancestors of key components of eukaryotic immunity, highlighting their importance in shaping the evolutionary trajectory of immune systems across the tree of life.

Polymorphisms in Anopheles gambiae immune genes associated with natural resistance to Plasmodium falciparum.

Many genes involved in the immune response of Anopheles gambiae, the main malaria vector in Africa, have been identified, but whether naturally occurring polymorphisms in these genes underlie variation in resistance to the human malaria parasite, Plasmodium falciparum, is currently unknown. Here we carried out a candidate gene association study to identify single nucleotide polymorphisms (SNPs) associated with natural resistance to P. falciparum. A. gambiae M form mosquitoes from Cameroon were experimentally challenged with three local wild P. falciparum isolates. Statistical associations were assessed between 157 SNPs selected from a set of 67 A. gambiae immune-related genes and the level of infection. Isolate-specific associations were accounted for by including the effect of the isolate in the analysis. Five SNPs were significantly associated to the infection phenotype, located within or upstream of AgMDL1, CEC1, Sp PPO activate, Sp SNAKElike, and TOLL6. Low overall and local linkage disequilibrium indicated high specificity in the loci found. Association between infection phenotype and two SNPs was isolate-specific, providing the first evidence of vector genotype by parasite isolate interactions at the molecular level. Four SNPs were associated to either oocyst presence or load, indicating that the genetic basis of infection prevalence and intensity may differ. The validity of the approach was verified by confirming the functional role of Sp SNAKElike in gene silencing assays. These results strongly support the role of genetic variation within or near these five A. gambiae immune genes, in concert with other genes, in natural resistance to P. falciparum. They emphasize the need to distinguish between infection prevalence and intensity and to account for the genetic specificity of vector-parasite interactions in dissecting the genetic basis of Anopheles resistance to human malaria.