A selective sweep in the Spike gene has driven SARS-CoV-2 human adaptation.

The coronavirus disease 2019 (COVID-19) pandemic underscores the need to better understand animal-to-human transmission of coronaviruses and adaptive evolution within new hosts. We scanned more than 182,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes for selective sweep signatures and found a distinct footprint of positive selection located around a non-synonymous change (A1114G; T372A) within the spike protein receptor-binding domain (RBD), predicted to remove glycosylation and increase binding to human ACE2 (hACE2), the cellular receptor. This change is present in all human SARS-CoV-2 sequences but not in closely related viruses from bats and pangolins. As predicted, T372A RBD bound hACE2 with higher affinity in experimental binding assays. We engineered the reversion mutant (A372T) and found that A372 (wild-type [WT]-SARS-CoV-2) enhanced replication in human lung cells relative to its putative ancestral variant (T372), an effect that was 20 times greater than the well-known D614G mutation. Our findings suggest that this mutation likely contributed to SARS-CoV-2 emergence from animal reservoirs or enabled sustained human-to-human transmission.

Humoral and T-cell-mediated responses to an insect-specific flavivirus-based Zika virus vaccine candidate.

Flaviviruses represent a significant global health threat and relatively few licensed vaccines exist to protect against them. Insect-specific flaviviruses (ISFVs) are incapable of replication in humans and have emerged as a novel and promising tool for flavivirus vaccine development. ISFV-based flavivirus vaccines have shown exceptional safety, immunogenicity, and efficacy, however, a detailed assessment of the correlates of protection and immune responses induced by these vaccines are still needed for vaccine optimization. Here, we explore the mechanisms of protective immunity induced by a previously created pre-clinical Zika virus (ZIKV) vaccine candidate, called Aripo/Zika (ARPV/ZIKV). In brief, immunocompromised IFN-αßR-/- mice passively immunized with ARPV/ZIKV immune sera experienced protection after lethal ZIKV challenge, although this protection was incomplete. ARPV/ZIKV-vaccinated IFN-αßR-/- mice depleted of CD4+ or CD8+ T-cells at the time of ZIKV challenge showed no morbidity or mortality. However, the adoptive transfer of ARPV/ZIKV-primed T-cells into recipient IFN-αßR-/- mice resulted in a two-day median increase in survival time compared to controls. Altogether, these results suggest that ARPV/ZIKV-induced protection is primarily mediated by neutralizing antibodies at the time of challenge and that T-cells may play a comparatively minor but cumulative role in the protection observed. Lastly, ARPV/ZIKV-vaccinated Tcra KO mice, which are deficient in T-cell responses, experienced significant mortality post-challenge. These results suggest that ARPV/ZIKV-induced cell-mediated responses are critical for development of protective immune responses at vaccination. Despite the strong focus on neutralizing antibody responses to novel flavivirus vaccine candidates, these results suggest that cell-mediated responses induced by ISFV-based vaccines remain important to overall protective responses.

The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Post-transcriptional control of the essential enzyme MurF by a small regulatory RNA in Brucella abortus.

Brucella abortus is a facultative, intracellular, zoonotic pathogen that resides inside macrophages during infection. This is a specialized niche where B. abortus encounters various stresses as it navigates through the macrophage. In order to survive this harsh environment, B. abortus utilizes post-transcriptional regulation of gene expression through the use of small regulatory RNAs (sRNAs). Here, we characterize a Brucella sRNAs called MavR (for MurF- and virulence-regulating sRNA), and we demonstrate that MavR is required for the full virulence of B. abortus in macrophages and in a mouse model of chronic infection. Transcriptomic and proteomic studies revealed that a major regulatory target of MavR is MurF. MurF is an essential protein that catalyzes the final cytoplasmic step in peptidoglycan (PG) synthesis; however, we did not detect any differences in the amount or chemical composition of PG in the ΔmavR mutant. A 6-nucleotide regulatory seed region within MavR was identified, and mutation of this seed region resulted in dysregulation of MurF production, as well as significant attenuation of infection in a mouse model. Overall, the present study underscores the importance of sRNA regulation in the physiology and virulence of Brucella.

The E2 glycoprotein holds key residues for Mayaro virus adaptation to the urban Aedes aegypti mosquito.

