Transcriptome sequencing and analysis of wild Amur Ide (Leuciscus waleckii) inhabiting an extreme alkaline-saline lake reveals insights into stress adaptation

Rh proteins and H+ transporters involved in ammonia excretion in Amur Ide (Leuciscus waleckii) under high alkali exposure

High alkaline environment can lead to respiratory alkalosis and ammonia toxification to freshwater fish. However, the Amur ide (Leuciscus waleckii), which inhabits an extremely alkaline lake in China with titratable alkalinity up to 53.57 mM (pH 9.6) has developed special physiological and molecular mechanisms to adapt to such an environment. Nevertheless, how the Amur ide can maintain acid-base balance and perform ammonia detoxification effectively remains unclear. Therefore, this study was designed to study the ammonia excretion rate (Tamm), total nitrogen accumulation in blood and tissues, including identification, expression, and localization of ammonia-related transporters in gills of both the alkali and freshwater forms of the Amur ide. The results showed that the freshwater form Amur ide does not have a perfect ammonia excretion mechanism exposed to high-alkaline condition. Nevertheless, the alkali form of Amur ide was able to excrete ammonia better than freshwater from Amur ide, which was facilitated by the ionocytes transporters (Rhbg, Rhcg1, Na+/H+ exchanger 2 (NHE2), and V-type H+ ATPase (VHA)) in the gills. Converting ammonia into urea served as an ammonia detoxication strategy to reduced endogenous ammonia accumulation under high-alkaline environment.

The Biochemical Composition and Quality of Adult Chinese Mitten Crab Eriocheir sinensis Reared in Carbonate-Alkalinity Water

Saline-alkaline aquaculture has become an important breakthrough in expanding the space available for aquaculture in China. However, the biochemical composition and quality of Eriocheir sinensis reared in carbonate alkalinity water are still unclear. Therefore, this study investigated the edible yield, coloration, and nutritional and flavor quality of Eriocheir sinensis. A significantly lower gonadosomatic index (GSI), meat yield (MY), and total edible yield (TEY) were detected in intensive pond (IP) samples than those in semi-intensive reed wetland (SIWR) (p < 0.05). Six color parameters in the hepatopancreas (p < 0.05) differed between IP and SIRW. The contents of crude protein and fat in the female hepatopancreas of IP were significantly higher than those in SIRW (p < 0.05). The concentrations of total monounsaturated fatty acids (∑MUFA), total essential fatty acids (∑EFA), and hypocholesterolaemic/hypercholesterolaemic ratio (h/H) in the female edible tissues checked were higher in IP than those in SIRW, with significant differences including ∑MUFA in the hepatopancreas and ovary, ∑EFA in the muscle, and h/H in the ovary (p < 0.05). Higher total free amino acid (∑FAA) contents of muscle were detected in SIRW than that in IP samples. Significantly higher K, Ca, Mg, Fe, and Zn contents in the ovary were detected in SIRW samples compared to IP (p < 0.05). Overall, Eriocheir sinensis reared in carbonate-alkalinity water is an important source of nutrients.

Effects of Saline-Alkaline Stress on Metabolome, Biochemical Parameters, and Histopathology in the Kidney of Crucian Carp (Carassius auratus)

The salinization of the water environment caused by human activities and global warming has increased which has brought great survival challenges to aquatic animals. Crucian carp (Carassius auratus) is an essential freshwater economic fish with superior adaptability to saline-alkali water. However, the physiological regulation mechanism of crucian carp adapting to saline-alkali stress remains still unclear. In this study, crucian carp were exposed to freshwater or 20, 40, and 60 mmol/L NaHCO₃ water environments for 30 days, the effects of saline-alkali stress on the kidney were evaluated by histopathology, biochemical assays and metabolomics analysis from renal function, antioxidant capacity and metabolites level. Our results showed different degrees of kidney damage at different exposure concentrations, which were characterized by glomerular atrophy and swelling, renal tubular degranulation, obstruction and degeneration, renal interstitial edema, renal cell proliferation and necrosis. Saline-alkali stress could change the levels of several physiological parameters with renal function and antioxidant capacity, including creatinine (CREA), urea nitrogen (BUN), superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GSH-Px) and malondialdehyde (MDA). In addition, metabolomics analysis showed that differential metabolites (DMs) were involved in various metabolic pathways, including phenylalanine, tyrosine, and tryptophan biosynthesis, aminoacyl-tRNA biosynthesis, purine metabolism, glycerophospholipid metabolism, sphingolipid metabolism, glycolysis/gluconeogenesis and the TCA cycle. In general, our study revealed that saline-alkaline stress could cause significant changes in renal function and metabolic profiles, and induce severe damage in the crucian carp kidney through destroying the anti-oxidant system and energy homeostasis, inhibiting protein and amino acid catabolism, as well as disordering purine metabolism and lipid metabolism. This study could contribute to a deeper understanding the adverse effects of saline-alkali stress on crucian carp kidney and the regulatory mechanism in the crucian carp of saline-alkali adaptation at the metabolic level.

Transcriptomic and proteomic analyses provide insights into the adaptive responses to the combined impact of salinity and alkalinity in Gymnocypris przewalskii

