Pollinator shifts despite hybridisation in the Cape's hyperdiverse heathers (Erica, Ericaceae).

Interrogating the ecological and geographic factors that influence population divergence dynamics can reveal why some groups of organisms diversify more prolifically than others. One such group is the heathers (Erica, Ericaceae), the largest plant genus in the Cape Floristic Region. We study Erica abietina, a highly variable species complex with four subspecies differing in geographic range, habitat and pollination syndrome. We test for population differentiation, hybridisation, introgression and pollinator-driven divergence using genotyping-by-sequencing on samples across the entire distribution. We find five variably distinct genetic groups, with one subspecies comprising two independent lineages that are geographically isolated and occur on different soil types. Phylogenetic analysis suggests two independent shifts between bird and insect pollination, with accompanying genetic divergence. However, for one pair of populations with different pollinators, we uncover several individuals of hybrid origin at a site of sympatry. These results suggest that floral differentiation driven by divergent selection acts in concert with geographic isolation to maintain reproductive isolation and promote speciation. Our investigations reveal a highly dynamic system whose diversity has been shaped by a variety of interacting forces. We suggest that such a system could be a model for much of the diversification of the Cape flora.

Chromosome-scale genome assembly of Rhododendron molle provides insights into its evolution and terpenoid biosynthesis.

BACKGROUND: Rhododendron molle (Ericaceae) is a traditional Chinese medicine, which has been used to treat rheumatism and relieve pain since ancient times. The characteristic grayanoids of this plant have been demonstrated to be the chemical basis for the analgesic activity. Moreover, unlike morphine, these diterpenoids are non-addictive. Grayanoids mainly distribute in the leaves, flowers, roots, and fruits of R. molle, with low content. Currently the research on the biosynthesis of grayanoids is hindered, partially due to lack of the genomic information. RESULTS: In the present study, a total of 744 Mb sequences were generated and assembled into 13 chromosomes. An ancient whole-genome duplication event (Ad-ß) was discovered that occurred around 70 million years ago. Tandem and segmental gene duplications led to specific gene expansions in the terpene synthase and cytochrome P450 (CYP450) gene families. Two diterpene synthases were demonstrated to be responsible for the biosynthesis of 16α-hydroxy-ent-kaurane, the key precursor for grayanoids. Phylogenetic analysis revealed a species-specific bloom of the CYP71AU subfamily, which may involve the candidate CYP450s responsible for the biosynthesis of grayanoids. Additionally, three putative terpene biosynthetic gene clusters were found. CONCLUSIONS: We reported the first genome assembly of R. molle and investigated the molecular basis underpinning terpenoids biosynthesis. Our work provides a foundation for elucidating the complete biosynthetic pathway of grayanoids and studying the terpenoids diversity in R. molle.

Mitochondrial genomic data are effective at placing mycoheterotrophic lineages in plant phylogeny.

Fully mycoheterotrophic plants can be difficult to place in plant phylogeny due to elevated substitution rates associated with photosynthesis loss. This potentially limits the effectiveness of downstream analyses of mycoheterotrophy that depend on accurate phylogenetic inference. Although mitochondrial genomic data sets are rarely used in plant phylogenetics, theory predicts that they should be resilient to long-branch artefacts, thanks to their generally slow evolution, coupled with limited rate elevation in heterotrophs. We examined the utility of mitochondrial genomes for resolving contentious higher-order placements of mycoheterotrophic lineages in two test cases: monocots (focusing on Dioscoreales) and Ericaceae. We find Thismiaceae to be distantly related to Burmanniaceae in the monocot order Dioscoreales, conflicting with current classification schemes based on few gene data sets. We confirm that the unusual Afrothismia is related to Taccaceae-Thismiaceae, with a corresponding independent loss of photosynthesis. In Ericaceae we recovered the first well supported relationships among its five major lineages: mycoheterotrophic Ericaceae are not monophyletic, as pyroloids are inferred to be sister to core Ericaceae, and monotropoids to arbutoids. Genes recovered from mitochondrial genomes collectively resolved previously ambiguous mycoheterotroph higher-order relationships. We propose that mitochondrial genomic data should be considered in standardised gene panels for inferring overall plant phylogeny.

