Chromosome-Level Reference Genome of the Ponza Grayling (Hipparchia sbordonii), an Italian Endemic and Endangered Butterfly.

Islands are crucial evolutionary hotspots, providing unique opportunities for differentiation of novel biodiversity and long-term segregation of endemic species. Islands are also fragile ecosystems, where biodiversity is more exposed to environmental and anthropogenic pressures than on continents. The Ponza grayling, Hipparchia sbordonii, is an endemic butterfly species that is currently found only in two tiny islands of the Pontine archipelago, off the coast of Italy, occupying an area smaller than 10 km2. It has been classified as Endangered (IUCN) because of the extremely limited area of occurrence, population fragmentation, and the recent demographic decline. Thanks to a combination of different assemblers of long and short genomic reads, bulk transcriptome RNAseq, and synteny analysis with phylogenetically close butterflies, we produced a highly contiguous, chromosome-scale annotated reference genome for the Ponza grayling, including 28 autosomes and the Z sexual chromosomes. The final assembly spanned 388.61 Gb with a contig N50 of 14.5 Mb and a BUSCO completeness score of 98.5%. Synteny analysis using four other butterfly species revealed high collinearity with Hipparchia semele and highlighted 10 intrachromosomal inversions longer than 10 kb, of which two appeared on the lineage leading to H. sbordonii. Our results show that a chromosome-scale reference genome is attainable also when chromatin conformation data may be impractical or present specific technical challenges. The high-quality genomic resource for H. sbordonii opens up new opportunities for the accurate assessment of genetic diversity and genetic load and for the investigations of the genomic novelties characterizing the evolutionary path of this endemic island species.

Unravelling the mystery of endemic versus translocated populations of the endangered Australian lungfish (Neoceratodus forsteri).

The Australian lungfish is a primitive and endangered representative of the subclass Dipnoi. The distribution of this species is limited to south-east Queensland, with some populations considered endemic and others possibly descending from translocations in the late nineteenth century shortly after European discovery. Attempts to resolve the historical distribution of this species have met with conflicting results based on descriptive genetic studies. Understanding if all populations are endemic or some are the result of, or influenced by, translocation events, has implications for conservation management. In this work, we analysed the genetic variation at three types of markers (mtDNA genomes, 11 STRs and 5196 nuclear SNPs) using the approximate Bayesian computation (ABC) algorithm to compare several demographic models. We postulated different contributions of Mary River and Burnett River gene pools into the Brisbane River and North Pine River populations, related to documented translocation events. We ran the analysis for each marker type separately, and we also estimated the posterior probabilities of the models combining the markers. Nuclear SNPs have the highest power to correctly identify the true model among the simulated datasets (where the model was known), but different marker types typically provided similar answers. The most supported demographic model able to explain the real dataset implies that an endemic gene pool is still present in the Brisbane and North Pine Rivers and coexists with the gene pools derived from past documented translocation events. These results support the view that ABC modelling can be useful to reconstruct complex historical translocation events with contemporary implications, and will inform ongoing conservation efforts for the endangered and iconic Australian lungfish.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes.

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758.

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

Fitness consequences and ancestry loss in the Apennine brown bear after a simulated genetic rescue intervention.

