Polyploidy is a major speciation process in vascular plants, and is postulated to be particularly important in shaping the diversity of extant ferns. However, limitations in the availability of bi-parental markers for ferns have greatly limited phylogenetic investigation of polyploidy in this group. With a large number of allopolyploid species, the genus Botrychium is a classic example in ferns where recurrent polyploidy is postulated to have driven frequent speciation events. Here, we use PacBio sequencing and the PURC bioinformatics pipeline to capture all homeologous or allelic copies of four long (∼1kb) low-copy nuclear regions from a sample of 45 specimens (25 diploids and 20 polyploids) representing 37 Botrychium taxa, and three outgroups. This sample includes most currently recognized Botrychium species in Europe and North America, and the majority of our specimens were genotyped with co-dominant nuclear allozymes to ensure species identification. We analyzed the sequence data using maximum likelihood (ML) and Bayesian inference (BI) concatenated-data (“gene tree”) approaches to explore the relationships among Botrychium species. Finally, we estimated divergence times among Botrychium lineages and inferred the multi-labeled polyploid species tree showing the origins of the polyploid taxa, and their relationships to each other and to their diploid progenitors. We found strong support for the monophyly of the major lineages within Botrychium and identified most of the parental donors of the polyploids; these results largely corroborate earlier morphological and allozyme-based investigations. Each polyploid had at least two distinct homeologs, indicating that all sampled polyploids are likely allopolyploids (rather than autopolyploids). Our divergence-time analyses revealed that these allopolyploid lineages originated recently—within the last two million years—and thus that the genus has undergone a recent radiation, correlated with multiple independent allopolyploidizations across the phylogeny. Also, we found strong parental biases in the formation of allopolyploids, with individual diploid species participating multiple times as either the maternal or paternal donor (but not both). Finally, we discuss the role of polyploidy in the evolutionary history of Botrychium and the interspecific reproductive barriers possibly involved in these parental biases.

The moonwort genus, Botrychium s. s., includes diploid and polyploid taxa that occur primarily in the northern hemisphere. Their evolutionary history, morphologically cryptic taxa and deep divergence of the family in the phylogeny of ferns has long fascinated pteridologists. Previous molecular studies did not include a complete taxonomic sampling of the taxa in the genus, nor multiple specimens from throughout the known geographical range of each taxon. Therefore, to investigate evolutionary relationships of the major clades of Botrychium s. s., we increased both taxonomic representativeness (multiple accessions per taxa), as well as phylogenetic resolution by including additional new chloroplast markers. To confirm identification and provide evidence from both maternal and paternal parentage of allopolyploids, we also included specimens that have been characterized by allozyme profiles determined by electrophoretic analysis of 20 nuclear enzyme loci for each taxon. We analyzed four chloroplast regions (matK intron, trnHGUG —psbA, andtrnLUAA —trnFGAA intergenic spacers, and rpL16 intron region) of 365 specimens fromAsia, Europe, North America, Oceania, and South America, sampling the geographical range of 34 of 35 accepted Botrychium s. s. taxa and thirteen putatively new taxa. We conducted a phylogenetic analysis of maternal lineages based on 2,385 aligned nucleotides using maximum likelihood and Bayesian inference to explore genetic diversity and phylogenetic relationships among taxa. We found strong support for themonophyly of three major clades: Lanceolatum, Lunaria, and Simplex-Campestre, and resolved 15 subclades. Our results suggest multiple origins for at least four polyploid taxa (B. boreale, B. michiganense, B. yaaxudakeit, and B. watertonense). The Simplex-Campestre clade had the largest number of species, despite having a similar total number of haplotypes as the Lunaria clade (62 and 59, respectively), which has the broadest worldwide distribution. In total, our new molecular phylogeny comprises 47 taxa, of which thirteen are discussed for possible taxonomic recognition.