Reproductive clonality of pathogens: a perspective on pathogenic viruses, bacteria, fungi, and parasitic protozoa

Diversity and dissemination of viruses in pathogenic protozoa

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region.

Forest health in the Anthropocene: the emergence of a novel tree disease is associated with poplar cultivation

Plant domestication and movement are large contributors to the success of new diseases. The introduction of new host species can result in accelerated evolutionary changes in pathogens, affecting long-established coevolutionary dynamics. This has been observed in poplars where severe epidemics of pathogens that were innocuous in their natural pathosystems occurred following host domestication. The North American fungus Sphaerulina musiva is responsible for endemic leaf spots on Populus deltoides. We show that the expansion of poplar cultivation resulted in the emergence of a new lineage of this pathogen that causes stem infections on a new host, P. balsamifera. This suggests a host shift since this is not a known host. Genome analysis of this emerging lineage reveals a mosaic pattern with islands of diversity separated by fixed genome regions, which is consistent with a homoploid hybridization event between two individuals that produced a hybrid swarm. Genome regions of extreme divergence and low diversity are enriched in genes involved in host-pathogen interactions. The specialization of this emerging lineage to a new host and its clonal propagation represents a serious threat to poplars and could affect both natural and planted forests. This work provides a clear example of the changes created by the intensification of tree cultivation that facilitate the emergence of specialized pathogens, jeopardizing the natural equilibrium between hosts and pathogens. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

Parasite hybridization promotes spreading of endosymbiotic viruses

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of endosymbiotic viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis parasites and their endosymbiotic Leishmania RNA virus. We show that parasite populations circulate in isolated pockets of suitable habitat and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites were geographically and ecologically dispersed, and commonly infected from a pool of genetically diverse viruses. Our results suggest that parasite hybridization, likely due to increased human migration and ecological perturbations, increased the frequency of endosymbiotic interactions known to play a key role in disease severity.

High-Throughput Microsatellite Markers Development for Genetic Characterization of Emerging Sporothrix Species

Sporotrichosis is the main subcutaneous mycosis worldwide transmitted by animal or plant vectors and often escalates to outbreaks or epidemics. The current cat-transmitted sporotrichosis driven by Sporothrix brasiliensis has become a significant public health issue in South America. Transmission dynamics remain enigmatic due to the lack of development of polymorphic markers for molecular epidemiological analysis. This study used a high-throughput mining strategy to characterize simple sequence repeat (SSR) markers from Sporothrix genomes. A total of 118,140-143,912 SSR loci were identified (82,841-98,369 unique markers), with a 3651.55-3804.65 SSR/Mb density and a majority of dinucleotides motifs (GC/CG). We developed a panel of 15 highly polymorphic SSR markers suitable for genotyping S. brasiliensis, S. schenckii, and S. globosa. PCR amplification revealed 240 alleles in 180 Sporothrix isolates with excellent polymorphic information content (PIC = 0.9101), expected heterozygosity (H = 0.9159), and discriminating power (D = 0.7127), supporting the effectiveness of SSR markers in uncovering cryptic genetic diversity. A systematic population genetic study estimated three clusters, corresponding to S. brasiliensis (population 1, n = 97), S. schenckii (population 2, n = 49), and S. globosa (population 3, n = 34), with a weak signature of mixed ancestry between populations 1 and 2 or 3 and 2. Partitioning of genetic variation via AMOVA revealed highly structured populations (ΦPT = 0.539; Nm = 0.213; p < 0.0001), with approximately equivalent genetic variability within (46%) and between (54%) populations. Analysis of SSR diversity supports Rio de Janeiro (RJ) as the center of origin for contemporary S. brasiliensis infections. The recent emergence of cat-transmitted sporotrichosis in northeastern Brazil indicates an RJ-Northeast migration resulting in founder effects during the introduction of diseased animals into sporotrichosis-free areas. Our results demonstrated high cross-species transferability, reproducibility, and informativeness of SSR genetic markers, helping dissect deep and fine-scale genetic structures and guiding decision making to mitigate the harmful effects of the expansion of cat-transmitted sporotrichosis.

Subspecific Nomenclature of Giardia duodenalis in the Light of a Compared Population Genomics of Pathogens

Genetic and genomic data have long recognized that the species Giardia duodenalis is subdivided into at least eight genetic clusters that have been named "assemblages" by specialists in the field. Some of these assemblages have been given the status of species, with Linnean binames. In the framework of the predominant clonal evolution model (PCE), we have shown that, from an evolutionary point of view, G. duodenalis assemblages are equatable to "near-clades", that is to say: clades whose discreteness is somewhat clouded by occasional genetic exchange, but remain discrete and stable in space and time. The implications of this evolutionary status for the species described within G. duodenalis are discussed in light of the most recent genetic and genomic studies. The pattern of this species' subspecific genetic variability and genetic clustering appears to be very similar to the ones of various parasitic, fungal and bacteria species. This underlines the relevance of a compared population genomics of pathogenic species allowed by the broad framework of the PCE model.

A unified global genotyping framework of dengue virus serotype-1 for a stratified coordinated surveillance strategy of dengue epidemics

Background

Dengue is the fastest spreading arboviral disease, posing great challenges on global public health. A reproduceable and comparable global genotyping framework for contextualizing spatiotemporal epidemiological data of dengue virus (DENV) is essential for research studies and collaborative surveillance.

Methods

Targeting DENV-1 spreading prominently in recent decades, by reconciling all qualified complete E gene sequences of 5003 DENV-1 strains with epidemiological information from 78 epidemic countries/areas ranging from 1944 to 2018, we established and characterized a unified global high-resolution genotyping framework using phylogenetics, population genetics, phylogeography, and phylodynamics.

Results

The defined framework was discriminated with three hierarchical layers of genotype, subgenotype and clade with respective mean pairwise distances 2–6%, 0.8–2%, and ≤ 0.8%. The global epidemic patterns of DENV-1 showed strong geographic constraints representing stratified spatial-genetic epidemic pairs of Continent-Genotype, Region-Subgenotype and Nation-Clade, thereby identifying 12 epidemic regions which prospectively facilitates the region-based coordination. The increasing cross-transmission trends were also demonstrated. The traditional endemic countries such as Thailand, Vietnam and Indonesia displayed as persisting dominant source centers, while the emerging epidemic countries such as China, Australia, and the USA, where dengue outbreaks were frequently triggered by importation, showed a growing trend of DENV-1 diffusion. The probably hidden epidemics were found especially in Africa and India. Then, our framework can be utilized in an accurate stratified coordinated surveillance based on the defined viral population compositions. Thereby it is prospectively valuable for further hampering the ongoing transition process of epidemic to endemic, addressing the issue of inadequate monitoring, and warning us to be concerned about the cross-national, cross-regional, and cross-continental diffusions of dengue, which can potentially trigger large epidemics.

Conclusions

The framework and its utilization in quantitatively assessing DENV-1 epidemics has laid a foundation and re-unveiled the urgency for establishing a stratified coordinated surveillance platform for blocking global spreading of dengue. This framework is also expected to bridge classical DENV-1 genotyping with genomic epidemiology and risk modeling. We will promote it to the public and update it periodically.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40249-022-01024-5.

Clonality, inbreeding, and hybridization in two extremotolerant black yeasts

BACKGROUND

The great diversity of lifestyles and survival strategies observed in fungi is reflected in the many ways in which they reproduce and recombine. Although a complete absence of recombination is rare, it has been reported for some species, among them 2 extremotolerant black yeasts from Dothideomycetes: Hortaea werneckii and Aureobasidium melanogenum. Therefore, the presence of diploid strains in these species cannot be explained as the product of conventional sexual reproduction.

RESULTS

Genome sequencing revealed that the ratio of diploid to haploid strains in both H. werneckii and A. melanogenum is about 2:1. Linkage disequilibrium between pairs of polymorphic loci and a high degree of concordance between the phylogenies of different genomic regions confirmed that both species are clonal. Heterozygosity of diploid strains is high, with several hybridizing genome pairs reaching the intergenomic distances typically seen between different fungal species. The origin of diploid strains collected worldwide can be traced to a handful of hybridization events that produced diploids, which were stable over long periods of time and distributed over large geographic areas.

CONCLUSIONS