Reduction in population size, with its predicted effects on population fitness, is the most alarming anthropogenic impact on endangered species. By introducing compatible individuals, genetic rescue (GR) is a promising but debated approach for reducing the genetic load unmasked by inbreeding and for restoring the fitness of declining populations. Although GR can improve genetic diversity and fitness, it can also produce loss of ancestry, hampering local adaptation, or replace with introduced variants the unique genetic pools evolved in endemic groups. We used forward genetic simulations based on empirical genomic data to assess fitness benefits and loss of ancestry risks of GR in the Apennine brown bear (Ursus arctos marsicanus). There are approximately 50 individuals of this isolated subspecies, and they have lower genetic diversity and higher inbreeding than other European brown bears, and GR has been suggested to reduce extinction risks. We compared 10 GR scenarios in which the number and genetic characteristics of migrants varied with a non-GR scenario of simple demographic increase due to nongenetic factors. The introduction of 5 individuals of higher fitness or lower levels of deleterious mutations than the target Apennine brown bear from a larger European brown bear population produced a rapid 10-20% increase in fitness in the subspecies and up to 22.4% loss of ancestry over 30 generations. Without a contemporary demographic increase, fitness started to decline again after a few generations. Doubling the population size without GR gradually increased fitness to a comparable level, but without losing ancestry, thus resulting in the best strategy for the Apennine brown bear conservation. Our results highlight the importance for management of endangered species of realistic forward simulations grounded in empirical whole-genome data.

Consecuencias en la aptitud y pérdida de ascendencia del oso pardo de los Apeninos después de un rescate genético simulado Resumen La reducción del tamaño poblacional, con los previsibles efectos sobre su aptitud, es el impacto antropogénico más alarmante sobre las especies amenazadas. Mediante la introducción de individuos compatibles, el rescate genético (RG) es una estrategia prometedora para reducir la carga genética revelada por la endogamia y restaurar la aptitud de las poblaciones en declive, aunque todavía se debate la eficiencia de esta. Aunque el RG puede mejorar la diversidad genética y la aptitud, también puede producir pérdida de ascendencia, lo que puede dificultar la adaptación local, o sustituir con variantes introducidas por los migrantes los acervos genéticos únicos que han evolucionado en grupos endémicos. En este trabajo realizamos simulaciones genéticas a futuro basadas en datos genómicos empíricos para evaluar los beneficios del RG en términos de aptitud y los riesgos de la pérdida de ascendencia en el oso pardo de los Apeninos (Ursus arctos marsicanus). Quedan aproximadamente 50 individuos de esta subespecie aislada que cuentan con una menor diversidad genética y un mayor nivel de endogamia comparado con otros osos pardos europeos y se ha sugerido que el RG podria reducir el riesgo de extinción de esta población. Comparamos 10 escenarios de RG en los que variaban el número y las características genéticas de los osos migrantes con un escenario sin RG con aumento demográfico causado por factores no genéticos. La introducción de 5 individuos procedentes de una población europea de oso pardo con mayor aptitud o niveles menores de mutaciones deletéreas que el oso pardo de los Apeninos produjo un rápido aumento de la aptitud del 10-20% en la subespecie y hasta un 22.4% de pérdida de ascendencia durante 30 generaciones. En las simulaciones sin un aumento demográfico, la aptitud empezó a disminuir de nuevo después de unas pocas generaciones. La duplicación del tamaño de la población sin RG aumentó gradualmente la aptitud hasta un nivel comparable al de algunos escenarios de RG, pero sin pérdida de ascendencia, por lo que parece ser la mejor estrategia para la conservación del oso pardo de los Apeninos. Nuestros resultados resaltan la importancia que tienen las simulaciones realistas a futuro basadas en datos empíricos del genoma completo para la gestión de especies amenazadas.

Tracing Eastern Wolf Origins From Whole-Genome Data in Context of Extensive Hybridization.

