The effects of recombination, mutation and selection on the evolution of the Rp1 resistance genes in grasses

Means, motive and opportunity for biological invasions: Genetic introgression in a fungal pathogen

Invasions by fungal plant pathogens pose a significant threat to the health of agricultural ecosystems. Despite limited standing genetic variation, many invasive fungal species can adapt and spread rapidly, resulting in significant losses to crop yields. Here, we report on the population genomics of Colletotrichum truncatum, a polyphagous pathogen that can infect more than 460 plant species, and an invasive pathogen of soybean in Brazil. We study the whole-genome sequences of 18 isolates representing 10 fields from two major regions of soybean production. We show that Brazilian C. truncatum is subdivided into three phylogenetically distinct lineages that exchange genetic variation through hybridization. Introgression affects 2%-30% of the nucleotides of genomes and varies widely between the lineages. We find that introgressed regions comprise secreted protein-encoding genes, suggesting possible co-evolutionary targets for selection in those regions. We highlight the inherent vulnerability of genetically uniform crops in the agro-ecological environment, particularly when faced with pathogens that can take full advantage of the opportunities offered by an increasingly globalized world. Finally, we discuss "the means, motive and opportunity" of fungal pathogens and how they can become invasive species of crops. We call for more population genomic studies because such analyses can help identify geographical areas and pathogens that pose a risk, thereby helping to inform control strategies to better protect crops in the future.

Evolutionary dissection of monkeypox virus: Positive Darwinian selection drives the adaptation of virus-host interaction proteins

The emerging and ongoing outbreak of human monkeypox (hMPX) in 2022 is a serious global threat. An understanding of the evolution of the monkeypox virus (MPXV) at the single-gene level may provide clues for exploring the unique aspects of the current outbreak: rapidly expanding and sustained human-to-human transmission. For the current investigation, alleles of 156 MPXV coding genes (which account for >95% of the genomic sequence) have been gathered from roughly 1,500 isolates, including those responsible for the previous outbreaks. Using a range of molecular evolution approaches, we demonstrated that intra-species homologous recombination has a negligible effect on MPXV evolution. Despite the fact that the majority of the MPXV genes (64.10%) were subjected to negative selection at the whole gene level, 10 MPXV coding genes (MPXVgp004, 010, 012, 014, 044, 098, 138, 178, 188, and 191) were found to have a total of 15 codons or amino acid sites that are known to evolve under positive Darwinian selection. Except for MPXVgp138, almost all of these genes encode proteins that interact with the host. Of these, five ankyrin proteins (MPXVgp004, 010, 012, 178, and 188) and one Bcl-2-like protein (MPXVgp014) are involved in poxviruses' host range determination. We discovered that the majority (80%) of positive amino acid substitutions emerged several decades ago, indicating that these sites have been under constant selection pressure and that more adaptable alleles have been circulating in the natural reservoir. This finding was also supported by the minimum spanning networks of the gene alleles. The three positive amino acid substitutions (T/A426V in MPXVgp010, A423D in MPXVgp012, and S105L in MPXVgp191) appeared in 2019 or 2022, indicating that they would be crucial for the virus' eventual adaptation to humans. Protein modeling suggests that positive amino acid substitutions may affect protein functions in a variety of ways. Further study should focus on revealing the biological effects of positive amino acid substitutions in the genes for viral adaptation to humans, virulence, transmission, and so on. Our study advances knowledge of MPXV's adaptive mechanism and provides insights for exploring factors that are responsible for the unique aspects of the current outbreak.

Virulence effector SidJ evolution in Legionella pneumophila is driven by positive selection and intragenic recombination