Gymnocypris przewalskii is the only high-land endemic teleost living in Qinghai Lake, the largest saline-alkaline lake in China. Its osmoregulatory physiology remains elusive due to a lack of precise identification of the response proteins. In the present study, DIA/SWATH was used to identify differentially expressed proteins (DEPs) under alkaline (pH = 10.1, carbonate buffer), saline (12‰, sodium chloride), and saline-alkaline [carbonate buffer (pH = 10.1) plus 11‰ sodium chloride] stresses. A total of 66,056 unique peptides representing 7,150 proteins and 230 DEPs [the false discovery rate (FDR) ≤ 0.05, fold change (FC) ≥ 1.5] were identified under different stresses. Comparative analyses of the proteome and transcriptome indicated that over 86% of DEPs did not show consistent trends with mRNA. In addition to consistent enrichment results under different stresses, the specific DEPs involved in saline-alkaline adaptation were primarily enriched in functions of homeostasis, hormone synthesis and reactions of defense response, complement activation and reproductive development. Meanwhile, a protein-protein interaction (PPI) network analysis of these specific DEPs indicated that the hub genes were ITGAX, MMP9, C3, F2, CD74, BTK, ANXA1, NCKAP1L, and CASP8. This study accurately isolated the genes that respond to stress, and the results could be helpful for understanding the physiological regulation mechanisms regarding salinity, alkalinity, and salinity-alkalinity interactions.

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.

Resequencing and SNP discovery of Amur ide (Leuciscus waleckii) provides insights into local adaptations to extreme environments

Amur ide (Leuciscus waleckii), a Cyprinid species, is broadly distributed in Northeast Asia. Different from its freshwater counterparts, the population in Lake Dali Nor has a strong alkalinity tolerance and can adapt to extremely alkali-saline water with bicarbonate over 50 mmol/L. To uncover the genetic basis of its alkaline adaptation, three populations, including one alkali form from Lake Dali Nor (DL), one freshwater form from its adjacent sister Lake Ganggeng Nor (GG), and one freshwater form from its historical origin, namely, the Songhua River (SH), were analyzed using genome resequencing technology. A total of 679.82 Gb clean data and 38,091,163 high-quality single-nucleotide polymorphism (SNP) loci were detected in the three populations. Nucleotide diversity and population structure analysis revealed that the DL and GG populations have lower nucleotide diversities and different genetic structures than those of the SH population. Selective sweeping showed 21 genes involved in osmoregulatory regulation (DLG1, VIPR1, AKT1, and GNAI1), inflammation and immune responses (DLG1, BRINP1, CTSL, TRAF6, AKT1, STAT3, GNAI1, SEC22b, and PSME4b), and cardiorespiratory development (TRAF6, PSME4b, STAT3, AKT1, and COL9A1) to be associated with alkaline adaption of the DL population. Interestingly, selective pressure (CodeML, MEME, and FEL) methods identified two functional codon sites of VIPR1 to be under positive selection in the DL population. The subsequent 3D protein modeling confirmed that these selected sites will incur changes in protein structure and function in the DL population. In brief, this study provides molecular evidence of population divergence and alkaline adaptation, which will be very useful for revealing the genetic basis of alkaline adaptation in Amur ide.

Analysis of Natural Selection of Immune Genes in Spinibarbus caldwelli by Transcriptome Sequencing

Spinibarbus caldwelli is an omnivorous cyprinid fish that is distributed widely in China. To investigate the adaptive evolution of S. caldwelli, the muscle transcriptome was sequenced by Illumina HiSeq 4000 platform. A total of 80,447,367 reads were generated by next-generation sequencing. Also, 211,386 unigenes were obtained by de novo assembly. Additionally, we calculated that the divergence time between S. caldwelli and Sinocyclocheilus grahami is 23.14 million years ago (Mya). And both of them diverged from Ctenopharyngodon idellus 46.95 Mya. Furthermore, 38 positive genes were identified by calculating Ka/Ks ratios from 9225 orthologs. Among them, several immune-related genes were identified as positively selected, such as POLR3B, PIK3C3, TOPORS, FASTKD3, CYPLP1A1, and UACA. Our results throw light on the nature of the natural selection of S. caldwelli and contribute to future immunological and transcriptome studies.

Characterization and analysis of the transcriptome in Gymnocypris selincuoensis on the Qinghai-Tibetan Plateau using single-molecule long-read sequencing and RNA-seq

The lakes on the Qinghai-Tibet Plateau (QTP) are the largest and highest lake group in the world. Gymnocypris selincuoensis is the only cyprinid fish living in lake Selincuo, the largest lake on QTP. However, its genetic resource is still blank, limiting studies on molecular and genetic analysis. In this study, the transcriptome of G. selincuoensis was first generated by using PacBio Iso-Seq and Illumina RNA-seq. A full-length (FL) transcriptome with 75,435 transcripts was obtained by Iso-Seq with N50 length of 3,870 bp. Among all transcripts, 75,016 were annotated to public databases, 64,710 contain complete open reading frames and 2,811 were long non-coding RNAs. Based on all- vs.-all BLAST, 2,069 alternative splicing events were detected, and 80% of them were validated by reverse transcription polymerase chain reaction (RT-PCR). Tissue gene expression atlas showed that the number of detected expressed transcripts ranged from 37,397 in brain to 19,914 in muscle, with 10,488 transcripts detected in all seven tissues. Comparative genomic analysis with other cyprinid fishes identified 77 orthologous genes with potential positive selection (Ka/Ks > 0.3). A total of 56,696 perfect simple sequence repeats were identified from FL transcripts. Our results provide valuable genetic resources for further studies on adaptive evolution, gene expression and population genetics in G. selincuoensis and other congeneric fishes.

Transcriptome changes for Nile tilapia (Oreochromis niloticus) in response to alkalinity stress

