Molecular phylogeny and historical biogeography of West Indian boid snakes (Chilabothrus)

Phylogenomic data resolve the historical biogeography and ecomorphs of Neotropical forest lizards (Squamata, Diploglossidae)

Few studies have been conducted on the biogeography and phylogenetic relationships of Neotropical forest lizards (Diploglossidae) because of incomplete taxon sampling, conflicting datasets, and low statistical support at phylogenetic nodes. Here, we enhance a recent nine-gene dataset with a genomic dataset of 3,232 loci and 642,775 aligned base pairs. The resulting phylogeny includes 30 diploglossid species, 10 of the 11 genera, and the three subfamilies. It shows significant support for all supra-specific taxa in either maximum likelihood or Bayesian analyses or both. With this well-supported phylogeny, we further investigate the historical biogeography of the group and how diploglossids reached the Caribbean islands. Our analyses indicate that Antillean diploglossid lizards originated from at least two overwater dispersals from South America. Our tests for the strength of convergent evolution between morphologically similar taxa support the recognition of a soil and a tree ecomorph. In addition, we propose grass, ground, rock, and swamp ecomorphs for species in this family based on ecological and morphological data and analyses.

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.

Genetic diversity in the US ex situ populations of the endangered Puerto Rican Boa, Chilabothrus inornatus

The Puerto Rican Boa (Chilabothrus inornatus) was placed on the US Endangered Species List in 1970. Progress has been made since to clarify the recovery status of this species, though the design of a new recovery plan must include information regarding genetic variation within and among populations of this species. While measures of genetic diversity in wild populations of this species are finally becoming available, relative genetic diversity represented in ex situ populations is unknown, which hampers efforts to develop an ex situ species management plan. Here, we provide an analysis of genetic diversity in US public and private collections (zoos and breeders) using mitochondrial sequence data and five highly polymorphic nuclear microsatellite loci. We analyzed 50 boas from the US ex situ population and determined overall genetic diversity and relatedness among these individuals. We then compared these data to mitochondrial and microsatellite data obtained from 176 individuals from wild populations across the native range of the species. We found little inbreeding and a large amount of retained genetic diversity in the US ex situ population of C. inornatus relative to wild populations. Genetic diversity in the ex situ population is similar to that found in wild populations. Ours is only the second explicit attempt to characterize genetic diversity at the molecular level in ex situ populations of boid snakes. We anticipate that these results will inform current breeding strategies as well as offer additional information that will facilitate the continuation of ex situ conservation breeding or management in boas.

Biogeography of the Caribbean Cyrtognatha spiders

Island systems provide excellent arenas to test evolutionary hypotheses pertaining to gene flow and diversification of dispersal-limited organisms. Here we focus on an orbweaver spider genus Cyrtognatha (Tetragnathidae) from the Caribbean, with the aims to reconstruct its evolutionary history, examine its biogeographic history in the archipelago, and to estimate the timing and route of Caribbean colonization. Specifically, we test if Cyrtognatha biogeographic history is consistent with an ancient vicariant scenario (the GAARlandia landbridge hypothesis) or overwater dispersal. We reconstructed a species level phylogeny based on one mitochondrial (COI) and one nuclear (28S) marker. We then used this topology to constrain a time-calibrated mtDNA phylogeny, for subsequent biogeographical analyses in BioGeoBEARS of over 100 originally sampled Cyrtognatha individuals, using models with and without a founder event parameter. Our results suggest a radiation of Caribbean Cyrtognatha, containing 11 to 14 species that are exclusively single island endemics. Although biogeographic reconstructions cannot refute a vicariant origin of the Caribbean clade, possibly an artifact of sparse outgroup availability, they indicate timing of colonization that is much too recent for GAARlandia to have played a role. Instead, an overwater colonization to the Caribbean in mid-Miocene better explains the data. From Hispaniola, Cyrtognatha subsequently dispersed to, and diversified on, the other islands of the Greater, and Lesser Antilles. Within the constraints of our island system and data, a model that omits the founder event parameter from biogeographic analysis is less suitable than the equivalent model with a founder event.

Genomic timetree and historical biogeography of Caribbean island ameiva lizards (Pholidoscelis: Teiidae)