Effector proteins translocated by the Dot/Icm type IV secretion system determine the virulence of Legionella pneumophila (L. pneumophila). Among these effectors, members of the SidE family (SidEs) regulate several cellular processes through a unique phosphoribosyl ubiquitination mechanism mediated by another effector, SidJ. Host-cell calmodulin (CaM) activates SidJ to glutamylate the SidEs of ubiquitin (Ub) ligases and to make a balanced Ub ligase activity. Given the central role of SidJ in this regulatory process, studying the nature of evolution of sidJ is important to understand the virulence of L. pneumophila and the interaction between the bacteria and its hosts. By studying sidJ from a large number of L. pneumophila strains and using various molecular evolution algorithms, we demonstrated that intragenic recombination drove the evolution of sidJ and contributed to sidJ diversification. Additionally, we showed that four codons of sidJ which are located in the N-terminal (NTD) (codons 58 and 200) and C-terminal (CTD) (codons 868 and 869) domains, but not in the kinase domain (KD) had been subjected to strong positive selection pressure, and variable mutation profiles of these codons were identified. Protein structural modeling of SidJ provided possible explanations for these mutations. Codons 868 and 869 mutations might engage in regulating the interactions of SidJ with CaM through hydrogen bonds and affect the CaM docking to SidJ. Mutation in codon 58 of SidJ might affect the distribution of main-chain atoms that are associated with the interaction with CaM. In contrast, mutations in codon 200 might influence the α-helix stability in the NTD. These mutations might be important to balance Ub ligase activity for different L. pneumophila hosts. This study first reported that intragenic recombination and positive Darwinian selection both shaped the genetic plasticity of sidJ, contributing to a deeper understanding of the adaptive mechanisms of this intracellular bacterium to different hosts.

The ETS-ETI cycle: evolutionary processes and metapopulation dynamics driving the diversification of pathogen effectors and host immune factors

The natural diversity of pathogen effectors and host immune components represents a snapshot of the underlying evolutionary processes driving the host-pathogen arms race. In plants, this arms race is manifested by an ongoing cycle of disease and resistance driven by pathogenic effectors that promote disease (effector-triggered susceptibility; ETS) and plant resistance proteins that recognize effector activity to trigger immunity (effector-triggered immunity; ETI). Here we discuss how this ongoing ETS-ETI cycle has shaped the natural diversity of both plant resistance proteins and pathogen effectors. We focus on the evolutionary forces that drive the diversification of the molecules that determine the outcome of plant-pathogen interactions and introduce the concept of metapopulation dynamics (i.e., the introduction of genetic variation from conspecific organisms in different populations) as an alternative mechanism that can introduce and maintain diversity in both host and pathogen populations.

A complex resistance locus in Solanum americanum recognizes a conserved Phytophthora effector

Late blight caused by Phytophthora infestans greatly constrains potato production. Many Resistance (R) genes were cloned from wild Solanum species and/or introduced into potato cultivars by breeding. However, individual R genes have been overcome by P. infestans evolution; durable resistance remains elusive. We positionally cloned a new R gene, Rpi-amr1, from Solanum americanum, that encodes an NRC helper-dependent CC-NLR protein. Rpi-amr1 confers resistance in potato to all 19 P. infestans isolates tested. Using association genomics and long-read RenSeq, we defined eight additional Rpi-amr1 alleles from different S. americanum and related species. Despite only ~90% identity between Rpi-amr1 proteins, all confer late blight resistance but differentially recognize Avramr1 orthologues and paralogues. We propose that Rpi-amr1 gene family diversity assists detection of diverse paralogues and alleles of the recognized effector, facilitating durable resistance against P. infestans.

Comparative genomics revealed adaptive admixture in Cryptosporidium hominis in Africa

