Phylogeography and paleoclimatic range dynamics explain variable outcomes to contact across a species' range.

Replicability of divergence after contact is a poorly characterized process, particularly in the contexts of phylogeography and postglacial range dynamics within species. Using contact zones located at the leading-, mid- and rear-edges of a species' range, we examined variation in outcomes to contact between divergent lineages of Campanula americana. We investigated whether contact zones vary in quantity and directionality of gene flow, how phylogeographic structure differs between contact zones, and how historic range dynamics may affect outcomes to contact. We found that all contact zones formed at similar times via primary contact yet detected significant admixture in only the rear-edge (RE) contact zone. In the northern leading-edge contact zone and the mid-range Virginia contact zone, gene flow was minimal and asymmetric. In the southern RE contact zone, gene flow was strong and symmetric. Asymmetric admixture in the leading-edge and Virginia contact zones matches the directionality of a known cosmopolitan cytonuclear incompatibility between lineages of C. americana. Our results emphasize the dependence of speciation processes on phylogeographic structure, evolutionary history and range dynamics.

Insight into the phylogeny and responses of species from the genus Sergia (Campanulaceae) to the climate changes predicted for the Mountains of Central Asia (a world biodiversity hotspot).

BACKGROUND: Together with other elevated areas, the Mountains of Central Asia are significantly threatened by ongoing climate change. The presence of refuges during the glaciations makes the region extremely rich in species, especially endemic ones. However, the limited potential for colonisation of other habitats makes rocky-related species with 'island-like' distribution, particularly vulnerable to climate change. To understand the processes underlying species response to climate warming, we assessed differences in ecological niches and phylogenetic relationship of two geographically disjunctive alpine species belonging to the genus Sergia. The taxa are considered Tertiary relicts, endemic to the Tian Shan and Pamir-Alai Mountains. To illustrate range dynamics and differences in occupied niches of Sergia species, we used Ecological Niche Modelling of current and future distribution. Whereas, to reconstruct the phylogenetic relationship within and between Sergia and other related Campanulaceae species from the region we used molecular data (ITS, cpDNA, DArTseq-derived SNPs). RESULTS: The results reveal that the genus Sergia is a polyphyletic group, and its representatives differ geographically, ecologically and genetically. Both S. regelii and S. sewerzowii constitute a common clade with Asyneuma group, however, S. sewerzowii is more closely related to Campanula alberti (a species that has never previously been considered closely related to the genus Asyneuma or Sergia) than to S. regelii. Sergia sewerzowii is adapted to lower elevations with higher temperatures, while S. regelii prefers higher elevations with lower temperatures. The future distribution models demonstrate a dramatic loss of S. regelii range with a shift to suitable habitats in higher elevations, while the potential range of S. sewerzowii increases and shifts to the north. CONCLUSIONS: This study shows that S. regelii and S. sewerzowii have a long and independent evolution history. Sergia regelii and S. sewerzowii significantly differ in realised niches. These differences are mirrored in the response of the studied endemics to future climate warming. As suitable habitats shrink, rapid changes in distribution can lead to species' range loss, which is also directly related to declines in genetic variability. The outcomes of this paper will help to more precisely assess the impact of climate changes on rocky-related plant species found in this world's biodiversity hotspot.

Maternal effects and inbreeding depression in post-translocation progeny of Campanula glomerata.

Evaluation of plant translocation success based on fitness-related quantitative traits combined with molecular markers may contribute to a finer assessment of inbreeding, selective and rescue processes, which might have long-term consequences for population dynamics and viability. We investigated fitness traits (seed germination, seedling viability, and juvenile growth and mortality) combined with 15 microsatellite loci of the first post-translocation seed progeny from two translocated populations of Campanula glomerata, an insect-pollinated, self-incompatible perennial herb. We examined whether inbreeding, heterosis through admixture, translocation site and maternal transplant seed source origin and lineage might affect seed quality and juvenile growth in controlled cultivation conditions. Flower production and seed germination of the transplants was higher in one of the two translocation sites, which might be related to differences in soil and vegetation composition and cover. Strong maternal effects related to seed source origin and lineage were found on progeny size, with the largest transplants producing the largest progeny. The differences in rosette diameter were maintained across the whole growth period measured. There was inbreeding depression (rather than heterosis) related to biparental inbreeding at the early progeny growth stage, also expressed through juvenile mortality. Our findings highlight that maternal transplant origin, especially when seed sources consisted of small, fragmented remnants, might have a selective value on fitness in the post-translocation generations. If maternal effects and inbreeding depression persist, they might affect global genetic diversity patterns in the long term. Further admixture in the next generations might buffer maternal and inbreeding effects or lead to outbreeding depression.

