Congruent population genetic structures and divergence histories in anther‐smut fungi and their host plants Silene italica and the Silene nutans species complex

Discordant population structure among rhizobium divided genomes and their legume hosts

Symbiosis often occurs between partners with distinct life history characteristics and dispersal mechanisms. Many bacterial symbionts have genomes comprising multiple replicons with distinct rates of evolution and horizontal transmission. Such differences might drive differences in population structure between hosts and symbionts and among the elements of the divided genomes of bacterial symbionts. These differences might, in turn, shape the evolution of symbiotic interactions and bacterial evolution. Here we use whole genome resequencing of a hierarchically structured sample of 191 strains of Sinorhizobium meliloti collected from 21 locations in southern Europe to characterize population structures of this bacterial symbiont, which forms a root nodule symbiosis with the host plant Medicago truncatula. S. meliloti genomes showed high local (within-site) variation and little isolation by distance. This was particularly true for the two symbiosis elements, pSymA and pSymB, which have population structures that are similar to each other, but distinct from both the bacterial chromosome and the host plant. Given limited recombination on the chromosome, compared to the symbiosis elements, distinct population structures may result from differences in effective gene flow. Alternatively, positive or purifying selection, with little recombination, may explain distinct geographical patterns at the chromosome. Discordant population structure between hosts and symbionts indicates that geographically and genetically distinct host populations in different parts of the range might interact with genetically similar symbionts, potentially minimizing local specialization.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Pollinator sharing between reproductively isolated genetic lineages of Silene nutans

High reciprocal pollination specialization leading to pollinator isolation can prevent interspecific pollen transfer and competition for pollinators. Sharing pollinators may induce mating costs, but it may also increase pollination services and pollen dispersal and offer more resources to pollinators, which may be important in case of habitat fragmentation leading to pollination disruption. We estimated pollen dispersal and pollinator isolation or sharing between two reproductively isolated genetic lineages of Silene nutans (Caryophyllaceae), which are rare and occur in parapatry in southern Belgium, forming two edaphic ecotypes. As inter-ecotypic crosses may lead to pollen wastage and inviable progeny, pollinator isolation might have evolved between ecotypes. Silene nutans is mainly pollinated by nocturnal moths, including nursery pollinators, which pollinate and lay their eggs in flowers, and whose caterpillars feed on flowers and seeds. Pollinator assemblages of the two ecotypes are largely unknown and inter-ecotypic pollen flows have never been investigated. Fluorescent powdered dyes were used as pollen analogues to quantify intra- and inter-ecotypic pollen transfers and seeds were germinated to detect chlorotic seedlings resulting from inter-ecotypic pollination. Nocturnal pollinators were observed using infrared cameras on the field, and seed-eating caterpillars were collected and reared to identify nursery pollinator species. No pollinator isolation was found: we detected long-distance (up to 5 km) inter-ecotypic dye transfers and chlorotic seedlings, indicating inter-ecotypic fertilization events. The rare moth Hadena albimacula, a nursery pollinator specialized on S. nutans, was found on both ecotypes, as well as adults visiting flowers (cameras recordings) as seed-eating caterpillars. However, S. nutans populations harbor different abundance and diversity of seed predator communities, including other rare nursery pollinators, suggesting a need for distinct conservation strategies. Our findings demonstrate the efficiency of moths, especially of nursery pollinators, to disperse pollen over long distances in natural landscapes, so to ensure gene flow and population sustainability of the host plant. Seed-predator specificities between the two reproductively isolated genetic lineages of S. nutans, and pollinator sharing instead of pollinator isolation when plants occur in parapatry, suggest that conservation of the host plant is also essential for sustaining (rare) pollinator and seed predator communities.

Unraveling coevolutionary dynamics using ecological genomics

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

Exploring Mitogenomes Diversity of Fusarium musae from Banana Fruits and Human Patients

Fusarium musae has recently been described as a cross-kingdom pathogen causing post-harvest disease in bananas and systemic and superficial infection in humans. The taxonomic identity of fungal cross-kingdom pathogens is essential for confirming the identification of the species on distant infected hosts. Understanding the level of variability within the species is essential to decipher the population homogeneity infecting human and plant hosts. In order to verify that F. musae strains isolated from fruits and patients are part of a common population and to estimate their overall diversity, we assembled, annotated and explored the diversity of the mitogenomes of 18 F. musae strains obtained from banana fruits and human patients. The mitogenomes showed a high level of similarity among strains with different hosts’ origins, with sizes ranging from 56,493 to 59,256 bp. All contained 27 tRNA genes and 14 protein-coding genes, rps3 protein, and small and large ribosomal subunits (rns and rnl). Variations in the number of endonucleases were detected. A comparison of mitochondrial endonucleases distribution with a diverse set of Fusarium mitogenomes allowed us to specifically discriminate F. musae from its sister species F. verticillioides and the other Fusarium species. Despite the diversity in F. musae mitochondria, strains from bananas and strains from human patients group together, indirectly confirming F. musae as a cross-kingdom pathogen.

