ABSTRACT
Mountains of the Balkan Peninsula are significant biodiversity hotspots with great species richness and a large proportion of narrow endemics. Processes that have driven the evolution of the rich Balkan mountain flora, however, are still insufficiently explored and understood. Here we focus on a group of Cardamine (Brassicaceae) perennials growing in wet, mainly mountainous habitats. It comprises several Mediterranean endemics, including those restricted to the Balkan Peninsula. We used target enrichment with genome skimming (Hyb-Seq) to infer their phylogenetic relationships, and, along with genomic in situ hybridization (GISH), to resolve the origin of tetraploid Cardamine barbaraeoides endemic to the Southern Pindos Mts. (Greece). We also explored the challenges of phylogenomic analyses of polyploid species and developed a new approach of allele sorting into homeologs that allows identifying subgenomes inherited from different progenitors. We obtained a robust phylogenetic reconstruction for diploids based on 1,168 low-copy nuclear genes, which suggested both allopatric and ecological speciation events. In addition, cases of plastid-nuclear discordance, in agreement with divergent nuclear ribosomal DNA (nrDNA) copy variants in some species, indicated traces of interspecific gene flow. Our results also support biogeographic links between the Balkan and Anatolian-Caucasus regions and illustrate the contribution of the latter region to high Balkan biodiversity. An allopolyploid origin was inferred for C. barbaraeoides, which highlights the role of mountains in the Balkan Peninsula both as refugia and melting pots favoring species contacts and polyploid evolution in response to Pleistocene climate-induced range dynamics. Overall, our study demonstrates the importance of a thorough phylogenomic approach when studying the evolution of recently diverged species complexes affected by reticulation events at both diploid and polyploid levels. We emphasize the significance of retrieving allelic and homeologous variation from nuclear genes, as well as multiple nrDNA copy variants from genome skim data.
ABSTRACT
BACKGROUND AND AIMS: Diploid representatives from the related polyploid complexes of Cardamine amara, C. pratensis and C. raphanifolia (Brassicaceae), were studied to elucidate phylogenetic relationships among the complexes and among the individual taxa included. METHODS: Two independent molecular data sets were used: nucleotide sequences from the internal transcribed spacers (ITS) of nrDNA, and amplified fragment length polymorphism (AFLP) markers. Seventeen diploid taxa from the studied groups were sampled. KEY RESULTS: Both ITS and AFLP analyses provided congruent results in inferred relationships, and revealed two main lineages. While the C. amara group, consisting of C. wiedemanniana and four subspecies of C. amara, was resolved as a well-supported monophyletic group, taxa from the C. pratensis and C. tenera groups (the latter representing diploid taxa of the complex of C. raphanifolia) all appeared together in a single clade/cluster with no support for the recognition of either of the groups. Intra-individual polymorphisms and patterns of nucleotide variation in the ITS region in C. uliginosa and C. tenera, together with the distribution of AFLP bands, indicate ancient hybridization and introgression among these Caucasian diploids. CONCLUSIONS: The lack of supported hierarchical structure suggests that extensive reticulate evolution between these groups, even at the diploid level, has occurred (although an alternative explanation, namely ancestral polymorphism in ITS data, cannot be completely excluded). Several implications for the investigation of the polyploid complexes of concern are drawn. When tracing origins of polyploid taxa, a much more complex scenario should be expected, taking into account all relatives as potential parents, irrespective of the group in which they are classified.