Adaptation to mosquito vectors suited for transmission in urban settings is a major driver in the emergence of arboviruses. To better anticipate future emergence events, it is crucial to assess their potential to adapt to new vector hosts. In this work, we used two different experimental evolution approaches to study the adaptation process of an emerging alphavirus, Mayaro virus (MAYV), to Ae. aegypti, an urban mosquito vector of many other arboviruses. We identified E2-T179N as a key mutation increasing MAYV replication in insect cells and enhancing transmission after escaping the midgut of live Ae. aegypti. In contrast, this mutation decreased viral replication and binding in human fibroblasts, a primary cellular target of MAYV in humans. We also showed that MAYV E2-T179N generates reduced viremia and displays less severe tissue pathology in vivo in a mouse model. We found evidence in mouse fibroblasts that MAYV E2-T179N is less dependent on the Mxra8 receptor for replication than WT MAYV. Similarly, exogenous expression of human apolipoprotein receptor 2 and Mxra8 enhanced WT MAYV replication compared to MAYV E2-T179N. When this mutation was introduced in the closely related chikungunya virus, which has caused major outbreaks globally in the past two decades, we observed increased replication in both human and insect cells, suggesting E2 position 179 is an important determinant of alphavirus host-adaptation, although in a virus-specific manner. Collectively, these results indicate that adaptation at the T179 residue in MAYV E2 may result in increased vector competence-but coming at the cost of optimal replication in humans-and may represent a first step towards a future emergence event.

Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel.

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

Seasonal changes in recombination characteristics in a natural population of Drosophila melanogaster.

Environmental seasonality is a potent evolutionary force, capable of maintaining polymorphism, promoting phenotypic plasticity and causing bet-hedging. In Drosophila, environmental seasonality has been reported to affect life-history traits, tolerance to abiotic stressors and immunity. Oscillations in frequencies of alleles underlying fitness-related traits were also documented alongside SNPs across the genome. Here, we test for seasonal changes in two recombination characteristics, crossover rate and crossover interference, in a natural D. melanogaster population from India using morphological markers of the three major chromosomes. We show that winter flies, collected after the dry season, have significantly higher desiccation tolerance than their autumn counterparts. This difference proved to hold also for hybrids with three independent marker stocks, suggesting its genetic rather than plastic nature. Significant between-season changes are documented for crossover rate (in 9 of 13 studied intervals) and crossover interference (in four of eight studied pairs of intervals); both single and double crossovers were usually more frequent in the winter cohort. The winter flies also display weaker plasticity of both recombination characteristics to desiccation. We ascribe the observed differences to indirect selection on recombination caused by directional selection on desiccation tolerance. Our findings suggest that changes in recombination characteristics can arise even after a short period of seasonal adaptation (~8-10 generations).

CAGm: a repository of germline microsatellite variations in the 1000 genomes project.

The human genome harbors an abundance of repetitive DNA; however, its function continues to be debated. Microsatellites-a class of short tandem repeat-are established as an important source of genetic variation. Array length variants are common among microsatellites and affect gene expression; but, efforts to understand the role and diversity of microsatellite variation has been hampered by several challenges. Without adequate depth, both long-read and short-read sequencing may not detect the variants present in a sample; additionally, large sample sizes are needed to reveal the degree of population-level polymorphism. To address these challenges we present the Comparative Analysis of Germline Microsatellites (CAGm): a database of germline microsatellites from 2529 individuals in the 1000 genomes project. A key novelty of CAGm is the ability to aggregate microsatellite variation by population, ethnicity (super population) and gender. The database provides advanced searching for microsatellites embedded in genes and functional elements. All data can be downloaded as Microsoft Excel spreadsheets. Two use-case scenarios are presented to demonstrate its utility: a mononucleotide (A) microsatellite at the BAT-26 locus and a dinucleotide (CA) microsatellite in the coding region of FGFRL1. CAGm is freely available at http://www.cagmdb.org/.

DNA sequencing of anatomy lab cadavers to provide hands-on precision medicine introduction to medical students.

