Interspecific hybridization and island colonization history, not rarity, most strongly affect the genetic diversity in a clade of Mascarene-endemic trees

Many factors shape the genetic diversity of island-endemic trees, with important implications for conservation. Oceanic island-endemic lineages undergo an initial founding bottleneck during the colonization process and subsequently accumulate diversity following colonization. Moreover, many island endemics occur in small populations and are further threatened by anthropogenic factors that cause population declines, making them susceptible to losses in genetic diversity through genetic drift, inbreeding, and bottlenecks. However, life-history traits commonly found in trees, such as outcrossing mechanisms, long lifespans, and a propensity for interspecific hybridization, may help buffer against losses of genetic variation. To assess the relative importance of colonization history, rarity, and distribution in shaping genetic diversity of island-endemic trees, we conducted a comparative population genomic analysis of 13 species of Diospyros (Ebenaceae) endemic to the Mascarene Islands that differ in island colonization history, distribution, population size, and IUCN threat status. We genotyped 328 individuals across the islands using 2b-RADseq, compared genetic diversity both among and within species, and assessed patterns of genetic structure. Genetic diversity did not vary significantly by IUCN status, but we found that species that co-occur with others on the same intermediate-aged island (Mauritius) had much greater genetic diversity than those that occur solitarily on an island (Réunion and Rodrigues), likely because of greater interspecific hybridization among species with overlapping distributions and processes related to time since island colonization. Results presented here were used to determine priority localities for in situ and ex situ conservation efforts to maximize the genetic diversity of each Mascarene Diospyros species.

Methods
leaf tissue was collected and stored in silica for DNA extraction. A voucher collection was made for each species at each sampling location; duplicates were deposited in the MO, CBMN, and/or MAU herbaria (Supplementary Table S1 in the manuscript).  DNA was extracted used a modified CTAB method (Doyle and Doyle 1987). 2bRAD libraries were prepared following Wang et al. 2012 using dual index barcoded adapters to multiplex 96 samples per library.  Samples were sequenced on an Illumina Hiseq 2000 using 1 × 50 bp reads.

reads were assembled using STACKS v. 2.2.