Gene Flow and Isolation in the Arid Nearctic Revealed by Genomic Analyses of Desert Spiny Lizards.

The opposing forces of gene flow and isolation are two major processes shaping genetic diversity. Understanding how these vary across space and time is necessary to identify the environmental features that promote diversification. The detection of considerable geographic structure in taxa from the arid Nearctic has prompted research into the drivers of isolation in the region. Several geographic features have been proposed as barriers to gene flow, including the Colorado River, Western Continental Divide, and a hypothetical Mid-Peninsular Seaway in Baja California. However, recent studies suggest that the role of barriers in genetic differentiation may have been overestimated when compared to other mechanisms of divergence. In this study, we infer historical and spatial patterns of connectivity and isolation in Desert Spiny Lizards (Sceloporus magister) and Baja Spiny Lizards (S. zosteromus), which together form a species complex composed of parapatric lineages with wide distributions in arid western North America. Our analyses incorporate mitochondrial sequences, genomic-scale data, and past and present climatic data to evaluate the nature and strength of barriers to gene flow in the region. Our approach relies on estimates of migration under the multispecies coalescent to understand the history of lineage divergence in the face of gene flow. Results show that the S. magister complex is geographically structured, but we also detect instances of gene flow. The Continental Divide is a strong barrier to gene flow, while the Colorado River is more permeable. Analyses yield conflicting results for the catalyst of differentiation of peninsular lineages in S. zosteromus. Our study shows how large-scale genomic data for thoroughly sampled species can shed new light on biogeography. Furthermore, our approach highlights the need for the combined analysis of multiple sources of evidence to adequately characterize the drivers of divergence.

Ontogenetic variation in the ecological stoichiometry of 10 fish species during early development.

The chemical composition and stoichiometry of vertebrate bodies changes greatly during ontogeny as phosphorus-rich bones form, but we know little about the variation among species during early development. Such variation is important because element ratios in animal bodies influence which element limits growth and how animals contribute to nutrient cycling. We quantified ontogenetic variation from embryos through 2-3 months of age in 10 species of fish in six different families, ranging in adult size from 73 to 720 mm in length. We measured whole-body concentrations (percentage of dry mass) and ratios of carbon (C), nitrogen (N), and phosphorus (P) as fish developed. We also quantified whole-body concentrations of calcium (Ca), because Ca should reflect bone development, and RNA, which can be a major pool of body P. To account for interspecific differences in adult size, we also examined how trends changed with relative size, defined as body length divided by adult length. Ontogenetic changes in body composition and ratios were relatively similar among species and were more similar when expressed as a function of relative size compared to age. Body P increased rapidly in all species (likely because of bone development) from embryos until individuals were ~5%-8% of adult size. Body N also increased, while body C, C:N, C:P, and N:P all decreased over this period. Body Ca increased with development but was more variable among species. Body RNA was low in embryos, increased rapidly in young larvae, then decreased as fish reached 5%-8% of adult size. After fish were about 5%-8% of adult size, changes in body composition were relatively slight for all elements and ratios. These results reveal a consistency in the dynamics of body stoichiometry during early ontogeny, presumably because of similar constraints on the allocation of elements to bones and other body pools. Because most changes occur when individuals are <1 month old (<10% of adult size for that species), early ontogenetic variation in body stoichiometry may be especially important for growth limitation of individuals and ecosystem-level nutrient cycling.

Signals of differential introgression in the genome of natural hybrids of Caribbean anoles.

