Genetic erosion in an endangered desert fish during a megadrought despite long-term supportive breeding.

Human water use combined with a recent megadrought have reduced river and stream flow through the southwest United States and led to periodic drying of formerly perennial river segments. Reductions in snowmelt runoff and increased extent of drying collectively threaten short-lived, obligate aquatic species, including the endangered Rio Grande silvery minnow (Hybognathus amarus). This species is subject to boom-and-bust population dynamics, under which large fluctuations in abundance are expected to lower estimates of effective population size and erode genetic diversity over time. Rates of diversity loss are also affected by additions of hatchery-origin fish used to supplement the wild population. We used demographic and genetic data from wild and hatchery individuals to examine the relationship of genetic diversity and effective population size to abundance over the last two decades. Genetic diversity was low during the early 2000s, but diversity and demographic metrics stabilized after the hatchery program was initiated and environmental conditions improved. Yet, from 2017 onward, allelic diversity declined (Cohen's d = 1.34) and remained low despite hatchery stocking and brief wild population recovery. Across the time series, single-sample estimates of effective population size based on linkage disequilibrium (LD Ne ) were positively associated (r = 0.53) with wild abundance and total abundance, but as the proportion of hatchery-origin spawners increased, LD Ne declined (r = -0.55). Megadrought limited wild spawner abundance and precluded refreshment of hatchery brood stocks with wild fish; hence, we predict a riverine population increasingly dominated by hatchery-origin individuals and accelerated loss of genetic diversity despite supplementation. We recommend an adaptive and accelerated management plan that integrates river flow management and hatchery operations to slow the pace of genetic diversity loss exacerbated by megadrought.

El uso humano del agua, combinado con una megasequía reciente, ha reducido el flujo de los ríos y arroyos en el suroeste de los Estados Unidos y ha provocado la seca periódica de segmentos de ríos que antes eran perennes. Las reducciones en la escorrentía del deshielo y el aumento de la sequía amenazan colectivamente a especies obligatoriamente acuáticas de vida corta, incluyendo la amenazada carpa chamizal (Hybognathus amarus). Esta especie está sujeta a una dinámica poblacional de explosión y colapso, bajo la cual se espera que grandes fluctuaciones en la abundancia reduzcan las estimaciones del tamaño efectivo de la población y erosionen la diversidad genética con el tiempo. Las tasas de pérdida de la diversidad también se ven afectadas por la adición de peces procedentes de criaderos usados para suplementar la población silvestre. Utilizamos datos demográficos y genéticos de individuos silvestres y de criaderos para examinar la relación entre la diversidad genética y el tamaño efectivo de la población con la abundancia durante las últimas dos décadas. La diversidad genética fue baja a principios de los 2000, pero las métricas de diversidad y demografía estabilizaron después de que se inició el programa de criadero y mejoraron las condiciones ambientales. Sin embargo, a partir de 2017, la diversidad alélica disminuyó (d de Cohen = 1,34) y se mantuvo baja a pesar de la suplementación con individuos de criaderos y la breve recuperación de la población silvestre. A lo largo del tiempo, las estimativas de muestras individuales del tamaño efectivo de la población basados en el desequilibrio de ligamiento (LD Ne) estaban asociadas positivamente (r = 0,53) con la abundancia silvestre y la abundancia total, pero a medida que la proporción de desovadores originados en criaderos aumentó, el LD Ne disminuyó (r = -0,55). La megasequía limitó la abundancia de desovadores silvestres e impidió el reabastecimiento de las poblaciones en cautiverio con peces silvestres; por lo tanto, predecimos una población ribereña cada vez más dominada por individuos procedentes de criaderos y una pérdida acelerada de diversidad genética a pesar de la suplementación. Recomendamos un plan de gestión adaptativo y acelerado que integre la gestión del caudal del río y las operaciones de criaderos para frenar el ritmo de pérdida de diversidad genética exacerbada por la megasequía.

The complete mitochondrial genomes of two imperiled species endemic to the Southwestern United States: Peppered Chub (Macrhybopsis tetranema) and Gila Trout (Oncorhynchus gilae).

Macrhybopsis tetranema and Oncorhynchus gilae are fish species endemic to the Southwestern United States. We present the complete mitochondrial genomes for these species. Each genome consisted of 13 protein-coding genes, two ribosomal (rRNA) genes, 22 transfer RNA (tRNA) genes, and the control region (D-loop). Mitogenome lengths were 16,916 base pairs (bp) for M. tetranema, and 16,976 bp for O. gilae. The GC content was 41% for M. tetranema and 46% for O. gilae. The relationships of M. tetranema and O. gilae were consistent with previous phylogenetic analyses.

Hybridization in late stages of speciation: Strong but incomplete genome-wide reproductive isolation and 'large Z-effect' in a moving hybrid zone.

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

Limited evidence for extensive genetic differentiation between X and Y chromosomes in Hybognathus amarus (Cypriniformes: Leuciscidae).

Sex determination systems and genetic sex differentiation across fishes are highly diverse but are unknown for most Cypriniformes, including Rio Grande silvery minnow (Hybognathus amarus). In this study, we aimed to detect and validate sex-linked markers to infer sex determination system and to demonstrate the utility of combining several methods for sex-linked marker detection in nonmodel organisms. To identify potential sex-linked markers, Nextera-tagmented reductively amplified DNA (nextRAD) libraries were generated from 66 females, 64 males, and 60 larvae of unknown sex. These data were combined with female and male de novo genomes from Nanopore long-read sequences. We identified five potential unique male nextRAD-tags and one potential unique male contig, suggesting an XY sex determination system. We also identified two single-nucleotide polymorphisms (SNPs) in the same contig with values of FST, allele frequencies, and heterozygosity conforming with expectations of an XY system. Through PCR we validated the marker containing the sex-linked SNPs and a single nextRAD-tag sex-associated marker but it was not male specific. Instead, more copies of this locus in the male genome were suggested by enhanced amplification in males. Results are consistent with an XY system with low differentiation between sex-determining regions. Further research is needed to confirm the level of differentiation between the sex chromosomes. Nonetheless, this study highlighted the power of combining reduced representation and whole-genome sequencing for identifying sex-linked markers, especially when reduced representation sequencing does not include extensive variation between sexes, either because such variation is not present or not captured.

Widespread introgression of MHC genes in Iberian Podarcis lizards.

Major histocompatibility complex (MHC) genes are crucial for the adaptive immune response of jawed vertebrates. Their variation, reaching extreme levels, is driven mainly by an arms race between hosts and pathogens. One hypothesised mechanism contributing to MHC polymorphism is adaptive introgression, the exchange of genetic variants between hybridising species favoured by selection, yet its effect on MHC variation is poorly understood. Detection of adaptive MHC introgression, though challenging, may be facilitated by the analysis of species complexes forming multiple hybrid zones. Here, we investigated MHC introgression in six hybrid zones formed by seven species of Podarcis lizards inhabiting the Iberian Peninsula. To differentiate adaptive introgression from neutral introgression, we compared the patterns of gene exchange in MHC and genome-wide markers. We found elevated sharing of MHC alleles in the proximity of contact beyond the areas of detectable genome-wide admixture in most hybrid zones and, in half of them, asymmetric MHC exchange. In general, the elevated MHC allele sharing between species pairs with abutting ranges compared to geographically isolated species pairs also supports the prevalence of introgression. Collectively, our results demonstrate widespread MHC introgression in the Iberian Podarcis complex and suggest its adaptiveness. Contrary to previous results from Triturus newts, we did not observe differences in the rate of introgression between MHC classes. Our work adds support to the emerging view of adaptive introgression as a key mechanism shaping MHC diversity. It also raises questions about the effect of elevated MHC variation and factors leading to the asymmetry of adaptive introgression.

Transitioning from microsatellites to SNP-based microhaplotypes in genetic monitoring programmes: Lessons from paired data spanning 20 years.

Many long-term genetic monitoring programmes began before next-generation sequencing became widely available. Older programmes can now transition to new marker systems usually consisting of 1000s of SNP loci, but there are still important questions about comparability, precision, and accuracy of key metrics estimated using SNPs. Ideally, transitioned programmes should capitalize on new information without sacrificing continuity of inference across the time series. We combined existing microsatellite-based genetic monitoring information with SNP-based microhaplotypes obtained from archived samples of Rio Grande silvery minnow (Hybognathus amarus) across a 20-year time series to evaluate point estimates and trajectories of key genetic metrics. Demographic and genetic monitoring bracketed multiple collapses of the wild population and included cases where captive-born repatriates comprised the majority of spawners in the wild. Even with smaller sample sizes, microhaplotypes yielded comparable and in some cases more precise estimates of variance genetic effective population size, multilocus heterozygosity and inbreeding compared to microsatellites because many more microhaplotype loci were available. Microhaplotypes also recorded shifts in allele frequencies associated with population bottlenecks. Trends in microhaplotype-based inbreeding metrics were associated with the fraction of hatchery-reared repatriates to the wild and should be incorporated into future genomic monitoring. Although differences in accuracy and precision of some metrics were observed between marker types, biological inferences and management recommendations were consistent.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae).

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species.

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.