BACKGROUND: Medical treatment informed by Precision Medicine is becoming a standard practice for many diseases, and patients are curious about the consequences of genomic variants in their genome. However, most medical students' understanding of Precision Medicine derives from classroom lectures. This format does little to foster an understanding for the potential and limitations of Precision Medicine. To close this gap, we implemented a hands-on Precision Medicine training program utilizing exome sequencing to prepare a clinical genetic report of cadavers studied in the anatomy lab. The program reinforces Precision Medicine related learning objectives for the Genetics curriculum. METHODS: Pre-embalmed blood samples and embalmed tissue were obtained from cadavers (donors) used in the anatomy lab. DNA was isolated and sequenced and illustrative genetic reports provided to the students. The reports were used to facilitate discussion with students on the implications of pathogenic genomic variants and the potential correlation of these variants in each "donor" with any anatomical anomalies identified during cadaver dissection. RESULTS: In 75% of cases, analysis of whole exome sequencing data identified a variant associated with increased risk for a disease/abnormal condition noted in the donor's cause of death or in the students' anatomical findings. This provided students with real-world examples of the potential relationship between genomic variants and disease risk. Our students also noted that diseases associated with 92% of the pathogenic variants identified were not related to the anatomical findings, demonstrating the limitations of Precision Medicine. CONCLUSION: With this study, we have established protocols and classroom procedures incorporating hands-on Precision Medicine training in the medical student curriculum and a template for other medical educators interested in enhancing their Precision Medicine training program. The program engaged students in discovering variants that were associated with the pathophysiology of the cadaver they were studying, which led to more exposure and understanding of the potential risks and benefits of genomic medicine.

The Komodo dragon (Varanus komodoensis) genome and identification of innate immunity genes and clusters.

BACKGROUND: We report the sequencing, assembly and analysis of the genome of the Komodo dragon (Varanus komodoensis), the largest extant lizard, with a focus on antimicrobial host-defense peptides. The Komodo dragon diet includes carrion, and a complex milieu of bacteria, including potentially pathogenic strains, has been detected in the saliva of wild dragons. They appear to be unaffected, suggesting that dragons have robust defenses against infection. While little information is available regarding the molecular biology of reptile immunity, it is believed that innate immunity, which employs antimicrobial host-defense peptides including defensins and cathelicidins, plays a more prominent role in reptile immunity than it does in mammals. . RESULTS: High molecular weight genomic DNA was extracted from Komodo dragon blood cells. Subsequent sequencing and assembly of the genome from the collected DNA yielded a genome size of 1.6 Gb with 45x coverage, and the identification of 17,213 predicted genes. Through further analyses of the genome, we identified genes and gene-clusters corresponding to antimicrobial host-defense peptide genes. Multiple ß-defensin-related gene clusters were identified, as well as a cluster of potential Komodo dragon ovodefensin genes located in close proximity to a cluster of Komodo dragon ß-defensin genes. In addition to these defensins, multiple cathelicidin-like genes were also identified in the genome. Overall, 66 ß-defensin genes, six ovodefensin genes and three cathelicidin genes were identified in the Komodo dragon genome. CONCLUSIONS: Genes with important roles in host-defense and innate immunity were identified in this newly sequenced Komodo dragon genome, suggesting that these organisms have a robust innate immune system. Specifically, multiple Komodo antimicrobial peptide genes were identified. Importantly, many of the antimicrobial peptide genes were found in gene clusters. We found that these innate immunity genes are conserved among reptiles, and the organization is similar to that seen in other avian and reptilian species. Having the genome of this important squamate will allow researchers to learn more about reptilian gene families and will be a valuable resource for researchers studying the evolution and biology of the endangered Komodo dragon.

Genomic signatures of experimental adaptive radiation in Drosophila.

Abiotic environmental factors play a fundamental role in determining the distribution, abundance and adaptive diversification of species. Empowered by new technologies enabling rapid and increasingly accurate examination of genomic variation in populations, researchers may gain new insights into the genomic background of adaptive radiation and stress resistance. We investigated genomic variation across generations of large-scale experimental selection regimes originating from a single founder population of Drosophila melanogaster, diverging in response to ecologically relevant environmental stressors: heat shock, heat knock down, cold shock, desiccation and starvation. When compared to the founder population, and to parallel unselected controls, there were more than 100,000 single nucleotide polymorphisms (SNPs) displaying consistent allelic changes in response to selective pressures across generations. These SNPs were found in both coding and noncoding sequences, with the highest density in promoter regions, and involved a broad range of functionalities, including molecular chaperoning by heat-shock proteins. The SNP patterns were highly stressor-specific despite considerable variation among line replicates within each selection regime, as reflected by a principal component analysis, and co-occurred with selective sweep regions. Only ~15% of SNPs with putatively adaptive changes were shared by at least two selective regimes, while less than 1% of SNPs diverged in opposite directions. Divergent stressors driving evolution in the experimental system of adaptive radiation left distinct genomic signatures, most pronounced in starvation and heat-shock selection regimes.

Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity.

Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.