Hybridization facilitates recombination between divergent genetic lineages and can be shaped by both neutral and selective processes. Upon hybridization, loci with no net fitness effects introgress randomly from parental species into the genomes of hybrid individuals. Conversely, alleles from one parental species at some loci may provide a selective advantage to hybrids, resulting in patterns of introgression that do not conform to random expectations. We investigated genomic patterns of differential introgression in natural hybrids of two species of Caribbean anoles, Anolis pulchellus and A. krugi in Puerto Rico. Hybrids exhibit A. pulchellus phenotypes but possess A. krugi mitochondrial DNA, originated from multiple, independent hybridization events, and appear to have replaced pure A. pulchellus across a large area in western Puerto Rico. Combining genome-wide SNP datasets with bioinformatic methods to identify signals of differential introgression in hybrids, we demonstrate that the genomes of hybrids are dominated by pulchellus-derived alleles and show only 10%-20% A. krugi ancestry. The majority of A. krugi loci in hybrids exhibit a signal of non-random differential introgression and include loci linked to genes involved in development and immune function. Three of these genes (delta like canonical notch ligand 1, jagged1 and notch receptor 1) affect cell differentiation and growth and interact with mitochondrial function. Our results suggest that differential non-random introgression for a subset of loci may be driven by selection favouring the inheritance of compatible mitochondrial and nuclear-encoded genes in hybrids.

A chromosome-level genome assembly and annotation of the desert horned lizard, Phrynosoma platyrhinos, provides insight into chromosomal rearrangements among reptiles.

BACKGROUND: The increasing number of chromosome-level genome assemblies has advanced our knowledge and understanding of macroevolutionary processes. Here, we introduce the genome of the desert horned lizard, Phrynosoma platyrhinos, an iguanid lizard occupying extreme desert conditions of the American southwest. We conduct analysis of the chromosomal structure and composition of this species and compare these features across genomes of 12 other reptiles (5 species of lizards, 3 snakes, 3 turtles, and 1 bird). FINDINGS: The desert horned lizard genome was sequenced using Illumina paired-end reads and assembled and scaffolded using Dovetail Genomics Hi-C and Chicago long-range contact data. The resulting genome assembly has a total length of 1,901.85 Mb, scaffold N50 length of 273.213 Mb, and includes 5,294 scaffolds. The chromosome-level assembly is composed of 6 macrochromosomes and 11 microchromosomes. A total of 20,764 genes were annotated in the assembly. GC content and gene density are higher for microchromosomes than macrochromosomes, while repeat element distributions show the opposite trend. Pathway analyses provide preliminary evidence that microchromosome and macrochromosome gene content are functionally distinct. Synteny analysis indicates that large microchromosome blocks are conserved among closely related species, whereas macrochromosomes show evidence of frequent fusion and fission events among reptiles, even between closely related species. CONCLUSIONS: Our results demonstrate dynamic karyotypic evolution across Reptilia, with frequent inferred splits, fusions, and rearrangements that have resulted in shuffling of chromosomal blocks between macrochromosomes and microchromosomes. Our analyses also provide new evidence for distinct gene content and chromosomal structure between microchromosomes and macrochromosomes within reptiles.

A photographic data set of reef and coral communities across Venezuela.

This dataset contains 2850 photographs of the seafloor in coral communities from Venezuela that were taken during 2017 and 2018. We used a hierarchical experimental design with four random factors representing four different spatial scales: (1) region (hundreds of kilometers), (2) localities (tens of kilometers), (2) reef sites (hundreds of meters) and (3) transects (a couple meters) across the Venezuelan coast. At each site, four 30-m transects were deployed parallel to the coastline, and 15 pictures were taken every other meter at each transect, containing an area of at least 80 × 90cm with enough resolution to identify benthic groups. This dataset covers spatial scales from a few meters to hundreds of kilometers; marine protected areas, and non-protected areas; coastal zones, continental and oceanic islands. These images have the potential to be further used for training researchers in benthic organisms identification, and training artificial intelligence classification algorithms. Also, they represent and updated baseline to perform spatial and temporal comparisons in Venezuela or further studies involving multiple spatial scales in the region.

The use of pseudo-multivariate standard error to improve the sampling design of coral monitoring programs.