Nile tilapia is an important economic fish in the world because of its fast growth, high meat yield and strong adaptability. It is more adaptable to high alkalinity than common freshwater fish and provides valuable material for developing alkaline-tolerant strains and understanding the adaptation mechanism of fish to extreme environmental stress. In this study, we employed high throughput RNA sequencing to reveal the tissues (gill, kidney and liver) transcriptome differences of O. niloticus at different carbonate alkalinities (FW, AW40 and AW60). A total of 1,369,381,790 raw reads were obtained, including 496,441,232 reads in FW group, 437,907,696 reads in AW40 and 435,032,862 reads in AW60. In addition, 484,555,626 reads in gill, 451,618,224 reads in kidney and 433,207,940 reads in liver. A large number of stress-regulated changes were detected comprehensively. We focused on 3 significantly change pathways (steroid biosynthesis, drug metabolism and protein digestion/absorption) and 17 DEGs (HMG-CoA reductase, UDP-glucuronosyltransferase, and carbonic anhydrase etc.) which were shared among compared groups (AW40 vs FW, AW60 vs FW, AW40 vs 60 AW60) in gill, kidney and liver, respectively. These pathways/genes are sensitive to alkalinity stress and crucial to the alkalinity adaptation of tilapia. Overall, we found a large number of candidate genes, which encode important regulators of stress tolerance and ultimately contribute to future alkaline-tolerant fish breeding. Among these genes, lipid metabolism (involving signal transduction), detoxification and immune related genes are more prominent to the response and adaptability of fish to alkalinity stress.

The complete mitochondrial genome of Amur ide (Leuciscus waleckii waleckii)

In this study, the complete mitochondrial genome of Leuciscus waleckii waleckii was sequenced and got a whole length of 16605 bp. This genome was contain 2 rRNA, 22 tRNA, 13 protein-coding genes, 1 control region (D-loop) and 1 replication origin. And the nucleotide composition of this mitochondrial genome is 27.72% for A, 26.28% for T, 27.23% for C and 18.77% for G. To clarify the phylogenetic relationship of the Leuciscus waleckii waleckii, we concluded the phylogenetic tree using 12 PCGs (except ND6) of mitochondrial genome in Leuciscus waleckii waleckii and 16 other cyprinid fish by Bayesian inference (BI) methods and maximum-likelihood (ML). And the result show that Leuciscus waleckii waleckii was close to other Leuciscus species, especially Leuciscus baicalensis.

Characterization, functional analysis, and expression levels of three carbonic anhydrases in response to pH and saline-alkaline stresses in the ridgetail white prawn Exopalaemon carinicauda

Carbonate alkalinity, salinity, and pH are three important stress factors for aquatic animals in saline-alkaline water. Carbonic anhydrases (CAs) catalyze the reversible reaction of CO2 reported to play an important role in the acid-base regulation in vertebrates. To explore the molecular mechanism of CAs efficacy in shrimp after their transfer into saline-alkaline water, the cDNAs of three CAs (EcCAc, EcCAg, and EcCAb) were cloned from Exopalaemon carinicauda. Sequence analysis showed that EcCAc and EcCAg both possessed a conserved α-CA domain and a proton acceptor site, and EcCAb contained a Pro-CA domain. Tissue expression analysis demonstrated that EcCAc and EcCAg were most abundantly in gills, and EcCAb was highly expressed in muscle. The cumulative mortalities remained below 25% under exposure to pH (pH 6 and pH 9), low salinity (5 ppt), or high carbonate alkalinity (5 and 10 mmol/L) after 72 h of exposure. However, mortalities increased up to 70% under extreme saline-alkaline stress (salinity 5 ppt, carbonate alkalinity 10 mmol/L, and pH 9) after 14 days of exposure. The EcCAc and EcCAg expressions in gills were significantly upregulated during the early period of pH and saline-alkaline stresses, while the EcCAb expressions showed no regular or large changes. The two-way ANOVA found significant interactions between salinity and carbonate alkalinity observed in EcCAc, EcCAg, and EcCAb expressions (p < 0.05). Furthermore, an RNA interference experiment resulted in increased mortality of EcCAc- and EcCAg-silenced prawns under saline-alkaline stress. EcCAc knockdown reduced expressions of Na+/H+ exchanger (EcNHE) and sodium bicarbonate cotransporter (EcNBC), and EcCAg knockdown reduced EcCAc, EcNHE, EcNBC, and V-type H+-ATPase (EcVTP) expressions. These results suggest EcCAc and EcCAg as important modulators in response to pH and saline-alkaline stresses in E. carinicauda.

Insights into the molecular basis of immunosuppression and increasing pathogen infection severity of ammonia toxicity by transcriptome analysis in pacific white shrimp Litopenaeus vannamei

The high concentration of ammonia resulting from intensive culture system and environmental pollution could cause disease occurrence in shrimp, but little information is available on its molecular mechanisms. In this study, we performed comparative transcriptome analysis among WSSV-infected shrimp under ammonia stress (LAV), WSSV-infected shrimp under normal water (LV), and normal shrimp under ammonia stress (LA) groups to identify the key genes and pathways involved in immunosuppression and increasing pathogen infection severity caused by ammonia toxicity in Litopenaeus vannamei. Totally, 526 significantly differential expressed genes (DEGs) were identified in LAV group compared to LV and LA groups, among which 270 genes were lost expressed and 67 genes uniquely expressed in the LAV group. According to the public functional reports for the annotated DEGs, they potentially involved in the following functions: (1) accelerating pathogen adhesion, invasion and multiplication; (2) reducing the ability for pathogen defense and immune response; (3) inhibiting positive regulation of apoptotic and antioxidant defense for host homeostasis; (4) inhibiting transcription and protein transport; (5) and increasing protein methylation and ubiquitination, etc. A total of 13 pathways were obtained mainly involving in this process, which mainly led to the following changes: (1) increasing the immunosuppression, anemia, endocrine dysfunction, neurotoxic effect and neuroinvasion, atherosclerosis and thrombogenesis, blood-brain barrier penetration, thyroid disorder, necrosis, inflammation, and circadian disturbance; (2) reducing the ability of vascular remodeling, angiogenesis, cell survival, migration, apoptosis, and lymph transferred to blood stream; (3) leading to cell hypertrophy, cellular shape changes, and mesangial matrix expansion. The present results firstly supplied molecular mechanisms for the ammonia toxicity inhibiting the immune system and increasing pathogen infection severity in shrimp, which is a prerequisite for better understanding the pathogenesis caused by ammonia toxicity.