Leaps and bounds: geographical and ecological distance constrained the colonisation of the Afrotemperate by Erica.

BACKGROUND: The coincidence of long distance dispersal (LDD) and biome shift is assumed to be the result of a multifaceted interplay between geographical distance and ecological suitability of source and sink areas. Here, we test the influence of these factors on the dispersal history of the flowering plant genus Erica (Ericaceae) across the Afrotemperate. We quantify similarity of Erica climate niches per biogeographic area using direct observations of species, and test various colonisation scenarios while estimating ancestral areas for the Erica clade using parametric biogeographic model testing. RESULTS: We infer that the overall dispersal history of Erica across the Afrotemperate is the result of infrequent colonisation limited by geographic proximity and niche similarity. However, the Drakensberg Mountains represent a colonisation sink, rather than acting as a "stepping stone" between more distant and ecologically dissimilar Cape and Tropical African regions. Strikingly, the most dramatic examples of species radiations in Erica were the result of single unique dispersals over longer distances between ecologically dissimilar areas, contradicting the rule of phylogenetic biome conservatism. CONCLUSIONS: These results highlight the roles of geographical and ecological distance in limiting LDD, but also the importance of rare biome shifts, in which a unique dispersal event fuels evolutionary radiation.

Evolutionary constraints on disparity of ericaceous pollen grains.

BACKGROUND AND AIMS: Flowering plants show a high diversity of pollen morphology, assumed to reflect not only variations in the underlying design, but also stress imposed by ecological conditions related to pollen survival and germination. Both components are expected to constrain the accumulation of pollen disparity. However, this assumption has rarely been tested using empirical data. METHODS: This study is designed to test this hypothesis by inferring the accumulation of pollen disparity in Ericaceae, a large eudicot family with recent, ongoing radiations, with focus on three functionally significant pollen characters using a dated phylogeny. KEY RESULTS: Multiple lines of evidence supported the hypothesis that pollen disparity in Ericaceae did not evolve steadily but rather pulsed over time, clearly decoupling from the relative constant rate pattern of species diversification inferred. In a 3-D pollen morphospace, most major clades appear to occupy distinct neighbouring regions, whereas the subfamily Epacridoideae overlaps extensively with other subfamilies. No evidence for correlations was found between dimension of pollen disparity and species diversity at either the subfamily or generic level. Furthermore, the distribution of species in present pollen morphospace showed a strong central tendency, with the core compartment containing a large number of species from species-rich genera. CONCLUSIONS: The recovered evidence fits well with the expectations of limitations on available pollen morphological disparity, and suggests that innovation of pollen germination traits may have little effect on species diversification.

RNA sequencing and anthocyanin synthesis-related genes expression analyses in white-fruited Vaccinium uliginosum.

BACKGROUND: Vaccinium uliginosum (Ericaceae) is an important wild berry having high economic value. The white-fruited V. uliginosum variety found in the wild lacks anthocyanin and bears silvery white fruits. Hence, it is a good resource for investigating the mechanism of fruit color development. This study aimed to verify the differences in the expression levels of some structural genes and transcription factors affecting the anthocyanin biosynthesis pathway by conducting high-throughput transcriptome sequencing and real-time PCR analysis by using the ripening fruits of V. uliginosum and the white-fruited variety. RESULTS: We annotated 42,837 unigenes. Of the 325 differentially expressed genes, 41 were up-regulated and 284 were down-regulated. Further, 11 structural genes of the flavonoid pathway were up-regulated, whereas two were down-regulated. Of the seven genes encoding transcription factors, five were up-regulated and two were down-regulated. The structural genes VuCHS, VuF3'H, VuFHT, VuDFR, VuANS, VuANR, and VuUFGT and the transcription factors VubHLH92, VuMYB6, VuMYBPA1, VuMYB11, and VuMYB12 were significantly down-regulated. However, the expression of only VuMYB6 and VuMYBPA1 rapidly increased during the last two stages of V. uliginosum when the fruit was ripening, consistent with anthocyanin accumulation. CONCLUSIONS: VuMYB6 was annotated as MYB1 by the BLAST tool. Thus, the white fruit color in the V. uliginosum variant can be attributed to the down-regulation of transcription factors VuMYB1 and VuMYBPA1, which leads to the down-regulation of structural genes associated with the anthocyanin synthesis pathway.