Cryptosporidiosis is a major cause of diarrhoeal illness among African children, and is associated with childhood mortality, malnutrition, cognitive development and growth retardation. Cryptosporidium hominis is the dominant pathogen in Africa, and genotyping at the glycoprotein 60 (gp60) gene has revealed a complex distribution of different subtypes across this continent. However, a comprehensive exploration of the metapopulation structure and evolution based on whole-genome data has yet to be performed. Here, we sequenced and analysed the genomes of 26 C. hominis isolates, representing different gp60 subtypes, collected at rural sites in Gabon, Ghana, Madagascar and Tanzania. Phylogenetic and cluster analyses based on single-nucleotide polymorphisms showed that isolates predominantly clustered by their country of origin, irrespective of their gp60 subtype. We found a significant isolation-by-distance signature that shows the importance of local transmission, but we also detected evidence of hybridization between isolates of different geographical regions. We identified 37 outlier genes with exceptionally high nucleotide diversity, and this group is significantly enriched for genes encoding extracellular proteins and signal peptides. Furthermore, these genes are found more often than expected in recombinant regions, and they show a distinct signature of positive or balancing selection. We conclude that: (1) the metapopulation structure of C. hominis can only be accurately captured by whole-genome analyses; (2) local anthroponotic transmission underpins the spread of this pathogen in Africa; (3) hybridization occurs between distinct geographical lineages; and (4) genetic introgression provides novel substrate for positive or balancing selection in genes involved in host-parasite coevolution.

Genetic mapping and survey of powdery mildew resistance in the wild Central Asian ancestor of cultivated grapevines in Central Asia

Cultivated grapevines (Vitis vinifera) lack resistance to powdery mildew (PM) with few exceptions. Resistance to this pathogen within V. vinifera has been reported in earlier studies and identified as the Ren1 locus in two Central Asian table grape accessions. Other PM-resistant cultivated varieties and accessions of the wild ancestor V. vinifera subsp. sylvestris were soon identified raising questions regarding the origin of the resistance. In this study, F1 breeding populations were developed with a PM susceptible V. vinifera subsp. vinifera breeding line and a PM-resistant subsp. sylvestris accession. Genotyping was carried out with five Ren1 locus linked SSR markers. A PM resistance locus explaining up to 96% of the phenotypic variation was identified in the same genomic position, where the Ren1 locus was previously reported. New SSR marker alleles linked with the resistance locus were identified. We report results of PM resistance in multiple accessions of subsp. sylvestris collected as seed lots or cuttings from five countries in the Caucasus and Central Asia. A total of 20 females from 11 seed lots and 19 males from nine seed lots collected from Georgia, Armenia, and Azerbaijan were resistant to PM. Three male and one female plant collected as cuttings from Afghanistan and Iran were also resistant to PM. Allelic analysis of markers linked with the Ren1 locus in conjunction with disease evaluation data found a high diversity of allelic haplotypes, which are only possible via recombination events occurring over a long time period. Sequence analysis of two alleles of the SSR marker that cosegregates with the resistance found SNPs that were present in the wild progenitor and in cultivated forms. Variable levels of PM resistance among the tested accessions were also observed. These lines of evidence suggest that the powdery mildew fungus may have been present in Asia for a longer time than currently thought, giving the wild progenitor V. vinifera subsp. sylvestris time to coevolve with and develop resistance to this pathogen.

Evidence for Adaptive Introgression of Disease Resistance Genes Among Closely Related Arabidopsis Species

The generation and maintenance of functional variation in the pathogen defense system of plants is central to the constant evolutionary battle between hosts and parasites. If a species is susceptible to a given pathogen, hybridization and subsequent introgression of a resistance allele from a related species can potentially be an important source of new immunity and is therefore expected to be selected for in a process referred to as adaptive introgression. Here, we survey sequence variation in 10 resistance (R-) genes and compare them with 37 reference genes in natural populations of the two closely related and interfertile species: Arabidopsis lyrata and A. halleri. The R-genes are highly polymorphic in both species and show clear signs of trans-species polymorphisms. We show that A. lyrata and A. halleri have had a history of limited introgression for the reference genes. For the R-genes, the introgression rate has been significantly higher than for the reference genes, resulting in fewer fixed differences between species and a higher sharing of identical haplotypes. We conclude that R-genes likely cross the species boundaries at a higher rate than reference genes and therefore also that some of the increased diversity and trans-specific polymorphisms in R-genes is due to adaptive introgression.

