A revised phylogeny of Boletaceae using whole genome sequences.

The porcini mushroom family Boletaceae is a diverse, widespread group of ectomycorrhizal (ECM) mushroom-forming fungi that so far has eluded intrafamilial phylogenetic resolution based on morphology and multilocus data sets. In this study, we present a genome-wide molecular data set of 1764 single-copy gene families from a global sampling of 418 Boletaceae specimens. The resulting phylogenetic analysis has strong statistical support for most branches of the tree, including the first statistically robust backbone. The enigmatic Phylloboletellus chloephorus from non-ECM Argentinian subtropical forests was recovered as a new subfamily sister to the core Boletaceae. Time-calibrated branch lengths estimate that the family first arose in the early to mid-Cretaceous and underwent a rapid radiation in the Eocene, possibly when the ECM nutritional mode arose with the emergence and diversification of ECM angiosperms. Biogeographic reconstructions reveal a complex history of vicariance and episodic long-distance dispersal correlated with historical geologic events, including Gondwanan origins and inferred vicariance associated with its disarticulation. Together, this study represents the most comprehensively sampled, data-rich molecular phylogeny of the Boletaceae to date, establishing a foundation for future robust inferences of biogeography in the group.

Contrasting continental patterns of adaptive population divergence in the holarctic ectomycorrhizal fungus Boletus edulis.

In the hyperdiverse fungi, the process of speciation is virtually unknown, including for the > 20 000 species of ectomycorrhizal mutualists. To understand this process, we investigated patterns of genome-wide differentiation in the ectomycorrhizal porcini mushroom, Boletus edulis, a globally distributed species complex with broad ecological amplitude. By whole-genome sequencing 160 individuals from across the Northern Hemisphere, we genotyped 792 923 single nucleotide polymorphisms to characterize patterns of genome-wide differentiation and to identify the adaptive processes shaping global population structure. We show that B. edulis exhibits contrasting patterns of genomic divergence between continents, with multiple lineages present across North America, while a single lineage dominates Europe. These geographical lineages are inferred to have diverged 1.62-2.66 million years ago, during a period of climatic upheaval and the onset of glaciation in the Pliocene-Pleistocene boundary. High levels of genomic differentiation were observed among lineages despite evidence of substantial and ongoing introgression. Genome scans, demographic inference, and ecological niche models suggest that genomic differentiation is maintained by environmental adaptation, not physical isolation. Our study uncovers striking patterns of genome-wide differentiation on a global scale and emphasizes the importance of local adaptation and ecologically mediated divergence, rather than prezygotic barriers such as allopatry or genomic incompatibility, in fungal population differentiation.

Lost in translation: Population genomics and long-read sequencing reveals relaxation of concerted evolution of the ribosomal DNA cistron.

Concerted evolution of the ribosomal DNA array has been studied in numerous eukaryotic taxa, yet is still poorly understood. rDNA genes are repeated dozens to hundreds of times in the eukaryotic genome (Eickbush and Eickbush, 2007) and it is believed that these arrays are homogenized through concerted evolution (Zimmer et al., 1980; Dover, 1993) preventing the accumulation of intragenomic, and intraspecific, variation. However, numerous studies have reported rampant intragenomic and intraspecific variation in the rDNA array (Ganley and Kobayashi, 2011; Naidoo et al., 2013; Hughes and Petersen, 2001; Lindner and Banik, 2011; Li et al., 2013; Lindner et al., 2013; Hughes et al., 2018), contradicting our current understanding of concerted evolution. The internal transcribed spacers (ITS) of the rDNA cistron are the most commonly used DNA barcoding region in Fungi (Schoch et al., 2012), and rely on concerted evolution to homogenize the rDNA array leading to a "barcode gap" (Puillandre et al., 2012). Here we show that in Boletus edulis Bull., ITS intragenomic variation persists at low allele frequencies throughout the rDNA array, this variation does not correlate with genomic relatedness between populations, and rDNA genes may not evolve in a strictly concerted fashion despite the presence of unequal recombination and gene conversion. Under normal assumptions, heterozygous positions found in ITS sequences represent hybridization between populations, yet through allelic mapping of the rDNA array we found numerous heterozygous alleles to be stochastically introgressed throughout, presenting a dishonest signal of gene flow. Moreover, despite the signal of gene flow in ITS, our organisms were highly inbred, indicating a disconnect between true gene flow and barcoding signals. In addition, we show that while the mechanisms of concerted evolution are ongoing in pseudo-heterozygous individuals, they are not fully homogenizing the ITS array. Concerted evolution of the rDNA array may insufficiently homogenize the ITS gene, allowing for misleading signals of gene flow to persist, vastly complicating the use of the ITS locus for DNA barcoding in Fungi.