Population Divergence Research Articles

Delimiting species, revealing cryptic diversity, and population divergence in Qinghai-Tibet Plateau weevils through DNA barcoding.

The Leptomias group represents one of the most diverse taxonomic group of weevils in the Qinghai-Tibet Plateau and its adjacent areas. Despite the potential of hidden diversity, relatively few comprehensive studies have been conducted on species diversity in this taxonomic group. In this study, we performed DNA barcoding analysis for species of the Leptomias group using a comprehensive DNA barcode dataset that included 476 sequences representing 54 morphospecies. Within the dataset, our laboratory contributed 474 sequences, and 390 sequences were newly generated for this study. The average Kimura 2-parameter distances among morphospecies and genera were 0.76% and 19.15%, respectively. In 94.4% of the species, the minimum interspecific distances exceeded the maximum intraspecific distances, indicating the presence of barcode gaps in most species of Leptomias group. The application of Automatic Barcode Gap Discovery, Assemble Species by Automatic Partitioning, Barcode Index Number, Bayesian Poisson tree processes, jMOTU, and Neighbor-joining tree methods revealed 45, 45, 63, 54, and 55 distinct clusters representing single species, respectively. Additionally, a total of four morphospecies, Leptomias kangmarensis, L. midlineatus, L. siahus, and L. sp.9RL, were found to be assigned to multiple subclade each, indicating the geographical divergences and the presence of cryptic diversity. Our findings of this study demonstrate that Qinghai-Tibet Plateau exhibits a higher species diversity of the Leptomias group, and it is imperative to investigate cryptic species within certain morphospecies using integrative taxonomic approaches in future studies. Moreover, the construction of a DNA barcode reference library presented herein establishes a robust foundational dataset to support forthcoming research on weevil taxonomy, phylogenetics, ecology, and evolution.

Variation of the 3'RR1 HS1.2 Enhancer and Its Genomic Context.

In humans, the HS1.2 enhancer in the Ig heavy-chain locus is modular, with length polymorphism. Previous studies have shown the following features for this variation: (i) strong population structuring; (ii) association with autoimmune diseases; and (iii) association with developmental changes in Ig expression. The HS1.2 region could then be considered as a contributor to inter-individual diversity in humoral response in adaptive immunity. We experimentally determined the HS1.2-length class genotype in 72 of the 1000 Genomes CEU cell lines and assigned the HS1.2 alleles to haplotypes defined by 18 landmark SNPs. We also sequenced the variable portion and ~200 bp of the flanking DNA of 34 HS1.2 alleles. Furthermore, we computationally explored the ability of different allelic arrangements to bind transcription factors. Non-random association between HS1.2 and Gm allotypes in the European population clearly emerged. We show a wealth of variation in the modular composition of HS1.2, with five SNPs further contributing to diversity. Longer alleles offer more potential sites for binding but, for same-length alleles, SNP variation creates/destroys potential binding sites. Altogether, the arrangements of modules and SNP alleles both inside and outside HS1.2 denote an organization of diversity far from randomness. In the context of the strong divergence of human populations for this genomic region and the reported disease associations, our results suggest that selective forces shaped the pattern of its diversity.

Long-term climatic stability drives accumulation and maintenance of divergent freshwater fish lineages in a temperate biodiversity hotspot

Anthropogenic climate change is forecast to drive regional climate disruption and instability across the globe. These impacts are likely to be exacerbated within biodiversity hotspots, both due to the greater potential for species loss but also to the possibility that endemic lineages might not have experienced significant climatic variation in the past, limiting their evolutionary potential to respond to rapid climate change. We assessed the role of climatic stability on the accumulation and persistence of lineages in an obligate freshwater fish group endemic to the southwest Western Australia (SWWA) biodiversity hotspot. Using 19,426 genomic (ddRAD-seq) markers and species distribution modelling, we explored the phylogeographic history of western (Nannoperca vittata) and little (Nannoperca pygmaea) pygmy perches, assessing population divergence and phylogenetic relationships, delimiting species and estimating changes in species distributions from the Pliocene to 2100. We identified two deep phylogroups comprising three divergent clusters, which showed no historical connectivity since the Pliocene. We conservatively suggest these represent three isolated species with additional intraspecific structure within one widespread species. All lineages showed long-term patterns of isolation and persistence owing to climatic stability but with significant range contractions likely under future climate change. Our results highlighted the role of climatic stability in allowing the persistence of isolated lineages in the SWWA. This biodiversity hotspot is under compounding threat from ongoing climate change and habitat modification, which may further threaten previously undetected cryptic diversity across the region.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits.

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Molecular characterization and virulence factor profiles of Staphylococcus aureus isolated from bovine mastitis in Brazilian herds

