Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments.

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family encodes phospholipases essential for the degradation of organelles in the lens of the eye. These phospholipases translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny and selection of plaat1 across bony fishes and tetrapods. We showed that plaat1 (probably ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual-acuity and blind mammals and fishes. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

Impacts of Quaternary climatic changes on the diversification of riverine cichlids in the lower Congo River.

Climatic and geomorphological changes during the Quaternary period impacted global patterns of speciation and diversification across a wide range of taxa, but few studies have examined these effects on African riverine fishes. The lower Congo River is an excellent natural laboratory for understanding complex speciation and population diversification processes as it is hydrologically extremely dynamic and recognized as a continental hotspot of diversity harboring many narrowly endemic species. A previous study using genome-wide SNP data highlighted the importance of dynamic hydrological regimes to the diversification and speciation in lower Congo River cichlids. However, historical climate and hydrological changes (e.g., reduced river discharge during extended dry periods) have likely also influenced ichthyofaunal diversification processes in this system. The lower Congo River offers a unique opportunity to study climate-driven changes in river discharge, given the massive volume of water from the entire Congo basin flowing through this short stretch of the river. Here, we, for the first time, investigate the impacts of paleoclimatic factors on ichthyofaunal diversification in this system by inferring divergence times and modeling patterns of gene flow in four endemic lamprologine cichlids, including the blind cichlid, Lamprologus lethops. Our results suggest that Quaternary climate changes associated with river discharge fluctuations may have impacted the diversification of species along the system. Our study, using reduced representation sequencing (2RADseq), indicates that the lower Congo River lamprologines emerged during the Early-Middle Pleistocene transition, characterized as one of the earth's major climatic transformation periods. Modeling results suggest that gene flow across populations and between species was not constant but occurred in temporally constrained pulses. We show that these results correlate with glacial-interglacial fluctuations. The current hyper-diverse fish assemblages of the lower Congo River riverscape likely reflect the synergistic effects of multiple drivers fueling complex evolutionary processes through time.

Estuarine Sediment Microbiomes from a Chronosequence of Restored Urban Salt Marshes.

Salt marshes play an important role in the global nutrient cycle. The sediments in these systems harbor diverse and complex bacterial communities possessing metabolic capacities that provide ecosystem services such as nutrient cycling and removal. On the East Coast of the USA, salt marshes have been experiencing degradation due to anthropogenic stressors. Salt marsh islands within Jamaica Bay, New York City (USA), are surrounded by a large highly urbanized watershed and have declined in area. Restoration efforts have been enacted to reduce further loss, but little is known about how microbial communities develop following restoration activities, or how processes such as nitrogen cycling are impacted. Sediment samples were collected at two sampling depths from five salt marsh islands to characterize the bacterial communities found in marsh sediment including a post-restoration chronosequence of 3-12 years. We used 16s rRNA amplicon sequencing to define alpha and beta diversity, taxonomic composition, and predicted metabolic profile of each sediment sample. We found significant differences in alpha diversity between sampling depths, and significant differences in beta diversity, taxonomic composition, and predicted metabolic capacity among the five sampling locations. The youngest restored site and the degraded natural sampling site exhibited the most distinct communities among the five sites. Our findings suggest that while the salt marsh islands are located in close proximity to each other, they harbor distinct bacterial communities that can be correlated with post-restoration age, marsh health, and other environmental factors such as availability of organic carbon. IMPORTANCE: Salt marshes play a critical role in the global nutrient cycle due to sediment bacteria and their metabolic capacities. Many East Coast salt marshes have experienced significant degradation over recent decades, thought largely to be due to anthropogenic stressors such as nitrogen loading, urban development, and sea-level rise. Salt marsh islands in Jamaica Bay (Queens/Brooklyn NY) are exposed to high water column nitrogen due to wastewater effluent. Several receding marsh islands have been subjected to restoration efforts to mitigate this loss. Little is known about the effect marsh restoration has on bacterial communities, their metabolic capacity, or how they develop post-restoration. Here, we describe the bacterial communities found in marsh islands including a post-restoration chronosequence of 3-12 years and one degraded marsh island that remains unrestored. We found distinct communities at marsh sites, despite their geographic proximity. Differences in diversity and community composition were consistent with changes in organic carbon availability that occur during marsh development, and may result in differences in ecosystem function among sites.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments.

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family are phospholipases essential for the degradation of organelles in the lens of the eye. They translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny, and selection of plaat1 across bony fishes and tetrapods. We show that plaat1 (likely ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual acuity and blind mammals and fish. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light-environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

Riverscape genomics of cichlid fishes in the lower Congo: Uncovering mechanisms of diversification in an extreme hydrological regime.

Freshwater fishes are notably diverse, given that freshwater habitat represents a tiny fraction of the earth's surface, but the mechanisms generating this diversity remain poorly understood. Rivers provide excellent models to understand how freshwater diversity is generated and maintained across heterogeneous habitats. In particular, the lower Congo River (LCR) consists of a dynamic hydroscape exhibiting extraordinary aquatic biodiversity, endemicity, morphological and ecological specialization. Previous studies have suggested that the numerous high-energy rapids throughout the LCR form physical barriers to gene flow, thus facilitating diversification and speciation, generating ichthyofaunal diversity. However, this hypothesis has not been fully explored using genome-wide SNPs for fish species distributed across the LCR. Here, we examined four lamprologine cichlids endemic to the LCR that are distributed along the river without range overlap. Using genome-wide SNP data, we tested the hypotheses that high-energy rapids serve as physical barriers to gene flow that generate genetic divergence at interspecific and intraspecific levels, and that gene flow occurs primarily in a downstream direction. Our results are consistent with the prediction that powerful rapids sometimes act as a barrier to gene flow but also suggest that, at certain temporal and spatial scales, they may provide multidirectional dispersal opportunities for riverine rheophilic cichlid fishes. These results highlight the complexity of diversification processes in rivers and the importance of assessing such processes across different riverscapes.

Trade-offs between reducing complex terminology and producing accurate interpretations from environmental DNA: Comment on "Environmental DNA: What's behind the term?" by Pawlowski et al., (2020).

In a recent paper, "Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring," Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.

Genomic Analysis of the Only Blind Cichlid Reveals Extensive Inactivation in Eye and Pigment Formation Genes.

Trait loss represents an intriguing evolutionary problem, particularly when it occurs across independent lineages. Fishes in light-poor environments often evolve "troglomorphic" traits, including reduction or loss of both pigment and eyes. Here, we investigate the genomic basis of trait loss in a blind and depigmented African cichlid, Lamprologus lethops, and explore evolutionary forces (selection and drift) that may have contributed to these losses. This species, the only known blind cichlid, is endemic to the lower Congo River. Available evidence suggests that it inhabits deep, low-light habitats. Using genome sequencing, we show that genes related to eye formation and pigmentation, as well as other traits associated with troglomorphism, accumulated inactivating mutations rapidly after speciation. A number of the genes affected in L. lethops are also implicated in troglomorphic phenotypes in Mexican cavefish (Astyanax mexicanus) and other species. Analysis of heterozygosity patterns across the genome indicates that L. lethops underwent a significant population bottleneck roughly 1 Ma, after which effective population sizes remained low. Branch-length tests on a subset of genes with inactivating mutations show little evidence of directional selection; however, low overall heterozygosity may reduce statistical power to detect such signals. Overall, genome-wide patterns suggest that accelerated genetic drift from a severe bottleneck, perhaps aided by directional selection for the loss of physiologically expensive traits, caused inactivating mutations to fix rapidly in this species.

First genome-wide analysis of the endangered, endemic lichen Cetradonia linearis reveals isolation by distance and strong population structure.

PREMISE OF THE STUDY: Lichenized fungi are evolutionarily diverse and ecologically important, but little is known about the processes that drive their diversification and genetic differentiation. Distributions are often assumed to be wholly shaped by ecological requirements rather than dispersal limitations. Furthermore, although asexual and sexual reproductive structures are observable, the lack of information about recombination rates makes inferences about reproductive strategies difficult. We investigated the population genomics of Cetradonia linearis, a federally endangered lichen in the southern Appalachians of eastern North America, to test the relative contributions of environmental and geographic distance in shaping genetic structure, and to characterize the mating system and genome-wide recombination. METHODS: Whole-genome shotgun sequencing was conducted to generate data for 32 individuals of C. linearis. A reference genome was assembled, and reads from all samples were aligned to generate a set of single-nucleotide polymorphisms for further analyses. KEY RESULTS: We found evidence for low rates of recombination and for isolation by distance, but not for isolation by environment. The species is putatively unisexual, given that only one mating-type locus was found. Hindcast species distribution models and the distribution of genetic diversity support C. linearis having a larger range during the Last Glacial Maximum in the southern portion of its current extent. CONCLUSIONS: Our findings contribute to the understanding of factors that shape genetic diversity in C. linearis and in fungi more broadly. Because all populations are highly genetically differentiated, the extirpation of any population would mean the loss of unique genetic diversity; therefore, our results support the continued conservation of this species.

Fishes in the desert: mitochondrial variation and phylogeography of Danakilia (Actinopterygii: Cichlidae) and Aphanius (Actinopterygii: Cyprinodontidae) in the Danakil Depression of northeastern Africa.

The Danakil Depression in northeastern Africa represents one of the harshest arid environments on Earth, yet two genera of fishes, Danakilia (Cichlidae) and Aphanius (Cyprinodontidae), share its sparse aquatic habitats. The evolutionary history of these fishes is investigated here in the context of genetic, geological and paleoenvironmental information. We collected samples from seven sites and assessed phylogeographic relationships using concatenated COI and cytb mtDNA genes. Danakilia morphospecies show low differentiation at mitochondrial markers, but variation is partitioned between a northern cluster containing D. dinicolai plus three undescribed riverine populations, and a southern cluster including two creek populations of D. franchettii separated by the hypersaline waters of Lake Afrera. Aphanius displayed four genetically distinct clades (A. stiassnyae in Lake Afrera; one distributed across the entire area; one in Lake Abaeded; and one in the Shukoray River), but without clear large-scale geographic structure. However, Danakil Aphanius are clearly differentiated from A. dispar sensu stricto from the Sinai Peninsula. Geological evidence suggests that after the Late Pleistocene closure of the Danakil-Red Sea connection, increased post-glacial groundwater availability caused the formation of a brackish paleo-lake flooding the entire region below the -50 m contour. Fish populations previously isolated in coastal oases during glaciation were able to mix in the paleo-lake. Subsequently, in a more arid phase starting ∼7300 BP, paleo-lake regression isolated fishes in separate drainages, triggering their still ongoing diversification.

Genomewide SNP data reveal cryptic phylogeographic structure and microallopatric divergence in a rapids-adapted clade of cichlids from the Congo River.

The lower Congo River is a freshwater biodiversity hot spot in Africa characterized by some of the world's largest rapids. However, little is known about the evolutionary forces shaping this diversity, which include numerous endemic fishes. We investigated phylogeographic relationships in Teleogramma, a small clade of rheophilic cichlids, in the context of regional geography and hydrology. Previous studies have been unable to resolve phylogenetic relationships within Teleogramma due to lack of variation in nuclear genes and discrete morphological characters among putative species. To sample more broadly across the genome, we analysed double-digest restriction-associated sequencing (ddRAD) data from 53 individuals across all described species in the genus. We also assessed body shape and mitochondrial variation within and between taxa. Phylogenetic analyses reveal previously unrecognized lineages and instances of microallopatric divergence across as little as ~1.5 km. Species ranges appear to correspond to geographic regions broadly separated by major hydrological and topographic barriers, indicating these features are likely important drivers of diversification. Mitonuclear discordance indicates one or more introgressive hybridization events, but no clear evidence of admixture is present in nuclear genomes, suggesting these events were likely ancient. A survey of female fin patterns hints that previously undetected lineage-specific patterning may be acting to reinforce species cohesion. These analyses highlight the importance of hydrological complexity in generating diversity in certain freshwater systems, as well as the utility of ddRAD-Seq data in understanding diversification processes operating both below and above the species level.

Melting barriers to faunal exchange across ocean basins.

Accelerated loss of sea ice in the Arctic is opening routes connecting the Atlantic and Pacific Oceans for longer periods each year. These changes may increase the ease and frequency with which marine birds and mammals move between the Pacific and Atlantic Ocean basins. Indeed, recent observations of birds and mammals suggest these movements have intensified in recent decades. Reconnection of the Pacific and Atlantic Ocean basins will present both challenges to marine ecosystem conservation and an unprecedented opportunity to examine the ecological and evolutionary consequences of interoceanic faunal exchange in real time. To understand these changes and implement effective conservation of marine ecosystems, we need to further develop modeling efforts to predict the rate of dispersal and consequences of faunal exchange. These predictions can be tested by closely monitoring wildlife dispersal through the Arctic Ocean and using modern methods to explore the ecological and evolutionary consequences of these movements.

Molecular phylogenetics reveals convergent evolution in lower Congo River spiny eels.

BACKGROUND: The lower Congo River (LCR) is a region of exceptional species diversity and endemism in the Congo basin, including numerous species of spiny eels (genus Mastacembelus). Four of these exhibit distinctive phenotypes characterized by greatly reduced optic globes deeply embedded into the head (cryptophthalmia) and reduced (or absent) melanin pigmentation, among other characteristics. A strikingly similar cryptophthalmic phenotype is also found in members of a number of unrelated fish families, strongly suggesting the possibility of convergent evolution. However, little is known about the evolutionary processes that shaped diversification in LCR Mastacembelus, their biogeographic origins, or when colonization of the LCR occurred. METHODS: We sequenced mitochondrial and nuclear genes from Mastacembelus species collected in the lower Congo River, and compared them with other African species and Asian representatives as outgroups. We analyzed the sequence data using Maximum Likelihood and Bayesian phylogenetic inference. RESULTS: Bayesian and Maximum Likelihood phylogenetic analyses, and Bayesian coalescent methods for species tree reconstruction, reveal that endemic LCR spiny eels derive from two independent origins, clearly demonstrating convergent evolution of the cryptophthalmic phenotype. Mastacembelus crassus, M. aviceps, and M. simbi form a clade, allied to species found in southern, eastern and central Africa. Unexpectedly, M. brichardi and brachyrhinus fall within a clade otherwise endemic to Lake Tanganikya (LT) ca. 1500 km east of the LCR. Divergence dating suggests the ages of these two clades of LCR endemics differ markedly. The age of the crassus group is estimated at ~4 Myr while colonization of the LCR by the brichardi-brachyrhinus progenitor was considerably more recent, dated at ~0.5 Myr. CONCLUSIONS: The phylogenetic framework of spiny eels presented here, the first to include LCR species, demonstrates that cryptophthalmia and associated traits evolved at least twice in Mastacembelus: once in M. brichardi and at least once in the M. crassus clade. Timing of diversification is broadly consistent with the onset of modern high-energy flow conditions in the LCR and with previous studies of endemic cichlids. The close genetic relationship between M. brichardi and M. brachyrhinus is particularly notable given the extreme difference in phenotype between these species, and additional work is needed to better understand the evolutionary history of diversification in this clade. The findings presented here demonstrate strong, multi-trait convergence in LCR spiny eels, suggesting that extreme selective pressures have shaped numerous phenotypic attributes of the endemic species of this region.

Climate impacts on transocean dispersal and habitat in gray whales from the Pleistocene to 2100.

Arctic animals face dramatic habitat alteration due to ongoing climate change. Understanding how such species have responded to past glacial cycles can help us forecast their response to today's changing climate. Gray whales are among those marine species likely to be strongly affected by Arctic climate change, but a thorough analysis of past climate impacts on this species has been complicated by lack of information about an extinct population in the Atlantic. While little is known about the history of Atlantic gray whales or their relationship to the extant Pacific population, the extirpation of the Atlantic population during historical times has been attributed to whaling. We used a combination of ancient and modern DNA, radiocarbon dating and predictive habitat modelling to better understand the distribution of gray whales during the Pleistocene and Holocene. Our results reveal that dispersal between the Pacific and Atlantic was climate dependent and occurred both during the Pleistocene prior to the last glacial period and the early Holocene immediately following the opening of the Bering Strait. Genetic diversity in the Atlantic declined over an extended interval that predates the period of intensive commercial whaling, indicating this decline may have been precipitated by Holocene climate or other ecological causes. These first genetic data for Atlantic gray whales, particularly when combined with predictive habitat models for the year 2100, suggest that two recent sightings of gray whales in the Atlantic may represent the beginning of the expansion of this species' habitat beyond its currently realized range.

Pre-whaling genetic diversity and population ecology in eastern Pacific gray whales: insights from ancient DNA and stable isotopes.

Commercial whaling decimated many whale populations, including the eastern Pacific gray whale, but little is known about how population dynamics or ecology differed prior to these removals. Of particular interest is the possibility of a large population decline prior to whaling, as such a decline could explain the ~5-fold difference between genetic estimates of prior abundance and estimates based on historical records. We analyzed genetic (mitochondrial control region) and isotopic information from modern and prehistoric gray whales using serial coalescent simulations and Bayesian skyline analyses to test for a pre-whaling decline and to examine prehistoric genetic diversity, population dynamics and ecology. Simulations demonstrate that significant genetic differences observed between ancient and modern samples could be caused by a large, recent population bottleneck, roughly concurrent with commercial whaling. Stable isotopes show minimal differences between modern and ancient gray whale foraging ecology. Using rejection-based Approximate Bayesian Computation, we estimate the size of the population bottleneck at its minimum abundance and the pre-bottleneck abundance. Our results agree with previous genetic studies suggesting the historical size of the eastern gray whale population was roughly three to five times its current size.

Comparing evolutionary patterns and variability in the mitochondrial control region and cytochrome B in three species of baleen whales.

The rapidly evolving mitochondrial control region remains an important source of information on phylogeography and demographic history for cetaceans and other vertebrates, despite great uncertainty in the rate of nucleotide substitution across both nucleotide positions and lineages. Patterns of variation in linked markers with slower rates of evolution can potentially be used to calibrate the rate of nucleotide substitution in the control region and to better understand the interplay of evolutionary and demographic forces across the mitochondrial genome above and below the species level. We have examined patterns of diversity within and between three baleen whale species (gray, humpback, and Antarctic minke whales) in order to determine how patterns of molecular evolution differ between cytochrome b and the control region. Our results show that cytochrome b is less variable than expected given the diversity in the control region for gray and humpback whales, even after functional differences are taken into account, but more variable than expected for minke whales. Differences in the frequency distributions of polymorphic sites and in best-fit models of nucleotide substitution indicate that these patterns may be the result of hypervariability in the control region in gray and humpback whales but, in minke whales, may result from a large, stable or expanding population size coupled with saturation at the control region. Using paired cytochrome b and control region data across individuals, we show that the average rate of nucleotide substitution in the control region may be on average 2.6 times higher than phylogenetically derived estimates in cetaceans. These results highlight the complexity of making inferences from control region data alone and suggest that applying simple rules of DNA sequence analyses across species may be difficult.

Mitochondrial and nuclear genetic variation across calving lagoons in Eastern North Pacific gray whales (Eschrichtius robustus).

Accurate knowledge of population structure in cetaceans is critical for preserving and managing breeding habitat, particularly when habitat is not uniformly protected. Most eastern gray whales return to their major breeding range each winter along the Pacific coast of Baja California, Mexico, concentrating in 3 major calving lagoons, but it is unknown whether genetic differences exist between lagoons. Previous photo-identification studies and genetic studies suggest that gray whales may return to their natal lagoons to breed, potentially resulting in the buildup of genetic differences. However, an earlier genetic study used only one genetic marker and did not include samples from Bahia Magdalena, a major calving lagoon not currently designated as a wildlife refuge. To expand on this previous study, we collected genetic data from the mitochondrial control region (442 bp) and 9 microsatellite markers from 112 individuals across all 3 major calving lagoons. Our data suggest that migration rates between calving lagoons are high but that a small but significant departure from panmixia exists between Bahia Magdalena and Laguna San Ignacio (Fisher's Exact test, P < 0.0001; F(ST) = 0.006, P = 0.025). Coalescent simulations show that the lack of extensive population structure may result from the disruption of structure due to whaling. Another possibility is that rates of migration have always been high (>10% per generation). In addition, microsatellite data showed evidence of a severe population bottleneck. Eastern gray whales are still recovering from the impacts of whaling on their breeding grounds, and these populations should be protected and monitored for future genetic changes.

DNA evidence for historic population size and past ecosystem impacts of gray whales.

Ecosystem restoration may require returning threatened populations of ecologically pivotal species to near their former abundances, but it is often difficult to estimate historic population size of species that have been heavily exploited. Eastern Pacific gray whales play a key ecological role in their Arctic feeding grounds and are widely thought to have returned to their prewhaling abundance. Recent mortality spikes might signal that the population has reached long-term carrying capacity, but an alternative is that this decline was due to shifting climatic conditions on Arctic feeding grounds. We used a genetic approach to estimate prewhaling abundance of gray whales and report DNA variability at 10 loci that is typical of a population of approximately 76,000-118,000 individuals, approximately three to five times more numerous than today's average census size of 22,000. Coalescent simulations indicate these estimates may include the entire Pacific metapopulation, suggesting that our average measurement of approximately 96,000 individuals was probably distributed between the eastern and currently endangered western Pacific populations. These levels of genetic variation suggest the eastern population is at most at 28-56% of its historical abundance and should be considered depleted. If used to inform management, this would halve acceptable human-caused mortality for this population from 417 to 208 per year. Potentially profound ecosystem impacts may have resulted from a decline from 96,000 gray whales to the current population. At previous levels, gray whales may have seasonally resuspended 700 million cubic meters of sediment, as much as 12 Yukon Rivers, and provided food to a million sea birds.