A transcriptome study on Macrobrachium nipponense hepatopancreas experimentally challenged with white spot syndrome virus (WSSV)

White spot syndrome virus (WSSV) is one of the most devastating pathogens of cultured shrimp, responsible for massive loss of its commercial products worldwide. The oriental river prawn Macrobrachium nipponense is an economically important species that is widely farmed in China and adult prawns can be infected by WSSV. However, the molecular mechanisms of the host pathogen interaction remain unknown. There is an urgent need to learn the host pathogen interaction between M. nipponense and WSSV which will be able to offer a solution in controlling the spread of WSSV. Next Generation Sequencing (NGS) was used in this study to determin the transcriptome differences by the comparison of control and WSSV-challenged moribund samples, control and WSSV-challenged survived samples of hepatopancreas in M. nipponense. A total of 64,049 predicted unigenes were obtained and classified into 63 functional groups. Approximately, 4,311 differential expression genes were identified with 3,308 genes were up-regulated when comparing the survived samples with the control. In the comparison of moribund samples with control, 1,960 differential expression genes were identified with 764 genes were up-regulated. In the contrast of two comparison libraries, 300 mutual DEGs with 95 up-regulated genes and 205 down-regulated genes. All the DEGs were performed GO and KEGG analysis, overall a total of 85 immune-related genes were obtained and these gene were groups into 13 functions and 4 KEGG pathways, such as protease inhibitors, heat shock proteins, oxidative stress, pathogen recognition immune receptors, PI3K/AKT/mTOR pathway, MAPK signaling pathway and Ubiquitin Proteasome Pathway. Ten genes that valuable in immune responses against WSSV were selected from those DEGs to furture discuss the response of host to WSSV. Results from this study contribute to a better understanding of the immune response of M. nipponense to WSSV, provide information for identifying novel genes in the absence of genome of M. nipponense. Furthermore, large number of transcripts obtained from this study could provide a strong basis for future genomic research on M. nipponense.

Transcriptome sequencing and comparative analysis of differentially-expressed isoforms in the roots of Halogeton glomeratus under salt stress

Although Halogeton glomeratus (H. glomeratus) has been confirmed to have a unique mechanism to regulate Na⁺ efflux from the cytoplasm and compartmentalize Na⁺ into leaf vacuoles, little is known about the salt tolerance mechanisms of roots under salinity stress. In the present study, transcripts were sequenced using the BGISEQ-500 sequencing platform (BGI, Wuhan, China). After quality control, approximately 24.08 million clean reads were obtained and the average mapping ratio to the reference gene was 70.00%. When comparing salt-treated samples with the control, a total of 550, 590, 1411 and 2063 DEIs were identified at 2, 6, 24 and 72h, respectively. Numerous differentially-expressed isoforms that play important roles in response and adaptation to salt condition are related to metabolic processes, cellular processes, single-organism processes, localization, biological regulation, responses to stimulus, binding, catalytic activity and transporter activity. Fifty-eight salt-induced isoforms were common to different stages of salt stress; most of these DEIs were related to signal transduction and transporters, which maybe the core isoforms regulating Na⁺ uptake and transport in the roots of H. glomeratus. The expression patterns of 18 DEIs that were detected by quantitative real-time polymerase chain reaction were consistent with their respective changes in transcript abundance as identified by RNA-Seq technology. The present study thoroughly explored potential isoforms involved in salt tolerance on H. glomeratus roots at five time points. Our results may serve as an important resource for the H. glomeratus research community, improving our understanding of salt tolerance in halophyte survival under high salinity stress.

De novo assembly and annotation of the whole transcriptome of Oratosquilla oratoria

As the representative species of Stomatopoda, Oratosquilla oratoria can tolerate the complex benthic environment. In this study, our goal was to develop the transcriptomic resource of O. oratoria that would support adaptation mechanism studies. We generated the whole transcriptome of O. oratoria from combined tissues (eyestalk, muscle, sexual and viscus) using Hiseq technology. A total of 51,305,284 high-quality clean reads were assembled to produce 59,054 non-redundant transcripts with a mean length of 987nt using Trinity and CORSET software. Among the predictable unigenes, a total of 32,451 unigenes were annotated based on protein databases. Finally, we predicted the coding sequences of 31,822 unigenes and obtained 19,057 SSRs in the present study. The present study will provide an important resource and foundational understanding for future genomic research of O. oratoria.

Genomic Basis of Adaptive Evolution: The Survival of Amur Ide (Leuciscus waleckii) in an Extremely Alkaline Environment

The Amur ide (Leuciscus waleckii) is a cyprinid fish that is widely distributed in Northeast Asia. The Lake Dali Nur population inhabits one of the most extreme aquatic environments on Earth, with an alkalinity up to 50 mmol/L (pH 9.6), thus providing an exceptional model with which to characterize the mechanisms of genomic evolution underlying adaptation to extreme environments. Here, we developed the reference genome assembly for L. waleckii from Lake Dali Nur. Intriguingly, we identified unusual expanded long terminal repeats (LTRs) with higher nucleotide substitution rates than in many other teleosts, suggesting their more recent insertion into the L. waleckii genome. We also identified expansions in genes encoding egg coat proteins and natriuretic peptide receptors, possibly underlying the adaptation to extreme environmental stress. We further sequenced the genomes of 10 additional individuals from freshwater and 18 from Lake Dali Nur populations, and we detected a total of 7.6 million SNPs from both populations. In a genome scan and comparison of these two populations, we identified a set of genomic regions under selective sweeps that harbor genes involved in ion homoeostasis, acid-base regulation, unfolded protein response, reactive oxygen species elimination, and urea excretion. Our findings provide comprehensive insight into the genomic mechanisms of teleost fish that underlie their adaptation to extreme alkaline environments.