Staphylococcus aureus is the main etiological agent of bovine mastitis worldwide and knowledge about its diversity and virulence factors is vital in controlling infections caused by this pathogen. The present study aimed to perform molecular characterization of a population of S. aureus (n=153 strains isolated from 1994 to 2014 in seven Brazilian states) by pulsed-field gel electrophoresis (PFGE) and evaluate their virulence profiles via polymerase chain reaction (PCR). PFGE identified 93 pulsotypes, with the isolates organized into 26 clusters and 20 unique pulsotypes. Predominant pulsotypes were observed, with variations according to the years of isolation and geographic origin of the isolates. Based on the PCR results for the genes encoding agglutination factors (ClfA and ClfB), binding proteins (fibronectin binding protein - FnBPA, elastin binding protein - Ebps, collagen binding protein - Cna), and toxins (Hla, Hlb and Luk-ED), 40 virulence profiles were detected. The frequency of virulence genes ranged from 58 to 98% (clfA:84.3%; clfB and hlb both 81%; hla:71.2%; fnBA:82.3%; Can:94.7%; ebps:58%; and lukED:98%). The existence of prevalent genotypes in some of the Brazilian states and the time period studied suggests that these genotypes are better adapted, with favorable characteristics in host/pathogen relationships. Genes of proven importance for S. aureus pathogenesis in bovine mastitis were widely distributed in genetically divergent populations, suggesting that most of these genes may be interesting candidates in the development of vaccines to control bovine mastitis in Brazil.

Divergent population dynamics: The case of the inner city of Riga

Population dynamics shape the spatial landscape of Europe. Although broadly studied, both suburban sprawl and inner-city reurbanisation have often been treated in isolation. Furthermore, these processes manifest distinctively across different local contexts. This research aimed to explore the population dynamics in the inner city of Riga from 2011 to 2021, contrasting them with the bordering municipalities of Riga, and to analyze how these dynamics align with the urban development model. Covering 58 neighborhoods and 7 bordering municipalities of Riga, the analysis revealed a shift in the trajectory of inner-city population dynamics during the latter half of the decade, which was finally characterized by growth. However, suburbanisation continued to outpace reurbanisation, indicating the coexistence of multiple urban development model stages.

Phylogeny, Genetic Diversity and Population Structure of Fritillaria cirrhosa and Its Relatives Based on Chloroplast Genome Data.

Fritillaria cirrhosa and its relatives have been utilized in traditional Chinese medicine for many years and are under priority protection in China. Despite their medicinal and protective value, research on their phylogeny, genetic diversity, and divergence remains limited. Here, we investigate the chloroplast genome variation architecture of 46 samples of F. cirrhosa and its relatives collected from various regions, encompassing the majority of wild populations across diverse geographical areas. The results indicate abundant variations in 46 accessions including 1659 single-nucleotide polymorphisms and 440 indels. Six variable markers (psbJ, ndhD, ycf1, ndhG, trnT-trnL, and rpl32-trnL) were identified. Phylogenetic and network analysis, population structure analysis, and principal component analysis showed that the 46 accessions formed five clades with significant divergence, which were related to their geographical distribution. The regions spanning from the southern Hengduan Mountains to the Qinghai-Tibet Plateau exhibited the highest levels of genetic diversity. F. cirrhosa and its relatives may have suffered a genetic bottleneck and have a relatively low genetic diversity level. Moreover, geographical barriers and discrete patches may have accelerated population divergence. The study offers novel perspectives on the phylogeny, genetic diversity, and population structure of F. cirrhosa and its relatives, information that can inform conservation and utilization strategies in the future.

Nascent transcription reveals regulatory changes in extremophile fishes inhabiting hydrogen sulfide-rich environments.