The First Complete Chloroplast Genome of Campanula carpatica: Genome Characterization and Phylogenetic Diversity.

Campanula carpatica is an ornamental flowering plant belonging to the family Campanulaceae. The complete chloroplast genome of C. carpatica was obtained using Illumina HiSeq X and Oxford Nanopore (Nanopore GridION) platforms. The chloroplast genome exhibited a typical circular structure with a total length of 169,341 bp, comprising a large single-copy region of 102,323 bp, a small single-copy region of 7744 bp, and a pair of inverted repeats (IRa/IRb) of 29,637 bp each. Out of a total 120 genes, 76 were protein-coding genes, 36 were transfer RNA genes, and eight were ribosomal RNA genes. The genomic characteristics of C. carpatica are similar to those of other Campanula species in terms of repetitive sequences, sequence divergence, and contraction/expansion events in the inverted repeat regions. A phylogenetic analysis of 63 shared genes in 16 plant species revealed that Campanula zangezura is the closest relative of C. carpatica. Phylogenetic analysis indicated that C. carpatica was within the Campanula clade, and C. pallida occupied the outermost position of that clade.

Increased resolution in the face of conflict: phylogenomics of the Neotropical bellflowers (Campanulaceae: Lobelioideae), a rapid plant radiation.

BACKGROUND AND AIMS: The centropogonid clade (Lobelioideae: Campanulaceae) is an Andean-centred rapid radiation characterized by repeated convergent evolution of morphological traits, including fruit type and pollination syndromes. While previous studies have resolved relationships of lineages with fleshy fruits into subclades, relationships among capsular species remain unresolved. This lack of resolution has impeded reclassification of non-monophyletic genera, whose current taxonomy relies heavily on traits that have undergone convergent evolution. METHODS: Targeted sequence capture using a probe-set recently developed for the centropogonid clade was used to obtain phylogenomic data from DNA extracted from both silica-dried and herbarium leaf tissue. These data were used to infer relationships among species using concatenated and partitioned species tree methods, and to quantify gene tree discordance. KEY RESULTS: While silica-dried leaf tissue resulted in longer assembled sequence data, the inclusion of herbarium samples improved taxonomic representation. Relationships among baccate lineages are similar to those inferred in previous studies, although they differ for lineages within and among capsular clades. We improve the phylogenetic resolution of Siphocampylus, which forms ten groups of closely related species which we informally name. Two subclades of Siphocampylus and two individual species are rogue taxa whose placement differs widely across analyses. Gene tree discordance (including cytonuclear discordance) is rampant. CONCLUSIONS: This first phylogenomic study of the centropogonid clade considerably improves our understanding of relationships in this rapid radiation. Differences across analyses and the possibility of additional lineage discoveries still hamper a solid and stable reclassification. Rapid morphological innovation corresponds with a high degree of phylogenomic complexity, including cytonuclear discordance, nuclear gene tree conflict and well-supported differences between analyses based on different nuclear loci. Together, these results point to a potential role of hemiplasy underlying repeated convergent evolution. This hallmark of rapid radiations is probably present in many other species-rich Andean plant radiations.

Development of highly polymorphic microsatellite markers for Campanula glomerata L. (Campanulaceae).