Characterization of Three Small Proteins in Brucella abortus Linked to Fucose Utilization.

Elucidating the function of proteins <50 amino acids in length is no small task. Nevertheless, small proteins can play vital roles in the lifestyle of bacteria and influence the virulence of pathogens; thus, the investigation of the small proteome is warranted. Recently, our group identified the Brucella abortus protein VtlR as a transcriptional activator of four genes, one of which is the well-studied small regulatory RNA AbcR2, while the other three genes encode hypothetical small proteins, two of which are highly conserved among the order Rhizobiales This study provides evidence that all three genes encode authentic small proteins and that all three are highly expressed under oxidative stress, low-pH, and stationary-phase growth conditions. Fractionation of the cells revealed that the proteins are localized to the membranes of B. abortus We demonstrate that the small proteins under the transcriptional control of VtlR are not accountable for attenuation observed with the B. abortusvtlR deletion strain. However, there is an association between VtlR-regulated genes and growth inhibition in the presence of the sugar l-fucose. Subsequent transcriptomic analyses revealed that B. abortus initiates the transcription of a locus encoding a putative sugar transport and utilization system when the bacteria are cultured in the presence of l-fucose. Altogether, our observations characterize the role of the VtlR-controlled small proteins BAB1_0914, BAB2_0512, and BAB2_0574 in the biology of B. abortus, particularly in the capacity of the bacteria to utilize l-fucose.IMPORTANCE Despite being one of the most common zoonoses worldwide, there is currently no human vaccine to combat brucellosis. Therefore, a better understanding of the pathogenesis and biology of Brucella spp., the causative agent of brucellosis, is essential for the discovery of novel therapeutics against these highly infectious bacteria. In this study, we further characterize the virulence-associated transcriptional regulator VtlR in Brucella abortus Our findings not only shed light on our current understanding of a virulence related genetic system in Brucella spp. but also increase our knowledge of small proteins in the field of bacteriology.

Unique divergence of the breast cancer 2 (BRCA2) gene in Neanderthals.

Unique divergence of the BRCA2, a tumor suppressor gene, in Neanderthals relative to other primates, including modern humans, is highlighted. This divergence with potentially pathogenic consequences raises a question about cancer susceptibility in the archaic species that was replaced by modern humans about 40,000 years ago.

Discovery of Novel Antimicrobial Peptides from Varanus komodoensis (Komodo Dragon) by Large-Scale Analyses and De-Novo-Assisted Sequencing Using Electron-Transfer Dissociation Mass Spectrometry.

Komodo dragons are the largest living lizards and are the apex predators in their environs. They endure numerous strains of pathogenic bacteria in their saliva and recover from wounds inflicted by other dragons, reflecting the inherent robustness of their innate immune defense. We have employed a custom bioprospecting approach combining partial de novo peptide sequencing with transcriptome assembly to identify cationic antimicrobial peptides from Komodo dragon plasma. Through these analyses, we identified 48 novel potential cationic antimicrobial peptides. All but one of the identified peptides were derived from histone proteins. The antimicrobial effectiveness of eight of these peptides was evaluated against Pseudomonas aeruginosa (ATCC 9027) and Staphylococcus aureus (ATCC 25923), with seven peptides exhibiting antimicrobial activity against both microbes and one only showing significant potency against P. aeruginosa. This study demonstrates the power and promise of our bioprospecting approach to cationic antimicrobial peptide discovery, and it reveals the presence of a plethora of novel histone-derived antimicrobial peptides in the plasma of the Komodo dragon. These findings may have broader implications regarding the role that intact histones and histone-derived peptides play in defending the host from infection. Data are available via ProteomeXChange with identifier PXD005043.

Nucleotide diversity inflation as a genome-wide response to experimental lifespan extension in Drosophila melanogaster.

