How genomics can help biodiversity conservation.

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

A high-quality reference genome for the critically endangered Aeolian wall lizard, Podarcis raffonei.

The Aeolian wall lizard, Podarcis raffonei, is an endangered species endemic to the Aeolian archipelago, Italy, where it is present only in 3 tiny islets and a narrow promontory of a larger island. Because of the extremely limited area of occupancy, severe population fragmentation and observed decline, it has been classified as Critically Endangered by the International Union for the Conservation of Nature (IUCN). Using Pacific Biosciences (PacBio) High Fidelity (HiFi) long-read sequencing, Bionano optical mapping and Arima chromatin conformation capture sequencing (Hi-C), we produced a high-quality, chromosome-scale reference genome for the Aeolian wall lizard, including Z and W sexual chromosomes. The final assembly spans 1.51 Gb across 28 scaffolds with a contig N50 of 61.4 Mb, a scaffold N50 of 93.6 Mb, and a BUSCO completeness score of 97.3%. This genome constitutes a valuable resource for the species to guide potential conservation efforts and more generally for the squamate reptiles that are underrepresented in terms of available high-quality genomic resources.

Global phylogeography of ridley sea turtles (Lepidochelys spp.): evolution, demography, connectivity, and conservation.

Globally distributed marine taxa are well suited for investigations of biogeographic impacts on genetic diversity, connectivity, and population demography. The sea turtle genus Lepidochelys includes the wide-ranging and abundant olive ridley (L. olivacea), and the geographically restricted and 'Critically Endangered' Kemp's ridley (L. kempii). To investigate their historical biogeography, we analyzed a large dataset of mitochondrial DNA (mtDNA) sequences from olive (n = 943) and Kemp's (n = 287) ridleys, and genotyped 15 nuclear microsatellite loci in a global sample of olive ridleys (n = 285). We found that the ridley species split ~ 7.5 million years ago, before the Panama Isthmus closure. The most ancient mitochondrial olive ridley lineage, located in the Indian Ocean, was dated to ~ 2.2 Mya. Both mitochondrial and nuclear markers revealed significant structure for olive ridleys between Atlantic (ATL), East Pacific (EP), and Indo-West Pacific (IWP) areas. However, the divergence of mtDNA clades was very recent (< 1 Mya) with low within- clade diversity, supporting a recurrent extinction-recolonization model for these ocean regions. All data showed that ATL and IWP groups were more closely related than those in the EP, with mtDNA data supporting recent recolonization of the ATL from the IWP. Individual olive ridley dispersal between the ATL, EP, and IN/IWP could be interpreted as more male- than female-biased, and genetic diversity was lowest in the Atlantic Ocean. All populations showed signs of recent expansion, and estimated time frames were concordant with their recent colonization history. Investigating species abundance and distribution changes over time is central to evolutionary biology, and this study provides a historical biogeographic context for marine vertebrate conservation and management. Supplementary Information: The online version contains supplementary material available at 10.1007/s10592-022-01465-3.

Population structure, genomic diversity and demographic history of Komodo dragons inferred from whole-genome sequencing.

Population and conservation genetics studies have greatly benefited from the development of new techniques and bioinformatic tools associated with next-generation sequencing. Analysis of extensive data sets from whole-genome sequencing of even a few individuals allows the detection of patterns of fine-scale population structure and detailed reconstruction of demographic dynamics through time. In this study, we investigated the population structure, genomic diversity and demographic history of the Komodo dragon (Varanus komodoensis), the world's largest lizard, by sequencing the whole genomes of 24 individuals from the five main Indonesian islands comprising the entire range of the species. Three main genomic groups were observed. The populations of the Island of Komodo and the northern coast of Flores, in particular, were identified as two distinct conservation units. Degrees of genomic divergence among island populations were interpreted as a result of changes in sea level affecting connectivity across islands. Demographic inference suggested that Komodo dragons probably experienced a relatively steep population decline over the last million years, reaching a relatively stable Ne during the Saalian glacial cycle (400-150 thousand years ago) followed by a rapid Ne decrease. Genomic diversity of Komodo dragons was similar to that found in endangered or already extinct reptile species. Overall, this study provides an example of how whole-genome analysis of a few individuals per population can help define population structure and intraspecific demographic dynamics. This is particularly important when applying population genomics data to conservation of rare or elusive endangered species.