BACKGROUND: Campanula glomerata L. (Campanulaceae) is a self-incompatible, insect-pollinated herb occurring in calcareous grasslands, and is declining and (critically) endangered in many parts of its European distribution range. It often exists as small and isolated populations. A recovery plan of C. glomerata has been implemented in southern Belgium, involving plant translocations. METHODS AND RESULTS: We developed microsatellite markers using an enriched genomic library and characterized 16 loci in 111 individuals from eight populations. These 16 loci were highly polymorphic, with 11 to 31 alleles per locus for a total of 329 alleles, and expected heterozygosity (He) ranging from 0.470 to 0.938. CONCLUSIONS: These highly polymorphic loci constitute a promising tool for detailed genetic analyses: assigning individuals to distinct multilocus genotypes will allow quantifying pollen dispersal, clonal propagation and sexual recruitment and identifying admixed seed progeny and their pollen donors. Evaluating the genetic status of existing populations and a genetic monitoring of the translocated populations will contribute to optimize success in restoring viable and evolutionary resilient populations.

Complete chloroplast genome sequence of Adenophora racemosa (Campanulaceae): Comparative analysis with congeneric species.

Adenophora racemosa, belonging to the Campanulaceae, is an important species because it is endemic to Korea. The goal of this study was to assemble and annotate the chloroplast genome of A. racemosa and compare it with published chloroplast genomes of congeneric species. The chloroplast genome was reconstructed using de novo assembly of paired-end reads generated by the Illumina MiSeq platform. The chloroplast genome size of A. racemosa was 169,344 bp. In total, 112 unique genes (78 protein-coding genes, 30 tRNAs, and 4 rRNAs) were identified. A Maximum likelihood (ML) tree based on 76 protein-coding genes divided the five Adenophora species into two clades, showing that A. racemosa is more closely related to Adenophora stricta than to Adenophora divaricata. The gene order and contents of the LSC region of A. racemosa were identical to those of A. divaricata and A. stricta, but the structure of the SSC and IRs was unique due to IR contraction. Nucleotide diversity (Pi) >0.05 was found in eleven regions among the three Adenophora species not included in sect. Remotiflorae and in six regions between two species (A. racemosa and A. stricta).

Strength in numbers? Cytotype frequency mediates effect of reproductive barriers in mixed-ploidy arrays.

When differentiated lineages come into contact, their fates depend on demographic and reproductive factors. These factors have been well-studied in taxa of the same ploidy, but less is known about sympatric lineages that differ in ploidy, particularly with respect to demographic factors. We assessed prezygotic, postzygotic, and total reproductive isolation in naturally pollinated arrays of diploid-tetraploid and tetraploid-hexaploid population mixes of Campanula rotundifolia by measuring pollinator transitions, seed yield, germination rate, and proportion of hybrid offspring. Four frequencies of each cytotype were tested, and pollinators consistently overvisited rare cytotypes. Seed yield and F1 hybrid production were greater in 4X-6X arrays than 2X-4X arrays, whereas germination rates were similar, creating two distinct patterns of reproductive isolation. In 2X-4X arrays, postzygotic isolation was near complete (3% hybrid offspring), and prezygotic isolation associated with pollinator preference is expected to facilitate the persistence of minority cytotypes. However, in 4X-6X arrays where postzygotic isolation permitted hybrid formation (44% hybrids), pollinator behavior drove patterns of reproductive isolation, with rare cytotypes being more isolated and greater gene flow expected from rare into common cytotypes. In polyploid complexes, both the specific cytotypes in contact and local cytotype frequency, likely reflecting spatial demography, will influence likelihood of gene exchange.

Comparative analysis of plastid genomes within the Campanulaceae and phylogenetic implications.

The conflicts exist between the phylogeny of Campanulaceae based on nuclear ITS sequence and plastid markers, particularly in the subdivision of Cyanantheae (Campanulaceae). Besides, various and complicated plastid genome structures can be found in species of the Campanulaceae. However, limited availability of genomic information largely hinders the studies of molecular evolution and phylogeny of Campanulaceae. We reported the complete plastid genomes of three Cyanantheae species, compared them to eight published Campanulaceae plastomes, and shed light on a deeper understanding of the applicability of plastomes. We found that there were obvious differences among gene order, GC content, gene compositions and IR junctions of LSC/IRa. Almost all protein-coding genes and amino acid sequences showed obvious codon preferences. We identified 14 genes with highly positively selected sites and branch-site model displayed 96 sites under potentially positive selection on the three lineages of phylogenetic tree. Phylogenetic analyses showed that Cyananthus was more closely related to Codonopsis compared with Cyclocodon and also clearly illustrated the relationship among the Cyanantheae species. We also found six coding regions having high nucleotide divergence value. Hotpot regions were considered to be useful molecular markers for resolving phylogenetic relationships and species authentication of Campanulaceae.