Early Duplication of a Single MHC IIB Locus Prior to the Passerine Radiations

A key characteristic of MHC genes is the persistence of allelic lineages over macroevolutionary periods, often through multiple speciation events. This phenomenon, known as trans-species polymorphism (TSP), is well documented in several major taxonomic groups, but has less frequently been observed in birds. The order Passeriformes is arguably the most successful terrestrial vertebrate order in terms of diversity of species and ecological range, but the reasons for this success remain unclear. Passerines exhibit the most highly duplicated MHC genes of any major vertebrate taxonomic group, which may generate increased immune response relative to other avian orders with fewer MHC loci. Here, we describe phylogenetic patterns of the MHC IIB in the passerine family Corvidae. Our results indicate wide-spread TSP within this family, with at least four supported MHC IIB allelic lineages that predate speciation by many millions of years. Markov chain Monte Carlo simulations indicate that divergence of these lineages occurred near the time of the divergence of the Passeriformes and other avian orders. We suggest that the current MHC diversity observed in passerines is due in part to the multiple duplication of a single MHC locus, DAB1, early in passerine evolution and that subsequent duplications of these paralogues have contributed to the enormous success of this order by increasing their ability to recognize and mount immune responses to novel pathogens.

The wheat homolog of putative nucleotide-binding site-leucine-rich repeat resistance gene TaRGA contributes to resistance against powdery mildew

Powdery mildew, one of the most destructive wheat diseases worldwide, is caused by Blumeria graminis f. sp. tritici (Bgt), a fungal species with a consistently high mutation rate that makes individual resistance (R) genes ineffective. Therefore, effective resistance-related gene cloning is vital for breeding and studying the resistance mechanisms of the disease. In this study, a putative nucleotide-binding site-leucine-rich repeat (NBS-LRR) R gene (TaRGA) was cloned using a homology-based cloning strategy and analyzed for its effect on powdery mildew disease and wheat defense responses. Real-time reverse transcription-PCR (RT-PCR) analyses revealed that a Bgt isolate 15 and salicylic acid stimulation significantly induced TaRGA in the resistant variety. Furthermore, the silencing of TaRGA in powdery mildew-resistant plants increased susceptibility to Bgt15 and prompted conidia propagation at the infection site. However, the expression of TaRGA in leaf segments after single-cell transient expression assay highly increased the defense responses to Bgt15 by enhancing callose deposition and phenolic autofluorogen accumulation at the pathogen invading sites. Meanwhile, the expression of pathogenesis-related genes decreased in the TaRGA-silenced plants and increased in the TaRGA-transient-overexpressing leaf segments. These results implied that the TaRGA gene positively regulates the defense response to powdery mildew disease in wheat.

HYBRIDCHECK: software for the rapid detection, visualization and dating of recombinant regions in genome sequence data

HYBRIDCHECK is a software package to visualize the recombination signal in large DNA sequence data set, and it can be used to analyse recombination, genetic introgression, hybridization and horizontal gene transfer. It can scan large (multiple kb) contigs and whole-genome sequences of three or more individuals. HYBRIDCHECK is written in the r software for OS X, Linux and Windows operating systems, and it has a simple graphical user interface. In addition, the r code can be readily incorporated in scripts and analysis pipelines. HYBRIDCHECK implements several ABBA-BABA tests and visualizes the effects of hybridization and the resulting mosaic-like genome structure in high-density graphics. The package also reports the following: (i) the breakpoint positions, (ii) the number of mutations in each introgressed block, (iii) the probability that the identified region is not caused by recombination and (iv) the estimated age of each recombination event. The divergence times between the donor and recombinant sequence are calculated using a JC, K80, F81, HKY or GTR correction, and the dating algorithm is exceedingly fast. By estimating the coalescence time of introgressed blocks, it is possible to distinguish between hybridization and incomplete lineage sorting. HYBRIDCHECK is libré software and it and its manual are free to download from http://ward9250.github.io/HybridCheck/.