Insights into Emydid Turtle Cytogenetics: The European Pond Turtle as a Model Species.

Our knowledge of Testudines evolution is limited by the lack of modern cytogenetic data. Compared to other reptiles, there is little information even on chromosome banding, let alone molecular cytogenetic data. Here, we provide detailed information on the karyotype of the European pond turtle Emys orbicularis, a model Emydidae, employing both chromosome banding and molecular cytogenetics. We provide a high-resolution G-banded karyotype and a map of rDNA genes and telomeric sequences using fluorescence in situ hybridization. We test hypotheses of sex-determining mechanisms in Emys by comparative genomic hybridization to determine if Emys has a cryptic sex-specific region. Our results provide valuable data to guide future efforts on genome sequencing and anchoring in Emydidae and for understanding karyotype evolution in Testudines.

Isolating Chromosomes of the Komodo Dragon: New Tools for Comparative Mapping and Sequence Assembly.

We developed new tools to build a high-quality chromosomal map of the Komodo dragon (Varanus komodoensis) available for cross-species phylogenomic analyses. First, we isolated chromosomes by flow sorting and determined the chromosome content of each flow karyotype peak by FISH. We then isolated additional Komodo dragon chromosomes by microdissection and amplified chromosome-specific DNA pools. The chromosome-specific DNA pools can be sequenced, assembled, and mapped by next-generation sequencing technology. The chromosome-specific paint probes can be used to investigate karyotype evolution through cross-species chromosome painting. Overall, the set of chromosome-specific DNA pools of V. komodoensis provides new tools for detailed phylogenomic analyses of Varanidae and squamates in general.

Host-microbiota interactions shed light on mortality events in the striped venus clam Chamelea gallina.

Mass mortalities due to disease outbreaks have recently affected a number of major taxa in marine ecosystems. Climate- and pollution-induced stress may compromise host immune defenses, increasing the risk of opportunistic diseases. Despite growing evidence that mass mortality events affecting marine species worldwide are strongly influenced by the interplay of numerous environmental factors, the reductionist approaches most frequently used to investigate these factors hindered the interpretation of these multifactorial pathologies. In this study, we propose a broader approach based on the combination of RNA-sequencing and 16S microbiota analyses to decipher the factors underlying mass mortality in the striped venus clam, Chamelea gallina, along the Adriatic coast. On one hand, gene expression profiling and functional analyses of microbial communities showed the over-expression of several genes and molecular pathways involved in xenobiotic metabolism, suggesting potential chemical contamination in mortality sites. On the other hand, the down-regulation of several genes involved in immune and stress response, and the over-representation of opportunistic pathogens such as Vibrio and Photobacterium spp. indicates that these microbial species may take advantage of compromised host immune pathways and defense mechanisms that are potentially affected by chemical exposure, resulting in periodic mortality events. We propose the application of our approach to interpret and anticipate the risks inherent in the combined effects of pollutants and microbes on marine animals in today's rapidly changing environment.

Conserved sex chromosomes and karyotype evolution in monitor lizards (Varanidae).