Pollen colour morphs take different paths to fitness.

Colour phenotypes are often involved in communication and are thus under selection by species interactions. However, selection may also act on colour through correlated traits or alternative functions of biochemical pigments. Such forms of selection are instrumental in maintaining petal colour diversity in plants. Pollen colour also varies markedly, but the maintenance of this variation is little understood. In Campanula americana, pollen ranges from white to dark purple, with darker morphs garnering more pollinator visits and exhibiting elevated pollen performance under heat stress. Here, we generate an F2 population segregating for pollen colour and measure correlations with floral traits, pollen attributes and plant-level traits related to fitness. We determine the pigment biochemistry of colour variants and evaluate maternal and paternal fitness of light and dark morphs by crossing within and between morphs. Pollen colour was largely uncorrelated with floral traits (petal colour, size, nectar traits) suggesting it can evolve independently. Darker pollen grains were larger and had higher anthocyanin content (cyanidin and peonidin) which may explain why they outperform light pollen under heat stress. Overall, pollen-related fitness metrics were greater for dark pollen, and dark pollen sires generated seeds with higher germination potential. Conversely, light pollen plants produce 61% more flowers than dark, and 18% more seeds per fruit, suggesting a seed production advantage. Results indicate that light and dark morphs may achieve fitness through different means-dark morphs appear to have a pollen advantage whereas light morphs have an ovule advantage-helping to explain the maintenance of pollen colour variation.

FAD2 Gene Radiation and Positive Selection Contributed to Polyacetylene Metabolism Evolution in Campanulids.

Polyacetylenes (PAs) are bioactive, specialized plant defense compounds produced by some species in the eudicot clade campanulids. Early steps of PA biosynthesis are catalyzed by Fatty Acid Desaturase 2 (FAD2). Canonical FAD2s catalyze desaturation, but divergent forms can catalyze hydroxylation, conjugation, acetylenation, and epoxygenation. These alternate reactions give rise to valuable unusual fatty acids, including the precursors to PAs. The extreme functional diversity of FAD2 enzymes and the origin of PA biosynthesis are poorly understood from an evolutionary perspective. We focus here on the evolution of the FAD2 gene family. We uncovered a core eudicot-wide gene duplication event giving rise to two lineages: FAD2-α and FAD2-ß. Independent neofunctionalizations in both lineages have resulted in functionally diverse FAD2-LIKEs involved in unusual fatty acid biosynthesis. We found significantly accelerated rates of molecular evolution in FAD2-LIKEs and use this metric to provide a list of uncharacterized candidates for further exploration of FAD2 functional diversity. FAD2-α has expanded extensively in Asterales and Apiales, two main clades of campanulids, by ancient gene duplications. Here, we detected positive selection in both Asterales and Apiales lineages, which may have enabled the evolution of PA metabolism in campanulids. Together, these findings also imply that yet uncharacterized FAD2-α copies are involved in later steps of PA biosynthesis. This work establishes a robust phylogenetic framework in which to interpret functional data and to direct future research into the origin and evolution of PA metabolism.

Effects of glaciation and whole genome duplication on the distribution of the Campanula rotundifolia polyploid complex.