Our results, based on the genomes of over 100 strains of 2 black yeasts, show that although they are clonal, they occasionally form stable and highly heterozygous diploid intraspecific hybrids. The mechanism of these apparently rare hybridization events, which are not followed by meiosis or haploidization, remains unknown. Both extremotolerant yeasts, H. werneckii and even more so A. melanogenum, a close relative of the intensely recombining and biotechnologically relevant Aureobasidium pullulans, provide an attractive model for studying the role of clonality and ploidy in extremotolerant fungi.

Characterization of fungal pathogens and germplasm screening for disease resistance in the main production area of the common bean in Argentina

The common bean (Phaseolus vulgaris L.) is the most important grain legume in the human diet, mainly in Africa and Latin America. Argentina is one of the five major producers of the common bean in the world, and the main cultivation areas are concentrated in the northwestern provinces of this country. Crop production of the common bean is often affected by biotic factors like some endemic fungal diseases, which exert a major economic impact on the region. The most important fungal diseases affecting the common bean in Argentina are white mold caused by Sclerotinia sclerotiorum, angular leaf spot caused by Pseudocercospora griseola, web blight and root rot caused by Rhizoctonia solani, which can cause production losses of up to 100% in the region. At the present, the most effective strategy for controlling these diseases is the use of genetic resistance. In this sense, population study and characterization of fungal pathogens are essential for developing cultivars with durable resistance. In this review we report diversity studies carried out on these three fungal pathogens affecting the common bean in northwestern Argentina, analyzing more than 200 isolates by means of molecular, morphological and pathogenic approaches. Also, the screening of physiological resistance in several common bean commercial lines and wild native germplasm is reviewed. This review contributes to the development of sustainable management strategies and cultural practices in bean production aimed to minimize yield losses due to fungal diseases in the common bean.

Microevolution and subspecific taxonomy of Trypanosoma cruzi

Trypanosoma cruzi, the agent of Chagas disease, is a highly polymorphic species, subdivided into 6 main evolutionary lineages or near-clades (formerly discrete typing units or DTUs). An additional near-clade (TC-bat) has recently been evidenced. This pattern is considered to be the result of predominant clonal evolution (PCE). PCE is compatible with occasional mating/hybridization, which do not break the prevalent pattern of clonal evolution, the main trait of it being the presence of Multigene Bifurcating Trees (MGBTs) at all evolutionary levels ("clonal frame"). The development of highly resolutive genetic (microsatellites*) and genomic (sequencing and multi-single nucleotide polymorphism {SNP}* typing) markers shows that PCE also operates at a microevolutionary* level within each of the near-clades ("Russian doll pattern"), in spite of occasional meiosis and hybridization events. Within each near-clade, one can evidence widespread clonal multilocus genotypes*, linkage disequilibrium*, Multigene Bifurcating Trees and lesser near-clades. The within near-clade population structure is like a miniature picture of that of the whole species, suggesting gradual rather than saltatory evolution. Additional data are required to evaluate the stability of these lesser near-clades in the long run and to evaluate the need for an adequate nomenclature for this microevolutionary level.

Genomic diversity and genetic variation of Leishmania panamensis within its endemic range

Species belonging to the Leishmania (Viannia) subgenus are important causative agents of cutaneous and mucocutaneous leishmaniasis in Central and South America. These parasites possess several distinctive biological features that are influenced by their genetics, population structure, and genome instability. To date, several studies have revealed varying degrees of genetic diversity within Leishmania species. Particularly, in species of the L. (Viannia) subgenus, a generalized high intraspecific genetic diversity has been reported, although, conflicting conclusions have been drawn using different molecular techniques. Despite being the most common Leishmania species circulating in Panama and Colombia, few studies have analyzed clinical samples of Leishmania panamensis using whole-genome sequencing, and their restricted number of samples has limited the information they can provide to understand the population structure of L. panamensis. Here, we used next generation sequencing (NGS) to explore the genetic diversity of L. panamensis within its endemic range, analyzing data from 43 isolates of Colombian and Panamanian origin. Our results show the occurrence of three well-defined geographically correlated groups, and suggests the possible occurrence of additional phylogeographic groups. Furthermore, these results support the existence of a mixed mode of reproduction in L. panamensis, with varying frequencies of events of genetic recombination occurring primarily within subpopulations of closely related strains. This study offers important insights into the population genetics and reproduction mode of L. panamensis, paving the way to better understand their population structure and the emergence and maintenance of key eco-epidemiological traits.

Microevolution of Trypanosoma cruzi reveals hybridization and clonal mechanisms driving rapid genome diversification

Protozoa and fungi are known to have extraordinarily diverse mechanisms of genetic exchange. However, the presence and epidemiological relevance of genetic exchange in Trypanosoma cruzi, the agent of Chagas disease, has been controversial and debated for many years. Field studies have identified both predominantly clonal and sexually recombining natural populations. Two of six natural T. cruzi lineages (TcV and TcVI) show hybrid mosaicism, using analysis of single-gene locus markers. The formation of hybrid strains in vitro has been achieved and this provides a framework to study the mechanisms and adaptive significance of genetic exchange. Using whole genome sequencing of a set of experimental hybrids strains, we have confirmed that hybrid formation initially results in tetraploid parasites. The hybrid progeny showed novel mutations that were not attributable to either (diploid) parent showing an increase in amino acid changes. In long-term culture, up to 800 generations, there was a variable but gradual erosion of progeny genomes towards triploidy, yet retention of elevated copy number was observed at several core housekeeping loci. Our findings indicate hybrid formation by fusion of diploid T. cruzi, followed by sporadic genome erosion, but with substantial potential for adaptive evolution, as has been described as a genetic feature of other organisms, such as some fungi.

Migration Drives the Replacement of Xanthomonas perforans Races in the Absence of Widely Deployed Resistance

Changes in Xanthomonas race and species composition causing bacterial spot of tomato have occurred throughout the world and are often associated with epidemics. Knowledge of bacterial population structure is key for resistance discovery and deployment. We surveyed Xanthomonas spp. composition from processing tomato fields in the Midwestern United States over a 4-year period between 2017 and 2020, compared these to strains collected previously, and found that X. perforans is currently the most prevalent species. We characterized 564 X. perforans isolates for sequence variation in avrXv3 to distinguish between race T3 and T4 and validated race designation using hypersensitive response (HR) assays for 106 isolates. Race T4 accounted for over 95% of X. perforans isolates collected in the Midwest between 2017 and 2020. Whole genome sequencing, Average Nucleotide Identity (ANI) analysis, core genome alignment and single nucleotide polymorphism (SNP) detection relative to a reference strain, and phylogenomic analysis suggest that the majority of Midwestern X. perforans strains collected between 2017 and 2020 were nearly identical, with greater than 99.99% ANI to X. perforans isolates collected from Collier County, Florida in 2012. These isolates shared a common SNP variant resulting an a premature stop codon in avrXv3. One sequenced isolate was identified with a deletion of avrXv3 and shared 99.99% ANI with a strain collected in Collier Co., Florida in 2006. A population shift to X. perforans T4 occurred in the absence of widely deployed resistance, with only 7% of tomato varieties tested having the resistant allele at the Xv3/Rx-4 locus. The persistence of nearly identical strains over multiple years suggests that migration led to the establishment of an endemic population. Our findings validate a genomics-based framework to track shifts in X. perforans populations due to migration, mutation, drift, or selection based on comparisons to 146 genomes.

Population genomics reveals historical and ongoing recombination in the Fusarium oxysporum species complex

The Fusarium oxysporum species complex (FOSC) is a group of closely related plant pathogens long-considered strictly clonal, as sexual stages have never been recorded. Several studies have questioned whether recombination occurs in FOSC, and if it occurs its nature and frequency are unknown. We analysed 410 assembled genomes to answer whether FOSC diversified by occasional sexual reproduction interspersed with numerous cycles of asexual reproduction akin to a model of predominant clonal evolution (PCE). We tested the hypothesis that sexual reproduction occurred in the evolutionary history of FOSC by examining the distribution of idiomorphs at the mating locus, phylogenetic conflict and independent measures of recombination from genome-wide SNPs and genes. A phylogenomic dataset of 40 single copy orthologs was used to define structure a priori within FOSC based on genealogical concordance. Recombination within FOSC was tested using the pairwise homoplasy index and divergence ages were estimated by molecular dating. We called SNPs from assembled genomes using a k-mer approach and tested for significant linkage disequilibrium as an indication of PCE. We clone-corrected and tested whether SNPs were randomly associated as an indication of recombination. Our analyses provide evidence for sexual or parasexual reproduction within, but not between, clades of FOSC that diversified from a most recent common ancestor about 500 000 years ago. There was no evidence of substructure based on geography or host that might indicate how clades diversified. Competing evolutionary hypotheses for FOSC are discussed in the context of our results.

