Comparative phylogenomics reveal complex evolution of life history strategies in a clade of bivalves with parasitic larvae (Bivalvia: Unionoida: Ambleminae)

Integrative phylogenetic, phylogeographic and morphological characterisation of the Unio crassus species complex reveals cryptic diversity with important conservation implications

The global decline of freshwater mussels and their crucial ecological services highlight the need to understand their phylogeny, phylogeography and patterns of genetic diversity to guide conservation efforts. Such knowledge is urgently needed for Unio crassus, a highly imperilled species originally widespread throughout Europe and southwest Asia. Recent studies have resurrected several species from synonymy based on mitochondrial data, revealing U. crassus to be a complex of cryptic species. To address long-standing taxonomic uncertainties hindering effective conservation, we integrate morphometric, phylogenetic, and phylogeographic analyses to examine species diversity within the U. crassus complex across its entire range. Phylogenetic analyses were performed using cytochrome c oxidase subunit I (815 specimens from 182 populations) and, for selected specimens, whole mitogenome sequences and Anchored Hybrid Enrichment (AHE) data on ∼ 600 nuclear loci. Mito-nuclear discordance was detected, consistent with mitochondrial DNA gene flow between some species during the Pliocene and Pleistocene. Fossil-calibrated phylogenies based on AHE data support a Mediterranean origin for the U. crassus complex in the Early Miocene. The results of our integrative approach support 12 species in the group: the previously recognised Unio bruguierianus, Unio carneus, Unio crassus, Unio damascensis, Unio ionicus, Unio sesirmensis, and Unio tumidiformis, and the reinstatement of five nominal taxa: Unio desectusstat. rev., Unio gontieriistat. rev., Unio mardinensisstat. rev., Unio nanusstat. rev., and Unio vicariusstat. rev. Morphometric analyses of shell contours reveal important morphospace overlaps among these species, highlighting cryptic, but geographically structured, diversity. The distribution, taxonomy, phylogeography, and conservation of each species are succinctly described.

Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America's Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia-a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species (Erimystax dissimilis) confirmed to serve as a parasitic larval host for P. sintoxia. Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

Community-wide correlations between species richness, abundance and population genomic diversity in a freshwater biodiversity hotspot

Understanding patterns of diversity across macro (e.g. species-level) and micro (e.g. molecular-level) scales can shed light on community function and stability by elucidating the abiotic and biotic drivers of diversity within ecological communities. We examined the relationships among taxonomic and genetic metrics of diversity in freshwater mussels (Bivalvia: Unionidae), an ecologically important and species-rich group in the southeastern United States. Using quantitative community surveys and reduced-representation genome sequencing across 22 sites in seven rivers and two river basins, we surveyed 68 mussel species and sequenced 23 of these species to characterize intrapopulation genetic variation. We tested for the presence of species diversity-abundance correlations (i.e. the more-individuals hypothesis, MIH), species-genetic diversity correlations (SGDCs) and abundance-genetic diversity correlations (AGDCs) across all sites to evaluate relationships between different metrics of diversity. Sites with greater cumulative multispecies density (a standardized metric of abundance) had a greater number of species, consistent with the MIH hypothesis. Intrapopulation genetic diversity was strongly associated with the density of most species, indicating the presence of AGDCs. However, there was no consistent evidence for SGDCs. Although sites with greater overall densities of mussels had greater species richness, sites with higher genetic diversity did not always exhibit positive correlations with species richness, suggesting that there are spatial and evolutionary scales at which the processes influencing community-level diversity and intraspecific diversity differ. Our work reveals the importance of local abundance as indicator (and possibly a driver) of intrapopulation genetic diversity.

A lizard is never late: Squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012 and 2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

A lizard is never late: squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012-2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

Oriental freshwater mussels arose in East Gondwana and arrived to Asia on the Indian Plate and Burma Terrane

Freshwater mussels cannot spread through oceanic barriers and represent a suitable model to test the continental drift patterns. Here, we reconstruct the diversification of Oriental freshwater mussels (Unionidae) and revise their taxonomy. We show that the Indian Subcontinent harbors a rather taxonomically poor fauna, containing 25 freshwater mussel species from one subfamily (Parreysiinae). This subfamily most likely originated in East Gondwana in the Jurassic and its representatives arrived to Asia on two Gondwanan fragments (Indian Plate and Burma Terrane). We propose that the Burma Terrane was connected with the Indian Plate through the Greater India up to the terminal Cretaceous. Later on, during the entire Paleogene epoch, these blocks have served as isolated evolutionary hotspots for freshwater mussels. The Burma Terrane collided with mainland Asia in the Late Eocene, leading to the origin of the Mekong's Indochinellini radiation. Our findings indicate that the Burma Terrane had played a major role as a Gondwanan "biotic ferry" alongside with the Indian Plate.