Regulating transcription allows organisms to respond to their environment, both within a single generation (plasticity) and across generations (adaptation). We examined transcriptional differences in gill tissues of fishes in the Poecilia mexicana species complex (family Poeciliidae), which have colonized toxic springs rich in hydrogen sulfide (H2S) in southern Mexico. There are gene expression differences between sulfidic and non-sulfidic populations, yet regulatory mechanisms mediating this gene expression variation remain poorly studied. We combined capped-small RNA sequencing (csRNA-seq), which captures actively transcribed (i.e. nascent) transcripts, and messenger RNA sequencing (mRNA-seq) to examine how variation in transcription, enhancer activity, and associated transcription factor binding sites may facilitate adaptation to extreme environments. csRNA-seq revealed thousands of differentially initiated transcripts between sulfidic and non-sulfidic populations, many of which are involved in H2S detoxification and response. Analyses of transcription factor binding sites in promoter and putative enhancer csRNA-seq peaks identified a suite of transcription factors likely involved in regulating H2S-specific shifts in gene expression, including several key transcription factors known to respond to hypoxia. Our findings uncover a complex interplay of regulatory processes that reflect the divergence of extremophile populations of P. mexicana from their non-sulfidic ancestors and suggest shared responses among evolutionarily independent lineages.

Fine-scale genetic structure and phenotypic divergence of a passerine bird population inhabiting a continuous Mediterranean woodland.

Genetic differentiation between populations inhabiting ecologically different habitats might appear because of limited dispersal and gene flow, which may lead to patterns of phenotypic divergence and local adaptation. In this study, we use dispersal, genotypic (24 microsatellite loci) and phenotypic (body size and clutch size) data to analyse patterns of genetic structuring and phenotypic divergence in a blue tit (Cyanistes caeruleus) population inhabiting a continuous and heterogeneous woodland along a valley. The two slopes of the valley differ in their forest formations and environmental conditions. Findings showed that most blue tits reproduced within their natal slope. Accordingly, microsatellite analyses revealed that populations of blue tits established in the two slopes show subtle genetic differentiation. The two genetic populations diverged in clutch size, exceeding the level of differentiation expected based on genetic drift, hence suggesting divergent selection (or other processes promoting divergence) on this life-history trait. Our findings reveal that restricted dispersal and spatial heterogeneity may lead to genetic differentiation among bird populations at a surprisingly small scale. In this respect, it is worth highlighting that such differentiation occurs for an organism with high dispersal capacity and within a continuous woodland. Moreover, we show that small-scale ecological differences, together with limited gene flow, can result in selection favouring different phenotypes even within the same continuum population.

Identification of common human infectious and potentially zoonotic novel genotypes of Enterocytozoon bieneusi in cavernicolous bats in Thailand.

Enterocytozoon bieneusi is a common cause of human microsporidiosis and can infect a variety of animal hosts worldwide. In Thailand, previous studies have shown that this parasite is common in domestic animals. However, information on the prevalence and genotypes of this parasite in other synanthropic wildlife, including bats, remains limited. Several pathogens have been previously detected in bats, suggesting that bats may serve as a reservoir for this parasite. In this study, a total of 105 bat guano samples were collected from six different sites throughout Thailand. Of these, 16 from Chonburi (eastern), Ratchaburi (western), and Chiang Rai (northern) provinces tested positive for E. bieneusi, representing an overall prevalence of 15.2%. Based on ITS1 sequence analysis, 12 genotypes were identified, including two known genotypes (D and type IV) frequently detected in humans and ten novel potentially zoonotic genotypes (TBAT01-TBAT10), all belonging to zoonotic group 1. Lyle's flying fox (Pteropus lylei), commonly found in Southeast Asia, was identified as the host in one sample that was also positive for E. bieneusi. Network analysis of E. bieneusi sequences detected in this study and those previously reported in Thailand also revealed intraspecific divergence and recent population expansion, possibly due to adaptive evolution associated with host range expansion. Our data revealed, for the first time, multiple E. bieneusi genotypes of zoonotic significance circulating in Thai bats and demonstrated that bat guano fertilizer may be a vehicle for disease transmission.

Genetic architecture and correlations between the gut microbiome and gut gene transcription in Chinook salmon (Oncorhynchus tshawytscha).