Analyses of Wheat Yellow Rust Populations Reveal Sexual Recombination and Seasonal Migration Pattern of Puccinia striiformis f. sp. tritici in Gangu, Northwestern China

Puccinia striiformis f. sp. tritici is the causal agent of wheat yellow rust with records of regular and severe epidemics in China. This study explored the population dynamics of the yellow rust pathogen in Gangu, northwestern China. In Gangu, the Weihe River runs from west to east and divides Gangu into three regions: North and South mountain, with the valley in between. To study the genetic structure of the pathogen in local populations, samples were collected over 3 years from the three regions at different altitudes both within and between the wheat cropping seasons. A total of 811 P. striiformis f. sp. tritici isolates were successfully genotyped using 16 simple sequence repeat markers. The results suggest that P. striiformis f. sp. tritici can survive year-round in Gangu. The P. striiformis f. sp. tritici populations migrated among the regions, and the migration pattern was not related to altitude. The oversummering populations in the North and South mountain regions were genetically different from each other; and the P. striiformis f. sp. tritici populations collected from the lower altitude in the valley had no relationship with any of the populations collected in the spring or fall, indicating that they too have a different origin. Signatures of random mating were found in the populations collected in both North and South mountain regions, but not in the valley populations.

Molecular Detection and Genotyping of Enterocytozoon bieneusi in Black Goats (Capra hircus) in Yunnan Province, Southwestern China

Simple Summary

Enterocytozoon bieneusi is one of the most common parasites in human and animals, and a threat to public health. So far, no data are available for E. bieneusi prevalence and genotypes in black goats in Yunnan Province, Southwestern China. Therefore, the objective of this study was to detect the prevalence and genotypes of E. bieneusi by examining 907 fecal samples collected from 5 counties in Yunnan Province. Ninety-three fecal samples (10.3%) were E. bieneusi-positive by PCR amplification. Four new genotypes and 11 known genotypes were identified, and all genotypes considered to be the zoonotic potential. Phylogenetic analysis showed that all of these genotypes were allocated into the zoonotic groups of E. bieneusi indicating its zoonotic potential. These results indicated that effective strategies and measures must be taken to prevent and control E. bieneusi transmission to other animals and humans.

Abstract

Enterocytozoon bieneusi is a fungus-like protist that can parasitize in the intestines of humans and various animals causing a threat to public health. However, there has been no data for E. bieneusi prevalence and genotypes in black goats in Yunnan Province, Southwestern China. In this study, 907 fecal samples were collected from black goats in 5 counties from Yunnan Province. The prevalence and genotypes of E. bieneusi were examined by nested PCR amplification targeting the nuclear internal transcribed spacer (ITS). Multi-locus sequence typing (MLST) was used to further examine the potential occurrence of genetic segregation. The overall prevalence of E. bieneusi in black goats in Yunnan Province was 10.3% (93/907). Statistical analysis revealed that E. bieneusi prevalence was significantly associated with the region, age and gender of black goats (p < 0.001). Four new genotypes (CYG-1, CYG-2, CYG-3, CYG-4) and 11 known genotypes (CHG1, CHG2, CHG3, CHG5, CHG28, J, D, BEB6, Wildboar3, CD6, SDD1) of E. bieneusi were identified. At the microsatellite and minisatellite loci, 15, 2, 17, and 33 sequences were obtained, respectively, forming one new multi-locus genotype (MLG27). Phylogenetic analysis revealed that all 15 genotypes were clustered into group 1 and group 2, with zoonotic potential. This is the first report of E. bieneusi prevalence and genotypes in black goats in Yunnan Province, China. Effective control strategies and measures should be taken to reduce the risk of E. bieneusi transmission between black goats, other animals, and humans.

Leishmania and the Model of Predominant Clonal Evolution

As it is the case for other pathogenic microorganisms, the respective impact of clonality and genetic exchange on Leishmania natural populations has been the object of lively debates since the early 1980s. The predominant clonal evolution (PCE) model states that genetic exchange in these parasites’ natural populations may have a high relevance on an evolutionary scale, but is not sufficient to erase a persistent phylogenetic signal and the existence of bifurcating trees. Recent data based on high-resolution markers and genomic polymorphisms fully confirm the PCE model down to a microevolutionary level.

The Population of Sclerotinia sclerotiorum in Brazil Is Structured by Mycelial Compatibility Groups

The genetic structure of the population of Sclerotinia sclerotiorum was analyzed using 238 individuals collected from different hosts. Individuals were characterized for microsatellite genotypes and mycelial compatibility groups (MCGs). A total of 22 MCGs and 64 multilocus lineages (MLLs) were identified. There was a close relationship between the MCGs and MLLs, but there was no association between MLLs and hosts or regions. At least 39 MCGs are present in Brazil, and 68.5% of the isolates were assigned to either MCG 1 or MCG 2. Eight new MCGs were found. Seven genetic groups were identified and associated with MCGs. Most genetic variation (70.0%) was because of differences among MCGs. High values of estimates of linkage disequilibrium among loci were more frequent in the total population (all MCGs). By contrast, there was evidence of random mating in subpopulations defined by MCGs 1 and 2. Additionally, there was evidence of outcrossing in the population of S. sclerotiorum in Brazil. The population was structured by MCGs; lineages originating from asexual reproduction or selfing prevail and are widely distributed in space, are persistent in time, and affect many hosts, but there is evidence of some degree of outcrossing, which may lead to a more genetically variable population in the future.

Trends in the molecular epidemiology and population genetics of emerging Sporothrix species