Transcriptome Analysis of the Hepatopancreas in the Pacific White Shrimp (Litopenaeus vannamei) under Acute Ammonia Stress

In the practical farming of Litopenaeus vannamei, the intensive culture system and environmental pollution usually results in a high concentration of ammonia, which usually brings large detrimental effects to shrimp, such as increasing the susceptibility to pathogens, reducing growth, decreasing osmoregulatory capacity, increasing the molting frequency, and even causing high mortality. However, little information is available on the molecular mechanisms of the detrimental effects of ammonia stress in shrimp. In this study, we performed comparative transcriptome analysis between ammonia-challenged and control groups from the same family of L. vannamei to identify the key genes and pathways response to ammonia stress. The comparative transcriptome analysis identified 136 significantly differentially expressed genes that have high homologies with the known proteins in aquatic species, among which 94 genes are reported potentially related to immune function, and the rest of the genes are involved in apoptosis, growth, molting, and osmoregulation. Fourteen GO terms and 6 KEGG pathways were identified to be significantly changed by ammonia stress. In these GO terms, 13 genes have been studied in aquatic species, and 11 of them were reported potentially involved in immune defense and two genes were related to molting. In the significantly changed KEGG pathways, all the 7 significantly changed genes have been reported in shrimp, and four of them were potentially involved in immune defense and the other three were related to molting, defending toxicity, and osmoregulation, respectively. In addition, majority of the significantly changed genes involved in nitrogen metabolisms that play an important role in reducing ammonia toxicity failed to perform the protection function. The present results have supplied molecular level support for the previous founding of the detrimental effects of ammonia stress in shrimp, which is a prerequisite for better understanding the molecular mechanism of the immunosuppression from ammonia stress.

A transcriptome study on Macrobrachium rosenbergii hepatopancreas experimentally challenged with white spot syndrome virus (WSSV)

The world production of shrimp such as the Malaysian giant freshwater prawn, Macrobrachium rosenbergii is seriously affected by the white spot syndrome virus (WSSV). There is an urgent need to understand the host pathogen interaction between M. rosenbergii and WSSV which will be able to provide a solution in controlling the spread of this infectious disease and lastly save the aquaculture industry. Now, using Next Generation Sequencing (NGS), we will be able to capture the response of the M. rosenbergii to the pathogen and have a better understanding of the host defence mechanism. Two cDNA libraries, one of WSSV-challenged M. rosenbergii and a normal control one, were sequenced using the Illumina HiSeq™ 2000 platform. After de novo assembly and clustering of the unigenes from both libraries, 63,584 standard unigenes were generated with a mean size of 698bp and an N50 of 1137bp. We successfully annotated 35.31% of all unigenes by using BLASTX program (E-value <10-5) against NCBI non-redundant (Nr), Swiss-Prot, Kyoto Encyclopedia of Genes and Genome pathway (KEGG) and Orthologous Groups of proteins (COG) databases. Gene Ontology (GO) assessment was conducted using BLAST2GO software. Differentially expressed genes (DEGs) by using the FPKM method showed 8443 host genes were significantly up-regulated whereas 5973 genes were significantly down-regulated. The differentially expressed immune related genes were grouped into 15 animal immune functions. The present study showed that WSSV infection has a significant impact on the transcriptome profile of M. rosenbergii's hepatopancreas, and further enhanced the knowledge of this host-virus interaction. Furthermore, the high number of transcripts generated in this study will provide a platform for future genomic research on freshwater prawns.

Comparative transcriptome analysis reveals molecular strategies of oriental river prawn Macrobrachium nipponense in response to acute and chronic nitrite stress

Macrobrachium nipponense is an economically and nutritionally important species threatened by ambient superfluous nitrite. De novo RNA-Seq was used to explore the molecular mechanism in M. nipponense exposed to the acute nitrite stress (26.05 mg/L nitrite-N) for 24 h and the chronic nitrite stress (1.38 mg/L nitrite-N) for 28 d A total of 175.13 million reads were obtained and assembled into 58,871 unigenes with an average length of 1028.7 bp and N50 of 1294 bp. Under the acute and chronic nitrite stress trials, 2824 and 2610 unigenes were significantly expressed. In GO analysis and KEGG pathway analysis, 30 pathways were significantly different between the two treatments while four pathways were in common and the markedly altered pathways were divided into four sections as immunity, metabolism, cell and others. The immunity section revealing the different depth of immunity provoked by nitrite stress contained the most pathways including the important pathways as phagosome, folate biosynthesis, glycerolipid metabolism, glycine, serine and threonine metabolism, amino sugar and nucleotide sugar metabolism under the acute nitrite stress, and lysosome, alanine, aspartate and glutamate metabolism, arginine and proline metabolism under the chronic nitrite stress. This is the first report of responses of M. nipponense under acute and chronic nitrite stress through de novo transcriptome sequencing on the transcriptome level. The results of transcriptome analysis improve our understanding on the underlying molecular mechanisms coping with nitrite stress in crustacean species.

Defining the role of omics in assessing ecosystem health: Perspectives from the Canadian environmental monitoring program