Phylogenomic inference in extremis: A case study with mycoheterotroph plastomes.

PREMISE OF THE STUDY: Phylogenomic studies employing large numbers of genes, including those based on plastid genomes (plastomes), are becoming common. Nonphotosynthetic plants such as mycoheterotrophs (which rely on root-associated fungi for essential nutrients, including carbon) tend to have highly elevated rates of plastome evolution, substantial genome reduction, or both. Mycoheterotroph plastomes therefore provide excellent test cases for investigating how extreme conditions impact phylogenomic inference. METHODS: We used parsimony and likelihood analysis of protein-coding gene sets from published and newly completed plastomes to infer the phylogenetic placement of taxa from the 10 angiosperm families in which mycoheterotrophy evolved. KEY RESULTS: Despite multiple very long branches that reflect elevated substitution rates, and frequently patchy gene recovery due to genome reduction, inferred phylogenetic placements of most mycoheterotrophic lineages in DNA-based likelihood analyses are both well supported and congruent with other studies. Amino-acid-based likelihood placements are broadly consistent with DNA-based inferences, but extremely rate-elevated taxa can have unexpected placements-albeit with weak support. In contrast, parsimony analysis is strongly misled by long-branch attraction among many distantly related mycoheterotrophic monocots. CONCLUSIONS: Mycoheterotrophic plastomes provide challenging cases for phylogenomic inference, as substitutional rates can be elevated and genome reduction can lead to sparse gene recovery. Nonetheless, diverse likelihood frameworks provide generally well-supported and mutually concordant phylogenetic placements of mycoheterotrophs, consistent with recent phylogenetic studies and angiosperm-wide classifications. Previous predictions of parallel photosynthesis loss within families are supported for Burmanniaceae, Ericaceae, Gentianaceae, and Orchidaceae. Burmanniaceae and Thismiaceae should not be combined as a single family in Dioscoreales.

On the brink: the highly reduced plastomes of nonphotosynthetic Ericaceae.

Ericaceae (the heather family) is a large and diverse group of plants that forms elaborate symbiotic relationships with mycorrhizal fungi, and includes several nonphotosynthetic lineages. Using an extensive sample of fully mycoheterotrophic (MH) species, we explored inter- and intraspecific variation as well as selective constraints acting on the plastomes of these unusual plants. The plastomes of seven MH genera were analysed in a phylogenetic context with two geographically disparate individuals sequenced for Allotropa, Monotropa, and Pityopus. The plastomes of nonphotosynthetic Ericaceae are highly reduced in size (c. 33-41 kbp) and content, having lost all photosynthesis-related genes, and are reduced to encoding housekeeping genes as well as a protease subunit (clpP)-like and acetyl-CoA carboxylase subunit D (accD)-like open reading frames. Despite an increase in the rate of their nucleotide substitutions, the remaining protein-coding genes are typically under purifying selection in full MHs. We also identified ribosomal proteins under relaxed or neutral selection. These plastomes also exhibit striking structural rearrangements. Intraspecific variation within MH Ericaceae ranges from a few differences (Allotropa) to extensive population divergences (Monotropa, Hypopitys), which indicates that cryptic speciation may be occurring in several lineages. The pattern of gene loss within fully MH Ericaceae plastomes suggests an advanced state of degradation.

Late Pleistocene origin of the entire circumarctic range of the arctic-alpine plant Kalmia procumbens.