PREMISE OF THE STUDY: Both intrinsic and extrinsic factors contribute to a species distribution. Among plants, the extrinsic effects of glaciation and intrinsic effects of whole genome duplication are powerful drivers of biogeographical patterns, but the interplay of these factors is poorly understood. Here, we investigate the roles glaciation and whole-genome duplication have played in the evolution of the widespread polyploid complex Campanula rotundifolia. METHODS: We assessed the cytotype of 37 populations that spanned the geographic and cytotypic range of the C. rotundifolia complex. We constructed a chloroplast phylogeny for these populations and used RAD-seq to create nuclear phylogenies and networks for a subset of 23 populations; and estimated divergence times of major clades using Bayesian estimation of substitution rates. KEY RESULTS: Campanula rotundifolia originated in south-central Europe and underwent range expansion throughout much of Europe and North America. Multiple genome duplications have occurred in C. rotundifolia-at least two tetraploid and three hexaploid formations. CONCLUSIONS: Nuclear and chloroplast phylogenies are largely congruent with a history of populations surviving glacial maxima in known Pleistocene refugia in Europe and North America. Divergent European clades are consistent with two disjunct glacial refugia within Europe. North America was colonized by hexaploids derived from Western European lineages. A glacial refugium in Midwestern North America likely facilitated post-glacial recolonization of North America and limited genetic divergence. These results implicate both glaciation and whole-genome duplication as contributing factors to the extant biogeography of C. rotundifolia.

A tale of two forests: ongoing aridification drives population decline and genetic diversity loss at continental scale in Afro-Macaronesian evergreen-forest archipelago endemics.

Background and Aims: Various studies and conservationist reports have warned about the contraction of the last subtropical Afro-Macaronesian forests. These relict vegetation zones have been restricted to a few oceanic and continental islands around the edges of Africa, due to aridification. Previous studies on relict species have generally focused on glacial effects on narrow endemics; however, little is known about the effects of aridification on the fates of previously widespread subtropical lineages. Methods: Nuclear microsatellites and ecological niche modelling were used to understand observed patterns of genetic diversity in two emblematic species, widely distributed in these ecosystems: Canarina eminii (a palaeoendemic of the eastern Afromontane forests) and Canarina canariensis (a palaeoendemic of the Canarian laurel forests). The software DIYABC was used to test alternative demographic scenarios and an ensemble method was employed to model potential distributions of the selected plants from the end of the deglaciation to the present. Key Results: All the populations assessed experienced a strong and recent population decline, revealing that locally widespread endemisms may also be alarmingly threatened. Conclusions: The detected extinction debt, as well as the extinction spiral to which these populations are subjected, demands urgent conservation measures for the unique, biodiversity-rich ecosystems that they inhabit.

Origins of East Asian Campanuloideae (Campanulaceae) diversity.

The Campanulaceae comprises approximately 2300 species that are distributed among five major lineages: Campanuloideae, Cyphioideae, Cyphocarpoideae, Lobelioideae, and Nemacladoideae. Of these, the Campanuloideae, a primarily Old World clade, has the largest diversity in East Asia. In this study, we reconstruct the phylogeny of East Asian Campanuloideae based on one nuclear gene (i.e., PPR70) and five plastid markers (i.e., atpB, matK, petD, rbcL, and trnL-trnF). We then use this phylogenetic framework to reconstruct the biogeographical history of the genus. Our molecular dataset includes 376 of the 1045 currently recognized species in the Campanuloideae. Of the 376 sampled species, 116 are from East Asia, representing ca. 60% of the East Asian Campanuloideae. Our PPR dataset included sequences for 156 accessions, representing 54 species, while our plastid dataset included sequences for 305 accessions, representing 354 species. Phylogenetic analyses recovered three large clades containing East Asian taxa: Campanulinae, Platycodinae, and Wahlenberginae. The historical assembly of Campanuloideae diversity in East Asia appears to have resulted from numerous, independent movements from Africa, Europe/W. Asia, and North America. Africa was inferred as the ancestral range for the Campanuloideae. Movement of the largest East Asian clade (Platycodinae) occurred at approximately 53.1 Ma (46.6-58.73 95% HPD) from Africa, with much of the current diversity found in East Asia having resulted from in situ diversification. Thirteen additional movements into East Asia, primarily from Europe/Western Asia, occurred subsequently. One dispersal event from western North America was also inferred. In contrast, only six movements out of East Asia were found. Our results suggest that East Asia has acted primarily as a sink for Campanuloideae diversity, with Europe, Western Asia, and Africa representing major source areas.

