The adaptive evolution of Leuciscus waleckii in Lake Dali Nur and convergent evolution of Cypriniformes fishes inhabiting extremely alkaline environments

Leuciscus waleckii is widely distributed in Northeast Asia and has high economic value. The population in Lake Dali Nur can adapt to extremely alkaline-saline water with bicarbonate over 50 mmol/L (pH 9.6), thus providing an exceptional model for exploring the mechanisms of adaptive evolution under extreme alkaline environments. Here, we assembled a high-quality chromosome-level reference genome for L. waleckii from Lake Dali Nur. Based on the resequencing of 85 individuals from divergent populations, the historical population size of L.waleckii in Lake Dali Nur dramatically expanded in a thousand years approximately 13,000 years ago, and experienced a cliff recession in the process of adapting to the alkaline environment of Lake Dali Nur approximately 6,000 years ago. Genome scans between freshwater and alkaline populations further revealed the significant selective sweep regions from Lake Dali Nur, which harbour a set of candidate genes involved in hypoxia tolerance, ion transport, acid-base regulation and nitrogen metabolism. 5 alkali population-specific nonsynonymous mutations were identified in CA15 gene copies. In addition, two sites with convergent amino acid mutation were detected in the RHCG-a gene among several alkali environment adapted Cypriniformes fish. Our findings provide comprehensive insight into the genomic mechanisms of L. waleckii and reveal their adaptative evolution under extreme alkaline environments.

The Chromosome-Level Genome of Triplophysa dalaica (Cypriniformes: Cobitidae) Provides Insights into Its Survival in Extremely Alkaline Environment

Lake Dali Nur, located in Inner Mongolia, North China, is alkaline, with Triplophysa dalaica one of the three fish species that not only survive, but thrive, in the lake. To investigate the presence of molecular mutations potentially responsible for this adaptation, the whole-genome sequence of the species was sequenced. A total of 126.5 and 106 Gb data, covering nearly 200× of the estimated genome, were generated using long-read sequencing and Hi-C technology, respectively. De novo assembly generated a genome totalled 607.91 Mb, with a contig N50 of 9.27 Mb. Nearly all whole-genome sequences were anchored and oriented onto 25 chromosomes, with telomeres for most chromosomes also being recovered. Repeats comprised approximately 35.01% of the whole genome. A total of 23,925 protein-coding genes were predicted, within which, 98.62% could be functionally annotated. Through comparisons of T. dalaica, T. tibetana, and T. siluroides gene models, a total of 898 genes were identified as likely being subjected to positive selection, with several of them potentially associated with alkaline adaptation, such as sodium bicarbonate cotransporter, SLC4A4. Demographic analyses suggested that the Dali population might have diverged from endemic freshwater Hai River populations, approximately 1 Ma. The high-quality T. dalaica genome, created in this study, not only aids in the analyses of alkaline adaptation, but may also assist in revealing the mysteries of the highly divergent genus Triplophysa in the future.

Identification and Analysis of Long Non-coding RNAs in Leuciscus waleckii Adapted to Highly Alkaline Conditions

Leuciscus waleckii is a freshwater fish that is known to inhabit the Dali Nor Lake, Inner Mongolia, China. The water in this lake has an HCO3 -/CO3 2- concentration of 54 mM (pH 9.6) and a salinity of 0.6‰. The physiological mechanisms that allow this fish to tolerate these saline/alkaline conditions have yet to be elucidated. Transcriptional component analysis has shown that the expression levels of a large number of genes involved in the pathways responsible for osmo-ionoregulation and arachidonic acid metabolism pathway expression change significantly (p < 0.05) during the regulation of acid-base balance under high alkaline stress. In this study, we investigated the role of long non-coding RNAs (lncRNAs) during adaptation to high alkaline conditions. Fish were challenged to an NaHCO3-adjusted alkalinity of 0 mM, 30 mM (pH 9.44 ± 0.08), and 50 mM (pH 9.55 ± 0.06) for 20 days in the laboratory. Gill and kidney tissues were then collected for high-throughput sequencing assays. A total of 159 million clean reads were obtained by high-throughput sequencing, and 41,248 lncRNA transcripts were identified. Of these, the mean number of exons and the mean length of the lncRNA transcripts were 4.8 and 2,079 bp, respectively. Based on the analysis of differential lncRNA transcript expression, a total of 5,244 and 6,571 lncRNA transcripts were found to be differentially expressed in the gills and kidneys, respectively. Results derived from Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the coding genes were correlated with the lncRNA expression profiles. GO analysis showed that many lncRNAs were enriched in the following processes: "transporter activity," "response to stimulus," and "binding." KEGG analysis further revealed that metabolic pathways were significantly enriched. A random selection of 16 lncRNA transcripts was tested by RT-qPCR; these results were consistent with our sequencing results. We found that a large number of genes, with the same expression profiles as those with differentially expressed lncRNAs, were associated with the regulation of acid-base balance, ion transport, and the excretion of ammonia and nitrogen. Collectively, our data indicate that lncRNA-regulated gene expression plays an important role in the process of adaptation to high alkaline conditions in L. waleckii.

