The genome of the king protea, Protea cynaroides.

The king protea (Protea cynaroides), an early-diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut-flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio-availability. Here, we present a high-quality chromosome-scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole-genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient-poor soils.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity.

PREMISE: Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives. METHODS: To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution. RESULTS: Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species. CONCLUSIONS: The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

Chromosome-level de novo genome assembly of Telopea speciosissima (New South Wales waratah) using long-reads, linked-reads and Hi-C.

Telopea speciosissima, the New South Wales waratah, is an Australian endemic woody shrub in the family Proteaceae. Waratahs have great potential as a model clade to better understand processes of speciation, introgression and adaptation, and are significant from a horticultural perspective. Here, we report the first chromosome-level genome for T. speciosissima. Combining Oxford Nanopore long-reads, 10x Genomics Chromium linked-reads and Hi-C data, the assembly spans 823 Mb (scaffold N50 of 69.0 Mb) with 97.8% of Embryophyta BUSCOs "Complete". We present a new method in Diploidocus (https://github.com/slimsuite/diploidocus) for classifying, curating and QC-filtering scaffolds, which combines read depths, k-mer frequencies and BUSCO predictions. We also present a new tool, DepthSizer (https://github.com/slimsuite/depthsizer), for genome size estimation from the read depth of single-copy orthologues and estimate the genome size to be approximately 900 Mb. The largest 11 scaffolds contained 94.1% of the assembly, conforming to the expected number of chromosomes (2n = 22). Genome annotation predicted 40,158 protein-coding genes, 351 rRNAs and 728 tRNAs. We investigated CYCLOIDEA (CYC) genes, which have a role in determination of floral symmetry, and confirm the presence of two copies in the genome. Read depth analysis of 180 "Duplicated" BUSCO genes using a new tool, DepthKopy (https://github.com/slimsuite/depthkopy), suggests almost all are real duplications, increasing confidence in the annotation and highlighting a possible need to revise the BUSCO set for this lineage. The chromosome-level T. speciosissima reference genome (Tspe_v1) provides an important new genomic resource of Proteaceae to support the conservation of flora in Australia and further afield.

High rates of evolution preceded shifts to sex-biased gene expression in Leucadendron, the most sexually dimorphic angiosperms.

Differences between males and females are usually more subtle in dioecious plants than animals, but strong sexual dimorphism has evolved convergently in the South African Cape plant genus Leucadendron. Such sexual dimorphism in leaf size is expected largely to be due to differential gene expression between the sexes. We compared patterns of gene expression in leaves among 10 Leucadendron species across the genus. Surprisingly, we found no positive association between sexual dimorphism in morphology and the number or the percentage of sex-biased genes (SBGs). Sex bias in most SBGs evolved recently and was species specific. We compared rates of evolutionary change in expression for genes that were sex biased in one species but unbiased in others and found that SBGs evolved faster in expression than unbiased genes. This greater rate of expression evolution of SBGs, also documented in animals, might suggest the possible role of sexual selection in the evolution of gene expression. However, our comparative analysis clearly indicates that the more rapid rate of expression evolution of SBGs predated the origin of bias, and shifts towards bias were depleted in signatures of adaptation. Our results are thus more consistent with the view that sex bias is simply freer to evolve in genes less subject to constraints in expression level.

Isolation of 27 polymorphic nuclear microsatellite markers for Roupala montana var. brasiliensis (Proteaceae).

Microsatellite primers pairs were developed for the Neotropical tree Roupala montana var. brasiliensis for use in studies on genetic diversity, mating system, and gene flow. Forty-two primer pairs were developed, resulting in 27 polymorphic loci, with two to 27 alleles per locus. The primer pairs were validated against 34 R. montana var. brasiliensis adult trees from four populations. The observed (H o) and expected (H e)heterozygosities ranged among loci from 0.061 to 0.930 (mean of 0.544) and from 0.116 to 0.950 (mean of 0.700), respectively. Null alleles were observed for ten loci. No genotypic linkage disequilibrium was detected in any pair of loci. This set of loci is suitable for population genetic studies of the species.

An exploration of assembly strategies and quality metrics on the accuracy of the rewarewa (Knightia excelsa) genome.

We used long read sequencing data generated from Knightia excelsa, a nectar-producing Proteaceae tree endemic to Aotearoa (New Zealand), to explore how sequencing data type, volume and workflows can impact final assembly accuracy and chromosome reconstruction. Establishing a high-quality genome for this species has specific cultural importance to Maori and commercial importance to honey producers in Aotearoa. Assemblies were produced by five long read assemblers using data subsampled based on read lengths, two polishing strategies and two Hi-C mapping methods. Our results from subsampling the data by read length showed that each assembler tested performed differently depending on the coverage and the read length of the data. Subsampling highlighted that input data with longer read lengths but perhaps lower coverage constructed more contiguous, kmers and gene-complete assemblies than short read length input data with higher coverage. The final genome assembly was constructed into 14 pseudochromosomes using an initial flye long read assembly, a racon/medaka/pilon combined polishing strategy, salsa2 and allhic scaffolding, juicebox curation, and Macadamia linkage map validation. We highlighted the importance of developing assembly workflows based on the volume and read length of sequencing data and established a robust set of quality metrics for generating high-quality assemblies. Scaffolding analyses highlighted that problems found in the initial assemblies could not be resolved accurately by Hi-C data and that assembly scaffolding was more successful when the underlying contig assembly was of higher accuracy. These findings provide insight into how quality assessment tools can be implemented throughout genome assembly pipelines to inform the de novo reconstruction of a high-quality genome assembly for nonmodel organisms.

Intraspecific trait variation influences physiological performance and fitness in the South Africa shrub genus Protea (Proteaceae).

BACKGROUND AND AIMS: Global plant trait datasets commonly identify trait relationships that are interpreted to reflect fundamental trade-offs associated with plant strategies, but often these trait relationships are not identified when evaluating them at smaller taxonomic and spatial scales. In this study we evaluate trait relationships measured on individual plants for five widespread Protea species in South Africa to determine whether broad-scale patterns of structural trait (e.g. leaf area) and physiological trait (e.g. photosynthetic rates) relationships can be detected within natural populations, and if these traits are themselves related to plant fitness. METHODS: We evaluated the variance structure (i.e. the proportional intraspecific trait variation relative to among-species variation) for nine structural traits and six physiological traits measured in wild populations. We used a multivariate path model to evaluate the relationships between structural traits and physiological traits, and the relationship between these traits and plant size and reproductive effort. KEY RESULTS: While intraspecific trait variation is relatively low for structural traits, it accounts for between 50 and 100 % of the variation in physiological traits. Furthermore, we identified few trait associations between any one structural trait and physiological trait, but multivariate regressions revealed clear associations between combinations of structural traits and physiological performance (R2 = 0.37-0.64), and almost all traits had detectable associations with plant fitness. CONCLUSIONS: Intraspecific variation in structural traits leads to predictable differences in individual-level physiological performance in a multivariate framework, even though the relationship of any particular structural trait to physiological performance may be weak or undetectable. Furthermore, intraspecific variation in both structural and physiological traits leads to differences in plant size and fitness. These results demonstrate the importance of considering measurements of multivariate phenotypes on individual plants when evaluating trait relationships and how trait variation influences predictions of ecological and evolutionary outcomes.

Climate explains population divergence in drought-induced plasticity of functional traits and gene expression in a South African Protea.

Long-term environmental variation often drives local adaptation and leads to trait differentiation across populations. Additionally, when traits change in an environment-dependent way through phenotypic plasticity, the genetic variation underlying plasticity will also be under selection. These processes could create a landscape of differentiation across populations in traits and their plasticity. Here, we performed a dry-down experiment under controlled conditions to measure responses in seedlings of a shrub species from the Cape Floristic Region, the common sugarbush (Protea repens). We measured morphological and physiological traits, and sequenced whole transcriptomes of leaf tissues from eight populations that represent both the climatic and the geographical distribution of this species. We found that there is substantial variation in how populations respond to drought, but we also observed common patterns such as reduced leaf size and leaf thickness, and up-regulation of stress-related and down-regulation of growth-related gene groups. Both high environmental heterogeneity and milder source site climates were associated with higher plasticity in various traits and co-expression gene networks. Associations between traits, trait plasticity, gene networks and the source site climate suggest that temperature may play a greater role in shaping these patterns when compared to precipitation, in line with recent changes in the region due to climate change. We also found that traits respond to climatic variation in an environment-dependent manner: some associations between traits and climate were apparent only under certain growing conditions. Together, our results uncover common responses of P. repens populations to drought, and climatic drivers of population differentiation in functional traits, gene expression and their plasticity.

Detailed Phenolic Characterization of Protea Pure and Hybrid Cultivars by Liquid Chromatography-Ion Mobility-High Resolution Mass Spectrometry (LC-IM-HR-MS).

In this study we report a detailed investigation of the polyphenol composition of Protea pure (P. cynaroides and P. neriifolia) and hybrid cultivars (Black beauty and Limelight). Aqueous methanol extracts of leaf and bract tissues were analyzed by ultrahigh pressure liquid chromatography hyphenated to photodiode array and ion mobility-high resolution mass spectrometric (UHPLC-PDA-IM-HR-MS) detection. A total of 67 metabolites were characterized based on their relative reversed phase (RP) retention, UV-vis spectra, low and high collision energy HR-MS data, and collisional cross section (CCS) values. These metabolites included 41 phenolic acid esters and 25 flavonoid derivatives, including 5 anthocyanins. In addition, an undescribed hydroxycinnamic acid-polygalatol ester, caffeoyl-O-polygalatol (1,5-anhydro-[6-O-caffeoyl]-sorbitol(glucitol)) was isolated and characterized by 1D and 2D NMR for the first time. This compound and its isomer are shown to be potential chemo-taxonomic markers.

Equilibrium and non-equilibrium phases in the radiation of Hakea and the drivers of diversity in Mediterranean-type ecosystems.

Mediterranean-type ecosystems (MTEs) contain exceptional plant diversity. Explanations for this diversity are usually classed as either "equilibrium," with elevated MTE diversity resulting from greater ecological carrying capacities, or "non-equilibrium," with MTEs having a greater accumulation of diversity over time than other types of ecosystems. These models have typically been considered as mutually exclusive. Here, we present a trait-based explanatory framework that incorporates both equilibrium and non-equilibrium dynamics. Using a large continental Australian plant radiation (Hakea) as a case study, we identify traits associated with niche partitioning in coexisting species (α-traits) and with environmental filtering (ß-traits), and reconstruct the mode and relative timing of diversification of these traits. Our results point to a radiation with an early non-equilibrium phase marked by divergence of ß-traits as Hakea diversified exponentially and expanded from the southwest Australian MTE into biomes across the Australian continent. This was followed from seven million years ago by an equilibrium phase, marked by diversification of α-traits and a slowdown in lineage diversification as MTE-niches became saturated. These results suggest that processes consistent with both equilibrium and non-equilibrium models have been important during different stages of the radiation of Hakea, and together they provide a richer explanation of present-day diversity patterns.

Environmental drivers and genomic architecture of trait differentiation in fire-adapted Banksia attenuata ecotypes.

Trait divergence between populations is considered an adaptive response to different environments, but to what extent this response is accompanied by genetic differentiation is less clear since it may be phenotypic plasticity. In this study, we analyzed phenotypic variation between two Banksia attenuata growth forms, lignotuberous (shrub) and epicormic resprouting (tree), in fire-prone environments to identify the environmental factors that have driven this phenotypic divergence. We linked genotype with phenotype and traced candidate genes using differential gene expression analysis. Fire intervals determined the phenotypic divergence between growth forms in B. attenuata. A genome-wide association study identified 69 single nucleotide polymorphisms, putatively associated with growth form, whereas no growth form- or phenotype-specific genotypes were identified. Genomic differentiation between the two growth forms was low (Fst = 0.024). Differential gene expression analysis identified 37 genes/transcripts that were differentially expressed in the two growth forms. A small heat-shock protein gene, associated with lignotuber presence, was differentially expressed in the two forms. We conclude that different fire regimes induce phenotypic polymorphism in B. attenuata, whereas phenotypic trait divergence involves the differential expression of a small fraction of genes that interact strongly with the disturbance regime. Thus, phenotypic plasticity among resprouters is the general strategy for surviving varying fire regimes.

How Does Evolution in Phosphorus-Impoverished Landscapes Impact Plant Nitrogen and Sulfur Assimilation?

Phosphorus (P) fertilisers, made from rock phosphate, are used to attain high crop yields. However, rock phosphate is a finite resource and excessive P fertilisers pollute our environment, stressing the need for more P-efficient crops. Some Proteaceae have evolved in extremely P-impoverished environments. One of their adaptations is to curtail the abundance of ribosomal RNA, and thus protein, and tightly control the acquisition and assimilation of nitrogen (N) and sulfur. This differs fundamentally from plants that evolved in environments where N limits plant productivity, but is likely common in many species that evolved in P-impoverished landscapes. Here, we scrutinise the relevance of these responses towards developing P-efficient crops, focusing on plant species where 'P is in the driver's seat'.

Anthropogenic fragmentation may not alter pre-existing patterns of genetic diversity and differentiation in perennial shrubs.

Many plant species have pollination and seed dispersal systems and evolutionary histories that have produced strong genetic structuring. These genetic patterns may be consistent with expectations following recent anthropogenic fragmentation, making it difficult to detect fragmentation effects if no prefragmentation genetic data are available. We used microsatellite markers to investigate whether severe habitat fragmentation may have affected the structure and diversity of populations of the endangered Australian bird-pollinated shrub Grevillea caleyi R.Br., by comparing current patterns of genetic structure and diversity with those of the closely related G. longifolia R.Br. that has a similar life history but has not experienced anthropogenic fragmentation. Grevillea caleyi and G. longifolia showed similar and substantial population subdivision at all spatial levels (global F'ST  = 0.615 and 0.454; Sp  = 0.039 and 0.066), marked isolation by distance and large heterozygous deficiencies. These characteristics suggest long-term effects of inbreeding in self-compatible species that have poor seed dispersal, limited connectivity via pollen flow and undergo population bottlenecks because of periodic fires. Highly structured allele size distributions, most notably in G. caleyi, imply historical processes of drift and mutation were important in isolated subpopulations. Genetic diversity did not vary with population size but was lower in more isolated populations for both species. Through this comparison, we reject the hypothesis that anthropogenic fragmentation has impacted substantially on the genetic composition or structure of G. caleyi populations. Our results suggest that highly self-compatible species with limited dispersal may be relatively resilient to the genetic changes predicted to follow habitat fragmentation.

Life-History Traits Evolved Jointly with Climatic Niche and Disturbance Regime in the Genus Leucadendron (Proteaceae).

Organisms have evolved a diversity of life-history strategies to cope with variation in their environment. Persistence as adults and/or seeds across recruitment events allows species to dampen the effects of environmental fluctuations. The evolution of life cycles with overlapping generations should thus permit the colonization of environments with uncertain recruitment. We tested this hypothesis in Leucadendron (Proteaceae), a genus with high functional diversity native to fire-prone habitats in the South African fynbos. We analyzed the joint evolution of life-history traits (adult survival and seed-bank strategies) and ecological niches (climate and fire regime), using comparative methods and accounting for various sources of uncertainty. In the fynbos, species with canopy seed banks that are unable to survive fire as adults display nonoverlapping generations. In contrast, resprouters with an underground seed bank may be less threatened by extreme climatic events and fire intervals, given their iteroparity and long-lasting seed bank. Life cycles with nonoverlapping generations indeed jointly evolved with niches with less exposure to frost but not with those with less exposure to drought. Canopy seed banks jointly evolved with niches with more predictable fire return, compared to underground seed banks. The evolution of extraordinary functional diversity among fynbos plants thus reflects, at least in part, the diversity of both climates and fire regimes in this region.

Isolation by distance and isolation by environment contribute to population differentiation in Protea repens (Proteaceae L.), a widespread South African species.

PREMISE OF THE STUDY: The Cape Floristic Region (CFR) of South Africa is renowned for its botanical diversity, but the evolutionary origins of this diversity remain controversial. Both neutral and adaptive processes have been implicated in driving diversification, but population-level studies of plants in the CFR are rare. Here, we investigate the limits to gene flow and potential environmental drivers of selection in Protea repens L. (Proteaceae L.), a widespread CFR species. METHODS: We sampled 19 populations across the range of P. repens and used genotyping by sequencing to identify 2066 polymorphic loci in 663 individuals. We used a Bayesian FST outlier analysis to identify single-nucleotide polymorphisms (SNPs) marking genomic regions that may be under selection; we used those SNPs to identify potential drivers of selection and excluded them from analyses of gene flow and genetic structure. RESULTS: A pattern of isolation by distance suggested limited gene flow between nearby populations. The populations of P. repens fell naturally into two or three groupings, which corresponded to an east-west split. Differences in rainfall seasonality contributed to diversification in highly divergent loci, as do barriers to gene flow that have been identified in other species. CONCLUSIONS: The strong pattern of isolation by distance is in contrast to the findings in the only other widespread species in the CFR that has been similarly studied, while the effects of rainfall seasonality are consistent with well-known patterns. Assessing the generality of these results will require investigations of other CFR species.

Characterization of Leaf Transcriptome in Banksia hookeriana.

Banksia is a significant element in vegetation of southwestern Australia, a biodiversity hotspot with global significance. In particular, Banksia hookeriana represents a species with significant economic and ecological importance in the region. For better conservation and management, we reported an overview of transcriptome of B. hookeriana using RNA-seq and de novo assembly. We have generated a total of 202.7 million reads (18.91 billion of nucleotides) from four leaf samples in four plants of B. hookeriana, and assembled 59,063 unigenes (average size=1098bp) through de novotranscriptome assembly. Among them, 39,686 unigenes were annotated against the Swiss-Prot, Clusters of Orthologous Groups (COG), and NCBI non-redundant (NR) protein databases. We showed that there was approximately one single nucleotide polymorphism (SNP) per 5.6-7.1kb in the transcriptome, and the ratio of transitional to transversional polymorphisms was approximately 1.82. We compared unigenes of B. hookeriana to those of Arabidopsis thaliana and Nelumbo nucifera through sequence homology, Gene Ontology (GO) annotation, and KEGG pathway analyses. The comparative analysis revealed that unigenes of B. hookeriana were closely related to those of N. nucifera. B. hookeriana, N. nucifera, and A. thaliana shared similar GO annotations but different distributions in KEGG pathways, indicating that B. hookeriana has adapted to dry-Mediterranean type shrublands via regulating expression of specific genes. In total 1927 potential simple sequence repeat (SSR) markers were discovered, which could be used in the genotype and genetic diversity studies of the Banksia genus. Our results provide valuable sequence resource for further study in Banksia.

Anchored phylogenomics improves the resolution of evolutionary relationships in the rapid radiation of Protea L.

PREMISE OF THE STUDY: Estimating phylogenetic relationships in relatively recent evolutionary radiations is challenging, especially if short branches associated with recent divergence result in multiple gene tree histories. We combine anchored enrichment next-generation sequencing with species tree analyses to produce a robust estimate of phylogenetic relationships in the genus Protea (Proteaceae), an iconic radiation in South Africa. METHODS: We sampled multiple individuals within 59 out of 112 species of Protea and 6 outgroup species for a total of 163 individuals, and obtained sequences for 498 low-copy, orthologous nuclear loci using anchored phylogenomics. We compare several approaches for building species trees, and explore gene tree-species tree discrepancies to determine whether poor phylogenetic resolution reflects a lack of informative sites, incomplete lineage sorting, or hybridization. KEY RESULTS: Phylogenetic estimates from species tree approaches are similar to one another and recover previously well-supported clades within Protea, in addition to providing well-supported phylogenetic hypotheses for previously poorly resolved intrageneric relationships. Individual gene trees are markedly different from one another and from species trees. Nonetheless, analyses indicate that differences among gene trees occur primarily concerning clades supported by short branches. CONCLUSIONS: Species tree methods using hundreds of nuclear loci provided strong support for many previously unresolved relationships in the radiation of the genus Protea. In cases where support for particular relationships remains low, these appear to arise from few informative sites and lack of information rather than strongly supported disagreement among gene trees.

Characterization of CYCLOIDEA-like genes in Proteaceae, a basal eudicot family with multiple shifts in floral symmetry.

BACKGROUND AND AIMS: The basal eudicot family Proteaceae (approx. 1700 species) shows considerable variation in floral symmetry but has received little attention in studies of evolutionary development at the genetic level. A framework for understanding the shifts in floral symmetry in Proteaceae is provided by reconstructing ancestral states on an upated phylogeny of the family, and homologues of CYCLOIDEA (CYC), a key gene for the control of floral symmetry in both monocots and eudicots, are characterized. METHODS: Perianth symmetry transitions were reconstructed on a new species-level tree using parsimony and maximum likelihood. CYC-like genes in 35 species (31 genera) of Proteaceae were sequenced and their phylogeny was reconstructed. Shifts in selection pressure following gene duplication were investigated using nested branch-site models of sequence evolution. Expression patterns of CYC homologues were characterized in three species of Grevillea with different types of floral symmetry. KEY RESULTS: Zygomorphy has evolved 10-18 times independently in Proteaceae from actinomorphic ancestors, with at least four reversals to actinomorphy. A single duplication of CYC-like genes occurred prior to the diversification of Proteaceae, with putative loss or divergence of the ProtCYC1 paralogue in more than half of the species sampled. No shifts in selection pressure were detected in the branches subtending the two ProtCYC paralogues. However, the amino acid sequence preceding the TCP domain is strongly divergent in Grevillea ProtCYC1 compared with other species. ProtCYC genes were expressed in developing flowers of both actinomorphic and zygomorphic Grevillea species, with late asymmetric expression in the perianth of the latter. CONCLUSION: Proteaceae is a remarkable family in terms of the number of transitions in floral symmetry. Furthermore, although CYC-like genes in Grevillea have unusual sequence characteristics, they display patterns of expression that make them good candidates for playing a role in the establishment of floral symmetry.

Pre-Gondwanan-breakup origin of Beauprea (Proteaceae) explains its historical presence in New Caledonia and New Zealand.

New Caledonia and New Zealand belong to the now largely submerged continent Zealandia. Their high levels of endemism and species richness are usually considered the result of transoceanic dispersal events followed by diversification after they re-emerged from the Pacific Ocean in the mid-Cenozoic. We explore the origin and evolutionary history of Beauprea (Proteaceae), which is now endemic to New Caledonia but was once spread throughout eastern Gondwana, including New Zealand. We review the extensive Beauprea-type pollen data in the fossil records and analyze the relationship of these fossil taxa to extant genera within Proteaceae. We further reconstruct the phylogenetic relations among nine extant species of Beauprea and estimate the age of the Beauprea clade. By incorporating extinct taxa into the Beauprea phylogenetic tree, we reconstruct the ancient distribution of this genus. Our analysis shows that Beauprea originated c. 88 Ma (million years ago) in Antarctica-Southeastern Australia and spread throughout Gondwana before its complete breakup. We propose that Beauprea, already existing as two lineages, was carried with Zealandia when it separated from the rest of Gondwana c. 82 Ma, thus supporting an autochthonous origin for Beauprea species now in New Caledonia and historically in New Zealand up to 1 Ma. We show that the presence of Beauprea through transoceanic dispersal is implausible. This means that neither New Caledonia nor New Zealand has been entirely submerged since the Upper Cretaceous; thus, possible vicariance and allopatry must be taken into account when considering the high levels of endemism and species richness of these island groups.

Evolutionary potential and adaptation of Banksia attenuata (Proteaceae) to climate and fire regime in southwestern Australia, a global biodiversity hotspot.

Substantial climate changes are evident across Australia, with declining rainfall and rising temperature in conjunction with frequent fires. Considerable species loss and range contractions have been predicted; however, our understanding of how genetic variation may promote adaptation in response to climate change remains uncertain. Here we characterized candidate genes associated with rainfall gradients, temperatures, and fire intervals through environmental association analysis. We found that overall population adaptive genetic variation was significantly affected by shortened fire intervals, whereas declining rainfall and rising temperature did not have a detectable influence. Candidate SNPs associated with rainfall and high temperature were diverse, whereas SNPs associated with specific fire intervals were mainly fixed in one allele. Gene annotation further revealed four genes with functions in stress tolerance, the regulation of stomatal opening and closure, energy use, and morphogenesis with adaptation to climate and fire intervals. B. attenuata may tolerate further changes in rainfall and temperature through evolutionary adaptations based on their adaptive genetic variation. However, the capacity to survive future climate change may be compromised by changes in the fire regime.