Evolution of diverse host infection mechanisms delineates an adaptive radiation of lampsiline freshwater mussels centered on their larval ecology

North American watersheds contain a high diversity of freshwater mussels (Unionoida). During the long-lived, benthic phase of their life cycle, up to 40 species can co-occur in a single riffle and there is typically little evidence for major differences in their feeding ecology or microhabitat partitioning. In contrast, their brief parasitic larval phase involves the infection of a wide diversity of fish hosts and female mussels have evolved a spectrum of adaptations for infecting host fish with their offspring. Many species use a passive broadcast strategy: placing high numbers of larvae in the water column and relying on chance encounters with potential hosts. Many other species, including most members of the Lampsilini, have a proactive strategy that entails the use of prey-mimetic lures to change the behavior of the hosts, i.e., eliciting a feeding response through which they become infected. Two main lure types are collectively produced: mantle tissue lures (on the female’s body) and brood lures, containing infective larvae, that are released into the external environment. In this study, we used a phylogenomic approach (ddRAD-seq) to place the diversity of infection strategies used by 54 North American lampsiline mussels into an evolutionary context. Ancestral state reconstruction recovered evidence for the early evolution of mantle lures in this clade, with brood lures and broadcast infection strategies both being independently derived twice. The most common infection strategy, occurring in our largest ingroup clade, is a mixed one in which mimetic mantle lures are apparently the predominant infection mechanism, but gravid females also release simple, non-mimetic brood lures at the end of the season. This mixed infection strategy clade shows some evidence of an increase in diversification rate and most members use centrarchids (Micropterus & Lepomis spp.) as their predominant fish hosts. Broad linkage between infection strategies and predominant fish host genera is also seen in other lampsiline clades: worm-like mantle lures of Toxolasma spp. with sunfish (Lepomis spp.); insect larvae-like brood lures (Ptychobranchus spp.), or mantle lures (Medionidus spp., Obovaria spp.), or mantle lures combined with host capture (Epioblasma spp.) with a spectrum of darter (Etheostoma & Percina spp.) and sculpin (Cottus spp.) hosts, and tethered brood lures (Hamiota spp.) with bass (Micropterus spp.). Our phylogenetic results confirm that discrete lampsiline mussel clades exhibit considerable specialization in the primary fish host clades their larvae parasitize, and in the host infection strategies they employ to do so. They are also consistent with the hypothesis that larval resource partitioning of fish hosts is an important factor in maintaining species diversity in mussel assemblages. We conclude that, taking their larval ecology and host-infection mechanisms into account, lampsiline mussels may be legitimately viewed as an adaptive radiation.

Establishing conservation units to promote recovery of two threatened freshwater mussel species (Bivalvia: Unionida: Potamilus)

Population genomics has significantly increased our ability to make inferences about microevolutionary processes and demographic histories, which have the potential to improve protection and recovery of imperiled species. Freshwater mussels (Bivalvia: Unionida) represent one of the most imperiled groups of organisms globally. Despite systemic decline of mussel abundance and diversity, studies evaluating spatiotemporal changes in distribution, demographic histories, and ecological factors that threaten long-term persistence of imperiled species remain lacking. In this study, we use genotype-by-sequencing (GBS) and mitochondrial sequence data (mtDNA) to define conservation units (CUs) for two highly imperiled freshwater mussel species, Potamilus amphichaenus and Potamilus streckersoni. We then synthesize our molecular findings with details from field collections spanning from 1901 to 2019 to further elucidate distributional trends, contemporary status, and other factors that may be contributing to population declines for our focal species. We collected GBS and mtDNA data for individuals of P. amphichaenus and P. streckersoni from freshwater mussel collections in the Brazos, Neches, Sabine, and Trinity drainages ranging from 2012 to 2019. Molecular analyses resolved disputing number of genetic clusters within P. amphichaenus and P. streckersoni; however, we find defensible support for four CUs, each corresponding to an independent river basin. Evaluations of historical and recent occurrence data illuminated a generally increasing trend of occurrence in each of the four CUs, which were correlated with recent increases in sampling effort. Taken together, these findings suggest that P. amphichaenus and P. streckersoni are likely rare throughout their respective ranges. Because of this, the establishment of CUs will facilitate evidence-based recovery planning and ensure potential captive propagation and translocation efforts are beneficial. Our synthesis represents a case study for conservation genomic assessments in freshwater mussels and provides a model for future studies aimed at recovery planning for these highly imperiled organisms.

A High-Quality Reference Genome for a Parasitic Bivalve with Doubly Uniparental Inheritance (Bivalvia: Unionida)