Scientific reviews and studies continue to describe omics technologies as the next generation of tools for environmental monitoring, while cautioning that there are limitations and obstacles to overcome. However, omics has not yet transitioned into national environmental monitoring programs designed to assess ecosystem health. Using the example of the Canadian Environmental Effects Monitoring (EEM) program, the authors describe the steps that would be required for omics technologies to be included in such an established program. These steps include baseline collection of omics endpoints across different species and sites to generate a range of what is biologically normal within a particular ecosystem. Natural individual variability in the omes is not adequately characterized and is often not measured in the field, but is a key component to an environmental monitoring program, to determine the critical effect size or action threshold for management. Omics endpoints must develop a level of standardization, consistency, and rigor that will allow interpretation of the relevance of changes across broader scales. To date, population-level consequences of routinely measured endpoints such as reduced gonad size or intersex in fish is not entirely clear, and the significance of genome-wide molecular, proteome, or metabolic changes on organism or population health is further removed from the levels of ecological change traditionally managed. The present review is not intended to dismiss the idea that omics will play a future role in large-scale environmental monitoring studies, but rather outlines the necessary actions for its inclusion in regulatory monitoring programs focused on assessing ecosystem health.

Comparative Transcriptome Analysis in the Hepatopancreas Tissue of Pacific White Shrimp Litopenaeus vannamei Fed Different Lipid Sources at Low Salinity

RNA-seq was used to compare the transcriptomic response of hepatopancreas in juvenile Litopenaeus vannamei fed three diets with different lipid sources, including beef tallow (BT), fish oil (FO), and an equal combination of soybean oil + BT + linseed oil (SBL) for 8 weeks at 3 practical salinity unit (psu). A total of 9622 isogenes were annotated in 316 KEGG pathways and 39, 42 and 32 pathways significantly changed in the paired comparisons of FO vs SBL, BT vs SBL, or FO vs BT, respectively. The pathways of glycerolipid metabolism, linoleic acid metabolism, arachidonic acid metabolism, glycerophospholipid metabolism, fatty acid biosynthesis, fatty acid elongation, fatty acid degradation, and biosynthesis of unsaturated fatty acid were significantly changed in all paired comparisons between dietary lipid sources, and the pathways of glycerolipid metabolism, linoleic acid metabolism, arachidonic acid metabolism and glycerophospholipid metabolism significantly changed in the FO vs SBL and BT vs SBL comparisons. These pathways are associated with energy metabolism and cell membrane structure. The results indicate that lipids sources affect the adaptation of L. vannamei to low salinity by providing extra energy or specific fatty acids to change gill membrane structure and control iron balance. The results of this study lay a foundation for further understanding lipid or fatty acid metabolism in L. vannamei at low salinity.

Transcriptome and Molecular Pathway Analysis of the Hepatopancreas in the Pacific White Shrimp Litopenaeus vannamei under Chronic Low-Salinity Stress

The Pacific white shrimp Litopenaeus vannamei is a euryhaline penaeid species that shows ontogenetic adaptations to salinity, with its larvae inhabiting oceanic environments and postlarvae and juveniles inhabiting estuaries and lagoons. Ontogenetic adaptations to salinity manifest in L. vannamei through strong hyper-osmoregulatory and hypo-osmoregulatory patterns and an ability to tolerate extremely low salinity levels. To understand this adaptive mechanism to salinity stress, RNA-seq was used to compare the transcriptomic response of L. vannamei to changes in salinity from 30 (control) to 3 practical salinity units (psu) for 8 weeks. In total, 26,034 genes were obtained from the hepatopancreas tissue of L. vannamei using the Illumina HiSeq 2000 system, and 855 genes showed significant changes in expression under salinity stress. Eighteen top Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly involved in physiological responses, particularly in lipid metabolism, including fatty-acid biosynthesis, arachidonic acid metabolism and glycosphingolipid and glycosaminoglycan metabolism. Lipids or fatty acids can reduce osmotic stress in L. vannamei by providing additional energy or changing the membrane structure to allow osmoregulation in relevant organs, such as the gills. Steroid hormone biosynthesis and the phosphonate and phosphinate metabolism pathways were also involved in the adaptation of L. vannamei to low salinity, and the differential expression patterns of 20 randomly selected genes were validated by quantitative real-time PCR (qPCR). This study is the first report on the long-term adaptive transcriptomic response of L. vannamei to low salinity, and the results will further our understanding of the mechanisms underlying osmoregulation in euryhaline crustaceans.

Transcriptional Profiling Reveals Differential Gene Expression of Amur Ide (Leuciscus waleckii) during Spawning Migration

Amur ide (Leuciscus waleckii), an important aquaculture species, inhabits neutral freshwater but can tolerate high salinity or alkalinity. As an extreme example, the population in Dali Nor lake inhabits alkalized soda water permanently, and migrates from alkaline water to neutral freshwater to spawn. In this study, we performed comparative transcriptome profiling study on the livers of Amur ide to interrogate the expression differences between the population that permanently inhabit freshwater in Ganggeng Nor lake (FW) and the spawning population that recently migrated from alkaline water into freshwater (SM). A total of 637,234,880 reads were generated, resulting in 53,440 assembled contigs that were used as reference sequences. Comparisons of these transcriptome files revealed 444 unigenes with significant differential expression (p-value ≤ 0.01, fold-change ≥ 2), including 246 genes that were up-regulated in SM and 198 genes that were up-regulated in FW. The gene ontology (GO) enrichment analysis and KEGG pathway analysis indicated that the mTOR signaling pathway, Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway, and oxidative phosphorylation were highly likely to affect physiological changes during spawning migration. Overall, this study demonstrates that transcriptome changes played a role in Amur ide spawning migration. These results provide a foundation for further analyses on the physiological and molecular mechanisms underlying Amur ide spawning migration.

RNA-seq analysis of Macrobrachium rosenbergii hepatopancreas in response to Vibrio parahaemolyticus infection

Background

The Malaysian giant freshwater prawn, Macrobrachium rosenbergii, is an economically important crustacean worldwide. However, production of this prawn is facing a serious threat from Vibriosis disease caused by Vibrio species such as Vibrio parahaemolyticus. Unfortunately, the mechanisms involved in the immune response of this species to bacterial infection are not fully understood. We therefore used a high-throughput deep sequencing technology to investigate the transcriptome and comparative expression profiles of the hepatopancreas from this freshwater prawn infected with V. parahaemolyticus to gain an increased understanding of the molecular mechanisms underlying the species’ immune response to this pathogenic bacteria.

