Inferred vs realized patterns of gene flow: an analysis of population structure in the Andros Island Rock Iguana

Phylogenomics and historical biogeography of West Indian Rock Iguanas (genus Cyclura)

The genus Cyclura includes nine extant species and six subspecies of West Indian Rock Iguanas and is one of the most imperiled genera of squamate reptiles globally. An understanding of species diversity, evolutionary relationships, diversification, and historical biogeography in this group is crucial for implementing sound long-term conservation strategies. We collected DNA samples from 1 to 10 individuals per taxon from all Cyclura taxa (n = 70 ingroup individuals), focusing where possible on incorporating individuals from different populations of each species. We also collected 1-2 individuals from each of seven outgroup species of iguanas (Iguana delicatissima; five Ctenosaura species) and Anolis sagrei (n = 12 total outgroup individuals). We used targeted genomic sequence capture to isolate and to sequence 1,872 loci comprising of 687,308 base pairs (bp) from each of the 82 individuals from across the nuclear genome. We extracted mitochondrial reads and assembled and annotated mitogenomes for all Cyclura taxa plus outgroup species. We present well-supported phylogenomic gene tree/species tree analyses for all extant species of Cyclura using ASTRAL-III, SVDQuartets, and StarBEAST2 methods, and discuss the taxonomic, biogeographic, and conservation implications of these data. We find a most recent common ancestor of the genus 9.91 million years ago. The earliest divergence within Cyclura separates C. pinguis from a clade comprising all other Cyclura. Within the latter group, a clade comprising C. carinata from the southern Lucayan Islands and C. ricordii from Hispaniola is the sister taxon to a clade comprising the other Cyclura. Among the other Cyclura, the species C. cornuta and C. stejnegeri (from Hispaniola and Isla Mona) form the sister taxon to a clade of species from Jamaica (C. collei), Cuba and Cayman Islands (C. nubila and C. lewisi), and the eastern (C. rileyi) and western (C. cychlura) Lucayan Islands. Cyclura cychlura and C. rileyi form a clade whose sister taxa are C. nubila and C. lewisi. Cyclura collei is the sister taxon to these four species combined.

Correlated population genetic structure in a three-tiered host-parasite system. The potential for coevolution and adaptive divergence

Three subspecies of Northern Bahamian Rock Iguanas, Cyclura cychlura, are currently recognized: C. c. cychlura, restricted to Andros Island, and C. c. figginsi and C. c. inornata, native to the Exuma Island chain. Populations on Andros are genetically distinct from Exuma Island populations, yet genetic divergence among populations in the Exumas is inconsistent with the two currently recognized subspecies from those islands. The potential consequences of this discrepancy might include the recognition of a single subspecies throughout the Exumas rather than two. That inference also ignores evidence that populations of C. cychlura are potentially adaptively divergent. We compared patterns of population relatedness in a three-tiered host-parasite system: C. cychlura iguanas, their ticks (genus Amblyomma, preferentially parasitizing these reptiles), and Rickettsia spp. endosymbionts (within tick ectoparasites). Our results indicate that while C. c. cychlura on Andros is consistently supported as a separate clade, patterns of relatedness among populations of C. c. figginsi and C. c. inornata within the Exuma Island chain are more complex. The distribution of the hosts, different tick species, and Rickettsia spp., supports the evolutionary independence of C. c. inornata. Further, these patterns are also consistent with two independent evolutionarily significant units within C. c. figginsi. Our findings suggest coevolutionary relationships between the reptile hosts, their ectoparasites, and rickettsial organisms, suggesting local adaptation. This work also speaks to the limitations of using neutral molecular markers from a single focal taxon as the sole currency for recognizing evolutionary novelty in populations of endangered species.

Heterozygosity-Fitness Correlations Reveal Inbreeding Depression in Neonatal Body Size in a Critically Endangered Rock Iguana

Inbreeding depression, though challenging to identify in nature, may play an important role in regulating the dynamics of small and isolated populations. Conversely, greater expression of genetic load can enhance opportunities for natural selection. Conditional expression concentrates these opportunities for selection and may lead to failure of detection. This study investigates the possibility for age-dependent expression of inbreeding depression in a critically endangered population of rock iguanas, Cyclura nubila caymanensis. We employ heterozygote-fitness correlations to examine the contributions of individual genetic factors to body size, a fitness-related trait. Nonsignificant reductions in homozygosity (up to 7%) were detected between neonates and individuals surviving past their first year, which may reflect natural absorption of inbreeding effects by this small, fecund population. The majority of variation in neonate body size was attributed to maternal or environmental effects (i.e., clutch identity and incubation length); however, heterozygosity across 22 microsatellite loci also contributed significantly and positively to model predictions. Conversely, effects of heterozygosity on fitness were not detectable when adults were examined, suggesting that inbreeding depression in body size may be age dependent in this taxon. Overall, these findings emphasize the importance of taking holistic, cross-generational approaches to genetic monitoring of endangered populations.

The Lesser Antillean Iguana (Iguana delicatissima) on St. Eustatius: Genetically Depauperate and Threatened by Ongoing Hybridization