Despite their long history with the basal split dating back to the Eocene, all species of monitor lizards (family Varanidae) studied so far share the same chromosome number of 2n = 40. However, there are differences in the morphology of the macrochromosome pairs 5-8. Further, sex determination, which revealed ZZ/ZW sex microchromosomes, was studied only in a few varanid species and only with techniques that did not test their homology. The aim of this study was to (i) test if cryptic interchromosomal rearrangements of larger chromosomal blocks occurred during the karyotype evolution of this group, (ii) contribute to the reconstruction of the varanid ancestral karyotype, and (iii) test homology of sex chromosomes among varanids. We investigated these issues by hybridizing flow sorted chromosome paints from Varanus komodoensis to metaphases of nine species of monitor lizards. The results show that differences in the morphology of the chromosome pairs 5-8 can be attributed to intrachromosomal rearrangements, which led to transitions between acrocentric and metacentric chromosomes in both directions. We also documented the first case of spontaneous triploidy among varanids in Varanus albigularis. The triploid individual was fully grown, which demonstrates that polyploidization is compatible with life in this lineage. We found that the W chromosome differs between species in size and heterochromatin content. The varanid Z chromosome is clearly conserved in all the analyzed species. Varanids, in addition to iguanas, caenophidian snakes, and lacertid lizards, are another squamate group with highly conserved sex chromosomes over a long evolutionary time.

Exploring mechanisms and origins of reduced dispersal in island Komodo dragons.