Result

A total of 59,122,940 raw reads were obtained from the control group, and 58,385,094 reads from the Vibrio-infected group. Via de novo assembly by Trinity assembler, 59,050 control unigenes and 73,946 Vibrio-infected group unigenes were obtained. By clustering unigenes from both libraries, a total of 64,411 standard unigenes were produced. The standard unigenes were annotated against the NCBI non-redundant, Swiss-Prot, Kyoto Encyclopaedia of Genes and Genome pathway (KEGG) and Orthologous Groups of Proteins (COG) databases, with 19,799 (30.73%), 16,832 (26.13%), 14,706 (22.83%) and 7,856 (12.19%) hits respectively, giving a final total of 22,455 significant hits (34.86% of all unigenes). A Gene Ontology (GO) analysis search using the Blast2GO program resulted in 6,007 unigenes (9.32%) being categorized into 55 functional groups. A differential gene expression analysis produced a total of 14,569 unigenes aberrantly expressed, with 11,446 unigenes significantly up-regulated and 3,103 unigenes significantly down-regulated. The differentially expressed immune genes fall under various processes of the animal immune system.

Conclusion

This study provided an insight into the antibacterial mechanism in M. rosenbergii and the role of differentially expressed immune genes in response to V. parahaemolyticus infection. Furthermore, this study has generated an abundant list of transcript from M.rosenbergii which will provide a fundamental basis for future genomics research in this field.

Electronic supplementary material

The online version of this article (doi:10.1186/s13099-015-0052-6) contains supplementary material, which is available to authorized users.

Transcriptomic profiling of the salt-stress response in the halophyte Halogeton glomeratus

BACKGROUND

Halogeton glomeratus (H. glomeratus) is an extreme halophyte that is widely distributed in arid regions, including foothills, the Gobi desert of northwest China, and the marginal loess of Central Asia. However, research on the salt-tolerant mechanisms and genes of this species are limited because of a lack of genomic sequences. In the present study, the transcriptome of H. glomeratus was analyzed using next-generation sequencing technology to identify genes involved in salt tolerance and better understand mechanisms of salt response in the halophyte H. glomeratus.

RESULTS

Illumina RNA-sequencing was performed in five sequencing libraries that were prepared from samples treated with 200 mM NaCl for 6, 12, 24, and 72 h and a control sample to investigate changes in the H. glomeratus transcriptome in response to salt stress. The de novo assembly of five transcriptomes identified 50,267 transcripts. Among these transcripts, 31,496 (62.66%) were annotated, including 44 Gene Ontology (GO) terms and 128 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Compared with transcriptomes from the control and NaCl-treated samples, there were 2,223, 5,643, 7,510 and 10,908 genes that were differentially expressed after exposure to NaCl for 6, 12, 24, and 72 h, respectively. One hundred and eighteen salt-induced genes were common to at least two stages of salt stress, and 291 up-regulated genes were common to various stages of salt stress. Numerous genes that are related to ion transport, reactive oxygen species scavenging, energy metabolism, hormone-response pathways, and responses to biotic and abiotic stress appear to play a significant role in adaptation to salinity conditions in this species. The detection of expression patterns of 18 salt-induced genes by quantitative real-time polymerase chain reaction were basically consistent with their changes in transcript abundance determined by RNA sequencing.

CONCLUSIONS

Our findings provide a genomic sequence resource for functional genetic assignments of an extreme halophyte, H. glomeratus. We believe that the transcriptome datasets will help elucidate the genetic basis of this species' response to a salt environment and develop stress-tolerant crops based on favorable wild genetic resources.

The positive regulatory roles of the TIFY10 proteins in plant responses to alkaline stress

The TIFY family is a novel plant-specific protein family, and is characterized by a conserved TIFY motif (TIFF/YXG). Our previous studies indicated the potential roles of TIFY10/11 proteins in plant responses to alkaline stress. In the current study, we focused on the regulatory roles and possible physiological and molecular basis of the TIFY10 proteins in plant responses to alkaline stress. We demonstrated the positive function of TIFY10s in alkaline responses by using the AtTIFY10a and AtTIFY10b knockout Arabidopsis, as evidenced by the relatively lower germination rates of attify10a and attify10b mutant seeds under alkaline stress. We also revealed that ectopic expression of GsTIFY10a in Medicago sativa promoted plant growth, and increased the NADP-ME activity, citric acid content and free proline content but decreased the MDA content of transgenic plants under alkaline stress. Furthermore, expression levels of the stress responsive genes including NADP-ME, CS, H⁺-ppase and P5CS were also up-regulated in GsTIFY10a transgenic plants under alkaline stress. Interestingly, GsTIFY10a overexpression increased the jasmonate content of the transgenic alfalfa. In addition, we showed that neither GsTIFY10a nor GsTIFY10e exhibited transcriptional activity in yeast cells. However, through Y2H and BiFc assays, we demonstrated that GsTIFY10a, not GsTIFY10e, could form homodimers in yeast cells and in living plant cells. As expected, we also demonstrated that GsTIFY10a and GsTIFY10e could heterodimerize with each other in both yeast and plant cells. Taken together, our results provided direct evidence supporting the positive regulatory roles of the TIFY10 proteins in plant responses to alkaline stress.

De novo assembly and characterization of the transcriptome in the desiccation-tolerant moss Syntrichia caninervis

BACKGROUND

Syntrichia caninervis is a desiccation-tolerant moss and the dominant bryophyte of the Biological Soil Crusts (BSCs) found in the Mojave and Gurbantunggut deserts. Next generation high throughput sequencing technologies offer an efficient and economic choice for characterizing non-model organism transcriptomes with little or no prior molecular information available.