Southeastern Canada is inhabited by an amalgam of hybridizing wolf-like canids, raising fundamental questions regarding their taxonomy, origins, and timing of hybridization events. Eastern wolves (Canis lycaon), specifically, have been the subject of significant controversy, being viewed as either a distinct taxonomic entity of conservation concern or a recent hybrid of coyotes (C. latrans) and grey wolves (C. lupus). Mitochondrial DNA analyses show some evidence of eastern wolves being North American evolved canids. In contrast, nuclear genome studies indicate eastern wolves are best described as a hybrid entity, but with unclear timing of hybridization events. To test hypotheses related to these competing findings we sequenced whole genomes of 25 individuals, representative of extant Canadian wolf-like canid types of known origin and levels of contemporary hybridization. Here we present data describing eastern wolves as a distinct taxonomic entity that evolved separately from grey wolves for the past ∼67,000 years with an admixture event with coyotes ∼37,000 years ago. We show that Great Lakes wolves originated as a product of admixture between grey wolves and eastern wolves after the last glaciation (∼8,000 years ago) while eastern coyotes originated as a product of admixture between "western" coyotes and eastern wolves during the last century. Eastern wolf nuclear genomes appear shaped by historical and contemporary gene flow with grey wolves and coyotes, yet evolutionary uniqueness remains among eastern wolves currently inhabiting a restricted range in southeastern Canada.

Selection and adaptive introgression guided the complex evolutionary history of the European common bean.

Domesticated crops have been disseminated by humans over vast geographic areas. Common bean (Phaseolus vulgaris L.) was introduced in Europe after 1492. Here, by combining whole-genome profiling, metabolic fingerprinting and phenotypic characterisation, we show that the first common bean cultigens successfully introduced into Europe were of Andean origin, after Francisco Pizarro's expedition to northern Peru in 1529. We reveal that hybridisation, selection and recombination have shaped the genomic diversity of the European common bean in parallel with political constraints. There is clear evidence of adaptive introgression into the Mesoamerican-derived European genotypes, with 44 Andean introgressed genomic segments shared by more than 90% of European accessions and distributed across all chromosomes except PvChr11. Genomic scans for signatures of selection highlight the role of genes relevant to flowering and environmental adaptation, suggesting that introgression has been crucial for the dissemination of this tropical crop to the temperate regions of Europe.

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.

A validated protocol for eDNA-based monitoring of within-species genetic diversity in a pond-breeding amphibian.

In light of the dramatic decline in amphibian biodiversity, new cost-efficient tools to rapidly monitor species abundance and population genetic diversity in space and time are urgently needed. It has been amply demonstrated that the use of environmental DNA (eDNA) for single-species detection and characterization of community composition can increase the precision of amphibian monitoring compared to traditional (observational) approaches. However, it has been suggested that the efficiency and accuracy of the eDNA approach could be further improved by more timely sampling; in addition, the quality of genetic diversity data derived from the same DNA has been confirmed in other vertebrate taxa, but not amphibians. Given the availability of previous tissue-based genetic data, here we use the common frog Rana temporaria Linnaeus, 1758 as our target species and an improved eDNA protocol to: (i) investigate differences in species detection between three developmental stages in various freshwater environments; and (ii) study the diversity of mitochondrial DNA (mtDNA) haplotypes detected in eDNA (water) samples, by amplifying a specific fragment of the COI gene (331 base pairs, bp) commonly used as a barcode. Our protocol proved to be a reliable tool for monitoring population genetic diversity of this species, and could be a valuable addition to amphibian conservation and wetland management.

Evidence of backcross inviability and mitochondrial DNA paternal leakage in sea turtle hybrids.

Hybridization is known to be part of many species' evolutionary history. Sea turtles have a fascinating hybridization system in which species separated by as much as 43 million years are still capable of hybridizing. Indeed, the largest nesting populations in Brazil of loggerheads (Caretta caretta) and hawksbills (Eretmochelys imbricata) have a high incidence of hybrids between these two species. A third species, olive ridleys (Lepidochelys olivacea), is also known to hybridize although at a smaller scale. Here, we used restriction site-associated DNA sequencing (RAD-Seq) markers, mitogenomes, and satellite-telemetry to investigate the patterns of hybridization and introgression in the Brazilian sea turtle population and their relationship with the migratory behaviours between feeding and nesting aggregations. We also explicitly test if the mixing of two divergent genomes in sea turtle hybrids causes mitochondrial paternal leakage. We developed a new species-specific PCR-assay capable of detecting mitochondrial DNA (mtDNA) inheritance from both parental species and performed ultra-deep sequencing to estimate the abundance of each mtDNA type. Our results show that all adult hybrids are first generation (F1) and most display a loggerhead migratory behaviour. We detected paternal leakage in F1 hybrids and different proportions of mitochondria from maternal and paternal species. Although previous studies showed no significant fitness decrease in hatchlings, our results support genetically-related hybrid breakdown possibly caused by cytonuclear incompatibility. Further research on hybrids from other populations in addition to Brazil and between different species will show if backcross inviability and mitochondrial paternal leakage is observed across sea turtle species.