The circumarctic ranges of arctic-alpine plants are thought to have been established in the late Pliocene/early Pleistocene, when the modern arctic tundra was formed in response to climate cooling. Previous findings of range-wide genetic structure in arctic-alpine plants have been thought to support this hypothesis, but few studies have explicitly addressed the temporal framework of the genetic structure. Here, we estimated the demographic history of the genetic structure in the circumarctic Kalmia procumbens using sequences of multiple nuclear loci and examined whether its genetic structure reflects prolonged isolation throughout the Pleistocene. Both Bayesian clustering and phylogenetic analyses revealed genetic distinction between alpine and arctic regions, whereas detailed groupings were somewhat discordant between the analyses. By assuming a population grouping based on the phylogenetic analyses, which likely reflects a deeper intraspecific divergence, we conducted model-based analyses and demonstrated that the intraspecific genetic divergence in K. procumbens likely originated during the last glacial period. Thus, there is no need to postulate range separation throughout the Pleistocene to explain the current genetic structure in this species. This study demonstrates that range-wide genetic structure in arctic-alpine plants does not necessarily result from the late Pliocene/early Pleistocene origin of their circumarctic ranges and emphasizes the importance of a temporal framework of the current genetic structure for understanding the biogeographic history of the arctic flora.

[Identification and expression analysis of receptor-like kinase gene ERECTA in mycoheterotrophic plant Monotropa hypopitys].

The precise spatial-temporal coordination of cell division and differentiation is necessary for the correct formation of tissues, organs, and the organism as a whole. This coordination has been implemented by the intercellular communication mediated by signaling molecules and receptors that selectively recognize them. Membrane receptor kinases of ERECTA family regulate inflorescence and flower structure, the formation of root epidermis and adaptation responses. The characterization of the ERECTA genes of flowering plant pinesap Monotropa hypopitys with unique development features can enrich the knowledge about the kinase ERECTA functions and conserved development processes with their participation. Transcriptomic and genomic search with the subsequent structural-phylogenetic analysis identified the mRNA of a gene of serine-threonine kinase receptor with leucine-rich repeats of MhyERL1, which is the only ortholog of the ERECTA family kinases of pinesap. A quantitative analysis of the MhyERL1 gene transcripts has revealed its expression in all analyzed pinesap tissues with maximum levels in the flowers. MhyERL1 is probably involved in defining the inflorescence and flower architecture, and the formation of the pinesap root epidermis. The cascades involving ERL1 are apparently conserved. The exception are pathways associated with the development of above-ground vegetative structures, and the immune response to fungal pathogens probably lost in the process of the pinesap adaptation to unfavorable environmental conditions.

Identification and expression analysis of chitinase genes in parasitic plant Monotropa hypopitys.

Genes encoding six chitinases, five of which belong to classes I (MhCHI3 and MhCHI4), IV (MhCHI1), V (MhCHI5), and VII (MhCHI2), were identified in the transcriptome of the parasitic mixoheterotrophic plant Monotropa hypopitys. The transcription level of MhCHI5 and MhCHI1 was low; however, in the leaves (bracts) and roots it was higher than in flowers. MhCHI4 transcripts were detected primarily in the flowers and were almost absent in the roots, whereas the expression level of MhCHI3 was relatively high in all organs but maximum in the leaves (bracts).

Identification and characterization of the flower meristem identity gene MhyLFY in mycoheterotrophic plant Monotropa hypopitys.

The gene encoding the transcription factor LEAFY was identified in the genome of the mycoheterotrophic plant, pinesap Monotropa hypopitys. In the transcriptomes of roots, bracts, and flowers of flowering pinesaps, the MhyLFY gene expression was absent. These data suggest the conservativeness of the LFY-dependent mechanism of flower meristem identity and flower formation in heterotrophic species with some differences associated to the specificity of development and the structure of such plants. The pinesap flowering under the control of the transcription factor MhyLFY may be initiated either in an embryonic inflorescence during spring dormancy release of adventitious root buds or in an inflorescence of a growing reproductive stem after photoperiodic induction.