Intercontinental dispersal and whole-genome duplication contribute to loss of self-incompatibility in a polyploid complex.

PREMISE OF THE STUDY: Angiosperm species often shift from self-incompatibility to self-compatibility following population bottlenecks. Across the range of a species, population bottlenecks may result from multiple factors, each of which may affect the geographic distribution and magnitude of mating-system shifts. We describe how intercontinental dispersal and genome duplication facilitate loss of self-incompatibility. METHODS: Self and outcross pollinations were performed on plants from 24 populations of the Campanula rotundifolia polyploid complex. Populations spanned the geographic distribution and three dominant cytotypes of the species (diploid, tetraploid, hexaploid). KEY RESULTS: Loss of self-incompatibility was associated with both intercontinental dispersal and genome duplication. European plants were largely self-incompatible, whereas North American plants were intermediately to fully self-compatible. Within both European and North American populations, loss of self-incompatibility increased as ploidy increased. Ploidy change and intercontinental dispersal both contributed to loss of self-incompatibility in North America, but range expansion did not affect self-incompatibility within Europe or North America. CONCLUSIONS: When species are subject to population bottlenecks arising through multiple factors, each factor can contribute to self-incompatibility loss. In a widespread polyploid complex, the loss of self-incompatibility can be predicted by the cumulative effects of whole-genome duplication and intercontinental dispersal.

The complete chloroplast genome sequences of three Adenophora species and comparative analysis with Campanuloid species (Campanulaceae).

We report the complete chloroplast genomes of three Adenophora species, and analyzed these compared them to five published Campanuloid plastomes. The total genome length of Adenophora divaricata, Adenophora erecta, and Adenophora stricta ranged from 159,759 to 176,331 bp. Among the eight Campanuloid species, many inversions were found to be only in the LSC region. IR contraction was also identified in the plastid genome of Adenophora stricta. Phylogenetic analyses based on 76 protein coding genes showed that Campanuloids are monophyletic, and are composed of two major groups: Campanula s. str. and Rapunculus. When we compared each homologous locus among the four Adenophora species, ten regions showed high nucleotide divergence value (>0.03). Among these, nine loci, excepting ycf3-rpoB, are considered to be useful molecular markers for phylogenetic studies and will be helpful to resolve phylogenetic relationships of Adenophora.

Does the seed bank contribute to the build-up of a genetic extinction debt in the grassland perennial Campanula rotundifolia?

Background and Aims: Habitat fragmentation threatens global biodiversity. Many plant species persist in habitat fragments via persistent life cycle stages such as seed banks, generating a species extinction debt. Here, seed banks are hypothesized to cause a temporal delay in the expected loss of genetic variation, which can be referred to as a genetic extinction debt, as a possible mechanism behind species extinction debts. Methods: Fragmented grassland populations of Campanula rotundifolia were examined for evidence of a genetic extinction debt, investigating if the seed bank contributed to the extinction debt build-up. The genetic make-up of 15 above- and below-ground populations was analysed in relation to historical and current levels of habitat fragmentation, both separately and combined. Key Results: Genetic diversity was highest in above-ground populations, though below-ground populations contained 8 % of unique alleles that were absent above-ground. Above-ground genetic diversity and composition were related to historical patch size and connectivity, but not current patch characteristics, suggesting the presence of a genetic extinction debt in the above-ground populations. No such relationships were found for the below-ground populations. Genetic diversity measures still showed a response to historical but not present landscape characteristics when combining genetic diversity of the above- and below-ground populations. Conclusions: The fragmented C. rotundifolia populations exhibited a genetic extinction debt. However, the role of the seed banks in the build-up of this extinction debt is probably small, since the limited, unique genetic diversity of the seed bank alone seems unable to counter the detrimental effects of habitat fragmentation on the population genetic structure of C. rotundifolia.

Embracing discordance: Phylogenomic analyses provide evidence for allopolyploidy leading to cryptic diversity in a Mediterranean Campanula (Campanulaceae) clade.