From a genomics perspective, bivalves (Mollusca: Bivalvia) have been poorly explored with the exception for those of high economic value. The bivalve order Unionida, or freshwater mussels, has been of interest in recent genomic studies due to their unique mitochondrial biology and peculiar life cycle. However, genomic studies have been hindered by the lack of a high-quality reference genome. Here, I present a genome assembly of Potamilus streckersoni using Pacific Bioscience single-molecule real-time long reads and 10X Genomics-linked read sequencing. Further, I use RNA sequencing from multiple tissue types and life stages to annotate the reference genome. The final assembly was far superior to any previously published freshwater mussel genome and was represented by 2,368 scaffolds (2,472 contigs) and 1,776,755,624 bp, with a scaffold N50 of 2,051,244 bp. A high proportion of the assembly was comprised of repetitive elements (51.03%), aligning with genomic characteristics of other bivalves. The functional annotation returned 52,407 gene models (41,065 protein, 11,342 tRNAs), which was concordant with the estimated number of genes in other freshwater mussel species. This genetic resource, along with future studies developing high-quality genome assemblies and annotations, will be integral toward unraveling the genomic bases of ecologically and evolutionarily important traits in this hyper-diverse group.

Gene family amplification facilitates adaptation in freshwater unionid bivalve Megalonaias nervosa

Freshwater unionid bivalves currently face severe anthropogenic challenges. Over 70% of species in the United States are threatened, endangered or extinct due to pollution, damming of waterways and overfishing. These species are notable for their unusual life history strategy, parasite-host co-evolution and biparental mitochondrial inheritance. Among this clade, the washboard mussel Megalonaias nervosa is one species that remains prevalent across the Southeastern United States, with robust population sizes. We have created a reference genome for M. nervosa to determine how genome content has evolved in the face of these widespread environmental challenges. We observe dynamic changes in genome content, with a burst of recent transposable element proliferation causing a 382 Mb expansion in genome content. Birth-death models suggest rapid expansions among gene families, with a mutation rate of 1.16 × 10^-8 duplications per gene per generation. Cytochrome P450 gene families have experienced exceptional recent amplification beyond expectations based on genome-wide birth-death processes. These genes are associated with increased rates of amino acid changes, a signature of selection driving evolution of detox genes. Fitting evolutionary models of adaptation from standing genetic variation, we can compare adaptive potential across species and mutation types. The large population size in M. nervosa suggests a 4.7-fold advantage in the ability to adapt from standing genetic variation compared with a low diversity endemic E. hopetonensis. Estimates suggest that gene family evolution may offer an exceptional substrate of genetic variation in M. nervosa, with P_sgv  = 0.185 compared with P_sgv  = 0.067 for single nucleotide changes. Hence, we suggest that gene family evolution is a source of 'hopeful monsters' within the genome that may facilitate adaptation when selective pressures shift. These results suggest that gene family expansion is a key driver of adaptive evolution in this key species of freshwater Unionidae that is currently facing widespread environmental challenges. This work has clear implications for conservation genomics on freshwater bivalves as well as evolutionary theory. This genome represents a first step to facilitate reverse ecological genomics in Unionidae and identify the genetic underpinnings of phenotypic diversity.

Novel Molecular Resources to Facilitate Future Genetics Research on Freshwater Mussels (Bivalvia: Unionidae)

Molecular data have been an integral tool in the resolution of the evolutionary relationships and systematics of freshwater mussels, despite the limited number of nuclear markers available for Sanger sequencing. To facilitate future studies, we evaluated the phylogenetic informativeness of loci from the recently published anchored hybrid enrichment (AHE) probe set Unioverse and developed novel Sanger primer sets to amplify two protein-coding nuclear loci with high net phylogenetic informativeness scores: fem-1 homolog C (FEM1) and UbiA prenyltransferase domain-containing protein 1 (UbiA). We report the methods used for marker development, along with the primer sequences and optimized PCR and thermal cycling conditions. To demonstrate the utility of these markers, we provide haplotype networks, DNA alignments, and summary statistics regarding the sequence variation for the two protein-coding nuclear loci (FEM1 and UbiA). Additionally, we compare the DNA sequence variation of FEM1 and UbiA to three loci commonly used in freshwater mussel genetic studies: the mitochondrial genes cytochrome c oxidase subunit 1 (CO1) and NADH dehydrogenase subunit 1 (ND1), and the nuclear internal transcribed spacer 1 (ITS1). All five loci distinguish among the three focal species (Potamilus fragilis, Potamilus inflatus, and Potamilus purpuratus), and the sequence variation was highest for ND1, followed by CO1, ITS1, UbiA, and FEM1, respectively. The newly developed Sanger PCR primers and methodologies for extracting additional loci from AHE probe sets have great potential to facilitate molecular investigations targeting supraspecific relationships in freshwater mussels, but may be of limited utility at shallow taxonomic scales.