The biodiversity hotspot as evolutionary hot-bed: spectacular radiation of Erica in the Cape Floristic Region.

BACKGROUND: The disproportionate species richness of the world's biodiversity hotspots could be explained by low extinction (the evolutionary "museum") and/or high speciation (the "hot-bed") models. We test these models using the largest of the species rich plant groups that characterise the botanically diverse Cape Floristic Region (CFR): the genus Erica L. We generate a novel phylogenetic hypothesis informed by nuclear and plastid DNA sequences of c. 60 % of the c. 800 Erica species (of which 690 are endemic to the CFR), and use this to estimate clade ages (using RELTIME; BEAST), net diversification rates (GEIGER), and shifts in rates of diversification in different areas (BAMM; MuSSE). RESULTS: The diversity of Erica species in the CFR is the result of a single radiation within the last c. 15 million years. Compared to ancestral lineages in the Palearctic, the rate of speciation accelerated across Africa and Madagascar, with a further burst of speciation within the CFR that also exceeds the net diversification rates of other Cape clades. CONCLUSIONS: Erica exemplifies the "hotbed" model of assemblage through recent speciation, implying that with the advent of the modern Cape a multitude of new niches opened and were successively occupied through local species diversification.

Deep-sequence profiling of miRNAs and their target prediction in Monotropa hypopitys.

Myco-heterotroph Monotropa hypopitys is a widely spread perennial herb used to study symbiotic interactions and physiological mechanisms underlying the development of non-photosynthetic plant. Here, we performed, for the first time, transcriptome-wide characterization of M. hypopitys miRNA profile using high throughput Illumina sequencing. As a result of small RNA library sequencing and bioinformatic analysis, we identified 55 members belonging to 40 families of known miRNAs and 17 putative novel miRNAs unique for M. hypopitys. Computational screening revealed 206 potential mRNA targets for known miRNAs and 31 potential mRNA targets for novel miRNAs. The predicted target genes were described in Gene Ontology terms and were found to be involved in a broad range of metabolic and regulatory pathways. The identification of novel M. hypopitys-specific miRNAs, some with few target genes and low abundances, suggests their recent evolutionary origin and participation in highly specialized regulatory mechanisms fundamental for non-photosynthetic biology of M. hypopitys. This global analysis of miRNAs and their potential targets in M. hypopitys provides a framework for further investigation of miRNA role in the evolution and establishment of non-photosynthetic myco-heterotrophs.

The loss of photosynthetic pathways in the plastid and nuclear genomes of the non-photosynthetic mycoheterotrophic eudicot Monotropa hypopitys.

BACKGROUND: Chloroplasts of most plants are responsible for photosynthesis and contain a conserved set of about 110 genes that encode components of housekeeping gene expression machinery and photosynthesis-related functions. Heterotrophic plants obtaining nutrients from other organisms and their plastid genomes represent model systems in which to study the effects of relaxed selective pressure on photosynthetic function. The most evident is a reduction in the size and gene content of the plastome, which correlates with the loss of genes encoding photosynthetic machinery which become unnecessary. Transition to a non-photosynthetic lifestyle is expected also to relax the selective pressure on photosynthetic machinery in the nuclear genome, however, the corresponding changes are less known. RESULTS: Here we report the complete sequence of the plastid genome of Monotropa hypopitys, an achlorophyllous obligately mycoheterotrophic plant belonging to the family Ericaceae. The plastome of M. hypopitys is greatly reduced in size (35,336 bp) and lacks the typical quadripartite structure with two single-copy regions and an inverted repeat. Only 45 genes remained presumably intact- those encoding ribosomal proteins, ribosomal and transfer RNA and housekeeping genes infA, matK, accD and clpP. The clpP and accD genes probably remain functional, although their sequences are highly diverged. The sets of genes for ribosomal protein and transfer RNA are incomplete relative to chloroplasts of a photosynthetic plant. Comparison of the plastid genomes of two subspecies-level isolates of M. hypopitys revealed major structural rearrangements associated with repeat-driven recombination and the presence of isolate-specific tRNA genes. Analysis of the M. hypopitys transcriptome by RNA-Seq showed the absence of expression of nuclear-encoded components of photosystem I and II reaction center proteins, components of cytochrome b6f complex, ATP synthase, ribulose bisphosphate carboxylase components, as well as chlorophyll from protoporphyrin IX biosynthesis pathway. CONCLUSIONS: With the complete loss of genes related to photosynthesis, NADH dehydrogenase, plastid-encoded RNA polymerase and ATP synthase, the M. hypopitys plastid genome is among the most functionally reduced ones characteristic of obligate non-photosynthetic parasitic species. Analysis of the M. hypopitys transcriptome revealed coordinated evolution of the nuclear and plastome genomes and the loss of photosynthesis-related functions in both genomes.