RESULTS

In this study, we employed next generation, high-throughput, Illumina RNA-Seq to analyze the poly-(A) + mRNA from hydrated, dehydrating and desiccated S. caninervis gametophores. Approximately 58.0 million paired-end short reads were obtained and 92,240 unigenes were assembled with an average size of 493 bp, N50 value of 662 bp and a total size of 45.48 Mbp. Sequence similarity searches against five public databases (NR, Swiss-Prot, COSMOSS, KEGG and COG) found 54,125 unigenes (58.7%) with significant similarity to an existing sequence (E-value ≤ 1e-5) and could be annotated. Gene Ontology (GO) annotation assigned 24,183 unigenes to the three GO terms: Biological Process, Cellular Component or Molecular Function. GO comparison between P. patens and S. caninervis demonstrated similar sequence enrichment across all three GO categories. 29,370 deduced polypeptide sequences were assigned Pfam domain information and categorized into 4,212 Pfam domains/families. Using the PlantTFDB, 778 unigenes were predicted to be involved in the regulation of transcription and were classified into 49 transcription factor families. Annotated unigenes were mapped to the KEGG pathways and further annotated using MapMan. Comparative genomics revealed that 44% of protein families are shared in common by S. caninervis, P. patens and Arabidopsis thaliana and that 80% are shared by both moss species.

CONCLUSIONS

This study is one of the first comprehensive transcriptome analyses of the moss S. caninervis. Our data extends our knowledge of bryophyte transcriptomes, provides an insight to plants adapted to the arid regions of central Asia, and continues the development of S. caninervis as a model for understanding the molecular aspects of desiccation-tolerance.

SNP discovery from transcriptome of the swimbladder of Takifugu rubripes

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies in many species. High-throughput sequencing of RNA was developed primarily to analyze global gene expression, while it is an efficient way to discover SNPs from the expressed genes. In this study, we conducted transcriptome sequencing of the swimbladder of Takifugu rubripes using Illumina HiSeq2000 platform to identify gene-associated SNPs in the swimbladder. A total of 30,312,181 unique-mapped-reads were obtained from 44,736,850 raw reads. A total of 62,270 putative SNPs were discovered, which were located in 11,306 expressed genes and 2,246 scaffolds. The average minor allele frequency (MAF) of the SNPs was 0.26. GO and KEGG pathway analysis were conducted to analyze the genes containing SNPs. Validation of selected SNPs revealed that 54% of SNPs (26/48) were true SNPs. The results suggest that RNA-Seq is an efficient and cost-effective approach to discover gene-associated SNPs. In this study, a large number of SNPs were identified and these data will be useful resources for population genetic study, evolution analysis, resource assessment, genetic linkage analysis and genome-wide association studies.

Transcriptome sequencing revealed the genes and pathways involved in salinity stress of Chinese mitten crab, Eriocheir sinensis

A total of 276.9 million reads were obtained and assembled into 206, 371 contigs with an average length of 614 bp and N50 of 1,470 bp. Comparison of digital gene expression between treatment and control group reveals 1,151 and 941 genes were significantly differentially expressed in crab gill and muscle, respectively. In gill and muscle, protein ubiquitination, ubiquinone biosynthesis, oxidative phosphorylation, and mitochondria dysfunction pathways were the top pathways differentially expressed following the challenge. EIF 2 signaling pathway and IGF-1 signaling pathway were the top ones among the signal-related pathways. Most of the amino acid metabolism pathways were found to be involved in this process. The expression patterns of 15 differentially expressed genes were validated by quantitative real-time RT-PCR (average correlation coefficient 0.80). This is the first report of expression analysis of genes and pathways involved in osmoregulation of Eriocheir sinensis through transcriptome sequencing. The findings of this study will further promote the understanding of the underlying molecular mechanism of salinity stress adaptation for crustacean species.

Gene expression changes leading extreme alkaline tolerance in Amur ide (Leuciscus waleckii) inhabiting soda lake

Background

Amur ide (Leuciscus waleckii) is an economically and ecologically important cyprinid species in Northern Asia. The Dali Nor population living in the soda lake Dali Nor can adapt the extremely high alkalinity, providing us a valuable material to understand the adaptation mechanism against extreme environmental stress in teleost.

Results

In this study, we generated high-throughput RNA-Seq data from three tissues gill, liver and kidney of L. waleckii living in the soda lake Dali Nor and the fresh water lake Ganggeng Nor, then performed parallel comparisons of three tissues. Our results showed that out of assembled 64,603 transcript contigs, 28,391 contigs had been assigned with a known function, corresponding to 20,371 unique protein accessions. We found 477, 2,761 and 3,376 differentially expressed genes (DEGs) in the gill, kidney, and liver, respectively, of Dali Nor population compared to Ganggeng Nor population with FDR ≤ 0.01and fold-change ≥ 2. Further analysis revealed that well-known functional categories of genes and signaling pathway, which are associated with stress response and extreme environment adaptation, have been significantly enriched, including the functional categories of “response to stimulus”, “transferase activity”, “transporter activity” and “oxidoreductase activity”, and signaling pathways of “mTOR signaling”, “EIF2 signaling”, “superpathway of cholesterol biosynthesis”. We also identified significantly DEGs encoding important modulators on stress adaptation and tolerance, including carbonic anhydrases, heat shock proteins, superoxide dismutase, glutathione S-transferases, aminopeptidase N, and aminotransferases.

Conclusions

Overall, this study demonstrated that transcriptome changes in L. waleckii played a role in adaptation to complicated environmental stress in the highly alkalized Dali Nor lake. The results set a foundation for further analyses on alkaline-responsive candidate genes, which help us understand teleost adaptation under extreme environmental stress and ultimately benefit future breeding for alkaline-tolerant fish strains.