From STRs to SNPs via ddRAD-seq: Geographic assignment of confiscated tortoises at reduced costs.

Assigning individuals to their source populations is crucial for conservation research, especially for endangered species threatened by illegal trade and translocations. Genetic assignment can be achieved with different types of molecular markers, but technical advantages and cost saving are recently promoting the shift from short tandem repeats (STRs) to single nucleotide polymorphisms (SNPs). Here, we designed, developed, and tested a small panel of SNPs for cost-effective geographic assignment of individuals with unknown origin of the endangered Mediterranean tortoise Testudo hermanni. We started by performing a ddRAD-seq experiment on 70 wild individuals of T. hermanni from 38 locations. Results obtained using 3182 SNPs are comparable to those previously obtained using STR markers in terms of genetic structure and power to identify the macro-area of origin. However, our SNPs revealed further insights into the substructure in Western populations, especially in Southern Italy. A small panel of highly informative SNPs was then selected and tested by genotyping 190 individuals using the KASP genotyping chemistry. All the samples from wild populations of known geographic origin were genetically re-assigned with high accuracy to the original population. This reduced SNPs panel represents an efficient molecular tool that enables individuals to be genotyped at low cost (less than €15 per sample) for geographical assignment and identification of hybrids. This information is crucial for the management in-situ of confiscated animals and their possible re-allocation in the wild. Our methodological pipeline can easily be extended to other species.

Genetic load: genomic estimates and applications in non-model animals.

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Identification of natural selection in genomic data with deep convolutional neural network.

BACKGROUND: With the increase in the size of genomic datasets describing variability in populations, extracting relevant information becomes increasingly useful as well as complex. Recently, computational methodologies such as Supervised Machine Learning and specifically Convolutional Neural Networks have been proposed to make inferences on demographic and adaptive processes using genomic data. Even though it was already shown to be powerful and efficient in different fields of investigation, Supervised Machine Learning has still to be explored as to unfold its enormous potential in evolutionary genomics. RESULTS: The paper proposes a method based on Supervised Machine Learning for classifying genomic data, represented as windows of genomic sequences from a sample of individuals belonging to the same population. A Convolutional Neural Network is used to test whether a genomic window shows the signature of natural selection. Training performed on simulated data show that the proposed model can accurately predict neutral and selection processes on portions of genomes taken from real populations with almost 90% accuracy.

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.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species.

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

Ancient genomes reveal early Andean farmers selected common beans while preserving diversity.