RAD-seq data point to a northern origin of the arctic-alpine genus Cassiope (Ericaceae).

Many arctic-alpine plants display a highly disjunct distribution between the Arctic/Boreal regions and the southern Asian mountains. Two main hypotheses have been proposed to explain the origin of this biogeographic pattern: (1) south-to-north migration in the late Pliocene/early Pleistocene, and (2) north-to-south migration during the Miocene. The genus Cassiope is disjunctly distributed between the Arctic/Boreal regions and the Himalayan-Hengduan Mountains (HHM) and was selected to test these hypotheses. We constructed a fossil-calibrated phylogeny of Ericaceae using two plastid regions to estimate the crown group age of Cassiope, and used sequence data from thousands of loci produced by restriction site associated DNA sequencing (RAD-seq) to reconstruct the phylogeny of Cassiope. We also performed Bayesian divergence time analysis and biogeographic analysis. The Cassiope crown group was estimated to have originated in the Miocene, which predates the onset of Northern hemisphere glaciation. All HHM species formed a clade together with one eastern Siberian species, and this clade was sister to all other Arctic/Boreal species. This topology implies a northern origin of Cassiope, which is confirmed by our biogeographic analysis. Our results thus suggest that the ancient north-to-south migration hypothesis is most consistent with the origin of Cassiope.

Importance of demographic history for phylogeographic inference on the arctic-alpine plant Phyllodoce caerulea in East Asia.

Arctic-alpine plants have enormous ranges in the Northern Hemisphere. Phylogeographic studies have provided insights into their glacial survival as well as their postglacial colonization history. However, our understanding of the population dynamics of disjunct alpine populations in temperate regions remains limited. During Pleistocene cold periods, alpine populations of arctic-alpine species in East Asia were either connected to an ice-free Beringia refugium or they persisted with prolonged isolation after their establishment. To estimate which of these scenarios is more likely, we elucidated the genetic structure of Phyllodoce caerulea (Ericaceae) in Beringia and northern Japan, East Asia. Sequence variation in multiple nuclear loci revealed that P. caerulea can be distinguished into northern and southern groups. A demographic analysis demonstrated that the north-south divergence did not predate the last glacial period and detected introgression from Phyllodoce aleutica, relative widely distributed in East Asia, exclusively into the southern group. Therefore, although there has been genetic divergence between northern Japan and Beringia in P. caerulea, the divergence is unlikely to have resulted from their prolonged geographic separation throughout several cycles of glacial and interglacial periods. Instead, our study suggests that the introgression contributed to the genetic divergence of P. caerulea and that the range of P. caerulea was plausibly connected between northern Japan and Beringia during the last glacial period. Overall, our study not only provides a biogeographic insight into alpine populations of arctic-alpine plants in East Asia but also emphasizes the importance of careful interpretation of genetic structure for inferring phylogeographic history.

Clonality increases with snow depth in the arctic dwarf shrub Empetrum hermaphroditum.