Population divergence through selection can drive local adaptation in natural populations which has implications for the effective restoration of declining and extirpated populations. However, adaptation to local environmental conditions is complicated when both the host and its associated microbiomes must respond via co-evolutionary change. Nevertheless, for adaptation to occur through selection, variation in both host and microbiome traits should include additive genetic effects. Here we focus on host immune function and quantify factors affecting variation in gut immune gene transcription and gut bacterial community composition in early life-stage Chinook salmon (Oncorhynchus tshawytscha). Specifically, we utilized a replicated factorial breeding design to determine the genetic architecture (sire, dam and sire-by-dam interaction) of gut immune gene transcription and microbiome composition. Furthermore, we explored correlations between host gut gene transcription and microbiota composition. Gene transcription was quantified using nanofluidic qPCR arrays (22 target genes) and microbiota composition using 16 S rRNA gene (V5-V6) amplicon sequencing. We discovered limited but significant genetic architecture in gut microbiota composition and transcriptional profiles. We also identified significant correlations between gut gene transcription and microbiota composition, highlighting potential mechanisms for functional interactions between the two. Overall, this study provides support for the co-evolution of host immune function and their gut microbiota in Chinook salmon, a species recognized as locally adapted. Thus, the inclusion of immune gene transcription profile and gut microbiome composition as factors in the development of conservation and commercial rearing practices may provide new and more effective approaches to captive rearing.

Living in a post‐industrial landscape: repeated patterns of genetic divergence in brown trout (Salmo trutta L.) across the British Isles

AbstractAimThe British Isles have been worked for millennia to extract metal ores to feed industrial development, leaving a legacy of mine water pollution that continues to impact freshwater communities in many regions. Brown trout (Salmo trutta L.) have long been observed to persist in these metal‐impacted systems as apex predators, with previous studies showing a small number of impacted populations to be highly genetically divergent. We sought to understand the scale of genetic diversity across regions and the repeatability of genetic divergence in trout populations affected by metal pollution.LocationWe examined four mine water‐impacted regions across the British Isles: west Wales, northeast England, southwest England and southeast Ireland.MethodsWe employed a panel of 95 SNP loci to screen 1236 individuals from 71 sites representing paired metal‐impacted and clean sites from across the four regions. From these, we obtained diversity statistics, assessed genetic structuring of populations and modelled historical demographic scenarios to understand which factors most credibly explain genetic variation in divergent populations.ResultsWe evidenced hierarchical population structure in the regions studied, in line with expectations from phylogeographic history. However, in a hierarchical analysis of genetic structuring the first level of differentiation was driven by the divergence of the metal‐impacted trout of Cornwall in southwest England. Within regions we observed reduced genetic diversity and repeated patterns of local genetic sub‐structuring between paired samples from metal‐impacted and relatively clean sites. Demographic history analyses suggested the timing of these splits to be relatively recent and to be associated with periods of peak mining activity.Main conclusionsOur findings demonstrate distinct patterns of genetic isolation and reduced diversity arising from legacy pollution in freshwater ecosystems, with impacts being most apparent where both chemical pollution and physical barriers are present. Management should focus on the amelioration of mine water wash‐out and the removal of barriers to fish movement to safeguard genetic diversity in impacted populations.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana.

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Fine-scale adaptive divergence and population genetic structure of Aedes aegypti in Metropolitan Manila, Philippines

BackgroundThe adaptive divergence of Aedes aegypti populations to heterogeneous environments can be a driving force behind the recent expansion of their habitat distribution and outbreaks of dengue disease in urbanized areas. In this study, we investigated the population genomics of Ae. aegypti at a regional scale in Metropolitan Manila, Philippines.MethodsWe used the Pool-Seq double digestion restriction-site association DNA sequencing (ddRAD-Seq) approach to generate a high number of single nucleotide polymorphisms (SNPs), with the aim to determine local adaptation and compare the population structure with 11 microsatellite markers. A total of 217 Ae. aegypti individuals from seven female and seven male populations collected from Metropolitan Manila were used in the assays.ResultsWe detected 65,473 SNPs across the populations, of which 76 were non-neutral SNPs. Of these non-neutral SNPs, the multivariate regression test associated 50 with eight landscape variables (e.g. open space, forest, etc.) and 29 with five climate variables (e.g. air temperature, humidity, etc.) (P-value range 0.005–0.045) in female and male populations separately. Male and female populations exhibited contrasting spatial divergence, with males exhibiting greater divergence than females, most likely reflecting the different dispersal abilities of male and female mosquitoes. In the comparative analysis of the same Ae. aegypti individuals, the pairwise FST values of 11 microsatellite markers were lower than those of the neutral SNPs, indicating that the neutral SNPs generated via pool ddRAD-Seq were more sensitive in terms of detecting genetic differences between populations at fine-spatial scales.ConclusionsOverall, our study demonstrates the utility of pool ddRAD-Seq for examining genetic differences in Ae. aegypti populations in areas at fine-spatial scales that could inform vector control programs such as Wolbachia-infected mosquito mass-release programs. This in turn would provide information on mosquito population dispersal patterns and the potential barriers to mosquito movement within and around the release area. In addition, the potential of environmental adaptability observed in Ae. aegypti could help population control efforts.Graphical

