Effective double-digest RAD sequencing and genotyping despite large genome size

Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers

Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.

Genomics and conservation: Guidance from training to analyses and applications

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

Genomic evidence for West Antarctic Ice Sheet collapse during the Last Interglacial

The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.

Population genetic structure and hybrid zone analyses for species delimitation in the Japanese toad (Bufo japonicus)

Hybridization following secondary contact may produce different outcomes depending on the extent to which genetic diversity and reproductive barriers have accumulated during isolation. The Japanese toad, Bufo japonicus, is distributed on the main islands of Japan. In the present study, we applied multiplexed inter-simple sequence repeat genotyping by sequencing to achieve the fine-scale resolution of the genetic cluster in B. j. japonicus and B. j. formosus. We also elucidated hybridization patterns and gene flow degrees across contact zones between the clusters identified. Using SNP data, we found four genetic clusters in B. j. japonicus and B. j. formosus and three contact zones of the cluster pairs among these four clusters. The two oldest diverged lineages, B. j. japonicus and B. j. formosus, formed a narrow contact zone consistent with species distinctiveness. Therefore, we recommend that these two subspecies be elevated to the species level. In contrast, the less diverged pairs of two clusters in B. j. japonicus and B. j. formosus, respectively, admixed over a hundred kilometers, suggesting that they have not yet developed strong reproductive isolation and need to be treated as conspecifics. These results will contribute to resolving taxonomic confusion in Japanese toads.

CRISPR-Cas9-based repeat depletion for high-throughput genotyping of complex plant genomes

High-throughput genotyping enables the large-scale analysis of genetic diversity in population genomics and genome-wide association studies that combine the genotypic and phenotypic characterization of large collections of accessions. Sequencing-based approaches for genotyping are progressively replacing traditional genotyping methods because of the lower ascertainment bias. However, genome-wide genotyping based on sequencing becomes expensive in species with large genomes and a high proportion of repetitive DNA. Here we describe the use of CRISPR-Cas9 technology to deplete repetitive elements in the 3.76-Gb genome of lentil (Lens culinaris), 84% consisting of repeats, thus concentrating the sequencing data on coding and regulatory regions (single-copy regions). We designed a custom set of 566,766 gRNAs targeting 2.9 Gbp of repeats and excluding repetitive regions overlapping annotated genes and putative regulatory elements based on ATAC-seq data. The novel depletion method removed ∼40% of reads mapping to repeats, increasing those mapping to single-copy regions by ∼2.6-fold. When analyzing 25 million fragments, this repeat-to-single-copy shift in the sequencing data increased the number of genotyped bases of ∼10-fold compared to nondepleted libraries. In the same condition, we were also able to identify ∼12-fold more genetic variants in the single-copy regions and increased the genotyping accuracy by rescuing thousands of heterozygous variants that otherwise would be missed because of low coverage. The method performed similarly regardless of the multiplexing level, type of library or genotypes, including different cultivars and a closely related species (L. orientalis). Our results showed that CRISPR-Cas9-driven repeat depletion focuses sequencing data on single-copy regions, thus improving high-density and genome-wide genotyping in large and repetitive genomes.

First Glimpse on Spring Starflower Domestication

The cultivation and domestication of plants are human-driven processes that change the biology and attributes of a plant. Ipheion uniflorum is a bulbous geophyte known as Spring Starflower whose cultivation dates back to the first half of the 19th century. At least seven cultivars have been developed from natural stands. However, comparative analyses of wild and cultivated materials are largely missing. In the present study, we provide a morphological evaluation and analyses of the cytological and genetic variability of I. uniflorum that reveal significant levels of differentiation and evidence of artificial selection in the Spring Starflower. Distinctive phenotypic characters in cultivated materials that are rarely found or lacking in wild plants and natural populations, such as pink or violet flowers, together with its reduced heterozygosity and starting genetic differentiation support the view of early mechanisms of domestication acting upon Spring Starflower plants. The probable geographic origin of the cultivated forms is discussed together with perspectives for plant breeding.

Hopping or Jumping on the Cliffs: The Unusual Phylogeographical and Demographic Structure of an Extremely Narrow Endemic Mediterranean Plant

Several past and recent climatic and geological events have greatly influenced the current distribution of coastal species around the Mediterranean Basin. As a consequence, the reconstruction of the distributional history of these species is challenging. In this study, we used both chloroplast and nuclear SNPs to assess the levels of genetic differentiation, contemporary/historical levels of gene flow, and demographic history for the three only known (one mainland and two insular) populations of Eokochia saxicola, a rare Mediterranean coastal rocky halophyte. Plastid genome analysis revealed very low intraspecific haplotype variation and partial admixture among Capri and Palinuro populations with at least two independent colonization events for the Strombolicchio islet. Nuclear SNPs variation consistently identified three distinct genetic clusters corresponding to our sampling localities. Furthermore, strong genetic isolation was confirmed by both historical and contemporary levels of migration among the three populations. The DIYABC analysis identified two introductions temporally separated from Palinuro to Capri (ca.25 Mya) and subsequently to Strombolicchio (ca.09 Mya) as the most likely hypothesis for the current distribution of E. saxicola. Regardless of their small population sizes, all study sites supported high-genetic diversity maintained by outcrossing and random mating between individuals owing largely to wind pollination, an exclusive trait among Mediterranean narrow endemics. In conclusion, the patterns observed confirm that some Mediterranean endemics are not necessarily "evolutionary dead-ends" but rather represent species that have extensive demographic stability and a strong evolutionary legacy.

Genome-Wide Evidence for Complex Hybridization and Demographic History in a Group of Cycas From China