PREMISE OF THE STUDY: Vegetative reproduction and spread through clonal growth plays an important role in arctic-alpine ecosystems with short cool growing seasons. Local variation in winter snow accumulation leads to discrete habitat types that may provide divergent conditions for sexual and vegetative reproduction. Therefore, we studied variation in clonal structure of a dominant, evergreen, dwarf shrub (Empetrum nigrum s.l. with the two taxa E. nigrum L. s.s. and E. hermaphroditum Hagerup) along a snow cover gradient and compared clonal diversity and spatial genetic structure between habitats. METHODS: We studied 374 individual shoots using 105 polymorphic AFLP markers and analyses based on hierarchical clustering, clonal diversity indices, and small-scale spatial genetic structure with pairwise kinship coefficient. We used two approaches to define a threshold of genotypic distance between two samples that are considered the same clone. Clonality was examined among three habitats (exposed ridges, sheltered depressions, birch forest) differing in snow conditions replicated in four study regions in Norway and Sweden. KEY RESULTS: Clonality of E. hermaphroditum differed between habitats with an increase in clonal diversity with decreasing snow depth. Small-scale spatial genetic structure increased with decreasing clonal diversity and increasing clone size. In three study regions, E. hermaphroditum was the only species, whereas in one region E. nigrum also occurred, largely confined to exposed ridges. CONCLUSIONS: Our results demonstrated that snow cover in conjunction with associated habitat conditions plays an important role for the mode of propagation of the dwarf shrub E. hermaphroditum.

As old as the mountains: the radiations of the Ericaceae.

Mountains are often more species-rich than lowlands. This could be the result of migration from lowlands to mountains, of a greater survival rate in mountains, or of a higher diversification rate in mountains. We investigated this question in the globally distributed family Ericaceae, which includes c. 4426 species ranging from sea level to > 5000 m. We predict that the interaction of low specific leaf area (SLA) and montane habitats is correlated with increased diversification rates. A molecular phylogeny of Ericaceae based on rbcL and matK sequence data was built and dated with 18 fossil calibrations and divergence time estimates. We identified radiations using bamm and correlates of diversification rate changes using binary-state speciation and extinction (BiSSE) and multiple-state speciation and extinction (MuSSE) analyses. Analyses revealed six largely montane radiations. Lineages in mountains diversified faster than nonmountain lineages (higher speciation rate, but no difference in extinction rate), and lineages with low SLA diversified faster than high-SLA lineages. Further, habitat and trait had a positive interactive effect on diversification. Our results suggest that the species richness in mountains is the result of increased speciation rather than reduced extinction or increased immigration. Increased speciation in Ericaceae was facilitated by low SLA.

Testing reticulate versus coalescent origins of Erica lusitanica using a species phylogeny of the northern heathers (Ericeae, Ericaceae).

Whilst most of the immense species richness of heathers (Calluna, Daboecia and Erica: Ericeae; Ericaceae) is endemic to Africa, particularly the Cape Floristic Region, the oldest lineages are found in the Northern Hemisphere. We present phylogenetic hypotheses for the major clades of Ericeae represented by multiple accessions of all northern Erica species and placeholder taxa for the large nested African/Madagascan clade. We identified consistent, strongly supported conflict between gene trees inferred from ITS and chloroplast DNA sequences with regard to the position of Erica lusitanica. We used coalescent simulations to test whether this conflict could be explained by coalescent stochasticity, as opposed to reticulation (e.g. hybridisation), given estimates of clade ages, generation time and effective population sizes (Ne). A standard approach, comparing overall differences between real and simulated trees, could not clearly reject coalescence. However, additional simulations showed that at the (higher) Ne necessary to explain conflict in E. lusitanica, further topological conflict would also be expected. Ancient hybridisation between ancestors of northern species is therefore a plausible scenario to explain the origin of E. lusitanica, and its morphological similarities to E. arborea. Assuming either process influences the results of species tree and further evolutionary inference. The coalescence scenario is equivocal with regard the standing hypothesis of stepping stone dispersal of Erica from Europe into Africa; whereas reticulate evolution in E. lusitanica would imply that the colonisation of Tropical East Africa by E. arborea instead occurred independently of dispersals within the rest of the African/Madagascan clade.