Genetic Analysis of the Emerging Citrus Yellow Vein Clearing Virus Reveals a Divergent Virus Population in American Isolates.

Citrus yellow vein clearing virus is a previously reported citrus virus from Asia with widespread distribution in China. In 2022, the California Department of Food and Agriculture conducted a multipest citrus survey targeting multiple citrus pathogens including citrus yellow vein clearing virus (CYVCV). In March 2022, a lemon tree with symptoms of vein clearing, chlorosis, and mottling in a private garden in the city of Tulare, California, tested positive for CYVCV, which triggered an intensive survey in the surrounding areas. A total of 3,019 plant samples, including citrus and noncitrus species, were collected and tested for CYVCV using conventional reverse transcription polymerase chain reaction, reverse transcription quantitative polymerase chain reaction, and Sanger sequencing. Five hundred eighty-six citrus trees tested positive for CYVCV, including eight citrus species not previously recorded infected under field conditions. Comparative genomic studies were conducted using 17 complete viral genomes. Sequence analysis revealed two major phylogenetic groups. Known Asian isolates and five California isolates from this study made up the first group, whereas all other CYVCV isolates from California formed a second group, distinct from all worldwide isolates. Overall, the CYVCV population shows rapid expansion and high differentiation indicating a population bottleneck typical of a recent introduction into a new geographic area.

Within population plastic responses to combined thermal-nutritional stress differ from those in response to single stressors, and are genetically independent across traits in both males and females.

Phenotypic plasticity helps animals to buffer the effects of increasing thermal and nutritional stress created by climate change. Plastic responses to single and combined stressors can vary among genetically diverged populations. However, less is known about how plasticity in response to combined stress varies among individuals within a population or whether such variation changes across life-history traits. This is important because individual variation within populations shapes population-level responses to environmental change. Here, we used isogenic lines of Drosophila melanogaster to assess plasticity of egg-to-adult viability and sex-specific body size for combinations of two temperatures (25°C or 28°C) and three diets (standard diet, low caloric diet, or low protein:carbohydrate ratio diet). Our results reveal substantial within-population genetic variation in plasticity for egg-to-adult viability and wing size in response to combined thermal-nutritional stress. This genetic variation in plasticity was a result of cross-environment genetic correlations that were often < 1 for both traits, as well as changes in the expression of genetic variation across environments for egg-to-adult viability. Cross-sex genetic correlations for body size were weaker when the sexes were reared in different conditions, suggesting that the genetic basis of traits may change with the environment. Further, our results suggest that plasticity in egg-to-adult viability is genetically independent from plasticity in body size. Importantly, plasticity in response to diet and temperature individually differed from plastic shifts in response to diet and temperature in combination. By quantifying plasticity and the expression of genetic variance in response to combined stress across traits, our study reveals the complexity of animal responses to environmental change, and the need for a more nuanced understanding of the potential for populations to adapt to ongoing climate change.

Shared xerophytic genes and their re-use in local adaptation to aridity in the desert plant Gymnocarpos przewalskii.

In order to thrive and survive, plant species need to combine stability in the long term and rapid response to environmental challenges in the short term. The former would be reflected by parallel or convergent adaptation across species, and the latter by pronounced local adaptation among populations of the same species. In the present study, we generated a high-quality genome and re-sequenced 177 individuals for Gymnocarpos przewalskii, an important desert plant species from North-West China, to detect local adaptation. We first focus on ancient adaptation to aridity at the molecular level by comparing the genomic data of 15 species that vary in their ability to withstand aridity. We found that a total of 118 genes were shared across xerophytic species but absent from non-xerophytic species. Of the 65 found in G. przewalskii, 63 were under purifying selection and two under positive selection. We then focused on local adaptation. Up to 20% of the G. przewalskii genome showed signatures of local adaptation to aridity during population divergence. Thirteen of the selected shared xerophytic genes were reused in local adaptation after population differentiation. Hence, only about 20% of the genes shared and specific to xerophytic species and associated with adaptation to aridity were later recruited for local adaptation in G. przewalskii.