Loss of dispersal typifies island biotas, but the selective processes driving this phenomenon remain contentious. This is because selection via, both indirect (e.g. relaxed selection or island syndromes) and direct (e.g. natural selection or spatial sorting) processes may be involved, and no study has yet convincingly distinguished between these alternatives. Here, we combined observational and experimental analyses of an island lizard, the Komodo dragon (Varanus komodoensis, the world's largest lizard), to provide evidence for the actions of multiple processes that could contribute to island dispersal loss. In the Komodo dragon, concordant results from telemetry, simulations, experimental translocations, mark-recapture, and gene flow studies indicated that despite impressive physical and sensory capabilities for long-distance movement, Komodo dragons exhibited near complete dispersal restriction: individuals rarely moved beyond the valleys they were born/captured in. Importantly, lizard site-fidelity was insensitive to common agents of dispersal evolution (i.e. indices of risk for inbreeding, kin and intraspecific competition, and low habitat quality) that consequently reduced survival of resident individuals. We suggest that direct selection restricts movement capacity (e.g. via benefits of spatial philopatry and increased costs of dispersal) alongside use of dispersal-compensating traits (e.g. intraspecific niche partitioning) to constrain dispersal in island species.

Genetic Pedigree Analysis of the Pilot Breeding Program for the Rediscovered Galapagos Giant Tortoise from Floreana Island.

An aim of many captive breeding programs is to increase population sizes for reintroduction and establishment of self-sustaining wild populations. Genetic analyses play a critical role in these programs: monitoring genetic variation, identifying the origin of individuals, and assigning parentage to track family sizes. Here, we use genetic pedigree analyses to examine 3 seasons of a pilot breeding program for the Floreana island Galapagos giant tortoise, C. niger, that had been declared extinct for ~150 years until individuals with mixed ancestry were recently discovered. We determined that 8 of 9 founding individuals were assigned parentage to at least 1 of 130 offspring produced, though there was considerable reproductive skew. In addition, we observed that genetic diversity of the progeny was lower than that of the founders. Despite the observed reproductive skew, we did not see evidence for assortative mating based on relatedness, but there was a trend toward reduced fitness when more related individuals bred. Finally, we found that the majority of progeny had ancestry assigned to the Floreana species (mean ± SE = 0.51 ± 0.02), though individual estimates varied. The success of these pilot seasons bodes well for a larger breeding program to help restore the previously extinct tortoise from Floreana island. Future efforts should continue to monitor for reproductive skew and assortative mating to maintain allelic diversity. We would also recommend forming smaller breeding groups and rotating individuals among them to prevent long-term reproductive skew among pairs.

Ecological allometries and niche use dynamics across Komodo dragon ontogeny.

Ontogenetic allometries in ecological habits and niche use are key responses by which individuals maximize lifetime fitness. Moreover, such allometries have significant implications for how individuals influence population and community dynamics. Here, we examined how body size variation in Komodo dragons (Varanus komodoensis) influenced ecological allometries in their: (1) prey size preference, (2) daily movement rates, (3) home range area, and (4) subsequent niche use across ontogeny. With increased body mass, Komodo dragons increased prey size with a dramatic switch from small (≤10 kg) to large prey (≥50 kg) in lizards heavier than 20 kg. Rates of foraging movement were described by a non-linear concave down response with lizard increasing hourly movement rates up until ∼20 kg body mass before decreasing daily movement suggesting reduced foraging effort in larger lizards. In contrast, home range area exhibited a sigmoid response with increased body mass. Intrapopulation ecological niche use and overlap were also strongly structured by body size. Thus, ontogenetic allometries suggest Komodo dragon's transition from a highly active foraging mode exploiting small prey through to a less active sit and wait feeding strategy focused on killing large ungulates. Further, our results suggest that as body size increases across ontogeny, the Komodo dragon exhibited marked ontogenetic niche shifts that enabled it to function as an entire vertebrate predator guild by exploiting prey across multiple trophic levels.

Lineage fusion in Galápagos giant tortoises.

Although many classic radiations on islands are thought to be the result of repeated lineage splitting, the role of past fusion is rarely known because during these events, purebreds are rapidly replaced by a swarm of admixed individuals. Here, we capture lineage fusion in action in a Galápagos giant tortoise species, Chelonoidis becki, from Wolf Volcano (Isabela Island). The long generation time of Galápagos tortoises and dense sampling (841 individuals) of genetic and demographic data were integral in detecting and characterizing this phenomenon. In C. becki, we identified two genetically distinct, morphologically cryptic lineages. Historical reconstructions show that they colonized Wolf Volcano from Santiago Island in two temporally separated events, the first estimated to have occurred ~199 000 years ago. Following arrival of the second wave of colonists, both lineages coexisted for approximately ~53 000 years. Within that time, they began fusing back together, as microsatellite data reveal widespread introgressive hybridization. Interestingly, greater mate selectivity seems to be exhibited by purebred females of one of the lineages. Forward-in-time simulations predict rapid extinction of the early arriving lineage. This study provides a rare example of reticulate evolution in action and underscores the power of population genetics for understanding the past, present and future consequences of evolutionary phenomena associated with lineage fusion.

Exploring the effects of seasonality and chemical pollution on the hepatopancreas transcriptome of the Manila clam.

The assessment of marine environmental health is a complex but fundamental task both for ecosystem conservation and food safety related to the human consumption of marine products. Manila clams inhabiting the Venice Lagoon constitute an excellent case study for evaluating the effects of complex mixtures of industrial and urban effluents on aquatic organisms. Clams were collected in different seasons at four locations within the Venice Lagoon. The sampling sites were characterized by a range of pollutant concentrations and included Porto Marghera, a highly polluted industrial area where clam harvesting for human consumption is strictly forbidden. Pooled soft tissues were subjected to mass spectroscopy analysis to measure the concentrations of PCDDs/PCDFs/PCBs-DL, PCBs, PBDEs, HCB and PAHs, and pooled digestive gland samples were used for gene expression profiling. While seasonal variation was found to be responsible for the largest proportion of transcriptional changes, significance analysis of microarrays quantitative correlation analysis identified 162 transcripts that were correlated with at least one class of chemicals measured in the samples from the four different sampling sites. Prediction Analysis of Microarrays (PAM) identified a minimal set of seven genes that correctly assigned samples collected in the restricted polluted area (Porto Marghera), independent of the season in which they were collected. An integrated approach combining transcriptomics and chemical analyses of the Manila clam provided a global picture of how Manila clams respond to complex mixtures of xenobiotics and their interplay with other biotic and abiotic factors. We were also able to identify gene expression signatures for different classes of chemicals and a set of robust biomarkers of exposure to these chemicals.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles.

BACKGROUND: Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles. RESULTS: The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species. CONCLUSIONS: Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

A chromosome-level reference genome and pangenome for barn swallow population genomics.

Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

Parthenogenesis in Komodo dragons.

Parthenogenesis, the production of offspring without fertilization by a male, is rare in vertebrate species, which usually reproduce after fusion of male and female gametes. Here we use genetic fingerprinting to identify parthenogenetic offspring produced by two female Komodo dragons (Varanus komodoensis) that had been kept at separate institutions and isolated from males; one of these females subsequently produced additional offspring sexually. This reproductive plasticity indicates that female Komodo dragons may switch between asexual and sexual reproduction, depending on the availability of a mate--a finding that has implications for the breeding of this threatened species in captivity. Most zoos keep only females, with males being moved between zoos for mating, but perhaps they should be kept together to avoid triggering parthenogenesis and thereby decreasing genetic diversity.

The era of reference genomes in conservation genomics.

Progress in genome sequencing now enables the large-scale generation of reference genomes. Various international initiatives aim to generate reference genomes representing global biodiversity. These genomes provide unique insights into genomic diversity and architecture, thereby enabling comprehensive analyses of population and functional genomics, and are expected to revolutionize conservation genomics.

Transcriptome sequencing and microarray development for the Manila clam, Ruditapes philippinarum: genomic tools for environmental monitoring.

BACKGROUND: The Manila clam, Ruditapes philippinarum, is one of the major aquaculture species in the world and a potential sentinel organism for monitoring the status of marine ecosystems. However, genomic resources for R. philippinarum are still extremely limited. Global analysis of gene expression profiles is increasingly used to evaluate the biological effects of various environmental stressors on aquatic animals under either artificial conditions or in the wild. Here, we report on the development of a transcriptomic platform for global gene expression profiling in the Manila clam. RESULTS: A normalized cDNA library representing a mixture of adult tissues was sequenced using a ultra high-throughput sequencing technology (Roche 454). A database consisting of 32,606 unique transcripts was constructed, 9,747 (30%) of which could be annotated by similarity. An oligo-DNA microarray platform was designed and applied to profile gene expression of digestive gland and gills. Functional annotation of differentially expressed genes between different tissues was performed by enrichment analysis. Expression of Natural Antisense Transcripts (NAT) analysis was also performed and bi-directional transcription appears a common phenomenon in the R. philippinarum transcriptome. A preliminary study on clam samples collected in a highly polluted area of the Venice Lagoon demonstrated the applicability of genomic tools to environmental monitoring. CONCLUSIONS: The transcriptomic platform developed for the Manila clam confirmed the high level of reproducibility of current microarray technology. Next-generation sequencing provided a good representation of the clam transcriptome. Despite the known limitations in transcript annotation and sequence coverage for non model species, sufficient information was obtained to identify a large set of genes potentially involved in cellular response to environmental stress.

Population genetic structure of Aldabra giant tortoises.

Evolution of population structure on islands is the result of physical processes linked to volcanism, orogenic events, changes in sea level, as well as habitat variation. We assessed patterns of genetic structure in the giant tortoise of the Aldabra atoll, where previous ecological studies suggested population subdivisions as a result of landscape discontinuity due to unsuitable habitat and island separation. Analysis of mitochondrial DNA (mtDNA) control region sequences and allelic variation at 8 microsatellite loci were conducted on tortoises sampled in 3 locations on the 2 major islands of Aldabra. We found no variation in mtDNA sequences. This pattern corroborated earlier work supporting the occurrence of a founding event during the last interglacial period and a further reduction in genetic variability during historical time. On the other hand, significant population structure recorded at nuclear loci suggested allopatric divergence possibly due to geographical barriers among islands and ecological partitions hindering tortoise movements within islands. This is the first attempt to study the population genetics of Aldabra tortoises, which are now at carrying capacity in an isolated terrestrial ecosystem where ecological factors appear to have a strong influence on population dynamics.