The phylogenetic relationships and biogeographic history of Caribbean island ameivas (Pholidoscelis) are not well-known because of incomplete sampling, conflicting datasets, and poor support for many clades. Here, we use phylogenomic and mitochondrial DNA datasets to reconstruct a well-supported phylogeny and assess historical colonization patterns in the group. We obtained sequence data from 316 nuclear loci and one mitochondrial marker for 16 of 19 extant species of the Caribbean endemic genus Pholidoscelis. Phylogenetic analyses were carried out using both concatenation and species tree approaches. To estimate divergence times, we used fossil teiids to calibrate a timetree which was used to elucidate the historical biogeography of these lizards. All phylogenetic analyses recovered four well-supported species groups (clades) recognized previously and supported novel relationships of those groups, including a (P. auberi + P. lineolatus) clade (western + central Caribbean), and a (P. exsul + P. plei) clade (eastern Caribbean). Divergence between Pholidoscelis and its sister clade was estimated to have occurred ~25 Ma, with subsequent diversification on Caribbean islands occurring over the last 11 Myr. Of the six models compared in the biogeographic analyses, the scenario which considered the distance among islands and allowed dispersal in all directions best fit the data. These reconstructions suggest that the ancestor of this group colonized either Hispaniola or Puerto Rico from Middle America. We provide a well-supported phylogeny of Pholidoscelis with novel relationships not reported in previous studies that were based on significantly smaller datasets. We propose that Pholidoscelis colonized the eastern Greater Antilles from Middle America based on our biogeographic analysis, phylogeny, and divergence time estimates. The closing of the Central American Seaway and subsequent formation of the modern Atlantic meridional overturning circulation may have promoted dispersal in this group.

Sex allocation and secondary sex ratio in Cuban boa (Chilabothrus angulifer): mother's body size affects the ratio between sons and daughters

Secondary sex ratios of animals with genetically determined sex may considerably deviate from equality. These deviations may be attributed to several proximate and ultimate factors. Sex ratio theory explains some of them as strategic decisions of mothers improving their fitness by selective investment in sons or daughters, e.g. local resource competition hypothesis (LRC) suggests that philopatric females tend to produce litters with male-biased sex ratios to avoid future competition with their daughters. Until now, only little attention has been paid to examine predictions of sex ratio theory in snakes possessing genetic sex determination and exhibiting large variance in allocation of maternal investment. Cuban boa is an endemic viviparous snake producing large-bodied newborns (∼200 g). Extremely high maternal investment in each offspring increases importance of sex allocation. In a captive colony, we collected breeding records of 42 mothers, 62 litters and 306 newborns and examined secondary sex ratios (SR) and sexual size dimorphism (SSD) of newborns. None of the examined morphometric traits of neonates appeared sexually dimorphic. The sex ratio was slightly male biased (174 males versus 132 females) and litter sex ratio significantly decreased with female snout-vent length. We interpret this relationship as an additional support for LRC as competition between mothers and daughters increases with similarity of body sizes between competing snakes.

Genetic divergence and diversity in the Mona and Virgin Islands Boas, Chilabothrus monensis (Epicrates monensis) (Serpentes: Boidae), West Indian snakes of special conservation concern

Habitat fragmentation reduces the extent and connectivity of suitable habitats, and can lead to changes in population genetic structure. Limited gene flow among isolated demes can result in increased genetic divergence among populations, and decreased genetic diversity within demes. We assessed patterns of genetic variation in the Caribbean boa Chilabothrus monensis (Epicrates monensis) using two mitochondrial and seven nuclear markers, and relying on the largest number of specimens of these snakes examined to date. Two disjunct subspecies of C. monensis are recognized: the threatened C. m. monensis, endemic to Mona Island, and the rare and endangered C. m. granti, which occurs on various islands of the Puerto Rican Bank. Mitochondrial and nuclear markers revealed unambiguous genetic differences between the taxa, and coalescent species delimitation methods indicated that these snakes likely are different evolutionary lineages, which we recognize at the species level, C. monensis and C. granti. All examined loci in C. monensis (sensu stricto) are monomorphic, which may indicate a recent bottleneck event. Each population of C. granti exclusively contains private mtDNA haplotypes, but five of the seven nuclear genes assayed are monomorphic, and nucleotide diversity is low in the two remaining markers. The faster pace of evolution of mtDNA possibly reflects the present-day isolation of populations of C. granti, whereas the slower substitution rate of nuDNA may instead mirror the relatively recent episodes of connectivity among the populations facilitated by the lower sea level during the Pleistocene. The small degree of overall genetic variation in C. granti suggests that demes of this snake could be managed as a single unit, a practice that would significantly increase their effective population size.

Continental monophyly of cichlid fishes and the phylogenetic position of Heterochromis multidens

The incredibly species-rich cichlid fish faunas of both the Neotropics and Africa are generally thought to be reciprocally monophyletic. However, the phylogenetic affinity of the African cichlid Heterochromis multidens is ambiguous, and this distinct lineage could make African cichlids paraphyletic. In past studies, Heterochromis has been variously suggested to be one of the earliest diverging lineages within either the Neotropical or the African cichlid radiations, and it has even been hypothesized to be the sister lineage to a clade containing all Neotropical and African cichlids. We examined the phylogenetic relationships among a representative sample of cichlids with a dataset of 29 nuclear loci to assess the support for the different hypotheses of the phylogenetic position of Heterochromis. Although individual gene trees in some instances supported alternative relationships, a majority of gene trees, integration of genes into species trees, and hypothesis testing of putative topologies all supported Heterochromis as belonging to the clade of African cichlids.