Divergent Aspergillus flavus corn population is composed of prolific conidium producers: Implications for saprophytic disease cycle

ABSTRACT The ascomycete fungus Aspergillus flavus infects and contaminates corn, peanuts, cottonseed, and tree nuts with toxic and carcinogenic aflatoxins. Subdivision between soil and host plant populations suggests that certain A. flavus strains are specialized to infect peanut, cotton, and corn despite having a broad host range. In this study, the ability of strains isolated from corn and/or soil in 11 Louisiana fields to produce conidia (field inoculum and male gamete) and sclerotia (resting bodies and female gamete) was assessed and compared with genotypic single-nucleotide polymorphism (SNP) differences between whole genomes. Corn strains produced upward of 47× more conidia than strains restricted to soil. Conversely, corn strains produced as much as 3000× fewer sclerotia than soil strains. Aspergillus flavus strains, typified by sclerotium diameter (small S-strains, <400 μm; large L-strains, >400 μm), belonged to separate clades. Several strains produced a mixture (M) of S and L sclerotia, and an intermediate number of conidia and sclerotia, compared with typical S-strains (minimal conidia, copious sclerotia) and L-strains (copious conidia, minimal sclerotia). They also belonged to a unique phylogenetic mixed (M) clade. Migration from soil to corn positively correlated with conidium production and negatively correlated with sclerotium production. Genetic differences correlated with differences in conidium and sclerotium production. Opposite skews in female (sclerotia) or male (conidia) gametic production by soil or corn strains, respectively, resulted in reduced effective breeding population sizes when comparing male:female gamete ratio with mating type distribution. Combining both soil and corn populations increased the effective breeding population, presumably due to contribution of male gametes from corn, which fertilize sclerotia on the soil surface. Incongruencies between aflatoxin clusters, strain morphotype designation, and whole genome phylogenies suggest a history of sexual reproduction within this Louisiana population, demonstrating the importance of conidium production, as infectious propagules and as fertilizers of the A. flavus soil population.

Evolvability predicts macroevolution under fluctuating selection.

Heritable variation is a prerequisite for evolutionary change, but the relevance of genetic constraints on macroevolutionary timescales is debated. By using two datasets on fossil and contemporary taxa, we show that evolutionary divergence among populations, and to a lesser extent among species, increases with microevolutionary evolvability. We evaluate and reject several hypotheses to explain this relationship and propose that an effect of evolvability on population and species divergence can be explained by the influence of genetic constraints on the ability of populations to track rapid, stationary environmental fluctuations.

Deep population structure linked to host vernalization requirement in the barley net blotch fungal pathogen.

Invasive fungal pathogens pose a substantial threat to widely cultivated crop species, owing to their capacity to adapt to new hosts and new environmental conditions. Gaining insights into the demographic history of these pathogens and unravelling the mechanisms driving coevolutionary processes are crucial for developing durably effective disease management programmes. Pyrenophora teres is a significant fungal pathogen of barley, consisting of two lineages, Ptt and Ptm, with global distributions and demographic histories reflecting barley domestication and spread. However, the factors influencing the population structure of P. teres remain poorly understood, despite the varietal and environmental heterogeneity of barley agrosystems. Here, we report on the population genomic structure of P. teres in France and globally. We used genotyping-by-sequencing to show that Ptt and Ptm can coexist in the same area in France, with Ptt predominating. Furthermore, we showed that differences in the vernalization requirement of barley varieties were associated with population differentiation within Ptt in France and at a global scale, with one population cluster found on spring barley and another population cluster found on winter barley. Our results demonstrate how cultivation conditions, possibly associated with genetic differences between host populations, can be associated with the maintenance of divergent invasive pathogen populations coexisting over large geographic areas. This study not only advances our understanding of the coevolutionary dynamics of the Pt-barley pathosystem but also prompts further research on the relative contributions of adaptation to the host versus adaptation to abiotic conditions in shaping Ptt populations.