This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Scaling macro- and microevolutionary dynamics in the Caesar’s mushrooms (Amanita sect. Caesareae)
by
Santiago Sánchez-Ramírez
A thesis submitted in conformity with the requirements for the degree of PhD in Ecology and Evolutionary Biology
Ecology and Evolutionary Biology University of Toronto
Ramírez, Felipe Ruán (and two of his students), and students Nallely Hernández and Mariano
Torres. I must thank many Amanita collectors that provided personal samples (US: Erin
Feldman; Canada: Renée Lebeuf, Greg Thorn; Europe: Mario Amalfi), and various Fungaria
personnel throughout North America.
Finally, the Mexican National Council of Science and Technology (CONACyT) deserves credit
for awarding me a full scholarship. I would also like to extend my gratitude to Dr. Daniel Piñero
for providing financial aid in times of need. The Centre for Global Change Science (CGCS) at
UofT also funded three months of fieldwork in Mexico and the eastern US. Both SciNet and
CAGEF provided access to high-performance computing infrastructure. Dr. Dax Torti at the
Donelly Sequencing Centre also provided advice and performed all wet lab procedures on exon-
targeted sequencing.
vi
Table of Contents ACKNOWLEDGMENTS!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!#$!
TABLE OF CONTENTS!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!$#!
LIST OF TABLES!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!%!
LIST OF FIGURES!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!%#!
LIST OF APPENDICES!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!%$!
1.1! BIOLOGICAL DIVERSITY AND EVOLUTIONARY SCALES!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!&!1.2! THE GENOMIC ERA!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!'!1.3! WHY FUNGI?!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!(!1.4! THE SYSTEM!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!)!1.5! QUESTIONS, OBJECTIVES, RESEARCH HYPOTHESES, AND SCIENTIFIC PRODUCTION!"""""""""""""""""""""!)!"#$#"! %&'()*+!,!#######################################################################################################################################################################!$!"#$#,! %&'()*+!-!#######################################################################################################################################################################!$!"#$#-! %&'()*+!.!#######################################################################################################################################################################!/!"#$#.! %&'()*+!$!#######################################################################################################################################################################!/!"#$#$! %&'()*+!/!#######################################################################################################################################################################!0!
! EVOLUTIOANRY HISTORY AND GLOBAL BIOGEOGRAPHY OF AMANITA SECT. 2
CAESAREAE!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!*!2.1! INTRODUCTION!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!*!2.2! MATERIAL AND METHODS!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!&&!,#,#"! 1'23+')3+4!(+3)3%315!'67!7')'!8')+9%*5!###################################################################################################!""!,#,#,! %'192+')936!5)+')*:9*5!'67!79;*+:*6%*<)98*!*5)98')936!######################################################################!",!,#,#-! '6%*5)+'1<'+*'!+*%365)+=%)936!'67!:*3795(*+5'1!837*15!##################################################################!"-!,#,#.! 79;*+59>9%')936!+')*!)*5)5!###################################################################################################################################!".!2.3! RESULTS!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!&+!
! WHOLE-GENOME SEQUENCING AND ANNOTATION OF AMANITA JACKSONII AND 4
PARTIAL GENOMIC SEQUENCES OF AMANITA BASII!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-&!4.1! INTRODUCTION!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-&!
viii
4.2! SEQUENCED STRAIN AND ACCESSION NUMBERS!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-'!4.3! MATERIAL AND METHODS!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-(!4.4! RESULTS AND DISCUSSION!""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-(!4.5! PARTIAL GENOMIC SEQUENCES FROM AMANITA BASII!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-)!4.6! DEVELOPMENT OF NOVEL GENOMIC LOCI!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!--!
! SPECIES DELIMITATION, COMPARATIVE PHYLOGEOGRAPHY AND GLACIAL 5
REFUGIA IN THE NORTH AMERICAN CAESAR’S MUSHROOM SPECIES COMPLEX!"""""""!-+!5.1! INTRODUCTION!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-+!5.2! MATERIAL AND METHODS!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!-,!$#,#"! 5'8(196:!######################################################################################################################################################################!$B!$#,#,! 1'23+')3+4!(+3)3%315!##########################################################################################################################################!$B!$#,#-! 7')'!8')+9%*5A!(&'596:!3>!&*)*+3C4:3)*5A!'67!76'!837*15!################################################################!$?!$#,#.! :*6*!)+**5A!%36%')*6')936A!'67!:*6*'13:9%'1!53+)96:!###########################################################################!$?!$#,#$! 2'4*59'6!5(*%9*5!7*1989)')936!############################################################################################################################!/D!$#,#/! 5(*%9*5!)+**!*5)98')936!#########################################################################################################################################!/"!$#,#0! 5(*%9*5!795)+92=)936!837*196:!###########################################################################################################################!/"!$#,#B! '6%*5)+'1!'+*'!+*%365)+=%)936!'67!79;*+59>9%')936!+')*5!#################################################################!/,!$#,#?! 5(*%9*5E!(3(=1')936!59C*!*5)98')936!>+38!)&*!8=1)9<5(*%9*5!%3'1*5%*6)!#########################################!/-!$#,#"D! +*:+*559365!'67!7*83:+'(&9%!&95)3+9*5!#####################################################################################################!/.!$#,#""! (3(=1')936!:*6*)9%!8*)+9%5!#############################################################################################################################!/$!5.3! RESULTS!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!+-!$#-#"! :*6*!)+**5!'67!%36%')*6')*7!7')'!#################################################################################################################!/$!$#-#,! 5(*%9*5!)+**!*5)98')936!'67!5(*%9*5!7*1989)')936!#####################################################################################!/B!$#-#-! 69%&*!837*196:A!293:*3:+'(&9%!+*%365)+=%)936!'67!79;*+59>9%')936!+')*5!##################################!/?!$#-#.! &95)3+9%'1!7*83:+'(&4!'67!(3(=1')936!:*6*)9%5!#####################################################################################!/?!5.4! DISCUSSION!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!/(!$#.#"! 5(*%9*5!7*1989)')936!'67!%+4()9%!79;*+59)4!##################################################################################################!0-!$#.#,! +*>=:9'1!8'%+3!'67!89%+3<*;31=)936'+4!746'89%5!##################################################################################!00!$#.#-! (3)*6)9'1!%';*')5!###################################################################################################################################################!B,!5.5! CONCLUSIONS!"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""!,(!
List of Tables )'21*!,#"!19G*19&337<+')93!)*5)!%38('+9536!3>!:*3795(*+5'1!837*15!96!!"!#$%!!5*%)#!&!'(!)'!'#!#############################!,D!)'21*!,#,!%38('+9536!3>!296'+4<5)')*!5(*%9')936!'67!*@)96%)936!H2955*I!837*15!2'5*7!36!'6!'6'14595!3>!;'+9'6%*!
)'21*!$#,#!POPULATION GENETIC SUMMARY STATISTICS PER SPECIES PER GENE.!#####################################################################!0D!TABLE 6.1. SUMMARIZED POLYMORPHISM DATA FOR 502 GENES IN THE A. JACKSONII COMPLEX. N, SAMPLE SIZE; S,
SEGREGATING SITES; !, NUCLEOTIDE DIVERSITY; ", WATTERSON’S THETA, DTAJ, TAJIMA’S D. IN THE TABLE
BELOW, POLYMORPHISM DATA IS SEPARATED BY SITE CLASS.!#############################################################################################!?$!
xi
List of Figures >9:=+*!,#"#!:+'(&9%'1!+*(+*5*6)')936!3>!)&*!:*3795(*+5'1!837*15!96!!+,5*%)#!&!'(!)'!'!96>*++*7!J9)&!1':+'6:*#!
>9:=+*!$#,#!MAXENT SPECIES DISTRIBUTION MODELING OF NORTH AMERICAN CAESAR’S MUSHROOMS IN PRESENT (A)
AND LAST GLACIAL MAXIMUM (B; C. 22 KA) CLIMATIC CONDITIONS. THE SCALE IN A INDICATES HABITAT
SUITABILITY SCORES. C AND D ARE DISTRIBUTION MAPS OF SAMPLES INCLUDED IN THIS STUDY COLORED BY
SPECIES GROUP. GREY SHADES IN C AND D REPRESENT FORESTED (TEMPERATE MIXED AND CONIFEROUS
xiii
FORESTS) AREAS ACCORDING TO WWF TERRESTRIAL ECOREGIONS (OLSON ET AL. 2001). THE DOTTED LINES IN
C AND D DELINEATE THE BOUNDARIES BETWEEN REFUGIAL AND EXPANSION AREAS.!################################################!0-!>9:=+*!$#-#!RECONSTRUCTED BAYESIAN SPECIES TREE FROM *BEAST SHOWING MEAN NODE HEIGHTS, AND
ANCESTRAL AREA RECONSTRUCTIONS FROM RASP (PIE CHARTS) AS WELL AS STOCHASTIC CHARACTER MAPPING
FROM DIVERSITREE. COLOR CODES REPRESENT SPECIES IN REFUGIAL (BLUE), WIDESPREAD (GREEN), OR IN
EXPANSION (RED) AREAS. TERMINALS IN GREY AND DOTTED LINEAGES WERE IGNORED AS THEY EITHER
REPRESENT EURASIAN OUTGROUPS OR WERE DEEMED OUTSIDE NORTH AMERICA (SANCHEZ-RAMIREZ ET AL.
2015). PLI = PLIOCENE, PLE = PLEISTOCENE. THE BOXPLOT ON THE LEFT SHOWS PER AREA MCMC ESTIMATES
OF SPECIATION (#), EXTINCTION (µ) AND DISPERSAL ($) RATES FROM MUSSE (DIVERSITREE).!########################!0.!>9:=+*!$#.#!(A) A PLOT OF SPECIES’ LATITUDES VERSUS THE MAGNITUDE OF POPULATION EXPANSION (M). THE DOTTED
LINE INDICATES M = 1. (B) PLOT OF LATITUDE VERSUS THE POPULATION SIZE AT PRESENT TIME FOR EXTANT
SPECIES AND BEFORE SPECIATION FOR ANCESTRAL SPECIES (!1). IN A AND B, THE THICK SOLID LINE
REPRESENTS THE BEST GAUSSIAN FIT TO THE DATA. (C) PLOT OF THE EXTENT OF RANGE SHIFT VERSUS M. IN C
THE SOLID LINE REPRESENTS THE BEST LINEAR FIT TO THE DATA (R2 = 0.23, P-VALUE = 0.002). IN ALL PLOTS,
GREY AND BLACK CIRCLES REPRESENT EXTANT AND ANCESTRAL SPECIES, RESPECTIVELY, AND THINNER SOLID
LINES REPRESENT VARIATION WITHIN THE DATA (E.G. POSTERIOR ESTIMATES AND SAMPLED LATITUDES). B,
AMANITA BASII; C, A. SP. ‘COCHISEANA’; J, A. JACKSONII; J2, A. SP-JACK2; J3, A. SP-JACK3; J5 A. SP-JACK5; J6 A. SP-
JACK6. * MEDIAN, † MEAN.!##############################################################################################################################################################!0?!>9:=+*!$#$#!PLOT OF DEMOGRAPHIC CHANGE AS A FUNCTION OF TIME FOR ANCESTRAL AND EXTANT SPECIES. THE
FIRST UPPER ROW SHOWS FOUR EXTANT SPECIES WITH FAST DEMOGRAPHIC CHANGE. THE SECOND ROW
INCLUDES FOUR EXTANT SPECIES WITH NEAR CONSTANT HISTORICAL DEMOGRAPHY. THE TWO LAST ROWS SHOW
EIGHT ANCESTRAL SPECIES WITH RELATIVELY FAST DEMOGRAPHIC CHANGE. IN THESE THE NUMBERS REPRESENT
NODES SHOWN IN FIG. 5.3.!###############################################################################################################################################################!B"!>9:=+*!/#"#!BOX-PLOTS OF TAJIMA’S D OVER 502 GENES SEPARATED BY SITE-CLASSES (NON-SYNONYMOUS, LIGHT
GREY; SYNONYMOUS, WHITE; INTRONS, DARK GREY) FOR THE A. JACKSONII COMPLEX.!##############################################!?0!>9:=+*!/#,#!BOX-PLOTS OF % RATIOS FOR DIVERGENCE (%D) AND POLYMORPHISM (%P) (UPPER PANEL), AND & VALUES
THE STANDARD MKT TEST (GREY) AND A MODIFIED VERSION THAT TAKES INTO ACCOUNT SITE-FREQUENCY
SPECTRA (WHITE) (MESSER AND PETROV 2013A).!################################################################################################################!"DD!>9:=+*!/#-#!DISTRIBUTIONS OF 'GSI AND VARGSI FOR ALL 502 GENES. 'GSI INDICATES HOW WELL SPECIES CLADES ARE
RESOLVED, AND VARGSI SHOWS THE EXTENT OF DEVIATION FROM THE MAXIMUM GSI VALUE POSSIBLE, 1.!########!"DD!>9:=+*!/#.#!SCATTER-PLOTS SHOWING THE RELATIONSHIP BETWEEN TAJIMA’S D AND GSI FOR EACH GENE, FOR EACH
SPECIES.!##############################################################################################################################################################################################!"D,!>9:=+*!/#$#!DEMOGRAPHIC TRENDS THROUGH TIME FOR EACH SPECIES BASED ON THE EBSP AND *BEAST MODELS.
MEAN POPULATION SIZE VALUES AND DOTTED LINES REPRESENT MEDIAN VALUES. IN *BEAST PLOTS, THE
DOTTED LINE REPRESENT MEAN VALUES. THE LOWER PHYLOGENETIC TREE IS THE MAXIMUM-CLADE-
CREDIBILITY SPECIES TREE. HORIZONTAL BARS REPRESENT DIVERGENCE 95% CI. NUMBERS OVER NODES ARE
xiv
CLADE POSTERIOR PROBABILITIES. VERTICAL DOTTED LINES POINT TO THE START OF THE PLEISTOCENE EPOCH
AND THE ONSET OF DEMOGRAPHIC EXPANSIONS.!##################################################################################################################!"D-!>9:=+*!/#/#!GENE ONTOLOGY TERMS WORDLES FOR GENES UNDER POSITIVE SELECTION (A: & > 0, B: PAML LTR)
COMBINED FOR ALL SPECIES. THE SIZE OF THE TERM INCLUDES THE EXTENT OF SELECTION AND HOW OFTEN THE
TERM WAS FOUND ACROSS GENES AND SPECIES.!###################################################################################################################!"D$!>9:=+*!/#0#!SPECIES DISTRIBUTION MAPS BASED ON GEOGRAPHIC DATA FROM SÁNCHEZ-RAMÍREZ ET AL. (2015) AND
MAXENT (PHILLIPS ET AL. 2006). THE DOTTED LINE MARKS AN HYPOTHESIZED LIMIT FOR GLACIAL REFUGIA.
>9:=+*!%B#$#!*BEAST SPECIES TREES USED FOR SPECIES DELIMITATION WITH CLADE POSTERIOR PROBABILITIES
(UPPER), BP&P SPECIES PROBABILITIES (LOWER AND BOLD), AND BRANCHING TIMES.!###########################################!"?.!>9:=+*!%B#/#!REGRESSION ANALYSES ON THE LATITUDE VERSUS POPULATION SIZE AND EXPANSION DATA. THE BLUE
LINE REPRESENTS THE LINEAR MODEL, WHILE THE GREEN AND RED REPRESENT THE GAUSSIAN AND
Figure 5.2. MaxEnt species distribution modeling of North American Caesar’s mushrooms in present (A) and last glacial maximum (B; c. 22 ka) climatic conditions. The scale in A
indicates habitat suitability scores. C and D are distribution maps of samples included in this study colored by species group. Grey shades in C and D represent forested (temperate
mixed and coniferous forests) areas according to WWF terrestrial ecoregions (Olson et al. 2001). The dotted lines in C and D delineate the boundaries between refugial and expansion
areas.
124.24) and Gaussian (AICm = 114.30, AIC! = 117.41) regressions (Appendix Fig. C8.6). For
both parameters, a Gaussian function was comparatively the best descriptor (Akaike weight for
m = 0.68 and for !1 = 0.88) (Fig. 5.4A,B). The least-squares estimation for m resulted in three
out of four parameters with a significant (alpha = 0.05) fit (n = 24.52, P-value = 2x10-15; s =
5.57, P-value = 0.02; a = 2.18, P-value = 5x10-3). A significant positive correlation was found
between the extent of range expansions and m (Fig. 5.4C, R2 = 0.23, P-value = 0.002). I detected
at least 12 population size changes with m > 1 (Fig. 5.4A, 5), of which four occurred in extant
species, and six other in ancestral species (Fig. 5.3). Moreover, six of these demographic changes
occurred in species found to be in exclusively refugial, five in widespread regions, and one in
exclusively expansion areas (Fig. 5.3-5). Nine out of 12 population size changes occurred during
the Pleistocene (Fig. 5.5 only shows seven, plots of A. sp-jack5 and A. sp-jack6 are not shown).
The average number of segregating sites in tef1, AJ0000-2, and AJ0103-132, were 7, 10, and 12,
respectively, with the highest values found in A. basii (32, 46, and 42, respectively). In the same
order, the average " per site for each gene was 0.016, 0.011, 0.013, being the highest in A. basii
for tef1 (0.181) and AJ0103-132 (0.072), and in A. sp-jack2 for AJ0000-2 (0.037). The average
Tajima’s D for each gene was -0.72, -0.29, -0.44, respectively. Values ranged from 2.03 to -2.41
across all genes. Amanita jacksonii had the lowest Tajima’s D for tef1 (-2.02) and AJ0103-132 (-
2.41), and A. basii for AJ0000-2 (-2.03). Details can be found in Table 5.2.
5.4 Discussion
5.4.1 Species delimitation and cryptic diversity
Cryptic species are defined as “two or more distinct species that are erroneously classified (and
hidden) under one species name” (Bickford et al. 2006). In most cases, the underestimation of
74
Figure 5.3. Reconstructed Bayesian species tree from *BEAST showing mean node heights,
and ancestral area reconstructions from RASP (pie charts) as well as stochastic character
mapping from DIVERSITREE. Color codes represent species in refugial (blue), widespread
(green), or in expansion (red) areas. Terminals in grey and dotted lineages were ignored as
they either represent Eurasian outgroups or were deemed outside North America
(Sanchez-Ramirez et al. 2015). Pli = Pliocene, Ple = Pleistocene. The boxplot on the left
shows per area MCMC estimates of speciation (!), extinction (µ) and dispersal (") rates
from MuSSE (DIVERSITREE).
75
species diversity is due to little or imperceptible morphological, ecological, or ethological
differences among cryptic taxa. Macrofungi themselves are cryptic in nature; remaining hidden
underground most their lifetime, only being conspicuous through fruiting bodies, on which most
taxonomic work is focused. Their fruiting bodies also have overall simple anatomy, which make
species boundaries based on morphology difficult to assess (Taylor et al. 2000, 2006). Since
mycologists started adopting molecular methods (Yang 2011) for species delimitation, many
cryptic species are now recognized (Vilgalys and Sun 1994; James et al. 1999; Matute et al.
2006; Arnold 2007; Le Gac et al. 2007; Carriconde et al. 2008; Geml et al. 2008; Crespo and
Lumbsch 2010; Sánchez-Ramírez et al. 2015a). Previous morphotaxonomic work (e.g. Guzmán
and Ramírez-Guillén 2001; Tulloss 2015) already pointed out suspicions of multiple cryptic
species in the A. caesarea-complex, but lacked genetic confirmation. This study not only
confirms part of this diversity genetically, it shows that cryptic lineages have different temporal
and spatial origins. Nineteen species were delimited from the combined information from
individual and concatenated gene trees, the gsi, and the multi-locus species delimitation based on
the multi-species coalescent. From these, only seven have been described, and four (A. sp.
“cochiseana”; A. sp.-AR01, A. sp.-F11, A. sp.-T31) other have been provisionally segregated
based on distinctive morphological features (Tulloss 2015). This means that between 63 and 42%
of species diversity was unknown for the NA Caesar’s mushrooms. In the Caesar’s mushroom
system, A. jacksonii, which was previously thought to have an ample distribution ranging from
Canada to Central America (Guzmán and Ramírez-Guillén 2001; Tulloss 2015), sets a good
example. Based on my results, I now recognize at least nine other cryptic species throughout the
range resembling A. jacksonii (Fig. 5.2-4). Some of these show complex speciation patterns with
no clear geographic separation (A. sp-AR01, A. sp-jack1, A. sp-jack2, A. sp-T31 and A. sp-W15
in SE USA), others have some degree of geographical structure (A. sp-jack3, A. sp-jack5, A. sp-
jack6 in southern Mexico). Except for A. sp-jack2, which can be found from Indiana (USA) to
Chiapas (Mexico), most species are range restricted. In the true A. jacksonii Pomerl., I delimit
both its genetic and geographic boundaries, showing that there is little geographic and ecological
overlap with other species.
A case that merits further exploration concerns the relationship between A. hayalyuy D. Arora
and G. H. Shepard and A. garabitoana Tulloss, Halling and G. M. Muell (Shepard et al. 2008;
Tulloss et al. 2011). These two names have been applied to material described from Chiapas
76
(Mexico) and Costa Rica, respectively. Both have similar morphological features and inhabit
similar habitats (e.g. pine-oak forests) (Shepard et al. 2008; Tulloss et al. 2011). Based on my
initial gene tree/gsi assessment, I considered both as a single entity, however I encourage further
analyses with larger genetic and geographic sampling in A. hayalyuy/garabitoana. Overall, part
of my up coming work is to formally publish the newly found species, and to match up genetic
entities to type specimens.
Incorporating the information of multiple unlinked loci can help improve species delimitations
(Carstens and Knowles 2007). In part, this has led to the use and development of multi-locus
approaches in phylogenetic and phylogeographic analyses (Knowles 2009; Fujita et al. 2012).
Some of the benefits of Bayesian species delimitation models based on the multi-species
coalescent are: (1) the assumption of independent gene histories and stochastic coalescent
processes, which take into account potential sources of gene tree incongruence such as
incomplete lineage sorting, (2) the multi-dimensional sampling of species delimitation models,
and (3) the quantification of uncertainty (Yang and Rannala 2010, 2014; Rannala and Yang
2003, 2013). Nonetheless, gene matrix concatenation (e.g. Rokas et al. 2003) is still one of the
most widely used approaches in phylogenetic inference, in spite of the adverse effects reported in
some groups (Degnan and Rosenberg 2006, 2009; Pollard et al. 2006; Carstens and Knowles
2007; Heled and Drummond 2010). The concatenation analysis was mostly consistent with
BP&P’s multi-locus delimitation (Fig. 5.1, Appendix Fig. C8.5). There were some minor
topological differences across gene trees, in particular at deeper nodes (Fig. 5.1), however, 14
out of 19 (74%) species were reciprocally monophyletic with posterior probability of > 0.99
(Table 5.1). Most inconsistencies were detected either at deeper nodes (deep coalescence) or
between species with shallow divergences (from 1.48 [1.08–2.13] to 0.41 [0.26–0.64] Ma) and/or
large population sizes (Table 5.1: A. basii and A. sp. “cochiseana” in clade ‘basii’; A. sp-jack2,
A. sp-jack5, A. sp-jack6 in clade ‘jacksonii’).
Incomplete lineage sorting within recently diverged species is expected, particularly if species
(and ancestral species) have large effective population sizes (Maddison and Knowles 2006). On
this front, I took advantage of the multi-species coalescent model in order to (1) mitigate the
effects of incomplete lineage sorting in phylogenetic reconstruction and divergence time
estimations, and to (2) estimate historical population sizes based on multi-locus data. The multi-
species coalescent model is considered a paradigmatic model, bridging between micro- and
77
macroevolutionary processes (Edwards 2009). Nonetheless, its demographic inferences are
rarely used in a hypothesis-testing context. Here, I intend to extend the range of its applications,
not only to delimit species and estimate species trees, but to integrate both macro- (speciation,
extinction, dispersal) and microevolutionary processes (population sizes) to better understand
diversity patterns in the context of Quaternary climate change.
5.4.2 Refugial macro and micro-evolutionary dynamics
Higher species richness in refugia is expected (Keppel et al. 2012). Most of the diversity in the
NA Caesar’s mushrooms was found in refugial areas, which together with species in widespread
areas make up about 80% of the total diversity (Fig. 5.1-3). This observation is consistent with
other studies that show high diversity patterns in refugia (Médail and Diadema 2009; Tzedakis
2009; Keppel et al. 2012). The nature of high refugial diversity is debatable, especially because
there are many factors, such as speciation, extinction, and dispersal which contribute to diversity
build-up (Birand et al. 2012). Some authors suggest that refugia harbour high diversity due to
ecological, geographic, and climatic conditions that accelerate diversification (Quian and
Ricklefs 2000; Jetz et al. 2004; Svenning and Skov 2007). MuSSE results show that high
speciation rates in refugial areas –but not high extinction in expansion areas or higher dispersal
to refugial areas– better explains current diversity patterns in the Caesar’s mushrooms (Fig. 5.3).
Particularly, most refugial species are found in Mexico and Central America, which are
important centers of diversification (Perry et al. 1998; Escalante et al. 2004; Zarza et al. 2008;
Bryson et al. 2011; De-Nova et al. 2012), and diverged relatively recently (Figs. 5.2 and 5.4)
suggesting a role of Pleistocene glaciations in promoting genetic isolation. A recent review by
Mastretta-Yanes et al. (2015) highlights many of the natural processes that promote divergence
and speciation in the Mexican highlands. In particular, they point out co-occurring climatic and
topographical changes in the Trans-Mexican Volcanic Belt during the Pleistocene, and their
implications for recent diversification. Another area with high relative diversity is SE USA,
which is considered a well-known refugium (Soltis et al. 2006). This view is also supported by
paleoclimatic modeling (Fig. 5.2A,B). In this region most species are sympatric, with the
exception of A. sp-F11, which is confined to temperate forests in Florida (Fig. 5.2). To show the
degree of overlap, I sampled four species of this complex (little morphological differences)
within few meters of each other in Duke Forest, North Carolina (J.-M. Moncalvo pers. obs.).
78
!
"
#
$
%
&' "! "' (! (' #! #'
!
"
#
$
%
! " # $ %
!
"
#
$
%
!)*+,-./012/)3.4/1567+89
-2-6:.512/)01;+)*!
&89!)*+,-./012/)3.4/1567+89
<.51567+=
+,5+/5)2>)?./4+)0@1>5=
A,5./5B/C+05?.:
D
D
E
E
C
C
D"
D"D'
D'
D(
D(
D$
D$
B
F
G
79
Figure 5.4. (A) A plot of species’ latitudes versus the magnitude of population expansion
(m). The dotted line indicates m = 1. (B) Plot of latitude versus the population size at
present time for extant species and before speciation for ancestral species (!1). In A and B,
the thick solid line represents the best Gaussian fit to the data. (C) Plot of the extent of
range shift versus m. In C the solid line represents the best linear fit to the data (R2 = 0.23,
P-value = 0.002). In all plots, grey and black circles represent extant and ancestral species,
respectively, and thinner solid lines represent variation within the data (e.g. posterior
estimates and sampled latitudes). b, Amanita basii; c, A. sp. ‘cochiseana’; j, A. jacksonii; j2,
A. sp-jack2; j3, A. sp-jack3; j5 A. sp-jack5; j6 A. sp-jack6. * median, † mean.
Divergence times in the species tree (Fig. 5.3) indicate that most species within SE USA
diverged before the Pleistocene.
Regression models did not support my prediction of a linear relationship of population size with
latitude. Nonetheless, we found an interesting Gaussian relationship on both metrics (m and !1)
showing that population sizes and expansions are larger at intermediate latitudes (Fig. 5.4A,B).
This observation suggests that historical and contemporary population sizes have been large
within a range between 20 and 35 degrees latitude. This finding gives some support to my “high
genetic diversity” hypothesis, showing that high genetic variation can be found within a refugial
range. We also note that ancestral species with this type of response are concentrated within this
range, which might be related to the high speciation rate found in refugial areas (Fig. 5.3), and/or
to high rates of speciation at worldwide temperate latitudes (Sánchez-Ramírez et al. 2015b). It is
interesting that predictions about Pleistocene refugia theory are mostly based on European
studies (Taberlet et al. 1998; Hewitt 1999, 2004a; Petit et al. 2003; Schmitt 2007). Notably
differences in topography (e.g. east-west versus north-south orogeny) can have implications for
different diversity patterns in both continents (Anderson et al. 2013). For instance, north-south
oriented mountain chains are more likely to cause refugia and endemisms along costal
environments (Brunsfeld et al. 2001). Also higher connectivity between temperate and tropical
environments (e.g. southern US and Mexico/Central America) make the extent of refugial limits
broader in North America, compared to Europe, where fragmentation is larger. Nonetheless,
patterns in the NA Caesar’s mushrooms suggest that high population sizes at lower latitudes (e.g.
southern end of the distribution, Central America) are probably less tenable because of less forest
coverage. Less frequent population size expansions towards higher latitudes might reflect range
80
81
Figure 5.5. Plot of demographic change as a function of time for ancestral and extant
species. The first upper row shows four extant species with fast demographic change. The
second row includes four extant species with near constant historical demography. The two
last rows show eight ancestral species with relatively fast demographic change. In these the
numbers represent nodes shown in Fig. 5.3.
shifts in species that move northward as glaciers retreated, without changes in population size. In
contrast, broader niches at mid-latitudes during glacial cycles (Fig. 5.2) may allow in part for
populations to both expand and conserve allele diversity (Fig. 5.4A,C). Additionally, we find that
the range shift extent is significantly correlated with population expansions (Fig. 5.4C). This
means that faster and farther movements northward or southward away from the ancestral range
are likely to influence population size changes in this group of fungi. This is a common
prediction of population genetic dynamics in response to Quaternary climate change (Hewitt
1996, 2000, 2004a,b).
Outliers in the data are likely to be driven by intrinsic environmental limitations leading to
species-specific local adaptation and colonization responses (Savolainen et al. 2007). To
exemplify, A. jacksonii has a broad range extending to the northern distribution limit (Fig. 5.2D),
and a demographic history that suggests a fast population expansion (Fig. 5.5). These
characteristics suggest that this species is a hardy colonizer. In contrast, other species with
similar latitudes, such as A. sp-W15, have narrower ranges (Fig. 5.2C,D) and lower and more
constant population sizes (Fig. 5.5). Similarly, A. sp-jack3, A. sp. ‘affin. hayalyuy’, and A. sp
‘affin. basii’ occur within the same latitude as A. basii, but do not show signs of population
expansions and have relatively smaller population size (Fig. 5.4A,B, 6). In widespread species,
only two (A. sp. ‘cochiseana’ and A. sp-jack2) species show signs of population expansions (Fig.
5.5), suggesting that they recolonized northern habitats from southern ‘source’ populations. In
both cases, genetic diversity statistics (Table 5.2) and population sizes (Fig. 5.4B) are not
comparatively low, which may suggest that population bottlenecks were not severe. According to
metapopulation theory, metapopulations that have sufficiently large ‘source’ populations and
high dispersal capabilities are likely to overcome high rates of extinction (Levin 1995), which
were likely during glacial cycles.
82
Several species (A. basii, A. sp-jack2, A. sp-jack5, A. sp-jack6) experienced population
expansions exclusively within refugial areas (Fig. 5.4A,B, 5), all of them occurring in Mexico
(Fig. 5.2C,D). This pattern may be attributed to the presence of microrefugia in this region (e.g.
Ornelas et al. 2010; Gutiérrez-Rodríguez et al. 2011; Ruiz-Sánchez and Ornelas 2014), coupled
with broader environmental niches during glacial periods (Fig. 5.2B), which may have facilitated
population expansions (Mastretta-Yanes et al. 2015). Moreover, I detected eight ancestral
demographic expansions, three occurring in widespread areas between the Miocene and
Pliocene, and three more occurring later within refugial areas during the Pleistocene (Fig. 5.3-5).
The former results are consistent with NA colonization patterns in Amanita sect. Caesareae
(Sánchez-Ramírez et al. 2015a), supporting the view of population expansion following the
colonization of a new region (Sakai et al. 2001); while the later events may also be attributed to
glacial-interglacial expansions. I find a higher number of population expansions occurring during
the Pleistocene (Fig. 5.5), which advocates to the importance of this epoch in shaping diversity
patterns in the NA Caesar’s mushrooms.
5.4.3 Potential caveats
Some results shown here should be interpreted with caution. (1) Population size (! ! !!!!)
estimates based on the coalescent often have several assumptions –for instance non-overlapping
generations, neutral evolution, infinite-site mutation models, random mating, and free
recombination between, but not within, loci–, which can be easily violated leading to biased
results. It has been shown that increasing the number of independent loci can improve the
accuracy of population size estimation, however, a large number (> 16) of unlink genetic loci is
often recommended (Heled and Drummond 2008). This is probably one of the reasons why
posterior population size estimations have large confidence intervals for some species (Fig. 4).
(2) The geographic limits of some of the newly discovered species are poorly known. This means
that their climatic demands are also poorly understood. Although I strive to include as much data
and as much uncertainty as possible for characterizing species distributions, we might later
discover, with more data, novel distributions or range shifts that might alter the discussed
interpretations. Finally, (4) certain flaws in BiSSE-type diversification models have been
recently pointed out. In particular, I am concerned about phylogenetic sudo-replication, the effect
of neutral traits, and low power in small-sized phylogenies (Maddison and FitzJohn 2015;
Rabosky and Goldberg 2015; Davis et al. 2013).
83
5.5 Conclusions
By applying multi-locus coalescent-based species tree models and broad geographic sampling, I
have performed a thorough assessment of the diversity in the Caesar’s mushrooms (A. affin.
caesarea) species complex in North America. With implementations of the multi-species
coalescent model in BP&P I delimited 19 species and inferred their divergence times and
historical population sizes across multiple single-copy genes. Results uncovered unexpected
levels of cryptic diversity (68%). High speciation rate in refugial areas and more abundant
population expansions during the Pleistocene are indicators of intricate dynamics during glacial
cycles, which evidence their role in explaining current and historical diversity patterns. Complex
population dynamics and high genetic diversity at mid-range (20-35 degrees latitude)
distributions point out Mexican highlands as a diversity pump. Also, I show that abrupt range
shifts are followed by population size expansions, which is expected under Quaternary climate
change. Finally, this study highlights the potential of the multi-species coalescent model in
comparative phylogeographic analyses and diversity assessments that include ancestral species.
84
Chapter 6
COMPARATIVE POPULATION GENOMICS AND MOLECULAR 6ADAPTATION IN THE AMANITA JACKSONII COMPLEX
6.1 Introduction
Understanding the nature of genomic variation can help illustrate how natural selection and
demography balance out throughout the history of a species. The neutral theory of molecular
evolution (Kimura 1986; Ohta 1992) postulates that most mutations contributing substantially to
the genetic pool of a population will be either neutral or nearly so. Behind this reasoning is the
assumption that changes that alter the amino acid composition of a protein will tend to be
deleterious and selected against, while those that are advantageous will be rare. In this
framework, one can test whether DNA variation at a given locus fits the expectations under the
neutral equilibrium model (Tajima 1989; Fu and Li 1993, 1999). These models primarily allow
for explicit testing of neutrality departures, however different kinds of evolutionary processes,
such as selective sweeps, recombination, and demographic changes, can mimic similar results
(Bachtrog and Andolfatto 2006; Ramírez-Soriano et al. 2008). Nonetheless, selection is likely to
affect specific regions through the fixation of linked neutral variation –a process known as
hitchhiking or selective sweeps–, while demographic changes are more likely to have genome-
wide effects (Andolfatto 2001; Nielsen 2001).
Natural selection at the molecular level is detectable in protein-coding DNA. This can be
accomplished by looking at the relationship between functional (amino acid-changing or non-
synonymous) and silent (neutral or synonymous) variation. A widely use test of selection is the
McDonald-Kreitman test (MKT, McDonald and Kreitman 1991), which is based on the
comparison between the ratio of functional !!!! and silent !!!! mutations within a species
(polymorphism) and the ratio of fixed !! and !! mutations between species (divergence). If both
polymorphism and divergence were solely driven by mutation and drift, one would expect their
ratio to equal one. This model assumes that all !! mutations are neutral and that !! mutations
are either strongly deleterious, neutral or strongly advantageous (Smith and Eyre-Walker 2002).
In this framework, advantageous mutations are more likely to contribute to divergence due to
85
their high fixation probability, while deleterious mutations are likely to be found at lower
frequencies due to the effects of purifying selection. On the other hand, the vast majority of
segragating mutations are likely to be neutral, as they will tend to persist in the population for
longer periods until they eventually go to fixation by random genetic drift (Sella et al. 2009).
However, the proportion of adaptive mutations can be under- or overestimated under certain
demographic scenarios combined with weak selection (Smith and Eyre-Walker 2002; Eyre-
Walker 2006; Messer and Petrov 2013a). For instance, in species with reduced population sizes,
slightly deleterious mutations are effectively neutral, thus they are at higher frequency in the
population resulting in an underestimation of the proportion of adaptive mutations. Similarly, a
recent population expansion will tend to decrease polymorphism, confounding sweep-driven
positive selection (Nielsen 2001; Nielsen et al. 2005; Wright and Gaut 2005). Furthermore,
another approach is to quantify selection directly as the ratio of substitution rates at !! and !!
sites (Kimura 1977; Huges and Nei 1988; Goldman and Yang 1994). In a neutral scenario,
!!!!! should approximate 1, with values above indicating diversifying (positive) selection, and
values below indicating purifying (negative) selection (Nielsen 1997; Yang and Nielsen 2000;
Nielsen 2001).
Detecting genome-wide patterns of selection and demography is now feasible through the
growing availability and accessibility of whole-genome data; the field now known as population
genomics (Charlesworth 2010). Most studies addressing questions such as genome-wide effects
of selection have usually focus on a number of model organisms (Hough et al. 2013; Ellegren
2014), most of which have good reference genomes and well-documented annotations, while for
the vast majority of species, empirical evidence is still lacking. Nonetheless, lowering
sequencing costs and more efficient bioinformatics tools are making these types of approaches
accessible to non-model organisms (Ekblom and Galindo 2011; Aguileta et al. 2010).
Fungi play multiple essential roles in the environment, mainly as decomposers of organic matter,
but also as pathogens, commensals, and mutualists. Not only due to their diverse ecology
(Selbmann et al. 2013), but also because of their genetic machinery (Anderson et al. 1992;
Schoustra et al. 2007), fungi are regarded as highly adaptable organisms. Most evidence comes
from experimental studies in yeasts (Suutari et al. 1990; Davies et al. 1995; Piper et al. 2001; Liu
2006; Dettman et al. 2007; Anderson et al. 2010; Gerstein et al. 2014), drug resistance
(Kontoyiannis and Lewis 2002; Anderson 2005), and the evolution of virulence in pathogens
86
(Bentrup and Russell 2001; Becher et al. 2010; Fisher et al. 2012). In spite of being great
candidates (mostly due to relatively smaller genomes and short generation times) to study
evolutionary genomics (Gladieux et al. 2014), adaptation in wild fungal populations and species
has been scarcely explored. However, recent work in few fungal species has evidenced adaptive
mechanisms molecularly and experimentally in the wild. For instance in Neurospora,
environmental factors such as temperature and latitude are contributing components to genome
divergence and adaptation in wild populations (Ellison et al. 2011). Moreover, Aguileta et al.
(2012) and Gladieux et al. (2013) highlight the importance of selection and molecular adaptation
during host switches and host specialization in Botrys and Microbotryum. Finally, Branco et al.
(2015) explores patterns of genome divergence in two populations of the ectomycorrhizal (EM)
–root-associated symbionts– fungus Suillus adapted coastal and montane environments. These
studies are few recent examples that evidence genome-wide adaptive mechanisms in wild fungal
populations.
The genus Amanita comprises 500-1000 species (Tulloss 2005), most of which live as EM
symbionts (Wolfe et al. 2012). The Caesar’s mushrooms (sect. Caesareae) are a clade of edible
EM Amanita distributed worldwide. Two studies by Sánchez-Ramírez et al. (2015a, accepted)
have shown that North America harbors high diversity in this group, which has derived, in part,
from recent continental radiations during the Plio-Pleistocene. With its eight species, the A.
jacksonii complex is an example of one of such radiations (Sánchez-Ramírez et al. 2015a,
accepted). I hypothesize that climatic dynamics during the Pleistocene may be responsible for
population size fluctuations and local adaptation to glacial refugia. Here, I explore the relative
contributions of natural selection and demography in explaining genome-wide diversity patterns
in the recently diverged A. jacksonii complex.
6.2 Material and Methods
6.2.1 Samples and genomic data
Forty-six dried specimens, representing eight species (Sánchez-Ramírez et al. accepted; Chapter
5) were processed from collections made in various locations in North America, all deposited in
fungal herbaria (Appendix Table D8.4). DNA was extracted from c. 30-50 mg of dried fungal
tissue (gill) and extracted using a modified standard CTAB/proteinase K/chorophrom:isoamilic
alcohol protocol (van der Nest et al. 2014; Chapter 4). In addition, whole-genome 454 sequences
87
were produced by the Duke Center for Genomic and Computational Biology for a North
American species outside the A. jacksonii complex –A. basii (Chapter 4). These sequences were
mapped to the A. jacksonii TRTC168611 draft genome (van der Nest et al. 2014), with the
purpose of building a whole-genome draft assembly of an outgroup species.
6.2.2 Identification of single-copy orthologs, probe design, and gene selection
Structural genes were predicted in A. jacksonii and A. basii using AUGUSTUS v3.0.3 (Stanke et al.
2006). All CDS in each species were reciprocally aligned against each other using command-line
(standalone) BLASTn, creating two files in tabular format (-outfmt 7). Next, a custom perl script
(rbh.pl, available at: https://sites.google.com/site/santiagosnchezrmirez/home/software/perl) was
used to find reciprocal single-copy hits, putting each pair in a single FASTA file. Each single-
copy pair was aligned using MUSCLE v3.6 (Edgar 2004). Probes for target hybridization were
designed based on conserved (identical) CDS regions, 60bp long, across both species. Only
genes with four or more non-overlapping regions were chosen. Two PCR primers (5’-
TAATACGACTCACTATAGGG-3’ and 5’-CTATAGTGTCACCTAAATC-3’) were added to
the 5’ and 3’ ends of the 60bp probe target, for a total of 100bp per oligonucleotide probe.
6.2.3 Sequence capture, library preparation, and sequencing
Custom oligos probes were synthesized by LCScience
(http://www.lcsciences.com/applications/genomics/oligomix/) using DNA microchip technology
(Gao et al. 2004) The synthesized probes were then PCR amplified as needed for downstream
applications. Probes were later transcribed into RNA using MegaScript T7 Kit (Invitrogen,
Carlsbad, CA) and Biotin-dUTP to produce labeled oligos. The DNA was sheared and a standard
Illumina TruSeq V2 DNA library kit was prepared, multiplexing all 46 samples. Sequence
captures (probe-target hybridization) were performed with 200ng of library, 500ng of
biotinylated probes, and Steptavidin beads (Invitrogen, Dynabeads M-270, CAT65305). The
capture procedure was repeated at least twice. The details of the protocols can be found in
Supplementary file 1. All target captures were directly sequenced in a single Illumina Hi-Seq
2500 lane.
88
6.2.4 Bioinformatics
Read alignments, sequence phasing, and quality control. Sequences were de-multiplexed into
individual FASTQ files and each one mapped onto the targets’ whole-gene sequences from A.
jacksonii using the Burrows-Wheeler Aligner (BWA v0.7.12) algorithm BWA-MEM (Li and
Durbin 2010). The resulting SAM files were then filtered with view in SAMtools v0.1.19 (Li et
al. 2009) for off target reads, low quality mapping (< 30), and PCR duplicates, and converted to
BAM files. The reference sequences were indexed and the BAM files sorted. I used
HAPCOMPASS (Aguiar and Istrail 2012), which uses a read-based graph algorithm, to generate
haploid sequences from BAM and VCF files. VCF files were produced by “piping” mpileup,
bcftools’ view, calling the “-cg” flags for genotype likelihood computing, and varFilter in
vcfutils.pl, only keeping SNPs with Fred quality ! 20. The SAM files produced by HapCompass
were then converted to BAM and sorted. Consensus FASTQ files were produced by “piping”
mpileup, bcftools’ view, and vcf2fq in vcfutils.pl, in a similar way as above. A custom perl script
(OrderFromSamtools.pl) was used to convert FASTQ files to FASTA format and to produce by-
gene alignments from by-sample files. A second perl script (getPhasedGenome.pl) parsed those
FASTA files to produce “a” and “b” allelic variants per haploid sequence. All missing data was
scored as Ns, while unsorted polymorphisms were either left as IUPAC ambiguity codes or
scored as Ns, depending on the analysis. Scripts are available at:
carbohydrate metabolic process [GO:0005975], hydrolase activity [GO:0016787]), and
homeostasis (calcium ion binding [GO:0005509]). Species-specific genes under selection varied
in their functional content and number (Appendix Fig. D8.14). A. sp-F11, confined to a unique
and humid environment (Florida), had comparatively more species-specific genes under positive
selection that also had unique GO enrichments (Appendix Fig. D8.14).
6.4 Discussion
Interpreting DNA variation from an evolutionary perspective has been a long-standing goal for
population geneticists, both doing empirical and theoretical work. The neutral equilibrium theory
105
Figure 6.6. Gene Ontology terms Wordles for genes under positive selection (A: ! > 0; B: M2a) combined for all species. The size of the term includes the extent of selection and how often the term was found across genes and species.
A
B
106
of molecular evolution (Kimura 1977, 1986) provides a null framework to test hypotheses about
natural selection, and other processes, including demography, acting on the genome. While the
core of population genomic research has develop around coalescent-based theory in a single –or
at least a pair– of species (Hough et al. 2013; Charlesworth 2010), there are clear benefits in
integrating comparative and phylogenetic data, in addition to divergence assessments, in multi-
species multi-locus genomic analyses (Cutter 2013). This study tackles common questions in
population genomics (i.e. how does selection and population size shape genome-wide
diversity?), but taking advantage of exhaustive species-level sampling, phylogeny, genealogy,
and time-scaled demographic inference.
6.4.1 The exon-targeted sequence-capture approach
Before and early into the “whole-genome sequencing” era, most population genomics studies
relied on PCR and Sanger sequencing of individual genes for data production. Some classic
studies included data ranging from tens to few hundreds of genes, in what were then called
“whole-genome” approaches. Exon-targeted sequencing takes advantage of the rapid and
massive data production of next-generation sequencing, while focusing on specific genes of
interest. Although any class of gene can be potentially targeted, I opted for single-copy genes,
attempting to avoid potential spurious DNA hybridizations with paralogues, which can bias
diversity estimates. However, targeting single-copy genes has the disadvantage of characterizing
genome-wide diversity based on a single “class” of genes, which may not necessarily represent a
heterogeneous sample from the genome. Nonetheless, the degree of variation and distribution of
the data within each species, in combination with my results, suggests that this sample of 502
genes maybe a good representation genome-wide variation. Furthermore, this method has been
proven successful when DNA material is only available in low quantities (Hancock-Hanser et al.
2013), or from degraded samples (Templeton et al. 2013). All of my samples came from dried
museum specimens and FTA plant saver cards (Dentinger et al. 2010), which would not reach
quantity/quality standards for other popular downstream genomic applications such as whole-
genome resequencing or reduced-library representation approaches (e.g. RAD-tag, GBS). In
particular, reduced-library representation approaches are practical for sampling genome-wide
SNPs, however the shortness of the fragments makes it impossible to conduct codon-based
divergence and polymorphism assessments. Therefore, exon-targeting methods are a practical
107
and efficient strategy for scoring many loci and complete structural genes in non-model
organisms.
6.4.2 Effects of selection and demography on genomic variation
One of the most striking results is the consistent excess of polymorphism compared to
divergence (Fig. 6.2), which translates into reduced adaptive amino acid substitutions. Studies in
flies (Smith and Eyre-Walker 2002; Bierne and Eyre-Walker 2004; Sella et al. 2009), plants
(Slotte et al. 2010), and humans (Fay et al. 2001) have suggested that the rate of adaptive
substitution can reach high proportions (20–40%), even in the presence of high levels of
selective constrain or purifying selection (Williamson et al. 2014; Slotte et al. 2010). The rate of
adaptive evolution is thought to be even higher in microorganisms, such as bacteria
(Charlesworth and Eyre-Walker 2006). In contrast, in these Amanita species, both MKT and
!!!!! likelihood-ratio tests indicate that the majority of genes approximate neutrality, with a
skew towards a negative alpha, and an overall small fraction of genes under adaptive evolution
(Fig. 6.2). Genome-wide diversity can arise from different evolutionary processes. Among them,
I will discuss the relative contributions of selection, population size changes, and gene flow in
explaining diversity patterns within and among members of the species complex.
Strong selective constraint is expected in protein coding DNA, largely because the majority of
amino acid changing mutations are likely to be deleterious (Ohta 1992; Andolfatto 2005). Both
polymorphism (Table 6.1, Appendix Fig. D8.8-9) and divergence (Fig. 6.2) data support overall
high levels of selective constraint in the A. jacksonii complex. One possible explanation for
values of ! ! ! and !! ! !! is purifying selection. If non-synonymous substitutions are not
being fixed between species, it is likely that they were deleterious early in the evolution of a
lineage, and are thus being purged from the population (Fay and Wu 2003). Mixed genome-wide
signals of positive and negative selection are not uncommon, and have been found in other
organisms such as plants (Williamson et al. 2014) and humans (Fay et al. 2001). An excess of
polymorphism relative to divergence can also be caused either by balancing selection (Verrelli et
al. 2002), a prevalence of soft selective sweeps (Messer and Petrov 2013b), or local adaptation in
subpopulations (Turner et al. 2010). In either case I would observe a high ratio of non-
synonymous to synonymous polymorphism, which would cause both ! ! ! and !! ! !!.
Genes under long-term balancing selection have longer than expected coalescence times due to
108
the retention of ancestral polymorphism (Edwards et al. 1997; Charlesworth 2009). Moreover,
high ! across multiple species is an indicator of positive selection on allelic diversity (Yang et
al. 2005; Aguilar and Garza 2007). Results show little support for long-term diversifying
selection. Only 10% of genes had a significant site class with ! ! !. In addition, genes with
strongly negative ! were not significantly affected by ancestral polymorphism (i.e. low gsi
scores; Appendix Fig. D8.12).
The fate of adaptive and slightly deleterious mutations is highly affected by !! (Ohta 1992;
Eyre-Walker and Keightley 2007; Charlesworth 2009). In small populations, adaptive mutations
can get lost by drift effects alone, while slightly deleterious mutations are effectively neutral and
can increase in frequency within the population. MKT-based approaches are generally
insensitive to the presence of slightly deleterious mutations (although see Eyre-Walker and
Keightley 2007 and Messer and Petrov 2013a). However, estimates of the effective population
size from silent diversity and the inference of demographic changes can help clarify if slightly
deleterious mutations might be common given a demographic history. Results from temporal
demographic models in BEAST indicate recent sharp increases in population size, particularly in
A. jacksonii, A. sp-jack2, A. sp-jack6, and A. sp-jack5 (Fig. 6.5). These changes are consistent
with skewed genome-wide frequency spectra (!Taj ! !, Fig. 6.1), also reflecting the steepness of
the population size change. From estimates of !! (Appendix Fig. D8.10) and demographic trends
(Fig. 6.5), I observe that A. sp-F11 has the largest !! and a relatively stable demography.
Amanita sp-F11 is likely to have had a stable population size throughout the glaciations confined
to a known Pleistocene refugium in the southeastern US (Fig. 6.7) (Soltis et al. 2006; Sánchez-
Ramírez et al. accepted). This explains not only the differences in !! and diversity between A.
sp-F11 and the rest of the complex, but it also allows us to use A. sp-F11 as reference species
with a constant population size. By doing so, I show that compared to A. sp-F11, species that
underwent steep demographic expansions (A. jacksonii, A. sp-jack2, A. sp-jack6, A. sp-jack5)
have lower !! (Appendix Fig. D8.10) in addition to an excess of low frequency variants (Fig.
6.1). When !! is low, slightly deleterious mutations are effectively neutral (Ohta 1992), and can
rise in frequency with demographic expansion (Nielsen 2001). This same situation can be
applied to species that did not undergo demographic expansion, but have low historical
population sizes such as A. sp-jack3, A sp-jack1, and A. sp-T31 (Fig. 6.5). Abundance of slightly
deleterious mutations can certainly underestimate MKT-based estimates of ! (Smith and Eyre-
109
Walker 2002). Since demography is more likely to affect the whole genome compared to soft
sweeps and local adaptation, I prefer the explanation that low !! caused by population
bottlenecks during the Pleistocene glaciations, followed by recent expansions, has recently
increased the frequency of slightly deleterious/effectively neutral mutations causing an excess of
polymorphism compared to divergence. This signal can be amplified by purifying selection,
further increasing the gap between divergence and polymorphism.
Adaptive responses can follow rapid environmental change or the colonization of a new
environment (Hermisson and Pennings 2005; Barrett and Schluter 2007; Hoffmann and Sgrò
2011). The Pleistocene is characterized by a series of abrupt temperature oscillations that have
deeply affected the genetic diversity of populations and species (Hewitt 1996, 2000, 2004a,b). If
I consider a situation where different alleles are selected during glacial and inter-glacial cycles, I
might expect to find a high degree of polymorphism. This is in fact a prediction during soft
selective sweeps under recurrent !! bottlenecks (Messer and Petrov 2013b). Although
demographic analysis did not support recurrent fluctuations of population size (Fig. 6.5), it is not
unlikely that these species might have suffered from recurrent fluctuations in population size
during glacial cycles. Under this view, I do not discard the effects of soft selective sweeps as a
potential cause for increased polymorphism.
Although I did not quantify gene flow directly, per gene gsi values provide an overview of the
extent of intraspecific gene sharing. Species that are more evolutionary distant (A. sp-F11, A. sp-
T31, A. sp-jack1) tend to share less alleles among species, than those that are more closely
related (A. jacksonii, A. sp-jack2, A. sp-jack3, A. sp-jack5, A. sp-jack6) (Figs. 6.4 and 6.5). In
particular, substantial allele sharing was observed (data not shown) between geographically close
species pairs (Fig. 6.7; A. jacksonii and A. sp-jack2; and A. sp-jack5 and A. sp-jack6). To less
extent A. sp-jack3 had some allele sharing with A. sp-jack2 and A. sp-jack5, which are also
geographically close, however, had more exclusive alleles. This is consistent with the view of A.
sp-jack3 being more range restricted and having a relatively small !!. Furthermore, it has been
shown that in subdivided populations different sampling strategies a can affect frequency-based
diversity estimators, such as ! and !Taj, to various degrees in the presence of migration (Städler
et al. 2009; Cutter et al. 2012). In particular, if high migration will tend to homogenize the
coalescence of samples making it converge to a similar !Taj value regardless of the sampling
110
Figure 6.7. Species distribution maps based on geographic data from Sánchez-Ramírez et
al. (accepted; Chapter 5) and MaxEnt (Phillips et al. 2006). The dotted line marks an
hypothesized limit for glacial refugia.
strategy. In contrast, low migration will tend to polarize and bias !Taj values towards high or low
values depending on the sampling (Cutter et al. 2012). I show that !Taj is not overall affected by
genes with low gsi. This means that the skew in !Taj estimates is not directly influenced by allele
111
sharing between species. Although the sampling strategy might still be a potential factor, I not
aware of any population structure within species that might bias within species diversity
estimates.
6.4.3 Functionality of adaptive genes and the EM habit
The genes that were assayed here were not selected for their specific functions. However,
comparing functional properties across species, among those genes under a given selective force,
may shed light on potential shared and local adaptive responses. Two of the most prevalent GO
terms with genes under positive directional selection (i.e. ! ! !) and positive diversifying
selection (! ! !) reflect ATP and protein binding functions. Although these terms are too
generic, probably due to the still insufficient characterization of fungal proteins, they suggest that
many of these genes are involved with energetic and metabolic cellular processes, as well as
protein-protein interactions. Overall, protein functions seem to be divided into three general
categories; those related to (1) gene transcription/expression activities, (2) signalling and
transport, and (3) metabolic oxidative processes (Fig. 6.6). Many of these functions coincide
with divergent regions in the genome of another EM fungus (S. brevipes) that are also thought to
be under selection (Branco et al. 2015). I also noted some functions related to carbohydrate
metabolism, specifically with catalytic and hydrolytic activities (Fig. 6.6, Appendix Fig. D8.14).
These may be compelling because all EM fungi depend on external carbohydrate sources that are
generally supplied by their hosts (Nehls et al. 2007). Throughout the evolution of the EM
symbiosis, many species have suffered from a convergent loss of decay mechanisms, while
retaining moderate gene repertoires with lignocellulitic capabilities (Kohler et al. 2015).
Adaptation in proteins related to carbohydrate metabolism, such as hydrolase activity, could
allow for a more efficient breakdown of complex polysaccharides during periods of stress and/or
low photosynthetic activity. Also worth mentioning, proteins with telomeric template reverse
transcriptase activity were found under diversifying selection (Fig. 6.6B). Studies suggest that
the dynamic nature of telomeric and subtelomeric regions are likely to contribute to rapid
evolution of fungi during biotic interactions (Aguileta et al. 2009).
I find that it is more common to find genes under adaptive divergence within species than genes
under diversifying selection across all species. This implies that different genes within species
have undergone adaptive evolution more commonly than a single gene across species. One of the
112
species that had the most exclusive set of genes under positive selection was A. sp-F11. To my
knowledge (Sánchez-Ramírez et al. 2015a, accepted) this species is restricted to a unique
environment in semi-tropical temperate forests in Florida (USA), with virtually no overlap with
other species in the group (Fig. 6.7). This observation could be an indication of local adaptation
to a Plio-Pleistocene refugium environment in this species.
6.5 Conclusions
I produced DNA sequence data for 502 structurally complete genes across multiple samples in
looks into genome-wide diversity patterns within a few closely related species. The “zooming” is
both at temporal and spatial scales. To my knowledge, no other study focusing on fungi has
investigated macro- and microevolutionary patterns throughout time, space, and genomes.
Figure 7.1. Temporal and spatial scales explored within the Caesar’s mushrooms.
0102030405060 0102030405060 01020 05
!"#$%&'() !"#$%&'(* !"#$%&'(+ !"#$%&'(,
Ma
114
Macroevolutionary investigations are generally complemented and strengthened by incorporating
data from fossil records. Many animals and plants have good fossil records that can make
macroevolutionary and biogeographic hypotheses robust. However, most fungal tissue decays
quickly and lacks rigid substances, making them poor fossilizers. By using time-scaled
phylogenies –from external calibrations with some of the few fungal fossils (Hibbett et al. 2007;
Smith et al. 2004)– I was able to infer ancestral distribution and dispersal patterns in the Caesar’s
mushrooms. Some of the results shown here may have implications for geographic and
evolutionary histories in other fungi and other organisms in general. For instance, botanists
studying biogeographic patterns of ectotrophic trees might find interesting that potential
symbiotic partners may have been found in Africa as early as the Eocene. Also that “Out-of-
Africa” and “Out-of-the-Paleotropics” dispersal patterns of EM species may have assisted
Boreotropical distributions during the Tertiary (Chapter 2). This pattern has also been suggested
for other groups of ectomycorrhizal fungi (Matheny et al. 2009; Du et al. 2012; Wilson et al.
2012). Moreover, the relationship between East Asia and North American floras is of interest to
many botanists (Quian and Ricklefs 2000; Xiang et al. 2000; Smith and Donoghue 2004).
Chapter 2 also documents that recurrent dispersal between these two continents was active
during the late Miocene and Pliocene, and probably assisted and/or followed established
ectotrophic tree communities.
Some attention has been given recently to global diversity patterns in fungi (Tedersoo et al.
2012, 2014). An interesting pattern that has arisen is that the diversity of EM fungi is higher in
subtropical and temperate regions, rather than following the general LDG pattern observed in
most other group of organisms. Although several hypotheses have been proposed (Tedersoo and
Nara 2010), I show that, at least in the Caesar’s mushrooms, this diversity bias might be a result
of high rates of speciation at temperate latitudes (Chapter 3). Evidently, more data from other
EM taxa/clades is needed to give more support to this pattern at a larger scale. Likewise, future
studies would benefit from the integration of fine-scale ecological and climatical data to larger
and better sampled phylogenies.
Genome data is part of the 21st century evolutionary biology (Losos et al. 2013). Comparative
genomics and population genomics provide the means to look deeper and more closely into the
evolutionary history of species and populations. In Chapter 4, I produced, assembled, and
annotated the genome of a strain of A. jacksonii, and obtained partial genomic sequences from a
115
strain of a closely related species, A. basii. I used these genomic data to developed novel single-
copy gene markers for multilocus phylogenetic analyzes (Chapter 5) and genome-wide diversity
exploration (Chapter 6). This is the first genome sequenced from an edible species of Amanita
and may provide insights later in the exploration of mushroom toxicity/edibility, and
complement our understanding about the evolution of ectomycorrhizal symbiosis (e.g. Kohler et
al. 2015) in the genus.
Fungal phylogeography is still in its early stages (Lumbsch et al. 2008). Part of the limitations
involve the lack of sufficient geographic sampling (Schmit and Mueller 2007), the failure to
recognize distinct evolutionary units (Taylor et al. 2000, 2006), and difficulties to characterize
diversity patterns in temporal and spatial contexts. Chapter 5 seeks to overcome some of these
limitations in the limitations of diversity in the North American Caesar’s mushrooms. As already
found in other fungi (Crespo and Lumbsch 2010), I uncover high levels of cryptic diversity in
North America. Most of the cryptic diversity occurs within refugia, probably as a result of high
speciation rates during glacial oscillations. Chapter 5 also reveals many population expansions
during the Pleistocene, most of which occur within refugial latitudes. This study evidences
responses to past climate change in fungi, which have been rarely documented.
Genome-wide DNA data can provide insights into demographic and adaptive processes that
shape diversity. In Chapter 6, a sample of 502 putative single-copy genes revealed differences in
DNA variation and frequency spectra among eight closely related species in the A. jacksonii
complex. Models for detecting selection and demographic trends indicate a strong influence of
negative selection and recent expansion events during the Pleistocene. In this Chapter, I also
identify potential genes involved in adaptive processes. Some of these may prove to be part of
fine-level symbiotic relationships or abiotic conditions that allow these fungi to thrive in the
North American temperate ecosystem.
Many new and interesting questions also emerged from this thesis. For instance, compared to
Asia and America, Europe apparently holds low diversity, in spite of being mostly temperate. Is
there really low diversity? Or is it a matter of sampling? To my knowledge there are many
geographic populations of “A. caesarea” in Europe and adjacent areas. For instance, there are
records of this species in highland forests in Iran, Georgia, Turkey, the Balkans, the Carpatians,
Italy, France, Spain, and Algeria, in addition a few small populations in Switzerland that are
116
protected. It would be interesting to know if these, in fact, represent a single or multiple cryptic
species, and if they diverged as a result of relictual isolation during the Quaternary. In North
America, A. basii and A. cochiseana are two sister species with a large effective population size
and a recent history of divergence. Although A. basii is generally found in southern Mexico, I
was able to identify one specimen of A. basii in New Mexico, U.S.A. With better geographic and
genomic sampling it would be interesting to see if this species pair forms a Mesoamerican
“species ring” fueled by Pleistocene dynamics. Many questions related to “speciation genomics”
also remain unanswered. For instance, how do sympatric species develop and maintain
reproductive isolation? Amanita calyptroderma and A. vernicoccora in western North America
appear to have decoupled fruiting phenology; the former fruits in the summer while the later
fruits in the spring. But, are there any environmental or genomic cues? What about other
sympatric species in southeastern U.S.A? Continuous sampling and genotyping in poorly known
areas is likely to provide better insights. Finally, further research on A. jacksonii may also prove
to be enlightening for studying the genomics of demographic expansions, selection, local
adaptation, and recombination in populations of ectomycorrhizal fungi.
117
References
Aanensenn D, Baguelinn M, Birrell P, Cauchemez S, Camacho A, Colijn C, Cori A, Didelot X, Eames K, Fraser C, Frost S, Hens N, Hugues J, Jombart T, Opatowski L, Ratmann O, Soubeyrand S, Suchard M, Wallinga J, Ypma R. 2015. OutBreakTools: Basic Tools for the Analysis of Disease Outbreaks. The Comprehensive R Archive Network. http://cran.r-project.org/web/packages/OutbreakTools/index.html
Agapow P-M, Bininda-Emonds ORP, Crandall KA et al. 2004. The impact of species concept on biodiversity studies. Quart. Rev. Biol. 79:161–179.
Aguiar D, Istrail S. 2012. HapCompass: a fast cycle basis algorithm for accurate haplotype assembly of sequence data. J. Comp. Biol. 19:577–590.
Aguilar A, Garza JC. 2007. Patterns of historical balancing selection on the salmonid major histocompatibility complex class II ! gene. J. Mol. Evol. 65:34–43.
Aguileta G, Hood M, Refrégier Ga, Giraud T. 2009. Genome Evolution in plant pathogenic and symbiotic fungi. Adv. Bot. Res. 49:151–193.
Aguileta G, Lengelle J, Chiapello H, Giraud T, Viaud M, Fournier E, Rodolphe F, Marthey S, Ducasse A, Gendrault A, Poulain J, Wincker P, Gout L. 2012. Genes under positive selection in a model plant pathogenic fungus, Botrytis. Infect. Genet. Evol. 12:987–996.
Aguileta G, Lengelle J, Marthey S, Chiapello H, Rodolphe F, Gendrault A, Yockteng R, Vercken E, Devier B, Fontaine MC, Wincker P, Dossat C, Cruaud C, Couloux A, Giraud T. 2010. Finding candidate genes under positive selection in non-model species: examples of genes involved in host specialization in pathogens. Mol. Ecol. 19:292–306.
Alexander IJ. 2006. Ectomycorrhizas – out of Africa? New Phytol. 172:589-591. Allen MB, Armstrong HA. 2008. Arabia–Eurasia collision and the forcing of mid-Cenozoic
global cooling. Palaeogeogr. Palaeoclimatol. Palaeoecol. 26:52–58. Alroy J. 2001. A multispecies overkill simulation of the end-Pleistocene megafaunal mass
extinction. Science 292:1893–1896. Alroy J. 2009. Speciation and extinction in the fossil record of North American mammals. In: R.
K. Butlin, J. R. Bridle and D. Schluter (eds.) Speciation and Patterns of Diversity, 301-323 pp. Cambridge University Press, Cambridge.
Álvarez-Lao DJ, García N. 2011. Geographical distribution of Pleistocene cold-adapted large mammal faunas in the Iberian Peninsula. Quatern. Int. 233:159–170.
Anderson D, Goudie A, Parker A. 2013. Global Environments Through the Quaternary. Oxford University Press, USA.
Anderson JB, Funt J, Thompson DA, Prabhu S, Socha A, Sirjusingh C, Dettman JR, Parreiras L, Guttman DS, Regev A, Kohn LM. 2010. Determinants of divergent adaptation and Dobzhansky-Muller interaction in experimental yeast populations. Curr. Biol. 20:1383–1388.
118
Anderson JB, Kohn LM, Leslie JF. 1992. Genetic mechanisms in fungal adaptation. In: The fungal community: its organization and role in the ecosystem. New York: Marcel Dekker. p, 73–98.
Anderson JB. 2005. Evolution of antifungal-drug resistance: mechanisms and pathogen fitness. Nat. Rev. Microbiol. 3:547–556.
Andolfatto P. 2001. Adaptive hitchhiking effects on genome variability. Curr. Opin. Genet. Dev. 11:635–641.
Andolfatto P. 2005. Adaptive evolution of non-coding DNA in Drosophila. Nature 437:1149–1152.
Anisimova M, Bielawski JP, Yang Z. 2001. Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution. Mol. Biol. Evol. 18:1585–1592.
Arnold AE. 2007. Understanding the diversity of foliar endophytic fungi: progress, challenges, and frontiers. Fungal Biol. Rev. 21:51–66.
Axelrod DI. 1966. Origin of deciduous and evergreen habits in temperate forests. Evolution 20:1!15.
Bachtrog D, Andolfatto P. 2006. Selection, recombination and demographic history in Drosophila miranda. Genetics 174:2045–2059.
Bacon CD, Baker WJ, Simmons MP. 2012. Miocene dispersal drives island radiations in the palm tribe Trachycarpeae (Arecaceae). Syst. Biol. 61:426–442.
Bahram M, Asef MR, Zarre S, Abbasi M, Reidl, S. 2006. Addition to the knowledge of Amanita (Agaricales, Plutaceae) from Iran. Rostaniha 7:107–119.
Barnosky AD, Koch PL, Feranec RS, Wing SL, Shabel AB. 2004. Assessing the causes of Late Pleistocene extinctions on the continents. Science, 306, 70–75.
Barraclough TG, Vogler AP. 2002. Recent diversification rates in North American tiger beetles estimated from a dated mtDNA phylogenetic tree. Mol. Biol. Evol., 19:1706–1716.
Barraclough T, Nee S. 2001. Phylogenetics and speciation. Trends Ecol. Evol. 16:39–399. Barrett RDH, Schluter D. 2008. Adaptation from standing genetic variation. Trends Ecol. Evol.
23:38–44.
Bas C. 1969. Morphology and subdivision of Amanita and a monograph of its section Lepidella. Persoonia 5:285–579.
Bas C. 1977. Species concepts in Amanita section Vaginatae. pp 79–103. In: The species concept in Hymenomycetes (ed. H. Clémençón). Bibliotheca Mycologica 61. J. Cramer, Vaduz, Liechtenstein.
Beatty GE, Provan JIM. 2010. Refugial persistence and postglacial recolonization of North America by the cold-tolerant herbaceous plant Orthilia secunda. Mol. Ecol. 19, 5009–5021.
Beaulieu J. M. and B. C. O’Meara. 2015. Extinction can be estimated from moderately sized molecular phylogenies. Evolution. doi: 10.1111/evo.12614.
Becher R, Hettwer U, Karlovsky P, Deising HB, Wirsel SGR. 2010. Adaptation of Fusarium graminearum to tebuconazole yielded descendants diverging for levels of fitness,
119
fungicide resistance, virulence, and mycotoxin production. Phytopathology. 100:444–453.
Beeli M. 1935) Amanita and Volvariella. Flore iconographique des champignons du Congo. Fasc. 1, pp. 11–27. Jardin Botanique de l’État, Brussels, Belgium.
Bennett KD, Tzedakis PC, Willis KJ. 1991. Quaternary refugia of north European trees. J. Biogeogr. 18:103–115.
Bentrup KZ, Russell DG. 2001. Mycobacterial persistence: adaptation to a changing environment. Trends Microbiol. 9:597–605.
Berbee ML, Taylor JW. 2010. Dating the molecular clock of fungi – how close are we? Fungal Biol. Rev. 24:1–16.
Bevan R, Lang B, Bryant D. 2005. Calculating the evolutionary rates of different genes: a fast, accurate estimator with applications to maximum likelihood phylogenetic analysis. Syst. Biol. 54:900–915.
Bibi F. 2011. Mio-Pliocene faunal exchanges and African biogeography: the record of fossil bovids. PLoS One 6:e16688.
Bickford D, Lohman DJ, Sodhi NS, Ng PK, Meier R, Winker K, Ingram KK, Das I. 2007. Cryptic species as a window on diversity and conservation. Trends Ecol. Evol. 22:148–155.
Bierne N, Eyre-Walker A. 2004. The genomic rate of adaptive amino acid substitution in Drosophila. Mol. Biol. Evol. 21:1350–1360.
Birand A, Vose A, Gavrilets S. 2012. Patterns of species ranges, speciation, and extinction. Am. Nat. 179:1–21.
Bivand RS, Pebesma E, Gómez-Rubio V. 2013. Applied spatial data analysis with R. Springer, New York.
Boa ER. 2004. Wild edible fungi: a global overview of their use and importance to people. In: Non-wood Forest Products, Vol. 17. Food and Agriculture Organization, Rome. 145 pp.
Boisvert S, Laviolette F, Corbeil J. 2010. Ray: simultaneous assembly of reads from a mix of high– throughput sequencing technologies. J. Comp. Biol. 17:1519–1533.
Braconnot P, Otto-Bliesner B, Harrison S, Joussaume S, Peterchmitt JY, Abe-Ouchi A, Crucifix M, Driesschaert E, Fichefet Th, Hewitt CD, Kageyama M, Kitoh A, Laîné A, Loutre M-F, Marti O, Merkel U, Ramstein G, Valdes P, Weber SL, Yu Y, Zhao Y. 2007. Results of PMIP2 coupled simulations of the Mid-Holocene and Last Glacial Maximum-Part 1: experiments and large-scale features. Clim. Past 3:261–277.
Branco S, Gladieux P, Ellison CE, Kuo A, LaButti K, Lipzen A, Grigoriev IV, Liao HL, Vilgalys R, Peay KG, Taylor JW, Bruns TD. 2015. Genetic isolation between two recently diverged populations of a symbiotic fungus. Mol. Ecol. 24:2747–2758.
Brown JH. 2014. Why are there so many species in the tropics? J. Biogeogr. 41:8-22. Brunsfeld SJ, Sullivan J, Soltis DE, Soltis PS. 2001. Comparative phylogeography of north-
western North America: a synthesis. In: Integrating Ecology and Evolution in a Spatial Context (eds Silvertown J, Antonovics J), pp. 319–339. Cambridge University Press, Oxford.
120
Bryson RW, Murphy RW, Lathrop A, Lazcano-Villareal D. 2011. Evolutionary drivers of phylogeographical diversity in the highlands of Mexico: a case study of the Crotalus triseriatus species group of montane rattlesnakes. J. Biogeogr. 38:697–710.
Buée M, Courty PE, Mignot D, Garbaye J. 2007. Soil niche effect on species diversity and catabolic activities in an ectomycorrhizal fungal community. Soil Biol. Biochem. 39:1947–1955.
Buermans HPJ, den Dunnen JT. 2014. Next generation sequencing technology: Advances and applications. Biochim. Biophys. Acta 1842:1932–1941.
Butler MA, King AA. 2004. Phylogenetic comparative analysis: a modeling approach for adaptive evolution. Am. Nat. 164:683–695.
Butlin R, The Marie Curie SPECIATION Network. 2012. What do we need to know about speciation? Trends Ecol. Evol. 27:27–39.
Butlin RK, Bridle JR, Schluter D. 2009. Speciation and patterns of biodiversity. 1-14 pp. In: Speciation and Patterns of Diversity (eds. RK Butlin, JR Bridle, D Schluter), Cambridge University Press, Cambridge.
Carbone I, Kohn L. 2004. Inferring process from pattern in fungal population genetics. Appl. Mycol. Biotechnol. 4:29–58.
Carriconde F, Gardes M, Jargeat P, Heilmann-Clausen J, Mouhamadou B, Gryta H. 2008. Population evidence of cryptic species and geographical structure in the cosmopolitan ectomycorrhizal fungus Tricholoma scalpturatum. Microb. Ecol. 56:513–24
Carstens BC, Knowles LL. 2007. Estimating species phylogeny from gene-tree probabilities despite incomplete lineage sorting: an example from Melanoplus grasshoppers. Syst. Biol. 56:400–411.
Charlesworth B. 2009. Fundamental concepts in genetics: effective population size and patterns of molecular evolution and variation. Nat. Rev. Genet. 10:195–205.
Charlesworth B. 2010. Molecular population genomics: a short history. Genet. Res. 92:397–411.
Charlesworth J, Eyre-Walker A. 2006. The rate of adaptive evolution in enteric bacteria. Mol. Biol. Evol. 23:1348–1356.
Churchill SE. 1998. Cold adaptation, heterochrony, and Neanderthals. Evol. Anthropol. 7:46–60. Coates AG, Obando JA. 1996. The geologic evolution of the Central American Isthmus. pp. 21–
56. Evolution and environment in tropical America (eds. J.B.C. Jackson, A.F. Budd and A.G. Coates), University of Chicago Press, Chicago.
Collar DC, JA Schulte II, BC O'Meara, JB Losos. 2010. Habitat use affects morphological diversification in dragon lizards. J. Evol. Biol. 23:1033–1049.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. 2005. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21:3674–3676.
Couvreur TLP, Pirie MD, Chatrou LW, Saunders RMK, Su YCF, Richardson JE, Erkens RHJ. 2011. Early evolutionary history of the flowering plant family Annonaceae: steady diversification and boreotropical geodispersal. J. Biogeogr. 38:664–680.
121
Coyne JA, Orr HA. 2004. Speciation. Sinauer Associates, Sunderland, MA. Crespo A, Lumbsch HT. 2010. Cryptic species in lichen-forming fungi. IMA Fungus 1:167–170.
Crisp MD, Cook LG. 2009. Explosive radiation or cryptic mass extinction? Interpreting signatures in molecular phylogenies. Evolution 63:2257–2265.
Cummings MP, Neel MC, Shaw KL. 2008. A genealogical approach to quantifying lineage divergence. Evolution 62:2411–2422.
Curtis TP, Wallbridge NC, Sloan WT. 2009. Theory, community assembly, diversity and evolution in the microbial world. 59-76 pp. In: Speciation and Patterns of Diversity (eds. R. K. Butlin, J. R. Bridle, D. Schluter), Cambridge University Press, Cambridge.
Cutter AD, Wang G-X, Ai H, Peng Y. 2012. Influence of finite-sites mutation, population subdivision and sampling schemes on patterns of nucleotide polymorphism for species with molecular hyperdiversity. Mol. Ecol. 21:1345–1359.
Cutter AD. 2013. Integrating phylogenetics, phylogeography and population genetics through genomes and evolutionary theory. Mol. Phylogenet. Evol. 69:1172–1185.
Darriba D, Taboada GL, Doallo R, Posada D. 2012. jModelTest 2: more models, new heuristics and parallel computing. Nature methods 9:772–772.
Davies JM, Lowry CV, Davies KJ. 1995. Transient adaptation to oxidative stress in yeast. Arch. Biochem. Biophys. 317:1–6.
Davis CC, Bell CD, Mathews S, Donoghue MJ. 2002. Laurasian migration explains Gondwanan disjunctions: evidence from Malpighiaceae. Proc. Natl. Acad. Sci. USA 99:6833–6837.
Davis MP, Midford PE, Maddison W. 2013. Exploring power and parameter estimation of the BiSSE method for analyzing species diversification. BMC Evol. Biol. 13:38.
de Queiroz K. 2007. Species concepts and species delimitation. Syst. Biol. 56:879–886. De-Nova JA, Medina R, Montero JC, Weeks A, Rosell JA, Olson ME, Eguiarte LE, Magallón S.
2012. Insights into the historical construction of species-rich Mesoamerican seasonally dry tropical forests: the diversification of Bursera (Burseraceae, Sapindales). New Phytol. 193:276–287.
Degnan JH, Rosenberg NA. 2006. Discordance of species trees with their most likely gene trees. PLoS Genet. 2:e68.
Degnan JH, Rosenberg NA. 2009. Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol. Evol. 24:332–340.
Dentinger BTM, Margaritescu S, Moncalvo J-M. 2010. Rapid and reliable high-throughput methods of DNA extraction for use in barcoding and molecular systematics of mushrooms. Mol. Ecol. Res. 10:628–633.
Dettman JR, Anderson JB, Kohn LM. 2010. Genome!wide investigation of reproductive isolation in experimental lineages and natural species of Neurospora: identifying candidate regions by microarray!based genotyping and mapping. Evolution 64:694–709.
Donoghue MJ, Smith SA. 2004. Patterns in the assembly of temperate forests around the Northern Hemisphere. Phil. Trans. R. Soc. B 359:1633–1644.
122
Donoghue MJ, Bell CD, Li J-H. 2001. Phylogenetic patterns in Northern Hemisphere plant geography. Int. J. Plant Sci. 162:S41–S52.
Dowle EJ, Morgan-Richards M, Trewick SA. 2013. Molecular evolution and the latitudinal biodiversity gradient. Heredity 110:501–510.
Drehmel D, Moncalvo J-M, Vilgalys R. 1999. Molecular phylogeny of Amanita based on large-subunit ribosomal DNA sequences: implications for taxonomy and character evolution. Mycologia 91:610–618.
Drummond AJ, Bouckaert RR. 2015. Bayesian evolutionary analysis with BEAST. Cambridge University Press. Available at: http://beast2.org/book.html
Drummond AJ, Suchard MA, Xie D, Rambaut A. 2012. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 29:1969–1973.
Du X-H, Zhao Q, O’Donnell K, Rooney AP, Yang ZL. 2012. Multigene molecular phylogenetics reveals true morels (Morchella) are especially species-rich in China. Fungal Genet. Biol. 49:455–469.
Edgar RC. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32:1792–1797.
Edwards SV. 2009. Is a new and general theory of molecular systematics emerging? Evolution 63:1–19.
Edwards SV, Chesnut K, Satta Y, Wakeland EK. 1997. Ancestral polymorphism of MHC class II genes in mice: implications for balancing selection and the mammalian molecular clock. Genetics 146:655–668.
Ekblom R, Galindo J. 2011. Application of next generation sequencing in molecular ecology of non-model organisms. Heredity 107:1–15.
Ellegren H. 2014. Genome sequencing and population genomics in non-model organisms. Trends Ecol. Evol. 29:51–63.
Ellison CE, Hall C, Kowbel D, Welch J, Brem RB, Glass NL, Taylor JW. 2011. Population genomics and local adaptation in wild isolates of a model microbial eukaryote. Proc. Natl. Acad. Sci. USA. 108:2831–2836.
Escalante T, Rodríguez G, Morrone JJ. 2004. The diversification of Nearctic mammals in the Mexican transition zone. Biological Journal of the Linnean Society 83:327–339.
Etienne RS, Haegman B. 2012. A conceptual and statistical framework for adaptive radiations with a key role for diversity dependence. Am. Nat. 180:75–89.
Etienne RS, Rosindell J. 2012. Prolonging the past counteracts the pull of the present: protracted speciation can explain observed slowdowns in diversification. Syst. Biol. 61:204–213.
Etienne RS, Haegeman B, Stadler T, Aze T, Pearson PN, Purvis A, Phillimore AB. 2012. Diversity-dependence brings molecular phylogenies closer to agreement with the fossil record. Proc. Royal Soc. B 279:1300–1309.
Etienne RS, Morlon H, Lambert A. 2014. Estimating the duration of speciation from phylogenies. Evolution 68:2430–2440.
123
Eyre-Walker A, Keightley PD. 2007. The distribution of fitness effects of new mutations. Nat. Rev. Genet. 8:610–618.
Eyre-Walker A. 2006. The genomic rate of adaptive evolution. Trends Ecol. Evol. 21:569–575. Fay JC, Wu C-I. 2003. Sequence divergence, functional constraint, and selection in protein
evolution. Annu. Rev. Genom. Human Genet. 4:213–235. Fay JC, Wyckoff GJ, Wu CI. 2001. Positive and negative selection on the human genome.
Emerging fungal threats to animal, plant and ecosystem health. Nature 484:186–194. Fitzjohn RG. 2010. Quantitative traits and diversification. Syst. Biol. 59:619–633.
FitzJohn RG. 2012. Diversitree: comparative phylogenetic analyses of diversification in R. Methods Ecol. Evol. 3:1084-1092.
FitzJohn RG, Maddison WP, Otto SP. 2009. Estimating trait-dependent speciation and extinction rates from incompletely resolved phylogenies. Syst. Biol. 58:595–611.
Flot JF. 2010. SeqPHASE: a web tool for interconverting PHASE input/output files and FASTA sequence alignments. Mol. Ecol. Res. 10:162–166.
Fu YX, Li WH. 1993. Statistical tests of neutrality of mutations. Genetics 133:693–709. Fu YX, Li WH. 1999. Coalescing into the 21st century: An overview and prospects of coalescent
theory. Theor. Popul. Biol. 56:1–10. Fujita MK, Leache AD, Burbrink FT, McGuire JA, Moritz C. 2012. Coalescent-based species
delimitation in an integrative taxonomy. Trends Ecol. Evol. 27:480–488. Galante TE, Horton TR, Swaney DP. 2011. 95% of basidiospores fall within 1 m of the cap: a
field- and modeling-based study. Mycologia 103:1175–1183. Gao X, Gulari E, Zhou X. 2004. In situ synthesis of oligonucleotide microarrays. Biopolymers
73:579–596. Garibay-Orijel R, Caballero J, Estrada-Torres A, Cifuentes J. 2007. Understanding cultural
significance, the edible mushrooms case. J. Ethnobiol. Ethnomed. 3:4. Gaston KJ. 2000. Global patterns of biodiversity. Nature 405:220–227.
Geml, J., A. Laursen, K. O'Neill, H. C. Nusbaum and D. L. Taylor. 2006. Beringian origins and cryptic speciation events in the fly agaric (Amanita muscaria). Mol. Ecol. 15:225–239.
Geml J, Tulloss RE, Laursen GA, Sazanova NA, Taylor DL. 2008. Evidence for strong inter- and intracontinental phylogeographic structure in Amanita muscaria, a wind-dispersed ectomycorrhizal basidiomycete. Mol. Phylogenet. Evol. 48:694–701.
Geml J, Timling I, Robinson CH, Lennon N, Nusbaum HC, Brochman C, Noordeloos ME, Taylor DL. 2012. An arctic community of symbiotic fungi assembled by long-distance dispersers: phylogenetic diversity of ectomycorrhizal basidiomycetes in Svalbard based on soil and sporocarp DNA. J. Biogeogr. 39:74–88.
124
Gerstein AC, Kuzmin A, Otto SP. 2014. Loss-of-heterozygosity facilitates passage through Haldane’s sieve for Saccharomyces cerevisiae undergoing adaptation. Nat. Comm. 5:3819.
Gillooly JF, Brown JH, West GB, Savage VM, Charnov EL. 2001. Effects of size and temperature on metabolic rate. Science 293:2248–2251.
Gladieux P, Devier B, Aguileta G, Cruaud C, Giraud T. 2013. Purifying selection after episodes of recurrent adaptive diversification in fungal pathogens. Infect. Genet. Evol. 17:123–131.
Gladieux P, Ropars J, Badouin H, Branca A, Aguileta G, de Vienne DM, Rodríguez de la Vega RC, Branco S, Giraud T. 2014. Fungal evolutionary genomics provides insight into the mechanisms of adaptive divergence in eukaryotes. Mol. Ecol. 23:753–773.
Glor. R. 2010. Phylogenetic insights on adaptive radiation. Annu. Rev. Ecol. Evol. Syst. 41:251–70.
Goldberg EE, Kohn JR, Lande R, Robertson KA, Smith SA, Igi! B. 2010. Species selection maintains self-incompatibility. Science 330:493–495.
Goldberg EE, Lancaster LT, Ree RH. 2011. Phylogenetic inference of reciprocal effects between geographic range evolution and diversification. Syst. Biol. 60:451-465.
Goldman N, Yang Z. 1994. A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol. Biol. Evol. 11:725–736.
Gräfe K, Frisch W, Villa IM, Meschede M. 2002. Geodynamic evolution of southern Costa Rica related to low-angle subduction of the Cocos Ridge: constraints from thermochronology. Tectonophysics 348:187–204.
Gutiérrez-Rodríguez C, Ornelas JF, Rodríguez-Gómez F. 2011. Chloroplast DNA phylogeography of a distylous shrub (Palicourea padifolia, Rubiaceae) reveals past fragmentation and demographic expansion in Mexican cloud forests. Mol. Phylogenet. Evol. 61:603–615.
Guzmán G, Ramírez-Guillén F. 2001. The Amanita caesarea-complex. Bibliotheca Mycologica 187. J. Cramer, Berlin, Germany.
Hallen HE, Luo H, Scott-Craig JS, Walton JD. 2007. Gene family encoding the major toxins of lethal Amanita mushrooms. Proc. Natl. Acad. Sci. USA 104:19097–19101.
Hancock-Hanser BL, Frey A, Leslie MS, Dutton PH, Archer FI, Morin PA. 2013. Targeted multiplex next-generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics. Mol. Ecol. Res. 13:254–268.
Harmon LJ, Schulte II JA, Larson A, Losos JB. 2003. Tempo and mode of evolutionary radiation in iguanian lizards. Science 301:961-964.
Heled J, Drummond AJ. 2008. Bayesian inference of population size history from multiple loci. BMC Evol. Biol. 8:289.
Heled J, Drummond AJ. 2010. Bayesian inference of species trees from multilocus data. Mol. Biol. Evol. 27:570–580.
Henkel TW, Meszaros R, Aime MC, Kennedy A. 2005. New Clavulina species from the Pakaraima Mountains of Guyana. Mycol. Progr. 4:343–350.
125
Hermisson J, Pennings PS. 2005. Soft Sweeps: molecular population genetics of adaptation from standing genetic variation. Genetics. 169:2335–2352.
Hernández CE, Rodríguez-Serrano E, Avaria-Llautureo J, Inostroza-Michael O, Morales-Pallero B, Boric-Bargetto D, Canales-Aguirre CB, Marquet PA, Meade A. 2013. Using phylogenetic information and the comparative method to evaluate hypotheses in macroecology. Methods Ecol. Evol. 4:401–415.
Hess J, Skrede I, Wolfe BE, LaButti K, Ohm RA, Grigoriev IV, Pringle A. 2014. Transposable element dynamics among asymbiotic and ectomycorrhizal Amanita fungi. Genome Biol. Evol. 6: 1564–1578.
Hewitt G. 2000. The genetic legacy of the Quaternary ice ages. Nature 405:907–913.
Hewitt GM. 1999. Post!glacial re!colonization of European biota. Biol. J. Linn. Soc. 68:87–112.
Hewitt GM. 1996. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol. J. Linn. Soc. 58:247–276.
Hewitt GM. 2004a. Genetic consequences of climatic oscillations in the Quaternary. Phil. Trans. R. Soc. B. 359:183–195.
Hewitt GM. 2004b. The structure of biodiversity – insights from molecular phylogeography. Front. Zool. 1:4.
Hey J. 2006. On the failure of modern species concepts. Trends Ecol. Evol. 21:447–450. Hibbett DS, Grimaldi D, Donoghue MJ. 1997. Fossil mushrooms from Miocene and Cretaceous
amber and the evolution of homobasidiomycetes. Am. J. Bot. 84:981–991. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. 2005. Very high resolution interpolated
climate surfaces for global land areas. Int. J. Climatol. 25:1965–1978. Hijmans RJ, Phillips S, Leathwick J, Elith J. 2011. R package ‘dismo’. Available online at:
http://cran.r-project.org/web/packages/dismo/index.html. Hijmans RJ, van Etten J. 2012. raster: Geographic analysis and modeling with raster data.
Version 2.0-12. Available at: http://CRAN.R-project.org/package=raster Hillebrand, H. 2004. On the generality of the latitudinal diversity gradient. Am. Nat. 163, 192–
211. Hillebrand H, Azovsky AI. 2001. Body size determines the strength of the latitudinal diversity
gradient. Ecography 24:251–256. Hinz UG, Fivaz J, Girod PA, Zyrd JP. 1997. The gene coding for the DOPA dioxygenase
involved in betalain biosynthesis in Amanita muscaria and its regulation. Mol. Gen. Genet. 256:1-6.
Hoffmann AA, Sgrò CM. 2011. Climate change and evolutionary adaptation. Nature 470:479–485.
Hosaka, K., Castellano, M.A. and Spatafora, J.W. 2008. Biogeography of Hysterangiales (Phallomycetidae, Basidiomycota). Mycol. Res. 112:448–462.
Hough J, Williamson RJ, Wright SI. 2013. Patterns of selection in plant genomes. Annu. Rev. Ecol. Evol. Syst. 44:31–49.
126
Huang S, Roy K, Jablonski D. 2014. Do past climate states influence diversity dynamics and the present!day latitudinal diversity gradient? Global Ecol. Biogeogr. 23:530–540.
Hughes AL, Nei M. 1988. Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature 335:167–170.
Ingram T, Mahler DL. 2013. SURFACE: detecting convergent evolution from comparative data by fitting Ornstein-Uhlenbeck models with stepwise AIC. Methods Ecol. Evol. 4:416–425.
Innan H. 2003. The coalescent and infinite-site model of a small multigene family. Genetics 163:803–810.
Irestedt M, Gelang M, Sangster G, Olsson U, Ericson PGP, Alström P. 2011. Neumann’s Warbler Hemitesia neumanni (Sylvioidea): the sole African member of a Palaeotropic Miocene avifauna. Ibis 153:78–86.
Jablonski D, Roy K, Valentine JW. 2006. Out of the tropics: evolutionary dynamics of the latitudinal diversity gradient. Nature 314:102–106.
James TY, Porter D, Hamrick JL, Vilgalys R. 1999. Evidence for limited intercontinental gene flow in the cosmopolitan mushroom, Schizophyllum commune. Evolution 53:1665–1677.
James TY, Moncalvo J-M, Li S, Vilgalys R. 2001. Polymorphism at the ribosomal DNA spacers and its relation to breeding structure of the widespread mushroom Schizophyllum commune. Genetics 157:149–161.
Jargeat P, Martos F, Carriconde F, Gryta H, Moreau PA, Gardes M. 2010. Phylogenetic species delimitation in ectomycorrhizal fungi and implications for barcoding: the case of the Tricholoma scalpturatum complex (Basidiomycota). Mol. Ecol. 19:5216–5230.
Jeandroz S, Murat C, Wang Y, Bonfante P, Le Tacon F. 2008. Molecular phylogeny and historical biogeography of the genus Tuber, the 'true truffles'. J. Biogeogr. 35:815–829.
Jetz W, Rahbek C, Colwell RK. 2004. The coincidence of rarity and richness and the potential signature of history in centres of endemism. Ecol. Lett. 7:1180–1191.
Justo A, Morgenstern I, Hallen-Adams HE, Hibbett DS. 2010. Convergent evolution of sequestrate forms in Amanita under Mediterranean climate conditions. Mycologia 102:675–688.
Kennedy PG, Matheny PB, Ryberg KM, Henkel TW, Uehling JK, Smith ME. 2012. Scaling up: examining the macroecology of ectomycorrhizal fungi. Mol. Ecol. 21:4151-4154.
Kennedy PG, Garibay-Orijel R, Higgins LM, Angeles-Arguiz R. 2011. Ectomycorrhizal fungi in Mexican Alnus forests support the host co-migration hypothesis and continental-scale patterns in phylogeography. Mycorrhiza 21:559–568.
Keppel G, Van Niel KP, Wardell-Johnson GW, Yates CJ, Byrne M, Mucina L, Schut AGT, Hopper SD, Franklin SE. 2012. Refugia: identifying and understanding safe havens for biodiversity under climate change. Global Ecol. Biogeogr. 21:393–404.
Kimura M. 1977. Preponderance of synonymous changes as evidence for the neutral theory of molecular evolution. Nature 267:275–276.
Kimura M. 1986. DNA and the neutral theory. Philos. Trans. R. Soc. Lond. B. 312:343–354.
127
Kindlmann P, Schödelbauerová I, Dixon AG. 2007. Inverse latitudinal diversity gradients in species diversity. pp. 276–257. In: Scaling Biodiversity (eds. D. Storch, P. A. Maquet and J. H. Brown), Cambridge University Press, Cambridge, UK.
Kirk PM, Cannon PF, Minter DW, Stalpers JA. 2008. Dictionary of the Fungi: CABI. Wallingford, UK.
Knowles LL. 2001. Did the Pleistocene glaciations promote divergence? Tests of explicit refugial models in montane grasshopprers. Mol. Ecol. 10, 691–701.
Kodandaramaiah U, Wahlberg N. 2007. Out-of-Africa origin and dispersal-mediated diversification of the butterfly genus Junonia (Nymphalidae: Nymphalinae). J. Evol. Biol. 20:2181–2191.
Kohler A, Kuo A, Nagy LG, Morin E, Barry KW, Buscot F, Canbäck B, Choi C, Cichocki N, Clum A, Colpaert J, Copeland A, Costa MD, Doré J, Floudas D, Gay G, Girlanda M, Henrissat B, Herrmann S, Hess J, Högberg N, Johansson T, Khouja HR, LaButti K, Lahrmann U, Levasseur A, Lindquist EA, Lipzen A, Marmeisse R, Martino E, Murat C, Ngan CY, Nehls U, Plett JM, Pringle A, Ohm RA, Perotto S, Peter M, Riley R, Rineau F, Ruytinx J, Salamov A, Shah F, Sun H, Tarkka M, Tritt A, Veneault-Fourrey C, Zuccaro A, Tunlid A, Grigoriev IV, Hibbett DS, Martin F. 2015. Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nat Genet 47:410–415.
Kontoyiannis DP, Lewis RE. 2002. Antifungal drug resistance of pathogenic fungi. Lancet 359:1135–1144.
Köstekci H, Yamaç M, Solak MH. 2005. Macrofungi of Türkmenbaba Mountain (Eski!ehir). Turk. J. Bot. 29:409–416.
Kryazhimskiy S, Plotkin JB. 2008. The population genetics of dN/dS. PLoS Genet. 4:e1000304. Lambert A, Morlon H, Etienne RS. 2015. The reconstructed tree in the lineage-based model of
protracted speciation. J. Math. Biol. 70:367–97. Lanfear R, Calcott B, Ho SYW, Guindon S. 2012. PartitionFinder: combined selection of
partitioning schemes and substitution models for phylogenetic analyses. Mol. Biol. Evol. 29:1695–1701.
Lavin M, Luckow M. 1993. Origins and relationships of tropical North America in the context of the boreotropics hypothesis. Am. J. Bot. 80:1–14.
Le Gac M, Hood ME, Giraud T. 2007. Evolution of reproductive isolation within a parasitic fungal species complex. Evolution 61:1781–1787.
Lemey P, Rambaut A, Drummond AJ, Suchard MA. 2009. Bayesian phylogeography finds its roots. PLoS Comp. Biol. 5:e1000520.
Lemmon A, Lemmon EM. 2008. A likelihood framework for estimating phylogeographic history on a continuous landscape. Syst. Biol. 57:544–561.
Lemmon EM, Lemmon AR. 2013. High-throughput genomic data in systematics and phylogenetics. Annu. Rev. Ecol. Evol. Syst. 44:99–121.
Levin DA. 1995. Metapopulations: an arena for local speciation. J. Evol. Biol. 8:635–644.
128
Li H, Durbin R. 2010. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26:589–595.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup. 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079.
Liao J, Li Z, Hiebeler DE, El-Bana M, Deckmyn G, Nijs I. 2013. Modelling plant population size and extinction thresholds from habitat loss and habitat fragmentation: effects of neighbouring competition and dispersal strategy. Ecol. Model. 268:9–17.
Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25:1451–1452.
Lieberman BS and Eldredge N. 2014. What is punctuated equilibrium? What is macroevolution? A response to Pennell et al. Trends Ecol. Evol. 29:185–186.
Liu L, Yu L, Kubatko L, Pearl D K, Edwards SV. 2009. Coalescent methods for estimating phylogenetic trees. Mol. Phylogenet. Evol. 53:320–328.
Liu ZL. 2006. Genomic adaptation of ethanologenic yeast to biomass conversion inhibitors. Appl. Microbiol. Biotechnol. 73:27–36.
Lumbsch HT, Buchanan PK, May TW, Mueller GM. 2008. Phylogeography and biogeography of fungi. Mycol. Res. 112:423–424.
Maddison WP, FitzJohn RG. 2015. The unsolved challenge to phylogenetic correlation tests for categorical characters. Syst. Biol. 64:127–136
Maddison WP, Midford PE, Otto SP. 2007. Estimating a binary character's effect on speciation and extinction. Syst. Biol. 56:701–710.
Mahler DL, Ingram T, Revell LJ, Losos JB. 2013. Exceptional convergence on the macroevolutionary landscape in island lizard radiations. Science 341:292–295.
Martin F, Cullen D, Hibbett D, Pisabarro A, Spatafora JW, Baker SE, Grigoriev IV. 2011. Sequencing the fungal tree of life. New Phytol. 190:818–821.
Martin F, Díez J, Dell B, Delaruelle C. 2002. Phylogeography of the ectomycorrhizal Pisolithus species as inferred from nuclear ribosomal DNA ITS sequences. New Phytol. 153:345–357.
Martiny JBH, Bohannan BJ, Brown JH, Colwell RK, Fuhrman JA, Green JL, Horner-Devine MC, Kane M, Adams Krumins J, Kuske CR, Morin PJ, Naeem S, Øvreås L, Reysenbach A-L, Smith VH, Staley JT. 2006. Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4:102–112.
Mastretta-Yanes A, Moreno-Letelier A, Piñero D, Jorgensen TH, Emerson BC (2015) Biodiversity in the Mexican highlands and the interaction of geology, geography and climate within the Trans-Mexican Volcanic Belt. J. Biogeogr. doi:10.1111/jbi.12546
Matheny PB, Aime MC, Bougler NL, Buyck B, Desjardin DE, Horak E, Kropp BR, Lodge DJ, Soytong K, Trappe JM, Hibbett DS. 2009. Out of the Paleotropics? Historical
129
biogeography and diversification of the cosmopolitan ectomycorrhizal mushroom family Inocybaceae. J. Biogeogr. 36:577–592.
Matheny PB, Wang Z, Binder M. et al. (2007) Contributions of rpb2 and tef1 to the phylogeny of mushrooms and allies (Basidiomycota, Fungi). Mol. Phylogenet. Evol. 43:430–451.
Matute DR, McEwen JG, Puccia R, Montes BA, San-Blas G, Bagagli E, Rauscher JT, Restrepo A, Morais F, Niño-Vega G, Taylor JW. 2006. Cryptic speciation and recombination in the fungus Paracoccidioides brasiliensis as revealed by gene genealogies. Mol. Biol. Evol. 23:65–73.
McDonald JH, Kreitman M. 1991. Adaptive protein evolution at the Adh locus in Drosophila. Nature 351:652–654.
McInnes L, David C, Orme L, Purvis A. 2011. Detecting shifts in diversity limits from molecular phylogenies: what can we know? Proc. R. Soc. B 278:3294–3302.
McPeek MA. 2008. The ecological dynamics of clade diversification and community assembly. Am. Nat. 172:270–284.
Médail F, Diadema K. 2009. Glacial refugia influence plant diversity patterns in the Mediterranean Basin. J. Biogeogr. 36:1333–1345.
Messer PW, Petrov DA. 2013a. Frequent adaptation and the McDonald-Kreitman test. Proc. Natl. Acad. Sci. USA. 110:8615–8620.
Messer PW, Petrov DA. 2013b. Population genomics of rapid adaptation by soft selective sweeps. Trends Ecol. Evol. 28:659–669.
Milne RI. 2006. Northern Hemisphere plant disjunctions: a window on Tertiary land bridges and climate change? Ann. Bot. 98:465–472.
Mittelbach GG, Schemske DW, Cornell HV, Allen AP, Brown JM, Bush MB, Harrison SP, Hurlbert AH, Knowlton N, Lessios HA, McCain CM, McCune AR, McDade LA, McPeek MA, Near TJ, Price TD, Ricklefs RE, Roy K, Sax DF, Schluter D, Sobel JM, Turelli M. 2007. Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography. Ecol. Lett. 10:315–331.
Moen D, Morlon H. 2014. Why does diversification slow down? Trends. Ecol. Evol. 29:190-197.
Moncalvo J-M, Drehmel D, Vilgalys R. 2000. Variation in modes and rates of evolution in nuclear and mitochondrial ribosomal DNA in the mushroom genus Amanita (Agaricales, Basidiomycota): phylogenetic implications. Mol. Phylogenet. Evol. 16:48–63.
Moncalvo J-M, Buchanan PK. 2008. Molecular evidence for long distance dispersal across the South Hemisphere in Ganoderma applanatum-australe species complex (Basidiomycota). Mycol. Res. 112:425–436.
Moreno-Hagelsieb G, Latimer K. 2008. Choosing BLAST options for better detection of orthologs as reciprocal best hits. Bioinformatics 24:319–324.
Morlon H, Potts MD, Plotkin JB. 2010. Inferring the dynamics of diversification: a coalescent approach. PLoS Biol. 8:e1000493.
Moyersoen B, Beever RE, Martin F. (2003) Genetic diversity of Pisolithus in New Zealand indicates multiple long-distance dispersal from Australia. New Phytol. 160:569–579.
130
Mulcahy DG. 2008. Phylogeography and species boundaries of the western North American Nightsnake (Hypsiglena torquata): Revisiting the subspecies concept. Mol. Phylogenet. Evol. 46:1095–1115.
Nee S. 2006. Birth-Death models in macroevolution. Annu. Rev. Ecol. Evol. Syst. 37:1–17.
Nee S, May RM, Harvey PH. 1994. The reconstructed evolutionary process. Phil. Trans. R. Soc. Lond. B 344:305–311.
Nehls U, Grunze N, Willmann M, Reich M, Küster H. 2007. Sugar for my honey: carbohydrate partitioning in ectomycorrhizal symbiosis. Phytochemistry 68:82–91.
Nei M. 1987. Molecular Evolutionary Genetics. Columbia University Press. Nei M, Li W-H. 1979. Mathematical model for studying genetic variation in terms of restriction
endonuclease. Proc. Natl. Acad. Sci. USA. 89:1477–1481. Nielsen R, Williamson S, Kim Y, Hubisz MJ, Clark AG, Bustamante C. 2005. Genomic scans
for selective sweeps using SNP data. Genome Res. 15:1566–1575. Nielsen R. 1997. The ratio of replacement to silent divergence and tests of neutrality. J. Evol.
Biol. 10:217–231. Nielsen R. 2001. Statistical tests of selective neutrality in the age of genomics. Heredity 86:641–
647. O'Meara BC. 2012. Evolutionary inferences from phylogenies: a review of methods. Annu. Rev.
Ecol. Evol. Syst. 43:267–85. Ohta T. 1992. The nearly neutral theory of molecular evolution. Annu. Rev. Ecol. Syst. 23:263–
D'amico JA, et al. 2001. Terrestrial Ecoregions of the World: A new map of life on Earth, A new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity. BioScience 51:933–938.
Ornelas JF, Ruiz-Sánchez E, Sosa V. 2010. Phylogeography of Podocarpus matudae (Podocarpaceae): pre-Quaternary relicts in northern Mesoamerican cloud forests. J. Biogeogr. 37:2384–2396.
Pagel M. 1994. Detecting correlated evolution on phylogenies: a general method for the comparative analysis of discrete characters. Proc. R. Soc. B 255:37–45.
Paradis E, Claude J, Strimmer K. 2004. APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20:289–290.
Parra G, Bradnam K, Korf I. 2007. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics 23:1061–1067.
Parra G, Bradnam K, Ning Z, Keane T, Korf I. 2009. Assessing the gene space in draft genomes. Nucleic Acids Res. 37:298–297.
Peay KG, Bruns TD, Kennedy PG, Bergemann SE, Garbelotto M. 2007. A strong species–area relationship for eukaryotic soil microbes: island size matters for ectomycorrhizal fungi. Ecol. Lett. 10:470–480.
131
Peay KG, Bidartondo MI, Arnold AE. 2010. Not every fungus is everywhere: scaling to the biogeography of fungal–plant interactions across roots, shoots and ecosystems. New Phytol. 185:878–882.
Pegler DN. 2002. Useful fungi of the world: Caesar's mushroom and the Christmas mushroom. Mycologist 16:140–41.
Perry JP, Graham A, Richardson M. 1998. The history of pines in Mexico and Central America. pp. 137–149. In: (ed. Richardson DM) Ecology and Biogeography of Pinus. Cambridge University Press, Cambridge, UK.
Petit RJ, Aguinagalde I, de Beaulieu JL, Bittkau C, Brewer S, Cheddadi R, Ennos R, Fineschi S, Grivet D, Lascoux M, Mohanty A, Müller-Starck G, Demesure-Musch B, Palmé A, Martín JP, Rendell S, Vendramin GG. 2003. Glacial refugia: hotspots but not melting pots of genetic diversity. Science 300:1563–1565.
Phillimore AB, Price TD. 2008. Density-dependent cladogenesis in birds. PLoS Biol. 6:e71. Phillips SJ, Anderson RP, Schapire RE. 2006. Maximum entropy modeling of species
geographic distributions. Ecol. Model. 190:231–259. Pianka ER. 1966. Latitudinal gradients in species diversity: a review of concepts. Am. Nat.
100:33–46. Pielou EC. 1991. After the Ice Age: the return of life to glaciated North America. The University
of Chicago Press, Chicago. Piper P, Calderon CO, Hatzixanthis K, Mollapour M. 2001. Weak acid adaptation: the stress
response that confers yeasts with resistance to organic acid food preservatives. Microbiology 147:2635–2642.
Pollard DA, Iyer VN, Moses AM, Eisen MB. 2006. Widespread discordance of gene trees with species tree in Drosophila: evidence for incomplete lineage sorting. PLoS Genet. 2:e173.
Põlme S, Bahram M, Yamanaka T, Nara K, Dai YC, Grebenc T, Kraigher H, Toivonen M, Wang P-H, Matsuda Y, Naadel T, Kennedy PG, Kõljalg U, Tedersoo L. 2013. Biogeography of ectomycorrhizal fungi associated with alders (Alnus spp.) in relation to biotic and abiotic variables at the global scale. New Phytol. 198:1239–1249.
Potter PE, Szatmari P. 2009. Global Miocene tectonics and the modern world. Earth-Sci. Rev. 96:279–295.
Pringle A, Bever JD, Gardes M, Parrent JL, Rillig MC, Klironomos JN. 2009. Mycorrhizal symbioses and plant invasions. Annu. Rev. Ecol. Evol. Syst. 40:699–715.
Pybus OG, Harvey PH. 2000. Testing macroevolutionary models using incomplete molecular phylogenies. Proc. R. Soc. B 267:2267–2272.
Pybus OG, Rambaut A, Holmes EC, Harvey PH. 2002. New inferences from tree shape: numbers of missing taxa and population growth rates. Syst. Biol. 51:881–888.
Pyron RA, Wiens JJ. 2013. Large-scale phylogenetic analyses reveal the causes of high tropical amphibian diversity. Proc. R. Soc. B 280:20131622.
132
Qian H, Ricklefs RE. 2000. Large-scale processes and the Asian bias in species diversity of temperate plants. Nature 407:180–182.
Queloz V, Sieber TN, Holdenrieder O, McDonald BA, Grünig CR. (2011) No biogeographical pattern for a root-associated fungal species complex. Global Ecol. Biogeogr. 20:160–169.
Quental TB, Marshall CR. 2010. Diversity dynamics: molecular phylogenies need the fossil record. Trends Ecol. Evol. 25:434–441.
R Core Team (2015) R: a language and environment for statistical computing. Version 3.0.1. R Foundation for Statistical Computing, Vienna, Austria.
Rabosky DL. 2009. Ecological limits on clade diversification in higher taxa. Am. Nat. 173:662–674.
Rabosky DL. 2010. Extinction rates should not be estimated from molecular phylogenies. Evolution 64:1816-1824.
Rabosky DL. 2013. Diversity-dependence, ecological speciation, and the role of competition in macroevolution. Annu. Rev. Ecol. Evol. Syst. 44:481-502.
Rabosky DL. 2014. Automatic detection of key innovations, rate shifts, and diversity-dependence on phylogenetic trees. PLoS One 9:e89543.
Rabosky DL, Goldberg EE. 2015. Model inadequacy and mistaken inferences of trait-dependent speciation. Syst. Biol. 64:340–355.
Rabosky DL, Lovette IJ. 2008. Explosive evolutionary radiations: decreasing speciation or increasing extinction through time? Evolution 62:1866-1875.
Rabosky DL, Slater GL, Alfaro ME. 2012. Clade age and species richness are decoupled across the eukaryotic tree of life. PLoS Biol. 10:e1001381.
Rambaut A. 2002. SeAl: sequence alignment editor. Version 2.0. Available at: http://tree.bio.ed.ac.uk/software/seal/
Rambaut A, Suchard M, Drummond A. 2015. Tracer. Version 1.6. Available at: http://tree.bio.ed.ac.uk/software/tracer
Ramírez-Soriano A, Ramos-Onsins SE, Rozas J, Calafell F, Navarro A. 2008. Statistical power analysis of neutrality tests under demographic expansions, contractions and bottlenecks with recombination. Genetics 179:555–567.
Rannala B, Yang Z. 2003. Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci. Genetics 164:1645–1656.
Rannala B, Yang Z. 2013. Improved reversible jump algorithms for Bayesian species delimitation. Genetics 194:245–253.
Ree RH, Smith SA. 2008. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst. Biol. 57:4–14.
Rehner SA. 2001. Primers for elongation factor 1-! (EF1-!). Assembling the Fungal Tree of Life. Available at: http://www.aftol.org/pdfs/EF1primer.pdf
Reid DA. 1980. A monograph of the Australian species of Amanita Persoon ex Hooker (Fungi). Aust. J. Bot. Suppl. 8:1–96.
133
Rohde K. 1992. Latitudinal gradients in species diversity: the search for the primary cause. Oikos 65:514–527.
Rokas A, Williams BL, King N, Carroll SB. 2003. Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature 425:798–804.
Rolland J, Condamine FL, Jiguet F, Morlon H. 2014. Faster speciation and reduced extinction in the tropics contribute to the mammalian latitudinal diversity gradient. PLoS Biol. 12:e1001775.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. 2012. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 61:539–542.
Ruiz-Sánchez E, Ornelas JF. 2014. Phylogeography of Liquidambar styraciflua (Altingiaceae) in Mesoamerica: survivors of a Neogene widespread temperate forest (or cloud forest) in North America? Ecol. Evol. 4:311–328.
Ryberg M, Matheny PB. 2011. Dealing with incomplete taxon sampling and diversification of a large clade of mushroom-forming fungi. Evolution 65:1862–1878.
Ryberg M, Matheny PB. 2012. Asynchronous origins of ectomycorrhizal clades of Agaricales. Proc. Royal Soc. B 279:2003–2011.
Sakai AK, Allendorf FW, Holt JS, Lodge DM, Molofsky J, With KA, Baughman S, Cabin RJ, Cohen JE, Ellstrand CN, McCauley DE, O’Neil P, Parker IM, Thompson JN, Weller SG. 2001. The population biology of invasive species. Annu. Rev. Ecol. Syst. 32:305–332.
Sánchez-Ramírez S, Tulloss RE, Amalfi M, Moncalvo J-M. 2015a. Palaeotropical origins, boreotropical distribution, and increasing rates of diversification in a clade of edible ectomycorrhizal fungi (Amanita sect. Caesareae). J. Biogeogr. 42:351–363.
Sánchez-Ramírez S, Etienne RS, Moncalvo J-M. 2015b. High speciation rate at temperate latitudes explains unusual diversity gradients in a clade of ectomycorrhizal fungi. Evolution 69:2196–2209.
Sánchez-Ramírez, S., R. E. Tulloss, J. Cifuentes-Blanco, L. Guzmán-Dávalos, R. Valenzuela, A. Estrada-Torres, F. Ruán-Soto, R. Díaz-Moreno, N. Hernández-Rico, M. Torres-Gómez, H. León and J.-M. Moncalvo. Accepted. In and out of refugia: historical patterns of diversity and demography in the North American Caesar’s mushroom species complex. Molecular Ecology.
Sanger TJ, Revell LJ, Gibson-Brown JJ, Losos JB. 2012. Repeated modification of early limb morphogenesis programmes underlies the convergence of relative limb length in Anolis lizards. Proc. R. Soc. B 279:739–748.
Sanmartín I, Ronquist F. 2004. Southern Hemisphere biogeography inferred by event-based models: plant versus animal patterns. Syst. Biol. 53:216–243.
Savolainen O, Pyhäjärvi T, Knürr T. 2007. Gene flow and local adaptation in trees. Annu. Rev. Ecol. Evol. Syst. 38:595–619.
Schemske 2009. Biotic interactions and speciation in the tropics. pp 219-239. In: (eds. RK Butlin, JR Bridle, D Schluter) Speciation and Patterns of Diversity, Cambridge University Press, Cambridge.
134
Schluter D. 2000. The ecology of adaptive radiation. Oxford University Press, New York, 288 pp.
Schmit JP, Mueller GM. 2007. An estimate of the lower limit of global fungal diversity. Biodivers. Conserv. 16:99–111.
Schmitt T. 2007. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Frontiers in Zoology 4:11.
Schoustra SE, Debets AJM, Slakhorst M, Hoekstra RF. 2007. Mitotic recombination accelerates adaptation in the fungus Aspergillus nidulans. PLoS Genet. 3:e68.
Selbmann L, Egidi E, Isola D, Onofri S, Zucconi L, de Hoog GS, Chinaglia S, Testa L, Tosi S, Balestrazzi A, Lantieri A, Compagno R, Tigini V, Varese GC. 2013. Biodiversity, evolution and adaptation of fungi in extreme environments. Plant Biosyst. 147:237–246.
Sella G, Petrov DA, Przeworski M, Andolfatto P. 2009. Pervasive natural selection in the Drosophila genome? PLoS Genet 5:e1000495.
Shepard GH, Arora D, Lampman A. 2008. The grace of the flood: classification and use of wild mushrooms among the highland Maya of Chiapas. Econ. Bot. 62:437–470.
Sherratt TN, Wilkinson DM. 2009. Big questions in ecology and evolution. Oxford Biology Series, Oxford University Press, Oxford, 297 pp.
Silvestro D, Cascales-Miñana B, Bacon C, Antonelli A. 2015. Revisiting the origin and diversification of vascular plants through a comprehensive Bayesian analysis of the fossil record. New Phytol. 207:424–436.
Silvestro D, Zizka G, Schulte K. 2014. Disentangling the effects of key innovations on the diversification of Bromelioideae (Bromeliaceae). Evolution 68:163–175.
Slotte T, Foxe JP, Hazzouri KM, Wright SI. 2010. Genome-wide evidence for efficient positive and purifying selection in Capsella grandiflora, a plant species with a large effective population size. Mol. Biol. Evol. 27:1813–1821.
Smith NGC, Eyre-Walker A. 2002. Adaptive protein evolution in Drosophila. Nature 415:1022–1024.
Smith SE, Read DJ. 2008. Mycorrhizal symbiosis. Third edition. Elsevier.
Smith SY, Currah RS, Stockey RA. 2004. Cretaceous and Eocene poroid hymenophores from Vancouver Island, British Columbia. Mycologia 96:180–186.
Soltis DE, Morris AB, McLachlan JS, Manos PS, Soltis PS. 2006. Comparative phylogeography of unglaciated eastern North America. Mol. Ecol. 15:4261–4293.
Städler T, Haubold B, Merino C, Stephan W, Pfaffelhuber P. 2009. The impact of sampling schemes on the site frequency spectrum in nonequilibrium subdivided populations. Genetics 182:205–216.
Stadler T. 2013. Recovering speciation and extinction dynamics based on phylogenies. J. Evol. Biol. 26:1203–1219.
135
Stadler T. 2011a. Inferring speciation and extinction processes from extant species data. Proc. Natl. Acad. Sci. USA 108:16145–16146.
Stadler T. 2011b. Mammalian phylogeny reveals recent diversification rate shifts. Proc. Natl. Acad. Sci. USA 108:6187–6192.
Stadler T. 2013a. Recovering speciation and extinction dynamics based on phylogenies. J. Evol. Biol. 26:1203–1219.
Stadler T. 2013b. How can we improve accuracy of macroevolutionary rate estimates? Syst. Biol. 62:321–329.
Stanke M, Keller O, Gunduz I, Hayes A, Waack S, Morgenstern B. 2006. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 34:W435–439.
Stanke M, Steinkamp R, Waack S, Morgenstern B. 2004. AUGUSTUS: a web server for gene finding in eukaryotes. Nucleic Acids Res. 32:W309-W31.
Stephens M, Donnelly P. 2003. A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am. J. Hum. Genet. 73:1162–1169.
Strack D, Vogt T, Schliemann W. 2003. Recent advances in betalain research. Phytochemistry 62:247–269.
Stuart AJ. 1991. Mammalian extinctions in the Late Pleistocene of northern Eurasia and North America. Biol. Rev. 66:453–562.
Suchard MA, Rambaut A. 2009. Many-core algorithms for statistical phylogenetics. Bioinformatics 25:1370–1376.
Suutari M, Liukkonen K, Laakso S. 1990. Temperature adaptation in yeasts: the role of fatty acids. J. Gen. Microbiol. 136:1469–1474.
Svenning JC, Skov F. 2007. Could the tree diversity pattern in Europe be generated by postglacial dispersal limitation? Ecol. Lett. 10:453–460.
Taberlet P, Cheddadi R. 2002. Ecology. Quaternary refugia and persistence of biodiversity. Science 297:2009–2010.
Taberlet P, Fumagalli L, Wust-Saucy AG, Cosson JF. 1998. Comparative phylogeography and postglacial colonization routes in Europe. Mol. Ecol. 7:453–464.
Tajima F. 1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585–595.
Talbot JM, Bruns TD, Taylor JW, Smith DP, Branco S, Glassman SI, Erlandson S, Vilgalys R, Liao H-L, Smith ME, Peay KG. 2014. Endemism and functional convergence across the North American soil mycobiome. Proc. Natl. Acad. Sci. USA 111:6341–6346.
Taylor AFS. and Alexander I. 2005. The ectomycorrhizal symbiosis: life in the real world. Mycologist 19:102-112.
Taylor JW, Jacobson DJ, Kroken S, Kasuga T, Geiser DM, Hibbett DS, Fisher MC (2000) Phylogenetic species recognition and species concepts in fungi. Fungal Genetics and Biology, 31, 21–32.
Taylor JW, Turner E, Townsend JP, Dettman JR, Jacobson D (2006) Eukaryotic microbes, species recognition and the geographic limits of species: examples from the kingdom
136
Fungi. Philosophical Transactions of the Royal Society B: Biological Sciences, 361, 1947–1963.
Taylor DL, Hollingsworth TN, McFarland JW, Lennon NJ, Nusbaum C, Ruess RW. 2014. A first comprehensive census of fungi in soil reveals both hyperdiversity and fine-scale niche partitioning. Ecol. Monogr. 84:3–20.
Taylor JW, Jacobson DJ, Kroken S, Kasuga T, Geiser MD, Hibbett DS, Fisher MC. 2000. Phylogenetic species recognition and species concepts in fungi. Fungal Genet. Biol. 31:21–32.
Tedersoo L, Nara K. 2010. General latitudinal gradient of biodiversity is reversed in ectomycorrhizal fungi. New Phytol. 185:351–354.
Tedersoo L, Smith ME. 2013. Lineages of ectomycorrhizal fungi revisited: Foraging strategies and novel lineages revealed by sequences from belowground. Fungal Biol. Rev. 27:83–99.
Tedersoo L, Sadam A, Zambrano M, Valencia R, Bahram M. 2010a. Low diversity and high host preference of ectomycorrhizal fungi in Western Amazonia, a neotropical biodiversity hotspot. ISME J. 4:465–471.
Tedersoo L, Bahram M, Toots M, Diédhiou AG, Henkel TW, Kjøller R, Morris MH, Nara K, Nouhra E, Peay KG, Põlme S, Ryberg M, Smith ME, Kõljalg U. 2012. Towards global patterns in the diversity and community structure of ectomycorrhizal fungi. Mol. Ecol. 21:4160–4170.
Tedersoo L, Bahram M, Põlme S, et al. 2014. Global diversity and geography of soil fungi. Science 346:1256688.
Templeton JEL, Brotherton PM, Llamas B, Soubrier J, Haak W, Cooper A, Austin JJ. 2013. DNA capture and next-generation sequencing can recover whole mitochondrial genomes from highly degraded samples for human identification. Investig Genet. 4:26.
Tiffney BH. 1985. Perspectives on the origin of the floristic similarity between Eastern Asia and eastern North America. J. Arnold Arbor. 66:73–94.
Tiffney BH, Manchester SR. 2001. The use of geological and paleontological evidence in evaluating plant phylogeographic hypotheses in the Northern Hemisphere Tertiary. Int. J. Plant Sci. 162:S3–S17.
Tulloss RE. 2015. Studies in the Amanitaceae. http://www.amanitaceae.org/ Accessed on: Feb 28, 2015.
Tulloss RE, Halling RE, Mueller GM. 2011. Studies in Amanita (Amanitaceae) of Central America. 1. Three new species from Costa Rica and Honduras. Mycotaxon 117:165–205.
Tulloss RE. 2005. Amanita-distribution in the Americas, with comparison to eastern and southern Asia and notes on spore character variation with latitude and ecology. Mycotaxon 93:189–231.
Tulloss, R.E. (1998) Provisional world key to species closely related to Amanita hemibapha with notes on the slender Caesar’s mushrooms of eastern North America. McIlvainea, 13, 46–53.
137
Tulloss RE, Yang Z-L. 2014. Section Caesareae. Studies in the Amanitaceae. Available at: http://amanitaceae.org/?section+Caesareae.
Tulloss RE, Ovrebo CL, Halling RE. 1992. Studies on Amanita (Amanitaceae) from Andean Colombia. Mem. NYBG 66:1–46.
Turner TL, Bourne EC, von Wettberg EJ, Hu TT, Nuzhdin SV. 2010. Population resequencing reveals local adaptation of Arabidopsis lyrata to serpentine soils. Nat. Genet. 42:260–263.
Tzedakis PC, Raynaud D, McManus JF, Berger A, Brovkin V, Kiefer T. 2009. Interglacial diversity. Nature Geosci. 2:751–755.
van der Nest MA, Beirn LA, Crouch JA, Demers JE, de Beer ZW, De Vos L, Gordon TR, Moncalvo J-M, Naidoo K, Sanchez-Ramirez S, Roodt D, Santana QC, Slinski SL, Stata M, Taerum SJ, Wilken PM, Wilson AM, Wingfield MJ, Wingfield BD. 2014. IMA Genome-F 3: Draft genomes of Amanita jacksonii, Ceratocystis albifundus, Fusarium circinatum, Huntiella omanensis, Leptographium procerum, Rutstroemia sydowiana, and Sclerotinia echinophila. IMA Fungus. 5:472–485.
Verrelli BC, McDonald JH, Argyropoulos G, Destro-Bisol G, Froment A, Drousiotou A, Lefranc G, Helal AN, Loiselet J, Tishkoff SA. 2002. Evidence for balancing selection from nucleotide sequence analyses of human G6PD. Am. J. Hum. Genet. 71:1112–1128.
Vetter J. 1998. Toxins of Amanita phalloides. Toxicon 36:13–24. Vilgalys R, Sun BL. 1994. Ancient and recent patterns of geographic speciation in the oyster
mushroom Pleurotus revealed by phylogenetic analysis of ribosomal DNA sequences. Proc. Natl. Acad. Sci. USA 91:4599–4603.
Vilgalys R, Hester M. 1990. Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptoccocus species. J. Bacteriol. 172:4238–4246.
Villanueva-Jiménez E, Villegas-Ríos M, Cifuentes-Blanco J, León-Avendaño H. 2006. Diversidad del género Amanita en dos áreas con diferentes condiciónes silvícolas en Ixtlán de Juárez, Oaxaca, México. Rev. Mex. Biodiv. 77:17–22.
Vo!í"ková J, Brabcová V, Cajthaml T, Baldrian P. 2014. Seasonal dynamics of fungal communities in a temperate oak forest soil. New Phytol. 201:269–278.
Waltari E, Hijmans RJ, Peterson AT, Nyári ÁS, Perkins SL, Guralnick RP. 2007. Locating Pleistocene refugia: comparing phylogeographic and ecological niche model predictions. PLoS One 2:e563.
Wang C-W. 1961. The forests of China with a survey of grassland and desert vegetation. Maria Moors Cabot Foundation, Harvard University, Cambridge, MA.
Watanabe M, Suzuki T, O’ishi R, Komuro Y, Watanabe S, Emori S, Takemura T, Chikira M, Ogura T, Sekiguchi M, Takata K, Yamazaki D, Yokohata T, Nozawa T, Hasumi H, Tatebe H, Kimoto M. 2010. Improved climate simulation by MIROC5: mean states, variability, and climate sensitivity. J Clim. 23:6312–6335.
Watterson GA. 1975. On the number of segregating sites in genetical models without recombination. Theor. Popul. Biol. 7:256–276.
138
Webb TIII, Bartlein PJ. 1992. Global changes during the last 3 million years: climatic controls and biotic responses. Annu. Rev. Ecol. Syst. 23:141–173.
Weir JT, Schluter D. 2004. Ice sheets promote speciation in boreal birds. Proc. R. Soc. B 271:1881–1887.
Weir JT, Schluter D. 2007. The latitudinal gradient in recent speciation and extinction rates of birds and mammals. Science 315:1574–76.
Weiß M, Yang Z-L, Oberwinker F. 1998. Molecular phylogenetic studies in the genus Amanita. Can. J. Bot. 76:1170–1179.
Wen J. 1999. Evolution of eastern Asian and eastern North American disjunct distributions in flowering plants. Annu. Rev. Ecol. Syst. 30:421–455.
Wiens JJ. 2007. Species delimitation: new approaches for discovering diversity. Syst. Biol. 56:875–878.
Wiens JJ. 2011. The causes of species richness patterns across space, time, and clades and the role of “ecological limits”. Quart. Rev. Biol. 86:75–96.
Wiens JJ, Donoghue MJ. 2004. Historical biogeography, ecology and species richness. Trends. Ecol. Evol. 19:639–644.
Williamson RJ, Josephs EB, Platts AE, Hazzouri KM, Haudry A, Blanchette M, Wright SI. 2014. Evidence for widespread positive and negative selection in coding and conserved noncoding regions of Capsella grandiflora. PLoS Genet. 10:e1004622.
Wilson AW, Binder M, Hibbett DS. 2012. Diversity and evolution of ectomycorrhizal host associations in the Sclerodermatineae (Boletales, Basidiomycota). New Phytol. 194:1079–1095.
Wolfe BE, Tulloss RE, Pringle A. 2012. The irreversible loss of a decomposition pathway marks the single origin of an ectomycorrhizal symbiosis. PLoS One. 7:e39597.
Wolfe JA. 1975. Some aspects of plant geography of the northern hemisphere during the Late Cretaceous and Tertiary. Ann. Miss. Bot. Gar., 62:264–279.
Wolfe JA. 1980. Tertiary climates and floristic relationships at high latitudes in the Northern Hemisphere. Palaeogeogr. Palaeoclimatol. Palaeoecol. 30:313–323.
Wood AE. 1997. Studies in the genus Amanita (Agaricales) in Australia. Australian Systematic Botany 10:723–854.
Wright SI, Gaut BS. 2005. Molecular population genetics and the search for adaptive evolution in plants. Mol. Biol. Evol. 22:506–519.
Xiang Q-Y, Soltis DE, Soltis PS, Manchester SR, Crawford DJ. 2000. Timing the eastern Asian–eastern North American floristic disjunction: molecular clock corroborates paleontological data. Mol. Phylogenet. Evol. 15:462–472.
Yang Z, Nielsen R. 2000. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol. Biol. Evol. 17:32–43.
Yang Z, Nielsen R. 2002. Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol. Biol. Evol. 19:908–917.
139
Yang Z, Rannala B. 2010. Bayesian species delimitation using multilocus sequence data. Proc. Natl. Acad. Sci. USA 107:9264–9269.
Yang Z, Rannala B. 2014. Unguided species delimitation using DNA sequence data from multiple loci. Mol. Biol. Evol. 31:3125–3135.
Yang Z, Wong WSW, Nielsen R. 2005. Bayes empirical bayes inference of amino acid sites under positive selection. Mol. Biol. Evol. 22:1107–1118.
Yang Z. 1998. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol. Biol. Evol. 15:568–573.
Yang Z. 2007. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Mol. Biol. Evol. 24:1586–1591.
Yang ZL. 2011. Molecular techniques revolutionize knowledge of basidiomycete evolution. Fungal Divers. 50:47–58.
Z-L. 1997. Die Amanita-Arten von Sudwestchina. Bibliotheca Mycologica 170. J. Cramer, Berlin, Germany.
Yu Y, Harris AJ, Blair C, He X. 2015. RASP (Reconstruct Ancestral State in Phylogenies): a tool for historical biogeography. Mol. Phylogenet. Evol. 87:46–49.
Yu Y, Harris AJ, He X-J. 2011. RASP (reconstruct ancestral state in phylogenies) 1.1. Sichuan University, Chengdu. Available at: http://mnh.scu.edu.cn/soft/blog/RASP.
Zarza E, Reynoso VH, Emerson BC. 2008. Diversification in the northern neotropics: mitochondrial and nuclear DNA phylogeography of the iguana Ctenosaura pectinata and related species. Mol. Ecol. 17:3259–3275.
Zdobnov EM, Apweiler R. 2001. InterProScan –an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 17:847–848.
Zhou L-L, Su YCF, Thomas DC, Saunders RMK. 2011. ‘Out-of-Africa’ dispersal of tropical floras during the Miocene climatic optimum: evidence from Uvaria (Annonaceae). J. Biogeogr. 39:322–335.
Zolan ME, Pukkila PJ. 1986. Inheritance of DNA methylation in Coprinus cinereus. Mol. Cell. Biol. 6:195–2.
density (HPD) intervals. White circles with letters indicate the mean age where fossil constrains were placed. Name clades were constrained as monophyletic. Only nodes with posterior probabilities greater than 0.8 are annotated. The density plot in on the top-left
corner shows age estimates from four runs on the most recent common ancestor (MRCA) of Amanita subg. Amanita.
157
!?#$%&'$(/3+.0
!;#$<=>&$./.*,
!5#$>5$.03
!5#$>5$.03*
!5#$>5$)*)
!5#$%&'$(,/+*
!5#$67$8797:$(/4*
!5#$%&'$(,)+1
!5#$%&'$(,)+(
!5#$@FCUV$304T*4
!5#$%&'$(,3+4
!-#$%&'$.(*+.0!-#$%&'$1)(+/!-#$%&'$(/,+*
!;#$<=>&$@CDEF7G$@(1
!;#$@CDEF7G$@()
!;#$<=>&$@CDEF7G$@,,!;#$%&'$130+3
!<#$%&'$0(3+1
!<#$%&'$.,1+*!<#$%&'$,14+..
!<#$%&'$,/0+(
!"#$%&'$0*1+3!"#$%&'$()/+(
!"#$%&'$()/+,
!"#$%&'$()/+4
!"#$%&'$()*+,
!"#$%&'$()*+/
!"#$%&'$()*+*
!?#$%&'$(.,+/!?#$%&'$(.4+/
!;#$<=>&$@CDEF7G$@(.
!;#$%&'$041+/
!;#$%&'$,*4$.0
!;#$%&'$,*)+.
!8#$%&'$(((+3
!6#$%&'$(,1+)
!?#$?N-$"NB?$/30!?#$%&'$1/4+.0
!?#$%&'$1/)+(
!?#$%&'$(,*+/
!2#$%&'$,/3+4
!?#$'%'=$./01)3
!?#$'%'=$./0(.,
!?#$'%'=$./0,11!?#$'%'=$./41/0
!?#$'%'=$.3..(,
!?#$'%'=$.3..3,!?#$'%'=$.3..4.
!?#$%&'$1/)+.
!-#$%&'$.0*+,
!-#$%&'$./,+.0
!-#$%&'$(./+)!-#$%&'$(/,+)
!-#$%&'$(*(+3
!-#$%&'$(*(+4!-#$%&'$(*(+)
!;#$<=>&$1.//0!;#$<=>&$1.3/1
!-#$<..(1.))
!;#$%&'$1/1+*!;#$A5H$=FCEID$/)*/
!5#$%&'$(,)+*
!5#$%&'$(,,+/
!&#$%&'$()(+.
!-#$%&'$(4(+*
!=#$"C::JDK$3).,
!=#$"C::JDK$3)./
!-#$%&'$.()+.!-#$%&'$.()+1
!-#$%&'$(3/+.
!-#$6&LM$,00!-#$6&LM$*/*4
!"#$%&'$0),+4!"#$%&'$1).+*!"#$%&'$()/+1
!"#$%&'$()/+3!"#$%&'$()/+)
!?#$;A**)04.T;A**)0.,
!-#$8=,1,/,3
!"#"$%
5<%$!5#
@&5$!?#
5B@$!=#
8&%$!6#
S8-$!&#
&B%$!<#
;5N$!2#
NOL$!"#
&O5$!-#
>&A$!;#
=%$!8#
SQPKUIQV$!H#
=IDPJDQ7R =IDPJDQ7R
158
Figure A8.2. Bayesian binary Markov chain monte carlo (BBM) results from rasp analysis
of Amanita sect. Caesareae (dataset 2) showing the most likely reconstructed areas (left)
and pie-charts with area posterior probabilities (right). The model used in rasp was F81+G.
!"#$,C=2;;;)
!"#$,C=2;;)*!"#$,C=2;;F2
!?#$%&'$(E2+E
!"#$'%'4$<)=;=F
!"#$'%'4$<)=;<=!"#$'%'4$<)=;<(
!>#$%&'$2)<+;
!"#$%&'$()3+)
!"#$%&'$()*+*!"#$'%'4$<)=(=*
!G#$'%'4$<)=(2)
!>#$%&'$=(=+;!>#$%&'$=F)+)
!D#$%&'$=*;+)
!?#$@&AB$2;(<;!?#$%&'$23(+3
!"#$B'4$HPK.7P0Q$<);R!4#$%&'$<(F+3
!>#$%&'$23E+<
!>#$%&'$(3;+2
!>#$4HB4$IJK.L5$;E=*"
!>#$%&'$=*F+<!>#$%&'$<2;+<=
!4#$%&'$;3F+;
!>#$%&'$(E(+2
!D#$%&'$2*(+(!D#$%&'$2*(+;
!"#$%&'$()=+<!"#$%&'$()<+3!"#$%&'$()F+)
!>#$%&'$2*F+<=!>#$%&'$2*F+2
!>#$%&'$(<)+<=
!D#$%&'$2*2+<
!>#$%&'$2FE+<=
!>#$%&'$(F<+F!>#$%&'$(3;+E
!"#$%&'$2)3+E
!,#$-.$/.0.1$(223!,#$-.$/.0.1$(3<;
!,#$-.$/.0.1$;(3E!,#$@,$EF<
!,#$%&'$2F<+(!,#$%&'$(;(+<=
!G#$&B=3<**=M&B=3<E3)
!C#$NO*2)(FE
!"#$DA**E=F3MDA**E=<(
4567869.: 4567869.:
159
Appendix B: Chapter 3
Figure B8.3. The set of credible diversification rate shifts from BAMM based on a Bayes
Factor criterion of 3. Each plot indicates the sampled frequency of each model and time-
varying speciation rate estimates. Colors are on a temperature scale (see Fig. 3.1).
f = 0.69 f = 0.078 f = 0.048
f = 0.029 f = 0.026 f = 0.025
f = 0.025 f = 0.02 f = 0.015
160
Figure B8.4. Boxplots of speciation and extinction rates in temperate and tropical subclades
based on the second best QuaSSE model (model 2 in Table 3.1) for the clade ‘Caesarea’
data set. In this model, speciation and extinction rates are constant functions of latitude.
!"#$%&'%()*+("%$&,
-&'#./,%)#&0#1.23+456
!"#$%&'%()*#7"#+&'#
89'%)$'%()*+("%$&,
89'%)$'%()*#7"#+&'#
:;::
:;:<
:;5:
:;5<
:;=:
161
Appendix C: Chapter 5
Table C8.3. Species, specimens, geographic locations, and GenBank accession numbers for DNA data used in Chapter 5.
Coordinates Accessions
Species Voucher Location latitude longitude tef1 AJ0000-2 AJ0132-103
Amanita vernicoccora RET 281-9 USA: Washington, Klickitat Co., North Fork Spring Creek Rd 45.807 -121.497 KP724683 KP724127 KP724403
194
Figure C8.5. *BEAST species trees used for species delimitation with clade posterior
probabilities (upper), BP&P species probabilities (lower and bold), and branching times.
!"#$%&''
!"#!"#$()(*'&+%,%%
!"#!"#$&''()*#$%&''%
!"#!"#-&''()*#*%.%/.0.1
!"#*%.%/.0.23%4%$'5)%,%
!"#%46%,&%,%
!"#!"+,-.
!"#(%/.754)8+49%
!"#!"+/&012
!"#!"+/&013
!"#!"+/&014
!"#!"+/&015
!"#:%(6&),''
!"#!"+/&01.
!"#!"+6..
!"#!"+74.
.*8!"#
8*24!"#
8*92!"#
8*:4!"#
8*;3!"#
8*::!"#
8*94#"$
.*8!"#
8*;9#"%%
8*:2!"#
.*8!"#
8*2;!"#
.*8!"#
8*9;!"#
8*9..*.5.*5<
5*2;
4*399*.39*35
.*49
.*9;
5*:.
4*:2
2*9:2*;43*4.
&'()**+,-.(/0 &1(-2)34**+,-.(/07(=>#?,&@
195
Figure C8.6. Regression analyses on the latitude versus population size and expansion data.
The blue line represents the linear model, while the green and red represent the Gaussian
and polynomial models, respectively.
!" #$ #" %$ %" &$ &"
$!
#%
&"
'()*+,-./01/202-3(,+0*/452(*6+0*
!" #$ #" %$ %" &$ &"
#&
78
9(,+,-.4
:02-3(,+0*/6+;4
196
Appendix A: Chapter 6
Table D8.4. Samples used in Chapter 6 for exon-targeted sequencing.
Species Voucher Location latitude longitude Amanita jacksonii 1Q Canada, Quebec, Boucherville 45.365 -73.423 Amanita jacksonii 2Q Canada, Quebec, Boucherville 45.591 -73.436 Amanita jacksonii 3Q Canada, Quebec, Ste. Anne de Bellevue 45.403 -73.950 Amanita jacksonii 4Q Canada, Quebec, Ste Anne de Bellevue 45.403 -73.950 Amanita jacksonii AW200 Canada, Ontario, Awenda Prov Park 44.844 -80.006 Amanita jacksonii AW201 Canada, Ontario, Awenda Prov. Park 44.844 -80.006 Amanita jacksonii AW203 Canada, Ontario, Awenda Prov Park 44.844 -80.006 Amanita jacksonii AW217 Canada, Ontario, Awenda Prov Park 44.844 -80.006 Amanita jacksonii AW228 Canada, Ontario, Awenda Prov Park 44.844 -80.006
Amanita jacksonii DEWV 1893 USA, West Virginia, Monongahela National Forest, Blue Bend 37.923 -80.267
Amanita jacksonii DEWV 2410 USA, West Virginia, Barbour, Audra State Park 39.041 -80.069
Amanita jacksonii DEWV 4715 USA, West Virginia, Upshur, French Creek, State Wildlife Center 38.886 -80.297
Amanita jacksonii DEWV 7136 USA, West Virginia, Blackwater Falls State Park 39.109 -79.492
Amanita jacksonii Feldman 20-07-13 USA, Pennsylvania, Mifflinburg, Raymond B. Winter State Park 40.992 -77.186
Amanita jacksonii RET 109-4 USA, Massashusets, Franklin 42.083 -71.397 Amanita jacksonii RET 154-10 USA, New York, Ulster 41.859 -74.312 Amanita jacksonii RET 315-9 USA, Connecticut, Mansfield, Storrs 41.808 -72.250 Amanita jacksonii RET 393-8 USA, Norh Carolina, Haywood 35.624 -82.993 Amanita jacksonii RET 489-5 USA, Connecticut, Middlesex Co, Devils Hopyard 41.476 -72.340 Amanita jacksonii RET 526-5 USA, Pennsylvania, Carbon Co, Hickory Run SP 41.036 -75.684 Amanita jacksonii RET 547-4 USA, Pennsylvania, Erie 42.086 -80.116 Amanita jacksonii RET 551-6 USA, Pennsylvania, Luzerne Co., Ricketts Glen Park 41.330 -76.291 Amanita jacksonii Thorn ON 2009 Canada, Ontario, Exeter 43.346 -81.480 Amanita sp-F11 RET 537-10 USA, Florida, Lee Co., Alva 26.797 -81.671
197
Species Voucher Location latitude longitude Amanita sp-F11 RET 537-5 USA, Florida, Lee Co., Alva 26.716 -81.610 Amanita sp-F11 RET 547-10 USA, Florida, Marion, Ocala Natl. Forest 29.241 -81.689 Amanita sp-F11 RET 548_8 USA, Florida, Marion, Ocala Natl. Forest 29.258 -81.779 Amanita sp-jack1 RV DUKE AJ20 USA, Norh Carolina, Durham 35.994 -78.899 Amanita sp-jack1 RV-DUKE-AJ10 USA, Norh Carolina, Durham, Duke Forest 36.021 -78.983 Amanita sp-jack1 TENN 57937 USA, North Carolina, Macon 36.439 -78.084 Amanita sp-jack1 TENN 66058 USA, Tennessee, Blount 35.671 -83.847 Amanita sp-jack2 Amanita 2004 Mexico, Chiapas, La Trinitaria, Lagunas de Montebello 16.757 -93.129 Amanita sp-jack2 F1132183 USA, Texas, Orange, Vidor 30.132 -94.015 Amanita sp-jack2 JM DUKE AJ18 USA, Norh Carolina, Durham 35.994 -78.899 Amanita sp-jack2 R. Doss 27708 USA, Arkansas 35.201 -91.832 Amanita sp-jack2 RET 354-8 USA, Texas, Newton, Bleakwood 30.692 -93.822 Amanita sp-jack2 RET 463_6 USA, Texas, Hardin , Big Thicket Natl. Preserve 30.474 -94.346 Amanita sp-jack2 RET 531-4 USA, Indiana, Lawrence, Lake Monroe 39.070 -86.445 Amanita sp-jack2 TENN 55460 USA, South Carolina, Oconee 34.749 -82.993 Amanita sp-jack2 TENN 61260 USA, Tennessee, Knoxville 35.961 -83.921 Amanita sp-jack3 ABM042 Mexico, Hidalgo, San Miguel El Cerezo 20.160 -98.728 Amanita sp-jack3 ABM044 Mexico, Hidalgo, San Miguel El Cerezo 20.171 -98.718 Amanita sp-jack3 ABM049 Mexico, Hidalgo, Acaxochitlan 20.158 -98.200 Amanita sp-jack3 ABM051 Mexico, Hidalgo, Acaxochitlan 20.172 -98.199 Amanita sp-jack3 RET 252-9 Mexico, Hidalgo, Zacaultipan 20.646 -98.655 Amanita sp-jack3 RS_364 Mexico, Chiapas, Chamula 16.787 -92.688 Amanita sp-jack3 V.V.R. 01 Mexico, Hidalgo, Zacaultipan, La Mojonera 20.643 -98.655 Amanita sp-jack5 FCME 12243 Mexico, Guerrero, Tixtla de Guerrero 17.556 -99.379 Amanita sp-jack5 FCME 5528 Mexico, Guerrero, Taxco 18.626 -99.672 Amanita sp-jack5 Martinez 19 Mexico, Oaxaca, Magdalena Mixtepec (Zaachila) 16.350 -96.603 Amanita sp-jack5 R y V 1090 Mexico, Oaxaca, San Pedro Tidaa 17.340 -97.379 Amanita sp-jack5 SSRH_03_B Mexico, Oaxaca, San Sebastian Rio Hondo 16.182 -96.464 Amanita sp-jack6 F. 241 Mexico, Jalisco, Sierra de Quila 20.278 -104.085
198
Species Voucher Location latitude longitude Amanita sp-jack6 FCME 17782 Mexico, Chihuahua 27.841 -107.585 Amanita sp-jack6 FCME 3859 Mexico, Nayarit, Tepic, Reserva Ecologica de San Juan 21.499 -104.931 Amanita sp-jack6 H. 191 Mexico, Jalisco, San Sebastian 20.366 -102.966 Amanita sp-jack6 H. 422 Mexico, Jalisco, Huaxtla, Rio Salado 20.727 -103.658
Amanita sp-jack6 SS MOR MX 2013 2 Mexico, Michoacan, Pino Real 19.650 -101.013
Amanita sp-T31 F1132188 USA, Texas, Angelina, Zavalla 31.159 -94.426 Amanita sp-T31 RET 453-1 USA, Missouri, St. Genevieve , Hawn St. Park 37.816 -90.253 Amanita sp-T31 RET 516-3 USA, Missouri, St. Genevieve , Hawn St. Park 37.816 -90.253 Amanita sp-T31 RET 552-4 USA, Missouri, St. Genevieve , Hawn St. Park 37.816 -90.253
199
Table D8.5. Genes used for demographic analyses, including a BLAST homolog, their functional annotations, as well as results from gsi and PAML analysis.
Linkage group Gene Best BLAST hit (Accession) Species Functional annotation GO terms ! vargsi "gsi LRT
scaffold-0 g0013 KIL67414.1 A. muscaria kinesin-like protein GO:0008017 0.186 0.000 8.000 Neu
scaffold-2 g0041 KIL69537.1 A. muscaria CMGC/DYRK domin containing protein GO:0004672 0.263 0.000 8.000 Pur
scaffold-2 g0082 KIL69598.1 A. muscaria protein prenyltransferase activity GO:0008318 0.234 0.000 8.000 Neu scaffold-8 g0162 KIL68638.1 A. muscaria hect-domain-containing protein GO:0004842 0.115 0.000 8.000 Neu scaffold-16 g0274 KIL62189.1 A. muscaria acyltransferase GO:0016747 0.344 0.000 8.000 Div
scaffold-31 g0464 KIL66181.1 A. muscaria proteophosphoglycan PPG4 protein GO:0008137 0.327 0.000 8.000 Neu
scaffold-31 g0472 KIL66211.1 A. muscaria ubiquitin-protein ligase GO:0004842 0.184 0.000 8.000 Pur scaffold-32 g0530 KIL71094.1 A. muscaria bzip-family transcription factor GO:0003700 0.263 0.000 8.000 Neu scaffold-81 g0828 KIL65846.1 A. muscaria chitin synthethase GO:0004100 0.099 0.001 7.924 Pur scaffold-102 g0981 KIL61570.1 A. muscaria phosphoketolase GO:0016832 0.197 0.000 8.000 Neu
scaffold-103 g1032 KIL65900.1 A. muscaria hypothetical protein, protein binding GO:0005515 0.404 0.000 8.000 Pur
scaffold-106 g1045 KIL69738.1 A. muscaria glycosyltransferase famlily 48 GO:0003843 0.115 0.000 8.000 Neu scaffold-106 g1074 KIL69785.1 A. muscaria STE/STE11 protein kinase GO:0004674 0.154 0.000 8.000 Pur
scaffold-106 g1111 KIL69836.1 A. muscaria kex protein, serine-type endopeptidase activity GO:0004252 0.248 0.000 8.000 Neu
scaffold-106 g1118 KIL69850.1 A. muscaria glutamate-tRNA ligase GO:0004818 0.342 0.000 8.000 Pur scaffold-106 g1153 KIL69899.1 A. muscaria vacuolar ATP synthase subunit D GO:0042626 0.093 0.000 8.000 Pur
scaffold-109 g1225 KIL70700.1 A. muscaria P-loop containing nucleoside triphosphate hydrolase GO:0004386 0.215 0.000 8.000 Neu
scaffold-116 g1367 KIL68906.1 A. muscaria hypothetical protein (palmitoyl-(protein) hydrolase activity) GO:0008474 0.386 0.000 8.000 Neu
scaffold-291 g1957 KIL63265.1 A. muscaria hypothetical protein NA 0.538 0.001 7.924 Neu
200
Linkage group Gene Best BLAST hit (Accession) Species Functional annotation GO terms ! vargsi "gsi LRT
scaffold-587 g3041 KIL67387.1 A. muscaria hypothetical protein (metal ion binding) GO:0046872 0.275 0.001 7.923 Pur
scaffold-663 g3502 KIL67827.1 A. muscaria metallopeptidase GO:0004222 0.162 0.000 8.000 Neu scaffold-666 g3586 KIL69438.1 A. muscaria ring finger ubiquitin ligase GO:0008270 0.230 0.000 8.000 Pur scaffold-666 g3629 KIL69372.1 A. muscaria carnitine acyl transferase GO:0016746 0.316 0.003 7.861 Neu scaffold-675 g3712 KIL69708.1 A. muscaria GTP binding GO:0005488 0.141 0.000 8.000 Neu
scaffold-685 g3819 KIL69245.1 A. muscaria von Willebrand ring finger domain-containing protein GO:0008270 0.126 0.003 7.861 Neu
scaffold-698 g4062 KIL68186.1 A. muscaria hypothetical protein NA 0.298 0.004 7.843 Pur scaffold-706 g4167 KIL71593.1 A. muscaria chitin synthethase GO:0004100 0.053 0.000 8.000 Pur scaffold-706 g4215 KIL71522.1 A. muscaria WD40 domain-containing protein GO:0005515 0.273 0.003 7.862 Neu
scaffold-723 g4374 KIL70861.1 A. muscaria hypothetical protein (GTPase activity) GO:0003924 0.326 0.003 7.803 Neu
scaffold-726 g4486 KJA15191.1 Hypholoma sublateritium kinase-like protein GO:0004674 0.291 0.000 8.000 Neu
scaffold-732 g4599 KIL62230.1 A. muscaria DNA mismatch repair protein GO:0030983 0.205 0.001 7.913 Neu
scaffold-742 g4668 KIL69515.1 A. muscaria pumilio RRM domin-containing protein GO:0000166 0.166 0.000 8.000 Pur
scaffold-757 g4744 KIL59086.1 A. muscaria AMP deaminase GO:0003876 0.188 0.004 7.843 Neu scaffold-831 g5023 KIL68278.1 A. muscaria growth-arrest protein GO:0007050 0.291 0.000 8.000 Neu
scaffold-871 g5212 KIL70662.1 A. muscaria chromatin remodeling complex protein NA 0.161 0.003 7.867 Neu
scaffold-8177 g7424 KIL70086.1 A. muscaria ubiquitin conjugating enzyme GO:0016881 0.131 0.000 8.000 Neu scaffold-9363 g7918 KIL70545.1 A. muscaria kinase-like protein GO:0004674 0.211 0.000 8.000 Neu
scaffold-9363 g7991 KIL68102.1 A. muscaria cellular morphogenesis-related protein NA 0.280 0.000 8.000 Neu
scaffold-9461 g8575 KIL70176.1 A. muscaria vacuolar sorting protein NA 0.318 0.000 8.000 Neu scaffold-9463 g8604 KIL71356.1 A. muscaria glycosyltransferase famlily 2 GO:0016758 0.091 0.000 8.000 Pur
201
Figure D8.7. Number of variable sites (grey) and percentage of missing data (black) versus
gene length.
!""" #""" $""" %""" &"""" &!""" &#"""
"!""
#""
$""
%""
&"""
&!""
&#""
'()(*+()',-*./01
2*34564/+(*76,(7
"&"
!"8"
#"9"
$"
:*;6776)'
202
Figure D8.8. The Watterson estimator, theta, separated by site classes.
!"#$%&'!()
*+**
*+*,
*+*-
*+*.
*+*/
0%1234566
3780%12-
3780%129
3780%12:
3780%12.
3780%12,
378;,,
378!.,
<453=5>=5?5$@45
203
Figure D8.9. Estimates of nucleotide diversity over all genes separated by class sites.
!"#$%&'()%*)(+%,-('.*/!0
1211
1213
1214
1215
1216
78#9-&:((
-;<78#94
-;<78#9=
-;<78#9>
-;<78#95
-;<78#93
-;<?33
-;<@53
!&:-.:A.:B:',&:
204
Figure D8.10. Estimates of effective population size.
!"#$%&'((
%)*!"#$+
%)*!"#$,
%)*!"#$-
%)*!"#$.
%)*!"#$/
%)*0//
%)*1./
2
/3/2,
+3/2,
.3/2,
43/2,
-3/2,
5667#8(97:)&);<"8(&':%(=7
205
Figure D8.11. Scatterplots of gsi values versus omega (polymorphism).
!"!
!"#
!"$
!"%
!"&
'"!
!"#$%&'((
!"#
%)*!"#$+
!"!
!"#
!"$
!"%
!"&
'"!
%)*!"#$, %)*!"#$-
!"!
!"#
!"$
!"%
!"&
'"!
%)*!"#$. %)*!"#$/
!"! !"( '"! '"( #"!
!"!
!"#
!"$
!"%
!"&
'"!
%)*0//
!)
!"! !"( '"! '"( #"!
%)*1./
206
Figure D8.12. Scatterplots of gsi values versus alpha.
!"!
!"#
!"$
!"%
!"&
'"!
!"#$%&'((
!"#
%)*!"#$+
!"!
!"#
!"$
!"%
!"&
'"!
%)*!"#$, %)*!"#$-
!"!
!"#
!"$
!"%
!"&
'"!
%)*!"#$. %)*!"#$/
('! (& (% ($ (# !
!"!
!"#
!"$
!"%
!"&
'"!
%)*0//
('! (& (% ($ (# !
%)*1./
!
207
Figure D8.13. Neutrality Index versus the scaled number of Stop codons.
0 2 4 6 8 10
05
1015
20
scaled number of Stop codons per gene
NI
208
A. jacksonii A. sp-jack2
A. sp-F11
A. sp-jack5
A. sp-jack6
A. sp-jack1
A. sp-jack3
A. sp-T31
209
Figure D8.14. Wordles of GO terms for each species. These represent genes that were found under positive selection (alpha > 0), but were species-specific.