Cycads represent one of the most ancestral living seed plants as well as one of the most threatened plant groups in the world. South China is a major center and potential origin of Cycas, the most rapidly diversified lineage of cycads. However, genomic-wide diversity of Cycas remains poorly understood due to the challenge of generating genomic markers associated with their inherent large genomes. Here, we perform a comprehensive conservation genomic study based on restriction-site associated DNA sequencing (RADseq) data in six representative species of Cycas in South China. Consistently low genetic diversity and strong genetic differentiation were detected across species. Both phylogenetic inference and genetic structure analysis via several methods revealed generally congruent groups among the six Cycas species. The analysis with ADMIXTURE showed low mixing of genetic composition among species, while individuals of C. dolichophylla exhibited substantial genetic admixture with C. bifida, C. changjiangensis, and C. balansae. Furthermore, the results from Treemix, f₄-statistic, and ABBA-BABA test were generally consistent and revealed the complex patterns of interspecific gene flow. Relatively strong signals of hybridization were detected between C. dolichophylla and C. szechuanensis, and the ancestor of C. taiwaniana and C. changjiangensis. Distinct patterns of demographic history were inferred for these species by Stairway Plot, and our results suggested that both climate fluctuation and frequent geological activities during the late Pleistocene exerted deep impacts on the population dynamics of these species in South China. Finally, we explore the practical implications of our findings for the development of conservation strategies in Cycas. The present study demonstrates the efficiency of RADseq for conservation genomic studies on non-model species with large and complex genomes. Given the great significance of cycads as a radical transition in the evolution of plant biodiversity, our study provides important insights into the mechanisms of diversification in such recently radiated living fossil taxa.

Genomic, karyological and morphological changes of South American garlics (Ipheion) provide insights into mechanisms of speciation in the Pampean region

Speciation proceeds through mechanisms that promote reproductive isolation and shape the extent of genetic variation in natural populations, and thus its study is essential to understand the evolutionary processes leading to increased biodiversity. Chromosomal rearrangements are known to facilitate reproductive isolation by hybrid sterility and favour speciation events. The genus Ipheion (Amaryllidaceae, Allioideae) is unique as its species exhibit a remarkable karyological variability but lack population-level genetic data. To unveil the diversification processes acting upon the formation of new lineages within Ipheion in the Pampas of South America, we combined morphology and karyology approaches with genotyping-by-sequencing. Our phylogenomic and population genomics results supported the taxonomic division of Ipheion into three morphological and genetically well-differentiated groups. The origin of Ipheion uniflorum was traced back to its current southern distribution area in the southern Pampean region (in Argentina), from where it had expanded to the north reaching Uruguay. Our results further suggested that chromosome rearrangements and ploidy shifts had triggered speciation events, first during the origin of I. uniflorum and later during its subsequent diversification into I. recurvifolium and I. tweedieanum, in both cases reinforced by extrinsic factors and biogeographical settings. The current study illustrates the analytical power of multidisciplinary approaches integrating phylo- and population genomics with classic analyses to reveal evolutionary processes in plants.

Ecological Observations on Hybrid Populations of European Plethodontid Salamanders, Genus Speleomantes

Speleomantes are the only plethodontid salamanders present in Europe. Multiple studies have been performed to investigate the trophic niche of the eight Speleomantes species, but none of these studies included hybrid populations. For the first time, we studied the trophic niche of five Speleomantes hybrid populations. Each population was surveyed twice in 2020, and stomach flushing was performed on each captured salamander; stomach flushing is a harmless technique that allows stomach contents to be inspected. We also assessed the potential divergence in size and body condition between natural and introduced hybrids, and their parental species. Previously collected data on Speleomantes were included to increase the robustness of these analyses. In only 33 out of 134 sampled hybrid Speleomantes we recognized 81 items belonging to 11 prey categories. The frequency of empty stomachs was higher in females and individuals from natural hybrid populations, whereas the largest number of prey was consumed by males. We compared the total length and body condition of 685 adult salamanders belonging to three types of hybrids and three parental (sub)species. Three group of salamanders (one hybrid and two parental species) showed significantly larger size, whereas no difference in body condition was observed. This study provided novel ecological information on Speleomantes hybrid populations. We also provided insights into the potential divergence between hybrids and parental species in terms of size and body condition. We discuss our findings, and formulate several hypotheses that should be tested in the future.

Comparative Analysis of SNP Discovery and Genotyping in Fagus sylvatica L. and Quercus robur L. Using RADseq, GBS, and ddRAD Methods

Next-generation sequencing of reduced representation genomic libraries (RRL) is capable of providing large numbers of genetic markers for population genetic studies at relatively low costs. However, one major concern of these types of markers is the precision of genotyping, which is related to the common problem of missing data, which appears to be particularly important in association and genomic selection studies. We evaluated three RRL approaches (GBS, RADseq, ddRAD) and different SNP identification methods (de novo or based on a reference genome) to find the best solutions for future population genomics studies in two economically and ecologically important broadleaved tree species, namely F. sylvatica and Q. robur. We found that the use of ddRAD method coupled with SNP calling based on reference genomes provided the largest numbers of markers (28 k and 36 k for beech and oak, respectively), given standard filtering criteria. Using technical replicates of samples, we demonstrated that more than 80% of SNP loci should be considered as reliable markers in GBS and ddRAD, but not in RADseq data. According to the reference genomes’ annotations, more than 30% of the identified ddRAD loci appeared to be related to genes. Our findings provide a solid support for using ddRAD-based SNPs for future population genomics studies in beech and oak.