All crops are the product of a domestication process that started less than 12,000 years ago from one or more wild populations1,2. Farmers selected desirable phenotypic traits (such as improved energy accumulation, palatability of seeds and reduced natural shattering3) while leading domesticated populations through several more or less gradual demographic contractions2,4. As a consequence, the erosion of wild genetic variation5 is typical of modern cultivars, making them highly susceptible to pathogens, pests and environmental change6,7. The loss of genetic diversity hampers further crop improvement programmes to increase food production in a changing world, posing serious threats to food security8,9. Using both ancient and modern seeds, we analysed the temporal dynamics of genetic variation and selection during the domestication process of the common bean (Phaseolus vulgaris) in the southern Andes. Here, we show that most domestic traits were selected for before 2,500 years ago, with no or only minor loss of whole-genome heterozygosity. In fact, most of the changes at coding genes and linked regions that differentiate wild and domestic genomes are already present in the ancient genomes analysed here, and all ancient domestic genomes dated between 600 and 2,500 years ago are highly variable (at least as variable as modern genomes from the wild). Single seeds from modern cultivars show reduced variation when compared with ancient seeds, indicating that intensive selection within cultivars in the past few centuries probably partitioned ancestral variation within different genetically homogenous cultivars. When cultivars from different Andean regions are pooled, the genomic variation of the pool is higher than that observed in the pool of ancient seeds from north and central western Argentina. Considering that most desirable phenotypic traits are probably controlled by multiple polymorphic genes10, a plausible explanation of this decoupling of selection and genetic erosion is that early farmers applied a relatively weak selection pressure2 by using many phenotypically similar but genetically diverse individuals as parents. Our results imply that selection strategies during the past few centuries, as compared with earlier times, more intensively reduced genetic variation within cultivars and produced further improvements by focusing on a few plants carrying the traits of interest, at the cost of marked genetic erosion within Andean landraces.

Mapping the geographic origin of captive and confiscated Hermann's tortoises: A genetic toolkit for conservation and forensic analyses.

The illegal trade has been threatening tortoise populations worldwide for decades. Nowadays, however, DNA typing and forensic genetic approaches allow us to investigate the geographic origin of confiscated animals and to relocate them into the wild, providing that suitable molecular tools and reference data are available. Here we assess the suitability of a small panel of microsatellite markers to investigate patterns of illegal translocations and to assist forensic genetic applications in the endangered Mediterranean land tortoise Testudo hermanni hermanni. Specific allelic ladders were created for each locus and tested on several reference samples. We used the microsatellite panel to (i) increase our understanding of the population genetic structure in wild populations with new data from previously unsampled geographic areas (overall 461 wild individuals from 28 sampling sites); (ii) detect the presence of non-native individuals in wild populations; and (iii) identify the most likely geographic area of origin of 458 confiscated individuals hosted in Italian seizure and recovery centers. Our analysis initially identified six major genetic clusters corresponding to different geographic macro-areas along the Mediterranean range. Long-distance migrants among wild populations, due to translocations, were found and removed from the reference database. Assignment tests allowed us to allocate approximately 70 % of confiscated individuals of unknown origin to one of the six Mediterranean macro-areas. Most of the assigned tortoises belonged to the genetic cluster corresponding to the area where the respective captivity center was located. However, we also found evidence of long-distance origins of confiscated individuals, especially in centers along the Adriatic coast and facing the Balkan regions, a well-known source of illegally traded individuals. Our results clearly show that the microsatellite panel and the reference dataset can play a beneficial role in reintroduction and repatriation projects when confiscated individuals need to be re-assigned to their respective macro-area of origin before release, and can assist future forensic genetic applications in detecting the illegal trade and possession of Testudo hermanni individuals.

Population Dynamics and Structural Effects at Short and Long Range Support the Hypothesis of the Selective Advantage of the G614 SARS-CoV-2 Spike Variant.

SARS-CoV-2 epidemics quickly propagated worldwide, sorting virus genomic variants in newly established propagules of infections. Stochasticity in transmission within and between countries or an actual selective advantage could explain the global high frequency reached by some genomic variants. Using statistical analyses, demographic reconstructions, and molecular dynamics simulations, we show that the globally invasive G614 spike variant 1) underwent a significant demographic expansion in most countries explained neither by stochastic effects nor by overrepresentation in clinical samples, 2) increases the spike S1/S2 furin-like site conformational plasticity (short-range effect), and 3) modifies the internal motion of the receptor-binding domain affecting its cross-connection with other functional domains (long-range effect). Our results support the hypothesis of a selective advantage at the basis of the spread of the G614 variant, which we suggest may be due to structural modification of the spike protein at the S1/S2 proteolytic site, and provide structural information to guide the design of variant-specific drugs.