Genomic Signature of Shifts in Selection and Alkaline Adaptation in Highland Fish

Understanding how organisms adapt to aquatic life at high altitude is fundamental in evolutionary biology. This objective has been addressed primarily related to hypoxia adaptation by recent comparative studies, whereas highland fish has also long suffered extreme alkaline environment, insight into the genomic basis of alkaline adaptation has rarely been provided. Here, we compared the genomes or transcriptomes of 15 fish species, including two alkaline tolerant highland fish species and their six alkaline intolerant relatives, three alkaline tolerant lowland fish species, and four alkaline intolerant species. We found putatively consistent patterns of molecular evolution in alkaline tolerant species in a large number of shared orthologs within highland and lowland fish taxa. Remarkably, we identified consistent signatures of accelerated evolution and positive selection in a set of shared genes associated with ion transport, apoptosis, immune response, and energy metabolisms in alkaline tolerant species within both highland and lowland fish taxa. This is one of the first comparative studies that began to elucidate the consistent genomic signature of alkaline adaptation shared by highland and lowland fish. This finding also highlights the adaptive molecular evolution changes that support fish adapting to extreme environments at high altitude.

Investigation of the Physiology of the Obligate Alkaliphilic Bacillus marmarensis GMBE 72T Considering Its Alkaline Adaptation Mechanism for Poly(3-hydroxybutyrate) Synthesis

The novel extreme obligate alkaliphilic Bacillus marmarensis DSM 21297 is known to produce polyhydroxybutyrate (PHB). However, the detailed mechanism of PHB synthesis in B. marmarensis is still unknown. Here, we investigated which metabolic pathways and metabolic enzymes are responsible for PHB synthesis in order to understand the regulatory pathway and optimize PHB synthesis in B. marmarensis. In accordance with the fact that beta-galactosidase, 3-hydroxyacyl-CoA dehydrogenase, and Enoyl-CoA hydratase together with acyl-CoA dehydrogenase and lipase were annotated in B. marmarensis according to the RAST server, we used glucose, lactose, and olive oil to understand the preferred metabolic pathway for the PHB synthesis. It was found that B. marmarensis produces PHB from glucose, lactose, and olive oil. However, the highest PHB titer and the highest amount of PHB synthesized per dry cell mass (Y_P/X) were achieved in the presence of lactose, as compared to glucose and olive oil. Additionally, in the absence of peptone, the amount of PHB synthesized is reduced for each carbon source. Interestingly, none of the carbon sources studied yielded an efficient PHB synthesis, and supplementation of the medium with potassium ions did not enhance PHB synthesis. According to these experimental results and the presence of annotated metabolic enzymes based on the RAST server, PHB accumulation in the cells of B. marmarensis could be improved by the level of the expression of 3-hydroxybutyryl-CoA dehydrogenase (1.1.1.157), which increases the production of NADPH. Additionally, the accumulation of 3-hydroxyacyl-CoA could enhance the production of PHB in B. marmarensis in the presence of fatty acids. To our knowledge, this is the first report investigating the regulatory system involved in the control of PHB metabolism of B. marmarensis.