The characteristics of coral reef sampling and monitoring are highly variable, with numbers of units and sampling effort varying from one study to another. Numerous works have been carried out to determine an appropriate effect size through statistical power; however, these were always from a univariate perspective. In this work, we used the pseudo multivariate dissimilarity-based standard error (MultSE) approach to assess the precision of sampling scleractinian coral assemblages in reefs of Venezuela between 2017 and 2018 when using different combinations of number of transects, quadrats and points. For this, the MultSE of 36 sites previously sampled was estimated, using four 30m-transects with 15 photo-quadrats each and 25 random points per quadrat. We obtained that the MultSE was highly variable between sites and is not correlated with the univariate standard error nor with the richness of species. Then, a subset of sites was re-annotated using 100 uniformly distributed points, which allowed the simulation of different numbers of transects per site, quadrats per transect and points per quadrat using resampling techniques. The magnitude of the MultSE stabilized by adding more transects, however, adding more quadrats or points does not improve the estimate. For this case study, the error was reduced by half when using 10 transects, 10 quadrats per transect and 25 points per quadrat. We recommend the use of MultSE in reef monitoring programs, in particular when conducting pilot surveys to optimize the estimation of the community structure.

On the importance of spatial scales on beta diversity of coral assemblages: a study from Venezuelan coral reefs.

Estimating variability across spatial scales has been a major issue in ecology because the description of patterns in space is extremely valuable to propose specific hypotheses to unveil key processes behind these patterns. This paper aims to estimate the variability of the coral assemblage structure at different spatial scales in order to determine which scales explain the largest variability on ß-diversity. For this, a fully-nested design including a series of hierarchical-random factors encompassing three spatial scales: (1) regions, (2) localities and (3) reefs sites across the Venezuelan territory. The variability among spatial scales was tested with a permutation-based analysis of variance (Permanova) based on Bray-Curtis index. Dispersion in species presence/absence across scales (i.e., ß-diversity) was tested with a PermDisp analysis based on Jaccard's index. We found the highest variability in the coral assemblage structure between sites within localities (Pseudo-F = 5.34; p-value = 0.001, CV = 35.10%). We also found that longitude (Canonical corr = 0.867, p = 0.001) is a better predictor of the coral assemblage structure in Venezuela, than latitude (Canonical corr = 0.552, p = 0.021). Largest changes in ß-diversity of corals occurred within sites (F = 2.764, df1= 35, df2 = 107, p = 0.045) and within localities (F = 4.438, df1= 6, df2 = 29, p = 0.026). Our results suggest that processes operating at spatial scales of hundreds of meters and hundreds of kilometers might both be critical to shape coral assemblage structure in Venezuela, whereas smaller scales (i.e., hundreds of meters) showed to be highly- important for the species turnover component of ß-diversity. This result highlights the importance of creating scale-adapted management actions in Venezuela and likely across the Caribbean region.

Systematic review and meta-analysis of 50 years of coral disease research visualized through the scope of network theory.

Coral disease research encompasses five decades of undeniable progress. Since the first descriptions of anomalous signs, we have come to understand multiple processes and environmental drivers that interact with coral pathologies. In order to gain a better insight into the knowledge we already have, we explored how key topics in coral disease research have been related to each other using network analysis. We reviewed 719 papers and conference proceedings published from 1965 to 2017. From each study, four elements determined our network nodes: (1) studied disease(s); (2) host genus; (3) marine ecoregion(s) associated with the study site; and (4) research objectives. Basic properties of this network confirmed that there is a set of specific topics comprising the majority of research. The top five diseases, genera, and ecoregions studied accounted for over 48% of the research effort in all cases. The community structure analysis identified 15 clusters of topics with different degrees of overlap among them. These clusters represent the typical sets of elements that appear together for a given study. Our results show that while some coral diseases have been studied considering multiple aspects, the overall trend is for most diseases to be understood under a limited range of approaches, e.g., bacterial assemblages have been considerably studied in Yellow and Black band diseases while immune response has been better examined for the aspergillosis-Gorgonia system. Thus, our challenge in the near future is to identify and resolve potential gaps in order to achieve a more comprehensive progress on coral disease research.