The Lesser Antillean Iguana (Iguana delicatissima) is an endangered species threatened by habitat loss and hybridization with non-native Green Iguanas (Iguana iguana). Iguana delicatissima has been extirpated on several islands, and the Green Iguana has invaded most islands with extant populations. Information is essential to protect this species from extinction. We collected data on 293 iguanas including 17 juveniles from St. Eustasius, one of the few remaining I. delicatissima strongholds. Genetic data were leveraged to test for hybridization presence with the Green Iguana using both mitochondrial and nuclear genes, including 16 microsatellite loci. The microsatellites were also analyzed to estimate genetic diversity, population structure, and effective population size. Using molecular and morphological data, we identified 286 I. delicatissima individuals captured during our first fieldwork effort, and 7 non-native iguanas captured during a second effort, showing hybridization occurs within this population. Comparing homologous microsatellites used in studies on Dominica and Chancel, the I. delicatissima population on St. Eustatius has extremely low genetic diversity (HO = 0.051; HE = 0.057), suggesting this population is genetically depauperate. Furthermore, there is significant evidence for inbreeding (FIS = 0.12) and weak spatial genetic structure (FST = 0.021, P = 0.002) within this population. Besides immediate threats including hybridization, this population's low genetic diversity, presence of physiological abnormalities and low recruitment could indicate presence of inbreeding depression that threatens its long-term survival. We conclude there is a continued region-wide threat to I. delicatissima and highlight the need for immediate conservation action to stop the continuing spread of Green Iguanas and to eliminate hybridization from St. Eustatius.

Genetic structure at three spatial scales is consistent with limited philopatry in Ricord's Rock Iguanas (Cyclura ricordii)

Cyclura ricordii is an endemic iguana from Hispaniola Island and is threatened on the IUCN Red List. The main threats are predation by introduced mammals, habitat destruction, and hunting pressure. The present study focused on two nesting sites from Pedernales Province in the Dominican Republic. The hypothesis that natal philopatry influences dispersal and nest-site selection was tested. Monitoring and sampling took place in 2012 and 2013. Polymorphic markers were used to evaluate whether natal philopatry limits dispersal at multiple spatial scales. Ripley's K revealed that nests were significantly clustered at multiple scales, when both nesting sites were considered and within each nesting site. This suggests a patchy, nonrandom distribution of nests within nest sites. Hierarchical AMOVA revealed that nest-site aggregations did not explain a significant portion of genetic variation within nesting sites. However, a small but positive correlation between geographic and genetic distance was detected using a Mantel's test. Hence, the relationship between geographic distance and genetic distance among hatchlings within nest sites, while detectable, was not strong enough to have a marked effect on fine-scale genetic structure. Spatial and genetic data combined determined that the nesting sites included nesting females from multiple locations, and the hypothesis of "natal philopatry" was not supported because females nesting in the same cluster were no more closely related to each other than to other females from the same nesting site. These findings imply that nesting aggregations are more likely associated with cryptic habitat variables contributing to optimal nesting conditions.

Age-dependent, negative heterozygosity-fitness correlations and local effects in an endangered Caribbean reptile, Iguana delicatissima

Inbreeding depression can have alarming impacts on threatened species with small population sizes. Assessing inbreeding has therefore become an important focus of conservation research. In this study, heterozygosity-fitness correlations (HFCs) were measured by genotyping 7 loci in 83 adult and 184 hatchling Lesser Antillean Iguanas, Iguana delicatissima, at a communal nesting site in Dominica to assess the role of inbreeding depression on hatchling fitness and recruitment to the adult population in this endangered species. We found insignificant correlations between multilocus heterozygosity and multiple fitness proxies in hatchlings and adults. Further, multilocus heterozygosity did not differ significantly between hatchlings and adults, which suggests that the survivorship of homozygous hatchlings does not differ markedly from that of their heterozygous counterparts. However, genotypes at two individual loci were correlated with hatching date, a finding consistent with the linkage between specific marker loci and segregating deleterious recessive alleles. These results provide only modest evidence that inbreeding depression influences the population dynamics of I. delicatissima on Dominica.

Genetic diversity and variation of Chinese fir from Fujian province and Taiwan, China, based on ISSR markers

Genetic diversity and variation among 11 populations of Chinese fir from Fujian province and Taiwan were assessed using inter-simple sequence repeat (ISSR) markers to reveal the evolutionary relationship in their distribution range in this report. Analysis of genetic parameters of the different populations showed that populations in Fujian province exhibited a greater level of genetic diversity than did the populations in Taiwan. Compared to Taiwan populations, significant limited gene flow were observed among Fujian populations. An UPGMA cluster analysis showed that the most individuals of Taiwan populations formed a single cluster, whereas 6 discrete clusters were formed by each population from Fujian. All populations were divided into 3 main groups and that all 5 populations from Taiwan were gathered into a subgroup combined with 2 populations, Dehua and Liancheng, formed one of the 3 main groups, which indicated relative stronger relatedness. It is supported by a genetic structure analysis. All those results are suggesting different levels of genetic diversity and variation of Chinese fir between Fujian and Taiwan, and indicating different patterns of evolutionary process and local environmental adaption.

Genetic diversity and structure in the Endangered Allen Cays Rock Iguana, Cyclura cychlura inornata

The Endangered Allen Cays Rock Iguana (Cyclura cychlura inornata) is endemic to the Allen Cays, a tiny cluster of islands in the Bahamas. Naturally occurring populations exist on only two cays (<4 ha each). However, populations of unknown origin were recently discovered on four additional cays. To investigate patterns of genetic variation among these populations, we analyzed nuclear and mitochondrial markers for 268 individuals. Analysis of three mitochondrial gene regions (2,328 bp) and data for eight nuclear microsatellite loci indicated low genetic diversity overall. Estimates of effective population sizes based on multilocus genotypes were also extremely low. Despite low diversity, significant population structuring and variation in genetic diversity measures were detected among cays. Genetic data confirm the source population for an experimentally translocated population while raising concerns regarding other, unauthorized, translocations. Reduced heterozygosity is consistent with a documented historical population decline due to overharvest. This study provides the first range-wide genetic analysis of this subspecies. We suggest strategies to maximize genetic diversity during ongoing recovery including additional translocations to establish assurance populations and additional protective measures for the two remaining natural populations.