BACKGROUND: Evolutionary theory predicts that antagonistically selected alleles, such as those with divergent pleiotropic effects in early and late life, may often reach intermediate population frequencies due to balancing selection, an elusive process when sought out empirically. Alternatively, genetic diversity may increase as a result of positive frequency-dependent selection and genetic purging in bottlenecked populations. RESULTS: While experimental evolution systems with directional phenotypic selection typically result in at least local heterozygosity loss, we report that selection for increased lifespan in Drosophila melanogaster leads to an extensive genome-wide increase of nucleotide diversity in the selected lines compared to replicate control lines, pronounced in regions with no or low recombination, such as chromosome 4 and centromere neighborhoods. These changes, particularly in coding sequences, are most consistent with the operation of balancing selection and the antagonistic pleiotropy theory of aging and life history traits that tend to be intercorrelated. Genes involved in antioxidant defenses, along with multiple lncRNAs, were among those most affected by balancing selection. Despite the overwhelming genetic diversification and the paucity of selective sweep regions, two genes with functions important for central nervous system and memory, Ptp10D and Ank2, evolved under positive selection in the longevity lines. CONCLUSIONS: Overall, the 'evolve-and-resequence' experimental approach proves successful in providing unique insights into the complex evolutionary dynamics of genomic regions responsible for longevity.

A Test for Gene Flow among Sympatric and Allopatric Hawaiian Picture-Winged Drosophila.

The Hawaiian Drosophila are one of the most species-rich endemic groups in Hawaii and a spectacular example of adaptive radiation. Drosophila silvestris and D. heteroneura are two closely related picture-winged Drosophila species that occur sympatrically on Hawaii Island and are known to hybridize in nature, yet exhibit highly divergent behavioral and morphological traits driven largely through sexual selection. Their closest-related allopatric species, D. planitibia from Maui, exhibits hybrid male sterility and reduced behavioral reproductive isolation when crossed experimentally with D. silvestris or D. heteroneura. A modified four-taxon test for gene flow was applied to recently obtained genomes of the three Hawaiian Drosophila species. The analysis indicates recent gene flow in sympatry, but also, although less extensive, between allopatric species. This study underscores the prevalence of gene flow, even in taxonomic groups considered classic examples of allopatric speciation on islands. The potential confounding effects of gene flow in phylogenetic and population genetics inference are discussed, as well as the implications for conservation.

Mapping Genomic Scaffolds to Chromosomes Using Laser Capture Microdissection in Application to Hawaiian Picture-Winged Drosophila.

Next-generation sequencing technologies have led to a decreased cost and an increased throughput in genome sequencing. Yet, many genome assemblies based on short sequencing reads have been assembled only to the scaffold level due to the lack of sufficient chromosome mapping information. Traditional ways of mapping scaffolds to chromosomes require a large amount of laboratory work and time to generate genetic and/or physical maps. To address this problem, we conducted a rapid technique which uses laser capture microdissection and enables mapping scaffolds of de novo genome assemblies directly to chromosomes in Hawaiian picture-winged Drosophila. We isolated and sequenced intact chromosome arms from larvae of D. differens. By mapping the reads of each chromosome to the recently assembled scaffolds from 3 Hawaiian picture-winged Drosophila species, at least 67% of the scaffolds were successfully assigned to chromosome arms. Even though the scaffolds are not ordered within a chromosome, the fast-generated chromosome information allows for chromosome-related analyses after genome assembling. We utilize this new information to test the faster-X evolution effect for the first time in these Hawaiian picture-winged Drosophila species.

Transcriptomic imprints of adaptation to fresh water: parallel evolution of osmoregulatory gene expression in the Alewife.

Comparative approaches in physiological genomics offer an opportunity to understand the functional importance of genes involved in niche exploitation. We used populations of Alewife (Alosa pseudoharengus) to explore the transcriptional mechanisms that underlie adaptation to fresh water. Ancestrally anadromous Alewives have recently formed multiple, independently derived, landlocked populations, which exhibit reduced tolerance of saltwater and enhanced tolerance of fresh water. Using RNA-seq, we compared transcriptional responses of an anadromous Alewife population to two landlocked populations after acclimation to fresh (0 ppt) and saltwater (35 ppt). Our results suggest that the gill transcriptome has evolved in primarily discordant ways between independent landlocked populations and their anadromous ancestor. By contrast, evolved shifts in the transcription of a small suite of well-characterized osmoregulatory genes exhibited a strong degree of parallelism. In particular, transcription of genes that regulate gill ion exchange has diverged in accordance with functional predictions: freshwater ion-uptake genes (most notably, the 'freshwater paralog' of Na+ /K+ -ATPase α-subunit) were more highly expressed in landlocked forms, whereas genes that regulate saltwater ion secretion (e.g. the 'saltwater paralog' of NKAα) exhibited a blunted response to saltwater. Parallel divergence of ion transport gene expression is associated with shifts in salinity tolerance limits among landlocked forms, suggesting that changes to the gill's transcriptional response to salinity facilitate freshwater adaptation.