The Mediterranean Basin harbors a remarkable amount of biodiversity, a high proportion of which is endemic to this region. Here, we present an in-depth study of an angiosperm species complex, in which cryptic taxonomic diversity has been hypothesized. Specifically, we focus on four currently recognized species in the Roucela complex, a well-supported clade in the Campanulaceae/Campanuloideae: Campanula creutzburgii, C. drabifolia, C. erinus, and C. simulans. This study takes a phylogenomic approach, utilizing near-complete plastomes and 130 nuclear loci, to uncover cryptic diversity and test hypotheses regarding hybridization and polyploidy within this clade. Genome size estimates recovered tetraploid and octoploid lineages within the currently recognized, widespread species C. erinus, showing an east-west geographic pattern. Though genomic data clearly differentiate these two cytotypes, we failed to discern morphological differences. The formation of a cryptic octoploid lineage, distributed across the eastern Mediterranean, is hypothesized to be the result of an allopolyploid event in which one parental morphology is retained. The tetraploid C. erinus and C. creutzburgii (also a tetraploid) are implicated as parental lineages. Our results highlight the utility of target-enrichment approaches for obtaining genomic datasets for thorough assessments of species diversity and the importance of carefully considering gene-tree discordance within such datasets.

Biparental chloroplast inheritance leads to rescue from cytonuclear incompatibility.

Although organelle inheritance is predominantly maternal across animals and plants, biparental chloroplast inheritance has arisen multiple times in the angiosperms. Biparental inheritance has the potential to impact the evolutionary dynamics of cytonuclear incompatibility, interactions between nuclear and organelle genomes that are proposed to be among the earliest types of genetic incompatibility to arise in speciation. We examine the interplay between biparental inheritance and cytonuclear incompatibility in Campanulastrum americanum, a plant species exhibiting both traits. We first determine patterns of chloroplast inheritance in genetically similar and divergent crosses, and then associate inheritance with hybrid survival across multiple generations. There is substantial biparental inheritance in C. americanum. The frequency of biparental inheritance is greater in divergent crosses and in the presence of cytonuclear incompatibility. Biparental inheritance helps to mitigate cytonuclear incompatibility, leading to increased fitness of F1 hybrids and recovery in the F2 generation. This study demonstrates the potential for biparental chloroplast inheritance to rescue cytonuclear compatibility, reducing cytonuclear incompatibility's contribution to reproductive isolation and potentially slowing speciation. The efficacy of rescue depended upon the strength of incompatibility, with a greater persistence of weak incompatibilities in later generations. These findings suggest that incompatible plastids may lead to selection for biparental inheritance.

Phylogenetic relationships of Burmeistera (Campanulaceae: Lobelioideae): Combining whole plastome with targeted loci data in a recent radiation.

The field of molecular systematics has benefited greatly with the advent of high-throughput sequencing (HTS), making large genomic datasets commonplace. However, a large number of targeted Sanger sequences produced by many studies over the last two decades are publicly available and should not be overlooked. In this study, we elucidate the phylogenetic relationships of the plant genus Burmeistera (Campanulaceae: Lobelioideae), while investigating how to best combine targeted Sanger loci with HTS data. We sequence, annotate, and analyze complete to nearly complete plastomes for a subset of the genus. We then combine these data with a much larger taxonomic dataset for which only Sanger sequences are available, making this the most comprehensively sampled study in the genus to date. We show that using a phylogeny inferred from the species with plastome data as a topological constraint for the larger dataset increases the resolution of our data and produces a more robust evolutionary hypothesis for the group. We then use the resulting phylogeny to study the evolution of morphological traits thought to be important in Burmeistera, and assess their usefulness in the current taxonomic classification of the genus. The main morphological character used to delimit subgeneric sections, the presence or absence of hairs on the apex of the two ventral anthers, shows a complex evolutionary history with many changes in the tree, suggesting that this character should not be used for taxonomic classification. Although it is too soon to propose a new subgeneric classification for Burmeistera, our results highlight some morphological traits shared by whole clades that could potentially be used in future taxonomic work.