Historical Differentiation and Recent Hybridization in Natural Populations of the Nematode-Trapping Fungus Arthrobotrys oligospora in China

Maintaining the effects of nematode-trapping fungi (NTF) agents in order to control plant-parasitic nematodes (PPNs) in different ecological environments has been a major challenge in biological control applications. To achieve such an objective, it is important to understand how populations of the biocontrol agent NTF are geographically and ecologically structured. A previous study reported evidence for ecological adaptation in the model NTF species Arthrobotrys oligospora. However, their large-scale geographic structure, patterns of gene flow, their potential phenotypic diversification, and host specialization remain largely unknown. In this study, we developed a new panel of 20 polymorphic short tandem repeat (STR) markers and analyzed 239 isolates of A. oligospora from 19 geographic populations in China. In addition, DNA sequences at six nuclear gene loci and strain mating types (MAT) were obtained for these strains. Our analyses suggest historical divergence within the A. oligospora population in China. The genetically differentiated populations also showed phenotypic differences that may be related to their ecological adaptations. Interestingly, our analyses identified evidence for recent dispersion and hybridization among the historically subdivided geographic populations in nature. Together, our results indicate a changing population structure of A. oligospora in China and that care must be taken in selecting the appropriate strains as biocontrol agents that can effectively reproduce in agriculture soil while maintaining their nematode-trapping ability.

Novel genomic approaches to study antagonistic coevolution between hosts and parasites

Host-parasite coevolution is ubiquitous, shaping genetic and phenotypic diversity and the evolutionary trajectory of interacting species. With the advances of high throughput sequencing technologies applicable to model and non-model organisms alike, it is now feasible to study in greater detail (a) the genetic underpinnings of coevolution, (b) the speed and type of dynamics at coevolving loci, and (c) the genomic consequences of coevolution. This review focuses on three recently developed approaches that leverage information from host and parasite full genome data simultaneously to pinpoint coevolving loci and draw inference on the coevolutionary history. First, co-genome-wide association study (co-GWAS) methods allow pinpointing the loci underlying host-parasite interactions. These methods focus on detecting associations between genetic variants and the outcome of experimental infection tests or on correlations between genomes of naturally infected hosts and their infecting parasites. Second, extensions to population genomics methods can detect genes under coevolution and infer the coevolutionary history, such as fitness costs. Third, correlations between host and parasite population size in time are indicative of coevolution, and polymorphism levels across independent spatially distributed populations of hosts and parasites can reveal coevolutionary loci and infer coevolutionary history. We describe the principles of these three approaches and discuss their advantages and limitations based on coevolutionary theory. We present recommendations for their application to various host (prokaryotes, fungi, plants, and animals) and parasite (viruses, bacteria, fungi, and macroparasites) species. We conclude by pointing out methodological and theoretical gaps to be filled to extract maximum information from full genome data and thereby to shed light on the molecular underpinnings of coevolution.

Comparison of microsatellite distribution patterns in twenty-nine beetle genomes

Simple sequence repeats (SSRs) represent an important source of genetic variation that provides a basis for adaptation to different environments in organisms. In this study, we examined the distribution patterns of SSRs in twenty-nine beetle genomes and carried out Gene Ontology (GO) analysis of CDSs embedded with perfect SSRs (P-SSRs). The results demonstrated that imperfect SSRs (I-SSRs) represented the most abundant SSR category in beetle genomes and in different genomic regions (CDS, exon, and intron regions). The numbers of P-SSRs, I-SSRs, compound SSRs, and variable number tandem repeats were positively correlated with beetle genome size, whereas neither the frequency nor the density of the SSRs was correlated with genome size. Moreover, our results demonstrated that common genomic features of P-SSRs within the same suborder or family of Coleoptera were rare. Mono-, di-, tri-, or tetranucleotide SSRs were the most abundant P-SSR categories in beetle genomes. The preferred predominant repeat motif among the mononucleotide P-SSRs was (A)n, but the most frequent repeat motifs for other length classes varied differentially among these genomes. Furthermore, the P-SSR type with the highest GC content differed in the beetle genomes and in different genomic regions. CV (coefficient of variability) analysis demonstrated that the repeat copy numbers of P-SSRs presented relatively higher variation in introns than in CDSs and exons. The GO terms of CDSs containing P-SSRs for molecular functions were mainly enriched in "binding" and "transcription". Our findings will be useful for studying the functional roles of microsatellite heterogeneity in beetle adaptation.