Scaling-up RADseq methods for large datasets of non-invasive samples: Lessons for library construction and data preprocessing.

Genetic non-invasive sampling (gNIS) is a critical tool for population genetics studies, supporting conservation efforts while imposing minimal impacts on wildlife. However, gNIS often presents variable levels of DNA degradation and non-endogenous contamination, which can incur considerable processing costs. Furthermore, the use of restriction-site-associated DNA sequencing methods (RADseq) for assessing thousands of genetic markers introduces the challenge of obtaining large sets of shared loci with similar coverage across multiple individuals. Here, we present an approach to handling large-scale gNIS-based datasets using data from the spotted hyena population inhabiting the Ngorongoro Crater in Tanzania. We generated 3RADseq data for more than a thousand individuals, mostly from faecal mucus samples collected non-invasively and varying in DNA degradation and contamination level. Using small-scale sequencing, we screened samples for endogenous DNA content, removed highly contaminated samples, confirmed overlap fragment length between libraries, and balanced individual representation in a sequencing pool. We evaluated the impact of (1) DNA degradation and contamination of non-invasive samples, (2) PCR duplicates and (3) different SNP filters on genotype accuracy based on Mendelian error estimated for parent-offspring trio datasets. Our results showed that when balanced for sequencing depth, contaminated samples presented similar genotype error rates to those of non-contaminated samples. We also showed that PCR duplicates and different SNP filters impact genotype accuracy. In summary, we showed the potential of using gNIS for large-scale genetic monitoring based on SNPs and demonstrated how to improve control over library preparation by using a weighted re-pooling strategy that considers the endogenous DNA content.

New Genetic Insights About Hybridization and Population Structure of Hawksbill and Loggerhead Turtles From Brazil.

An extremely high incidence of hybridization among sea turtles is found along the Brazilian coast. This atypical phenomenon and its impact on sea turtle conservation can be elucidated through research focused on the evolutionary history of sea turtles. We assessed high-quality multilocus haplotypes of 143 samples of the 5 species of sea turtles that occur along the Brazilian coast to investigate the hybridization process and the population structure of hawksbill (Eretmochelys imbricata) and loggerhead turtles (Caretta caretta). The multilocus data were initially used to characterize interspecific hybrids. Introgression (F2 hybrids) was only confirmed in hatchlings of F1 hybrid females (hawksbill × loggerhead), indicating that introgression was either previously overestimated and F2 hybrids may not survive to adulthood, or the first-generation hybrid females nesting in Brazil were born as recent as few decades ago. Phylogenetic analyses using nuclear markers recovered the mtDNA-based Indo-Pacific and Atlantic lineages for hawksbill turtles, demonstrating a deep genetic divergence dating from the early Pliocene. In addition, loggerhead turtles that share a common feeding area and belong to distinct Indo-Pacific and Atlantic mtDNA clades present no clear genetic differentiation at the nuclear level. Finally, our results indicate that hawksbill and loggerhead rookeries along the Brazilian coast are likely connected by male-mediated gene flow.

Global phylogeography of the critically endangered hawksbill turtle (Eretmochelys imbricata).

The hawksbill turtle is a broadly distributed, highly migratory and critically endangered sea turtle species. The paucity of studies restricts the comprehension of its behavior and life history. In this work, we performed a global phylogeographic analysis using a compilation of previously published mitochondrial haplotype data to understand the dynamics and diversity of hawksbill populations worldwide. Our results revealed a complex demographic pattern associated to hawksbill phylogeography since the Pliocene. Isolation by distance is not enough to explain distinct demographic units of hawksbill turtles, which are also influenced by other factors as oceanic currents, coral reef distribution and nesting timing. The foraging aggregations are typically mixed stocks of individuals originating from multiple nesting areas, but there is also a trend of foragers coming from nearby natal beaches. Phylogenetic analysis indicates two highly divergent major lineages split between Atlantic and Indo-Pacific rookeries, but there is also a more recent Atlantic Ocean colonization from the Indo-Pacific Ocean. Long-distance dispersal events are likely responsible for homogenization between distant populations within oceans. Our findings provided new insights about population connectivity, identified gaps that should be prioritized in future research and highlighted the need for international efforts aiming at hawksbill's conservation.