Sporothrix (Ophiostomatales) comprises species that are pathogenic to humans and other mammals as well as environmental fungi. Developments in molecular phylogeny have changed our perceptions about the epidemiology, host-association, and virulence of Sporothrix. The classical agent of sporotrichosis, Sporothrix schenckii, now comprises several species nested in a clinical clade with S. brasiliensis, S. globosa, and S. luriei. To gain a more precise view of outbreaks dynamics, structure, and origin of genetic variation within and among populations of Sporothrix, we applied three sets of discriminatory AFLP markers (#3 EcoRI-GA/MseI-TT, #5 EcoRI-GA/MseI-AG, and #6 EcoRI-TA/MseI-AA) and mating-type analysis to a large collection of human, animal and environmental isolates spanning the major endemic areas. A total of 451 polymorphic loci were amplified in vitro from 188 samples, and revealed high polymorphism information content (PIC = 0.1765-0.2253), marker index (MI = 0.0001-0.0002), effective multiplex ratio (E = 15.1720-23.5591), resolving power (Rp = 26.1075-40.2795), discriminating power (D = 0.9766-0.9879), expected heterozygosity (H = 0.1957-0.2588), and mean heterozygosity (Havp  = 0.000007-0.000009), demonstrating the effectiveness of AFLP markers to speciate Sporothrix. Analysis using the program structure indicated three genetic clusters matching S. brasiliensis (population 1), S. schenckii (population 2), and S. globosa (population 3), with the presence of patterns of admixture amongst all populations. AMOVA revealed highly structured clusters (PhiPT = 0.458-0.484, P < 0.0001), with roughly equivalent genetic variability within (46-48 %) and between (52-54 %) populations. Heterothallism was the exclusive mating strategy, and the distributions of MAT1-1 or MAT1-2 idiomorphs were not significantly skewed (1:1 ratio) for S. schenckii (χ2 = 2.522; P = 0.1122), supporting random mating. In contrast, skewed distributions were found for S. globosa (χ2 = 9.529; P = 0.0020) with a predominance of MAT1-1 isolates, and regional differences were highlighted for S. brasiliensis with the overwhelming occurrence of MAT1-2 in Rio de Janeiro (χ2 = 14.222; P = 0.0002) and Pernambuco (χ2 = 7.364; P = 0.0067), in comparison to a higher prevalence of MAT1-1 in the Rio Grande do Sul (χ2 = 7.364; P = 0.0067). Epidemiological trends reveal the geographic expansion of cat-transmitted sporotrichosis due to S. brasiliensis via founder effect. These data support Rio de Janeiro as the centre of origin that has led to the spread of this disease to other regions in Brazil. Our ability to reconstruct the source, spread, and evolution of the ongoing outbreaks from molecular data provides high-quality information for decision-making aimed at mitigating the progression of the disease. Other uses include surveillance, rapid diagnosis, case connectivity, and guiding access to appropriate antifungal treatment.

Towards an ecosystem model of infectious disease

Increasingly intimate associations between human society and the natural environment are driving the emergence of novel pathogens, with devastating consequences for humans and animals alike. Prior to emergence, these pathogens exist within complex ecological systems that are characterized by trophic interactions between parasites, their hosts and the environment. Predicting how disturbance to these ecological systems places people and animals at risk from emerging pathogens-and the best ways to manage this-remains a significant challenge. Predictive systems ecology models are powerful tools for the reconstruction of ecosystem function but have yet to be considered for modelling infectious disease. Part of this stems from a mistaken tendency to forget about the role that pathogens play in structuring the abundance and interactions of the free-living species favoured by systems ecologists. Here, we explore how developing and applying these more complete systems ecology models at a landscape scale would greatly enhance our understanding of the reciprocal interactions between parasites, pathogens and the environment, placing zoonoses in an ecological context, while identifying key variables and simplifying assumptions that underly pathogen host switching and animal-to-human spillover risk. As well as transforming our understanding of disease ecology, this would also allow us to better direct resources in preparation for future pandemics.

Reproduction in Trypanosomatids: Past and Present

Simple Summary

The reproduction of trypanosomatids is a fundamental issue for host–parasite interaction, and its biological importance lies in knowing how these species acquire new defense mechanisms against the countermeasures imposed by the host, which is consistent with the theory of the endless race or the Red Queen hypothesis for the existence of meiotic sex. Moreover, the way these species re-produce may also be at the origin of novel and more virulent clades and is relevant from a thera-peutic or vaccination point of view, as sex may contribute to increased tolerance and even to the rapid acquisition of drug resistance mechanisms. Kinetoplastids are single-celled organisms, many of them being responsible for important parasitic diseases, globally termed neglected diseases, which are endemic in low-income countries. Leishmaniasis, African (sleeping sickness) and American trypanosomiasis (Chagas disease) caused by trypanosomatids are among the most ne-glected tropical scourges related to poverty and poor health systems. The reproduction of these microorganisms has long been considered to be clonal due to population genetic observations. However, there is increasing evidence of true sex and genetic exchange events under laboratory conditions. We would like to highlight the importance of this topic in the field of host/parasite in-terplay, virulence, and drug resistance.

Abstract

Diseases caused by trypanosomatids (Sleeping sickness, Chagas disease, and leishmaniasis) are a serious public health concern in low-income endemic countries. These diseases are produced by single-celled parasites with a diploid genome (although aneuploidy is frequent) organized in pairs of non-condensable chromosomes. To explain the way they reproduce through the analysis of natural populations, the theory of strict clonal propagation of these microorganisms was taken as a rule at the beginning of the studies, since it partially justified their genomic stability. However, numerous experimental works provide evidence of sexual reproduction, thus explaining certain naturally occurring events that link the number of meiosis per mitosis and the frequency of mating. Recent techniques have demonstrated genetic exchange between individuals of the same species under laboratory conditions, as well as the expression of meiosis specific genes. The current debate focuses on the frequency of genomic recombination events and its impact on the natural parasite population structure. This paper reviews the results and techniques used to demonstrate the existence of sex in trypanosomatids, the inheritance of kinetoplast DNA (maxi- and minicircles), the impact of genetic exchange in these parasites, and how it can contribute to the phenotypic diversity of natural populations.

Evolutionary genomic and bacteria GWAS analysis of Mycobacterium avium subsp. paratuberculosis and dairy cattle Johne's disease phenotypes

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle. The Regional Dairy Quality Management Alliance (RDQMA) project is a multistate research program involving MAP isolates taken from three intensively studied commercial dairy farms in the northeastern United States, which emphasized longitudinal data collection of both MAP isolates and animal health in three regional dairy herds for a period of about 7 years. This paper reports the results of a pan-GWAS analysis involving 318 MAP isolates and dairy cow Johne's disease phenotypes, taken from these three farms. Based on our highly curated accessory gene count the pan-GWAS analysis identified several MAP genes associated with bovine Johne's disease phenotypes scored from these three farms, with some of the genes having functions suggestive of possible cause/effect relationships to these phenotypes. This paper reports a pan-genomic comparative analysis between MAP and Mycobacterium tuberculosis, assessing functional Gene Ontology category enrichments between these taxa. Finally, we also provide a population genomic perspective on the effectiveness of herd isolation, involving closed dairy farms, in preventing MAP inter-farm cross infection on a micro-geographic scale.IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle, and enormous economic consequences for the dairy industry. Understanding which genes in this bacterium are correlated with key disease phenotypes can lead to functional experiments targeting these genes and ultimately lead to improved control strategies. This study represents a rare example of a prolonged longitudinal study of dairy cattle where the disease was measured and the bacteria were isolated from the same cows. The genome sequences of over 300 MAP isolates were analyzed for genes that were correlated with a wide range of Johne's disease phenotypes. A number of genes were identified that were significantly associated with several aspects of the disease and suggestive of further experimental follow-up.

Hidden Diversity within Common Protozoan Parasites as Revealed by a Novel Genomotyping Scheme

Giardia duodenalis (syn. Giardia lamblia, Giardia intestinalis) is the causative agent of giardiasis, one of the most common diarrheal infections in humans. Evolutionary relationships among G. duodenalis genotypes (or subtypes) of assemblage B, one of two genetic assemblages causing the majority of human infections, remain unclear due to poor phylogenetic resolution of current typing methods. In this study, we devised a methodology to identify new markers for a streamlined multilocus sequence typing (MLST) scheme based on comparisons of all core genes against the phylogeny of whole-genome sequences (WGS). Our analysis identified three markers with resolution comparable to that of WGS data. Using newly designed PCR primers for our novel MLST loci, we typed an additional 68 strains of assemblage B. Analyses of these strains and previously determined genome sequences showed that genomes of this assemblage can be assigned to 16 clonal complexes, each with unique gene content that is apparently tuned to differential virulence and ecology. Obtaining new genomes of Giardia spp. and other eukaryotic microbial pathogens remains challenging due to difficulties in culturing the parasites in the laboratory. Hence, the methods described here are expected to be widely applicable to other pathogens of interest and advance our understanding of their ecology and evolution.IMPORTANCE Giardia duodenalis assemblage B is a major waterborne pathogen and the most commonly identified genotype causing human giardiasis worldwide. The lack of morphological characters for classification requires the use of molecular techniques for strain differentiation; however, the absence of scalable and affordable next-generation sequencing (NGS)-based typing methods has prevented meaningful advancements in high-resolution molecular typing for further understanding of the evolution and epidemiology of assemblage B. Prior studies have reported high sequence diversity but low phylogenetic resolution at standard loci in assemblage B, highlighting the necessity of identifying new markers for accurate and robust molecular typing. Data from comparative analyses of available genomes in this study identified three loci that together form a novel high-resolution typing scheme with high concordance to whole-genome-based phylogenomics and which should aid in future public health endeavors related to this parasite. In addition, data from newly characterized strains suggest evidence of biogeographic and ecologic endemism.

The discernible and hidden effects of clonality on the genotypic and genetic states of populations: Improving our estimation of clonal rates

Partial clonality is widespread across the tree of life, but most population genetic models are designed for exclusively clonal or sexual organisms. This gap hampers our understanding of the influence of clonality on evolutionary trajectories and the interpretation of population genetic data. We performed forward simulations of diploid populations at increasing rates of clonality (c), analysed their relationships with genotypic (clonal richness, R, and distribution of clonal sizes, Pareto β) and genetic (F_IS and linkage disequilibrium) indices, and tested predictions of c from population genetic data through supervised machine learning. Two complementary behaviours emerged from the probability distributions of genotypic and genetic indices with increasing c. While the impact of c on R and Pareto β was easily described by simple mathematical equations, its effects on genetic indices were noticeable only at the highest levels (c > 0.95). Consequently, genotypic indices allowed reliable estimates of c, while genetic descriptors led to poorer performances when c < 0.95. These results provide clear baseline expectations for genotypic and genetic diversity and dynamics under partial clonality. Worryingly, however, the use of realistic sample sizes to acquire empirical data systematically led to gross underestimates (often of one to two orders of magnitude) of c, suggesting that many interpretations hitherto proposed in the literature, mostly based on genotypic richness, should be reappraised. We propose future avenues to derive realistic confidence intervals for c and show that, although still approximate, a supervised learning method would greatly improve the estimation of c from population genetic data.

Models in parasite and pathogen evolution: Genomic analysis reveals predominant clonality and progressive evolution at all evolutionary scales in parasitic protozoa, yeasts and bacteria

The predominant clonal evolution (PCE) model of pathogenic microorganisms postulates that the impact of genetic recombination in those pathogens' natural populations is not enough to erase a persistent phylogenetic signal at all evolutionary scales from microevolution till geological times in the whole ecogeographical range of the species considered. We have tested this model with a set of representative parasitic protozoa, yeasts and bacteria in the light of the most recent genomic data. All surveyed species, including those that were considered as highly recombining, exhibit similar PCE patterns above and under the species level, from macro- to micro-evolutionary scales (Russian doll pattern), suggesting gradual evolution. To our knowledge, it is the first time that such a strong common evolutionary feature among very diverse pathogens has been evidenced. The implications of this model for basic biology and applied research are exposed. These implications include our knowledge on the pathogens' reproductive mode, their population structure, the possibility to type strain and to follow up epidemics (molecular epidemiology) and to revisit pathogens' taxonomy through a flexible use of the phylogenetic species concept (Cracraft, 1983).

Trypanosoma Cruzi Genome: Organization, Multi-Gene Families, Transcription, and Biological Implications

Chagas disease caused by the parasite Trypanosoma cruzi affects millions of people. Although its first genome dates from 2005, its complexity hindered a complete assembly and annotation. However, the new sequencing methods have improved genome annotation of some strains elucidating the broad genetic diversity and complexity of this parasite. Here, we reviewed the genomic structure and regulation, the genetic diversity, and the analysis of the principal multi-gene families of the recent genomes for several strains. The telomeric and sub-telomeric regions are sites with high recombination events, the genome displays two different compartments, the core and the disruptive, and the genome plasticity seems to play a key role in the survival and the infection process. Trypanosoma cruzi (T. cruzi) genome is composed mainly of multi-gene families as the trans-sialidases, mucins, and mucin-associated surface proteins. Trans-sialidases are the most abundant genes in the genome and show an important role in the effectiveness of the infection and the parasite survival. Mucins and MASPs are also important glycosylated proteins of the surface of the parasite that play a major biological role in both insect and mammal-dwelling stages. Altogether, these studies confirm the complexity of T. cruzi genome revealing relevant concepts to better understand Chagas disease.

Staphylococcus aureus Host Tropism and Its Implications for Murine Infection Models

Staphylococcus aureus (S. aureus) is a pathobiont of humans as well as a multitude of animal species. The high prevalence of multi-resistant and more virulent strains of S. aureus necessitates the development of new prevention and treatment strategies for S. aureus infection. Major advances towards understanding the pathogenesis of S. aureus diseases have been made using conventional mouse models, i.e., by infecting naïve laboratory mice with human-adapted S.aureus strains. However, the failure to transfer certain results obtained in these murine systems to humans highlights the limitations of such models. Indeed, numerous S. aureus vaccine candidates showed promising results in conventional mouse models but failed to offer protection in human clinical trials. These limitations arise not only from the widely discussed physiological differences between mice and humans, but also from the lack of attention that is paid to the specific interactions of S. aureus with its respective host. For instance, animal-derived S. aureus lineages show a high degree of host tropism and carry a repertoire of host-specific virulence and immune evasion factors. Mouse-adapted S.aureus strains, humanized mice, and microbiome-optimized mice are promising approaches to overcome these limitations and could improve transferability of animal experiments to human trials in the future.

Genetic analysis reveals long-standing population differentiation and high diversity in the rust pathogen Melampsora lini

A priority for research on infectious disease is to understand how epidemiological and evolutionary processes interact to influence pathogen population dynamics and disease outcomes. However, little is understood about how population adaptation changes across time, how sexual vs. asexual reproduction contribute to the spread of pathogens in wild populations and how diversity measured with neutral and selectively important markers correlates across years. Here, we report results from a long-term study of epidemiological and genetic dynamics within several natural populations of the Linum marginale-Melampsora lini plant-pathogen interaction. Using pathogen isolates collected from three populations of wild flax (L. marginale) spanning 16 annual epidemics, we probe links between pathogen population dynamics, phenotypic variation for infectivity and genomic polymorphism. Pathogen genotyping was performed using 1567 genome-wide SNP loci and sequence data from two infectivity loci (AvrP123, AvrP4). Pathogen isolates were phenotyped for infectivity using a differential set. Patterns of epidemic development were assessed by conducting surveys of infection prevalence in one population (Kiandra) annually. Bayesian clustering analyses revealed host population and ecotype as key predictors of pathogen genetic structure. Despite strong fluctuations in pathogen population size and severe annual bottlenecks, analysis of molecular variance revealed that pathogen population differentiation was relatively stable over time. Annually, varying levels of clonal spread (0–44.8%) contributed to epidemics. However, within populations, temporal genetic composition was dynamic with rapid turnover of pathogen genotypes, despite the dominance of only four infectivity phenotypes across the entire study period. Furthermore, in the presence of strong fluctuations in population size and migration, spatial selection may maintain pathogen populations that, despite being phenotypically stable, are genetically highly dynamic.

Melampsora lini is a rust fungus that infects native flax, Linum marginale in south-eastern Australia where its epidemiology and evolution have been intensively studied since 1987. Over that time, substantial diversity in the pathotypic structure of M. lini has been demonstrated but an understanding of how genetic diversity in pathogen populations is maintained through space and time is lacking. Here we integrated phenotypic, genotypic and epidemiological datasets spanning 16 annual epidemics across three host populations to examine long-term pathogen genetic dynamics. The results show that host ecotype is the dominant selective force in the face of strong bottlenecks and annual patterns of genetic turnover. Results from previous studies indicate that in this geographic region, M. lini lacks the capacity to reproduce sexually–we thus expected to find limited genetic diversity and evidence for strong clonality influencing genetic dynamics within growing seasons. However, the breadth of genomic coverage provided by the SNP markers revealed high levels of genotypic variation within M. lini populations. This discovery contrasts with observed phenotypic dynamics as the epidemics of this pathogen were largely dominated by four pathotypes across the study period. Based on a detailed assessment and comparison of pathotypic and genotypic patterns, our study increases the understanding of how genetic diversity is generated and maintained through space and time within wild pathogen populations. The implications for the management of resistance to pathogens in agricultural or conservation contexts are significant: the appearance of clonality may be hiding high levels of pathogen diversity and recombination. Understanding how this diversity is generated could provide new and unique ways to mitigate or suppress the emergence of infectious strains, allowing to efficiently combat harmful diseases.

Founder events influence structures of Aspergillus flavus populations

In warm regions, agricultural fields are occupied by complex Aspergillus flavus communities composed of isolates in many vegetative compatibility groups (VCGs) with varying abilities to produce highly toxic, carcinogenic aflatoxins. Aflatoxin contamination is reduced with biocontrol products that enable atoxigenic isolates from atoxigenic VCGs to dominate the population. Shifts in VCG frequencies similar to those caused by the introduction of biocontrol isolates were detected in Sonora, Mexico, where biocontrol is not currently practiced. The shifts were attributed to founder events. Although VCGs reproduce clonally, significant diversity exists within VCGs. Simple sequence repeat (SSR) fingerprinting revealed that increased frequencies of VCG YV150 involved a single haplotype. This is consistent with a founder event. Additionally, great diversity was detected among 82 YV150 isolates collected over 20 years across Mexico and the United States. Thirty-six YV150 haplotypes were separated into two populations by Structure and SplitsTree analyses. Sixty-five percent of isolates had MAT1-1 and belonged to one population. The remaining had MAT1-2 and belonged to the second population. SSR alleles varied within populations, but recombination between populations was not detected despite co-occurrence at some locations. Results suggest that YV150 isolates with opposite mating-type have either strongly restrained or lost sexual reproduction among themselves.

Founder events influence structures of Aspergillus flavus populations

In warm regions, agricultural fields are occupied by complex Aspergillus flavus communities composed of isolates in many vegetative compatibility groups (VCGs) with varying abilities to produce highly toxic, carcinogenic aflatoxins. Aflatoxin contamination is reduced with biocontrol products that enable atoxigenic isolates from atoxigenic VCGs to dominate the population. Shifts in VCG frequencies similar to those caused by the introduction of biocontrol isolates were detected in Sonora, Mexico, where biocontrol is not currently practiced. The shifts were attributed to founder events. Although VCGs reproduce clonally, significant diversity exists within VCGs. Simple sequence repeat (SSR) fingerprinting revealed that increased frequencies of VCG YV150 involved a single haplotype. This is consistent with a founder event. Additionally, great diversity was detected among 82 YV150 isolates collected over 20 years across Mexico and the United States. Thirty-six YV150 haplotypes were separated into two populations by Structure and SplitsTree analyses. Sixty-five percent of isolates had MAT1-1 and belonged to one population. The remaining had MAT1-2 and belonged to the second population. SSR alleles varied within populations, but recombination between populations was not detected despite co-occurrence at some locations. Results suggest that YV150 isolates with opposite mating-type have either strongly restrained or lost sexual reproduction among themselves.

Comparative genomics and pangenome-oriented studies reveal high homogeneity of the agronomically relevant enterobacterial plant pathogen Dickeya solani

Background

Dickeya solani is an important plant pathogenic bacterium causing severe losses in European potato production. This species draws a lot of attention due to its remarkable virulence, great devastating potential and easier spread in contrast to other Dickeya spp. In view of a high need for extensive studies on economically important soft rot Pectobacteriaceae, we performed a comparative genomics analysis on D. solani strains to search for genetic foundations that would explain the differences in the observed virulence levels within the D. solani population.

Results

High quality assemblies of 8 de novo sequenced D. solani genomes have been obtained. Whole-sequence comparison, ANIb, ANIm, Tetra and pangenome-oriented analyses performed on these genomes and the sequences of 14 additional strains revealed an exceptionally high level of homogeneity among the studied genetic material of D. solani strains. With the use of 22 genomes, the pangenome of D. solani, comprising 84.7% core, 7.2% accessory and 8.1% unique genes, has been almost completely determined, suggesting the presence of a nearly closed pangenome structure. Attribution of the genes included in the D. solani pangenome fractions to functional COG categories showed that higher percentages of accessory and unique pangenome parts in contrast to the core section are encountered in phage/mobile elements- and transcription- associated groups with the genome of RNS 05.1.2A strain having the most significant impact. Also, the first D. solani large-scale genome-wide phylogeny computed on concatenated core gene alignments is herein reported.

Conclusions

The almost closed status of D. solani pangenome achieved in this work points to the fact that the unique gene pool of this species should no longer expand. Such a feature is characteristic of taxa whose representatives either occupy isolated ecological niches or lack efficient mechanisms for gene exchange and recombination, which seems rational concerning a strictly pathogenic species with clonal population structure. Finally, no obvious correlations between the geographical origin of D. solani strains and their phylogeny were found, which might reflect the specificity of the international seed potato market.

Low lineage diversity and increased virulence of group C Streptococcus dysgalactiae subsp. equisimilis

Introduction. In some species, the population structure of pathogenic bacteria is clonal. However, the mechanisms that determine the predominance and persistence of specific bacterial lineages of group C Streptococcus remain poorly understood. In Brazil, a previous study revealed the predominance of two main lineages of Streptococcus dysgalactiae subsp. equisimilis (SDSE).Aim. The aim of this study was to assess the virulence and fitness advantages that might explain the predominance of these SDSE lineages for a long period of time.Methodology. emm typing was determined by DNA sequencing. Adhesion and invasion tests were performed using human bronchial epithelial cells (16HBE14o-). Biofilm formation was tested on glass surfaces and the presence of virulence genes was assessed by PCR. Additionally, virulence was studied using Caenorhabditis elegans models and competitive fitness was analysed in murine models.Results. The predominant lineages A and B were mostly typed as emm stC839 and stC6979, respectively. Notably, these lineages exhibited a superior ability to adhere and invade airway cells. Furthermore, the dominant lineages were more prone to induce aversive olfactory learning and more likely to kill C. elegans. In the competitive fitness assays, they also showed increased adaptability. Consistent with the increased virulence observed in the ex vivo and in vivo models, the predominant lineages A and B showed a higher number of virulence-associated genes and a superior ability to accumulate biofilm.Conclusion. These results suggest strongly that this predominance did not occur randomly but rather was due to adaptive mechanisms that culminated in increased colonization and other bacterial properties that might confer increased bacteria-host adaptability to cause disease.

Genomics and High-Resolution Typing Confirm Predominant Clonal Evolution Down to a Microevolutionary Scale in Trypanosoma cruzi

Trypanosoma cruzi, the agent of Chagas disease, is a paradigmatic case of the predominant clonal evolution (PCE) model, which states that the impact of genetic recombination in pathogens' natural populations is not sufficient to suppress a persistent phylogenetic signal at all evolutionary scales. In spite of indications for occasional recombination and meiosis, recent genomics and high-resolution typing data in T. cruzi reject the counterproposal that PCE does not operate at lower evolutionary scales, within the evolutionary units (=near-clades) that subdivide the species. Evolutionary patterns in the agent of Chagas disease at micro- and macroevolutionary scales are strikingly similar ("Russian doll pattern"), suggesting gradual, rather than saltatory evolution.

Gene duplications, divergence and recombination shape adaptive evolution of the fish ectoparasite Gyrodactylus bullatarudis

Determining the molecular basis of parasite adaptation to its host is an important component in understanding host-parasite coevolution and the epidemiology of parasitic infections. Here, we investigate short- and long-term adaptive evolution in the eukaryotic parasite Gyrodactylus bullatarudis infecting Caribbean guppies (Poecilia reticulata), by comparing the reference genome of Tobagonian G. bullatarudis with other Platyhelminthes, and by analysing resequenced samples from local Trinidadian populations. At the macroevolutionary timescale, we observed duplication of G-protein and serine proteases genes, which are probably important in host-parasite arms races. Serine protease also showed strong evidence of ongoing, diversifying selection at the microevolutionary timescale. Furthermore, our analyses revealed that a hybridization event, involving two divergent genomes, followed by recombination has dramatically affected the genetic composition of Trinidadian populations. The recombinant genotypes invaded Trinidad and replaced local parasites in all populations. We localized more than 300 genes in regions fixed in local populations for variants of different origin, possibly due to diversifying selection pressure from local host populations. In addition, around 70 genes were localized in regions identified as heterozygous in some, but not all, individuals. This pattern is consistent with a very recent spread of recombinant parasites. Overall, our results are consistent with the idea that recombination between divergent genomes can result in particularly successful parasites.

Evidence of hybridization, mitochondrial introgression and biparental inheritance of the kDNA minicircles in Trypanosoma cruzi I

Background

Genetic exchange in Trypanosoma cruzi is controversial not only in relation to its frequency, but also to its mechanism. Parasexual genetic exchange has been proposed based on laboratory hybrids, but population genomics strongly suggests meiosis in T. cruzi. In addition, mitochondrial introgression has been reported several times in natural isolates although its mechanism is not fully understood yet. Moreover, hybrid T. cruzi DTUs (TcV and TcVI) have inherited at least part of the kinetoplastic DNA (kDNA = mitochondrial DNA) from both parents.

Methodology/Principal findings

In order to address such topics, we sequenced and analyzed fourteen nuclear DNA fragments and three kDNA maxicircle genes in three TcI stocks which are natural clones potentially involved in events of genetic exchange. We also deep-sequenced (a total of 6,146,686 paired-end reads) the minicircle hypervariable region (mHVR) of the kDNA in such three strains. In addition, we analyzed the DNA content by flow cytometry to address cell ploidy. We observed that most polymorphic sites in nuclear loci showed a hybrid pattern in one cloned strain and the other two cloned strains were compatible as parental strains (or nearly related to the true parents). The three clones had almost the same ploidy and the DNA content was similar to the reference strain Sylvio (a nearly diploid strain). Despite maxicircle genes evolve faster than nuclear housekeeping ones, we detected no polymorphisms in the sequence of three maxicircle genes showing mito-nuclear discordance. Lastly, the hybrid stock shared 66% of its mHVR clusters with one putative parent and 47% with the other one; in contrast, the putative parental stocks shared less than 30% of the mHVR clusters between them.

Conclusions/significance

The results suggest a reductive division, a natural hybridization, biparental inheritance of the minicircles in the hybrid and maxicircle introgression. The models including such phenomena and explaining the relationships between these three clones are discussed.

Chagas disease, an important public health problem in Latin America, is caused by the parasite Trypanosoma cruzi. Despite being a widely studied parasite, several questions on the biology of genetic exchange remain unanswered. Population genomic studies have inferred meiosis in T. cruzi, but this cellular division mechanism has not been observed in laboratory yet. In addition, previous results suggest that mitochondrial DNA (called kDNA) may be inherited from both parents in hybrids. Here, we analyzed a hybrid strain and its potential parents to address the mechanisms of genetic exchange at nuclear and mitochondrial levels. We observed that the hybrid strain had heterozygous patterns and DNA content compatible with a meiosis event. Also, we observed that the evolutionary histories of nuclear DNA and kDNA maxicircles were discordant and that the three strains shared identical DNA sequences. Mitochondrial introgression of maxicircle DNA from one genotype to another may explain this observation. In addition, we demonstrated that the hybrid strain shared kDNA minicircles with both parental strains. Our results suggest that hybridization implied meiosis and biparental inheritance of the kDNA. Further research is required to address such phenomena in detail.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

Longitudinal dynamics of herd-level Mycobacterium bovis MLVA type surveillance in cattle in Northern Ireland 2003-2016

Investigating genetically-structured diversity in pathogen populations over time is important to better understand disease maintenance and spread. Herd-level surveillance of Mycobacterium bovis genotypes (multi-locus VNTR analysis types, MLVA types) from all culture-confirmed bovine tuberculosis (TB) herd cases was undertaken in Northern Ireland (NI), generating an unparalleled, longitudinal, population-level 14-year survey for this pathogen. Across this population, 295 genetically-distinct M. bovis MLVA types were identified in the 19,717 M. bovis isolates surveyed. Of these, the most frequent was MLVA type 002 (23.0%); 151 MLVA types were represented more than once, in groups ranging from 2 to 4438 isolates. Only 23 MLVA types were isolated in all 14 years. Investigating inter-annual frequency of M. bovis MLVA types, examples of statistically-significant expansions (MLVA types 002, 004, 006, 009 and 027), contractions (MLVA types 001, 007 and 011) and maintenance (MLVA types 003 and 005) were disclosed, during a period of fluctuating bovine TB herd-level incidence at the NI scale. The fixed period frequency distribution of MLVA types remained highly right-skewed. Novel VNTR copy number variant MLVA types (N = 242; an average of 17 per annum) were identified throughout the survey. The MLVA type distribution in the landscape was not random; MLVA types showed statistically-significant geographical localization and strong spatial associations with Divisional Veterinary Office (DVO) regions. There was also evidence of differential risk of particular MLVA types across breeds (Holstein/Friesian vs. other), age-class, and sex and some evidence of an association between the number of animals testing positive for bovine TB during the disclosing test and particular MLVA types, although there was substantial variation.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

Fifty Aureobasidium pullulans genomes reveal a recombining polyextremotolerant generalist

The black yeast Aureobasidium pullulans is a textbook example of a generalistic and ubiquitous fungus thriving in a wide variety of environments. To investigate whether A. pullulans is a true generalist, or alternatively, whether part of its versatility can be attributed to intraspecific specialization masked by cryptic diversification undetectable by traditional phylogenetic analyses, we sequenced and analysed the genomes of 50 strains of A. pullulans from different habitats and geographic locations. No population structure was observed in the sequenced strains. Decay of linkage disequilibrium over shorter physical distances (<100 bp) than in many sexually reproducing fungi indicates a high level of recombination in the species. A homothallic mating locus was found in all of the sequenced genomes. Aureobasidium pullulans appears to have a homogeneous population genetics structure, which is best explained by good dispersal and high levels of recombination. This means that A. pullulans is a true generalist that can inhabit different habitats without substantial specialization to any of these habitats at the genomic level. Furthermore, in the future, the high level of A. pullulans recombination can be exploited for the identification of genomic loci that are involved in the many biotechnologically useful traits of this black yeast.

Population Genomics of an Obligately Halophilic Basidiomycete Wallemia ichthyophaga

Background

Wallemia ichthyophaga is a highly specialized basidiomycetous fungus. It is one of the most halophilic fungi ever described, only able to grow at low water activity. This specialization is thought to explain why it is only rarely isolated from nature.

Results

Genomes of 21 W. ichthyophaga strains were sequenced with PE150 reads on BGISEQ500 platform. The genomes shared high similarity with the reference genome of the species, they were all smaller than 10 Mbp and had a low number of predicted genes. Groups of strains isolated in the same location encompassed clones as well as very divergent strains. There was little concordance between phylogenies of predicted genes. Linkage disequilibrium of pairs of polymorphic loci decayed relatively quickly as a function of distance between the loci (LD decay distance 1270 bp). For the first time a putative mating-type locus was identified in the genomes of W. ichthyophaga.

Conclusion

Based on the comparison of W. ichthyophaga genomes it appears that some phylogenetic lineages of the species can persist in the same location over at least several years. Apart from this, the differences between the strains do not reflect the isolation habitat or geographic location. Together with results supporting the existence of (sexual) recombination in W. ichthyophaga, the presented results indicate that strains of W. ichthyophaga can form a single recombining population even between different habitats and over large geographical distances.

Deciphering how plant pathogenic bacteria disperse and meet: Molecular epidemiology of Xanthomonas citri pv. citri at microgeographic scales in a tropical area of Asiatic citrus canker endemicity

Although some plant pathogenic bacteria represent a significant threat to agriculture, the determinants of their ecological success and evolutionary potential are still poorly understood. Refining our understanding of bacterial strain circulation at small spatial scales and the biological significance and evolutionary consequences of co-infections are key questions. The study of bacterial population biology can be challenging, because it requires high-resolution markers that can be genotyped with a high throughput. Here, we overcame this difficulty for Xanthomonas citri pv. citri, a genetically monomorphic bacterium causing Asiatic citrus canker (ACC). Using a genotyping method that did not require cultivating the bacterium or purifying DNA, we deciphered the pathogen's spatial genetic structure at several microgeographic scales, down to single lesion, in a situation of ACC endemicity. In a grove where copper was recurrently applied for ACC management, copper-susceptible and copper-resistant X. citri pv. citri coexisted and the bacterial population structured as three genetic clusters, suggesting a polyclonal contamination. The range of spatial dependency, estimated for the two largest clusters, was four times greater for the cluster predominantly composed of copper-resistant bacteria. Consistently, the evenness value calculated for this cluster was indicative of increased transmission. Linkage disequilibrium was high even at a tree scale, probably due to a combination of clonality and admixture. Approximately 1% of samples exhibited within-lesion multilocus polymorphism, explained at least in part by polyclonal infections. Canker lesions, which are of major biological significance as an inoculum source, may also represent a preferred niche for horizontal gene transfer. This study points out the potential of genotyping data for estimating the range of spatial dependency of plant bacterial pathogens, an important parameter for guiding disease management strategies.

Novel clinical and dual infection by Histoplasma capsulatum genotypes in HIV patients from Northeastern, Brazil

Histoplasmosis is a worldwide-distributed deep mycosis that affects healthy and immunocompromised hosts. Severe and disseminated disease is especially common in HIV-infected patients. At least 11 phylogenetic species are recognized and the majority of diversity is found in Latin America. The northeastern region of Brazil has one of the highest HIV/AIDS prevalence in Latin America and Ceará State has one of the highest death rates due to histoplasmosis in the world, where the mortality rate varies between 33-42%. The phylogenetic distribution and population genetic structure of 51 clinical isolates from Northeast Brazil was studied. For that morphological characteristics, exoantigens profile, and fungal mating types were evaluated. The genotypes were deduced by a MSLT in order to define local population structure of this fungal pathogen. In addition, the relationships of H. capsulatum genotypes with clinically relevant phenotypes and clinical aspects were investigated. The results suggest two cryptic species, herein named population Northeast BR1 and population Northeast BR2. These populations are recombining, exhibit a high level of haplotype diversity, and contain different ratios of mating types MAT1-1 and MAT1-2. However, differences in phenotypes or clinical aspects were not observed within these new cryptic species. A HIV patient can be co-infected by two or more genotypes from Northeast BR1 and/or Northeast BR2, which may have significant impact on disease progression due to the impaired immune response. We hypothesize that co-infections could be the result of multiple exposure events and may indicate higher risk of disseminated histoplasmosis, especially in HIV infected patients.

Genotyping genetically heterogeneous Cyclospora cayetanensis infections to complement epidemiological case linkage

Sexually reproducing pathogens such as Cyclospora cayetanensis often produce genetically heterogeneous infections where the number of unique sequence types detected at any given locus varies depending on which locus is sequenced. The genotypes assigned to these infections quickly become complex when additional loci are analysed. This genetic heterogeneity confounds the utility of traditional sequence-typing and phylogenetic approaches for aiding epidemiological trace-back, and requires new methods to address this complexity. Here, we describe an ensemble of two similarity-based classification algorithms, including a Bayesian and heuristic component that infer the relatedness of C. cayetanensis infections. The ensemble requires a set of haplotypes as input and assigns arbitrary distances to specimen pairs reflecting their most likely relationships. The approach was applied to data generated from a test cohort of 88 human fecal specimens containing C. cayetanensis, including 30 from patients whose infections were associated with epidemiologically defined outbreak clusters of cyclosporiasis. The ensemble assigned specimens to plausible clusters of genetically related infections despite their complex haplotype composition. These relationships were corroborated by a significant number of epidemiological linkages (P < 0.0001) suggesting the ensemble's utility for aiding epidemiological trace-back investigations of cyclosporiasis.

Genomic Evidence of Recombination in the Basidiomycete Wallemia mellicola

One of the most commonly encountered species in the small basidiomycetous sub-phylum Wallemiomycotina is Wallemia mellicola, a xerotolerant fungus with a widespread distribution. To investigate the population characteristics of the species, whole genomes of twenty-five strains were sequenced. Apart from identification of four strains of clonal origin, the distances between the genomes failed to reflect either the isolation habitat of the strains or their geographical origin. Strains from different parts of the world appeared to represent a relatively homogenous and widespread population. The lack of concordance between individual gene phylogenies and the decay of linkage disequilibrium indicated that W. mellicola is at least occasionally recombining. Two versions of a putative mating-type locus have been found in all sequenced genomes, each present in approximately half of the strains. W. mellicola thus appears to be capable of (sexual) recombination and shows no signs of allopatric speciation or specialization to specific habitats.

Sexual reproduction in a natural Trypanosoma cruzi population

Background

Sexual reproduction provides an evolutionary advantageous mechanism that combines favorable mutations that have arisen in separate lineages into the same individual. This advantage is especially pronounced in microparasites as allelic reassortment among individuals caused by sexual reproduction promotes allelic diversity at immune evasion genes within individuals which is often essential to evade host immune systems. Despite these advantages, many eukaryotic microparasites exhibit highly-clonal population structures suggesting that genetic exchange through sexual reproduction is rare. Evidence supporting clonality is particularly convincing in the causative agent of Chagas disease, Trypanosoma cruzi, despite equally convincing evidence of the capacity to engage in sexual reproduction.

Methodology/ Principle Findings

In the present study, we investigated two hypotheses that can reconcile the apparent contradiction between the observed clonal population structure and the capacity to engage in sexual reproduction by analyzing the genome sequences of 123 T. cruzi isolates from a natural population in Arequipa, Peru. The distribution of polymorphic markers within and among isolates provides clear evidence of the occurrence of sexual reproduction. Large genetic segments are rearranged among chromosomes due to crossing over during meiosis leading to a decay in the genetic linkage among polymorphic markers compared to the expectations from a purely asexually-reproducing population. Nevertheless, the population structure appears clonal due to a high level of inbreeding during sexual reproduction which increases homozygosity, and thus reduces diversity, within each inbreeding lineage.

Conclusions/ Significance

These results effectively reconcile the apparent contradiction by demonstrating that the clonal population structure is derived not from infrequent sex in natural populations but from high levels of inbreeding. We discuss epidemiological consequences of this reproductive strategy on genome evolution, population structure, and phenotypic diversity of this medically important parasite.

The rearrangement of alleles among individuals in a population during sexual reproduction maintains high allelic diversity within individuals in a population at polymorphic genes. Allelic diversity within individuals can be particularly important for parasites as it enhances their ability to evade host immune systems. Despite the potential benefits of sexual reproduction for parasites, natural populations of the protozoan parasite—and causative agent of human Chagas disease—Trypanosoma cruzi, exhibit clonal population structures indicative of asexual reproduction. This is particularly surprising as T. cruzi has the capacity for sexual reproduction. Here, we resolve this apparent contradiction by sequencing whole genomes of 123 T. cruzi isolates from a natural population in Arequipa, Peru. Evidence of past sexual reproduction and allelic rearrangements are common in this T. cruzi population. However, the majority of sexual reproduction events occur between close relatives resulting in an apparent clonal population structure. Sexual reproduction with distant relatives in areas with greater strain diversity has the potential to affect public health by increasing diversity in immune evasion genes within individuals and enhancing within-host survival, rapidly diversifying antigens that could affect the sensitivity of serological diagnostics, and by generating diversity in pathogenicity or drug resistance.

Host Specificity of Enterocytozoon bieneusi and Public Health Implications

Enterocytozoon bieneusi is the most common cause of human microsporidiosis and it also infects a wide range of mammals and birds worldwide. The role of animals in the transmission of this parasite to humans and its public health importance remain poorly elucidated. This review summarizes all E. bieneusi genotypes identified thus far based on sequence analysis of the ribosomal internal transcribed spacer (ITS) from specimens obtained from humans, domestic and wild animals, and water sources; it examines genotypes, host and geographical distribution, analyzes inter- and intragenotype group host specificity, and interprets the public health significance of genotype groups and major zoonotic genotypes, with the goal of improving our understanding of host specificity in E. bieneusi and its implications for interspecies and zoonotic transmission.

Multilocus Sequence Typing and Population Genetic Analysis of Enterocytozoon bieneusi: Host Specificity and Its Impacts on Public Health

Microsporidia comprise a large class of unicellular eukaryotic pathogens that are medically and agriculturally important, but poorly understood. There have been nearly 1,500 microsporidian species described thus far, which are variable in biology, genetics, genomics, and host specificity. Among those, Enterocytozoon bieneusi is the well-known species responsible for the most recorded cases of human microsporidian affections. The pathogen can colonize a broad range of mammals and birds and most of the animals surveyed share some genotypes with humans, posing a threat to public health. Based on DNA sequence analysis of the ribosomal internal transcribed spacer (ITS) and phylogenetic analysis, several hundreds of E. bieneusi genotypes have been defined and clustered into different genetic groups with varied levels of host specificity. However, single locus-based typing using ITS might have insufficient resolution to discriminate among E. bieneusi isolates with complex genetic or hereditary characteristics and to assess the elusive reproduction or transmission modes of the organism, highlighting the need for exploration and application of multilocus sequence typing (MLST) and population genetic tools. The present review begins with a primer on microsporidia and major microsporidian species, briefly introduces the recent advances on E. bieneusi ITS genotyping and phylogeny, summarizes recent MLST and population genetic data, analyzes the inter- and intragroup host specificity at the MLST level, and interprets the public health implications of host specificity in zoonotic or cross-species transmission of this ubiquitous fungus.

Novel insights on colonization routes and evolutionary potential of Colletotrichum kahawae, a severe pathogen of Coffea arabica

Pathogenic fungi are emerging at an increasing rate on a wide range of host plants, leading to tremendous threats to the global economy and food safety. Several plant pathogens have been considered to be invasive species, rendering large-scale population genomic analyses crucial to better understand their demographic history and evolutionary potential. Colletotrichum kahawae (Ck) is a highly aggressive and specialized pathogen, causing coffee berry disease in Arabica coffee in Africa. This pathogen leads to severe production losses and its dissemination out of Africa is greatly feared. To address this issue, a population genomic approach using thousands of single nucleotide polymorphisms (SNPs) spaced throughout the genome was used to unveil its demographic history and evolutionary potential. The current study confirms that Ck is a true clonal pathogen, perfectly adapted to green coffee berries, with three completely differentiated populations (Angolan, Cameroonian and East African). Two independent clonal lineages were found within the Angolan population as opposed to the remaining single clonal populations. The most probable colonization scenario suggests that this pathogen emerged in Angola and immediately dispersed to East Africa, where these two populations began to differentiate, followed by the introduction in Cameroon from an Angolan population. However, the differentiation between the two Angolan clonal lineages masks the mechanism for the emergence of the Cameroonian population. Our results suggest that Ck is completely differentiated from the ancestral lineage, has a low evolutionary potential and a low dispersion ability, with human transport the most likely scenario for its potential dispersion, which makes the fulfilment of the quarantine measures and management practices implemented crucial.