- Open Access
Genome-wide alternative splicing profile in the posterior kidney of brown trout (Salmo trutta) during proliferative kidney disease
BMC Genomics volume 23, Article number: 446 (2022)
The cnidarian myxozoan parasite Tetracapsuloides bryosalmonae causes chronic proliferative kidney disease (PKD) in salmonids. This parasite is a serious threat to wild and cultured salmonids. T. bryosalmonae undergoes intra-luminal sporogonic development in the kidney of brown trout (Salmo trutta) and the viable spores are released via urine. We investigated the alternative splicing pattern in the posterior kidney of brown trout during PKD.
RNA-seq data were generated from the posterior kidney of brown trout collected at 12 weeks post-exposure to T. bryosalmonae. Subsequently, this data was mapped to the brown trout genome. About 153 significant differently expressed alternatively spliced (DEAS) genes, (delta PSI = 5%, FDR P-value < 0.05) were identified from 19,722 alternatively spliced events. Among the DEAS genes, the least and most abundant alternative splicing types were alternative 5′ splice site (5.23%) and exon skipping (70.59%), respectively. The DEAS genes were significantly enriched for sodium-potassium transporter activity and ion homeostasis (ahcyl1, atp1a3a, atp1a1a.1, and atp1a1a.5). The protein-protein interaction network analysis enriched two local network clusters namely cation transporting ATPase C-terminus and Sodium/potassium ATPase beta chain cluster, and mixed inclusion of Ion homeostasis and EF-hand domain cluster. Furthermore, the human disease-related salmonella infection pathway was significantly enriched in the protein-protein interaction network.
This study provides the first baseline information about alternative splicing in brown trout during PKD. The generated data lay a foundation for further functional molecular studies in PKD - brown trout infection model. The information generated from the present study can help to develop therapeutic strategies for PKD in the future.
Proliferative kidney disease (PKD) is an important disease of salmonids that has both economic and ecological significance. PKD is quite prevalent across Europe and North America that affects both wild salmonid population and domestic salmonid aquaculture . PKD is caused by an evolutionarily degenerated cnidarian myxozoan endoparasite Tetracapsuloides bryosalmonae . The two hosts parasitic life-cycle of T. bryosalmonae involve salmonid and bryozoans as intermediate and definitive hosts respectively [3, 4]. Climate change driven warmer water temperature favours the spread of T. bryosalmonae and increases the severity of PKD in the affected fish . T. bryosalmonae infection along with warmer water temperature is responsible for the diminishing brown trout (Salmo trutta) population in the streams of the Alps . A massive fish kill event, mainly mountain whitefish (Prosopium williamsoni) in the Yellowstone River, Montana, USA was caused due to T. bryosalmonae that led to the subsequent temporary lockdown of the river for public access followed by a huge loss in the associated recreational industry .
The portal entry of T. bryosalmonae into the salmonid host is through gills and eventually, the parasite reaches various internal organs like the kidney, spleen, and liver . However, the kidney is the primary organ for infection where the parasite undergoes extrasporogonic development and is subsequently released via urine . The released parasite can readily infect the bryozoan host. The development of T. bryosalmonae in the affected fish causes a chronic inflammatory reaction leading to the swelling of the kidney and spleen . Interestingly, the European strain of T. bryosalmonae can infect and cause PKD in rainbow trout (Oncorhynchus mykiss), brown trout, and brook trout (Salvelinus fontinalis). However, the mature spores of the parasites can be released only by native brown trout and brook trout, whereas the exotic rainbow trout remains as a dead-end host . Moreover, brown trout can serve as a reservoir host in which the parasite can persist and get released even after five years post-recovery from the clinical signs of PKD .
Due to the PKD associated huge economic losses in the rainbow trout aquaculture industry, most of the laboratory experiments to understand the host-pathogen interaction of this disease was done in the dead-end rainbow trout host . In rainbow trout, PKD is manifested by dysregulated proliferation of B lymphocytes, hyper-secretion of immunoglobulins, imbalance in the cytokines of T-helper-like cells and up-regulation of SOCS genes [12,13,14,15,16,17,18]. However, few experiments were also done in brown trout to understand the host response during PKD at the molecular level suggests increased B cell response and higher Th1- like cytokine production [17, 19,20,21,22].
Alternative splicing is a post-transcriptional mechanism that gives rise to distinct multiple mRNAs from a single gene that may lead to proteins with a distinct function, stability, and sub-cellular localization . Alternative splicing is observed in more than 90% of human genes . Specific cis- and trans-regulatory elements control the inclusions of different exons spliced in the mature mRNA. Based on the splicing pattern alternative splicing can be classified into five basic types, (i) alternative 3′ splice site (A3SS), (ii) alternative 5′ splice site (A5SS), (iii) exon skipping (SE), (iv) intron retention (RI), and (v) mutually exclusive exons (MXE) . Furthermore, abnormal mRNA transcripts having premature termination codons (PTCs) are eliminated by non-sense mediated decay . Interestingly, about one-third of the spliced transcripts shall have such PTCs and are targeted for non-sense mediated decay .
RNA splicing increases the diversity of mRNA translated from the genome and thereby increases proteome diversity in eukaryotes. The protein isoforms formed due to alternative spliced gene variants may have different structural and functional properties . Studies suggest that any alternation in the splicing pattern can affect the normal physiological state of the organism. In humans, chronic diseases like cancer and arthritis are associated with alternatively spliced genes. Furthermore, pathogens can alter the splicing pattern of the host during host-pathogen interaction . The mechanism of alternative splicing is highly regulated and the dynamics of the splicing event changes under different physiological conditions . Moreover, alternative splicing is known to regulate various molecular functions of an organism including immune response [29, 30] and disease resistance [31, 32]. Only a few studies are available in fish to understand the alternative splicing mechanism during host-pathogen interaction. Alternative splicing plays an important role in channel catfish (Ictalurus punctatus) during heat stress to increase the tolerance of the fish . Furthermore, genes involved in the host defense mechanism were differentially expressed alternatively spliced (DEAS) in channel catfish during bacterial infections such as Edwardsiella ictaluri  and Flavobacterium columnare .
Recently, next-generation sequencing technology has made a substantial contribution in sequencing the chromosome level assembly of brown trout genome . RNA-seq based transcriptomics has opened up avenues to explore differential gene expression, post-transcriptional modifications, and novel transcripts during host-pathogen interaction . Furthermore, generating molecular information about the pathogen shall be valuable in developing efficient prophylactic measures . RNA-seq has been employed to understand and generate the transcriptome resources of T. bryosalmonae [38, 39] and its hosts [21, 40,41,42]. However, there are no studies to explore the differentially expressed alternatively spliced transcripts in salmonids during proliferative kidney disease. In this study, we utilized the RNA-seq data generated in our previous experiment from the PKD-affected brown trout , to gain insights on the alternative splicing events during host-parasite interaction.
Mapping of RNA-seq data
Relatively, a high number of interstitial pre-sporogenic stages of T. bryosalmonae in the posterior kidney of brown trout were observed using immunohistochemistry at 12 weeks post exposure (wpe) [21, 22]. Hence, the samples collected at 12 wpe were used for RNA-seq analysis. A total of 421.8 million raw 100 bp single-end RNA-seq reads were generated from the posterior kidney of brown trout (6 exposed and 6 controls) using Illumina platform. Subsequently, 420.6 million high quality reads were obtained after quality control and trimming. About 381.01 million reads (90.58%) mapped to the brown trout genome. A detailed summary of RNA-seq data and its mapping with the brown trout genome is provided in Supplementary File S1.
Alternative splicing landscape in PKD affected brown trout
A total of 19,722 alternative spliced genes were differentially expressed in the posterior kidney of brown trout during PKD. Significant DEAS genes were filtered based on differential percent spliced-in (PSI) value and FDR P-value. About 153 alternative spliced genes were filtered and identified as significant DEAS genes (delta PSI = 5%, FDR P-value < 0.05) (Table 1). Moreover, out of these 153 DEAS genes, about 77 and 76 DEAS genes had increased and decreased inclusion levels, respectively (Table 1 and Supplementary File S2). All the five basic types of alternatively spliced events A3SS, A5SS, SE, RI, and MXE (Fig. 1A) were observed among the 153 DEAS genes. The most abundant DEAS event was SE (70.59%) followed by A3SS (9.8%), RI (7.19%), MXE (7.19%), and A5SS (5.23%) (Fig. 1B). The distribution of all the five basic types of alternatively spliced events in different chromosomes (chr) of brown trout during PKD is illustrated in Fig. 1C. The significant DEAS genes were distributed in all most every chromosome of brown trout except chr11, chr16, and chr32. Moreover, the highest DEAS gene distribution was observed in chr17, chr24, and chr34 with each chromosome having 8 events. The splicing plot of filtered 153 DEAS genes computed based on mean PSI values and categorised under all the five different splicing types is demonstrated in Fig. 2. Hierarchical clustering of individual PSI values of the 153 significant DEAS genes clustered the exposed and control groups separately (Fig. 3).
Functional enrichment analysis of DEAS genes
The orthologous Ensembl gene identifiers of zebrafish were identified for the significant DEAS genes of brown trout and used for functional enrichment analysis (Supplementary File S3). Overrepresentation or functional enrichment of DEAS genes in the PKD affected brown trout is illustrated in Fig. 4. The DEAS genes were significantly (adjusted P-value < 0.05) enriched for biological processes such as cytoskeleton organization (GO:0007010) and establishment or maintenance of transmembrane electrochemical gradient (GO:0010248). Furthermore, the DEAS genes were enriched for molecular function such as P-type sodium/potassium-exchanging transporter activity (GO:0005391), P-type sodium transporter activity (GO:0008554) and P-type potassium transmembrane transporter activity (GO:0008556). In the Reactome pathaway database the DEAS genes significantly enriched for Ion homeostasis (REAC:R-DRE-5578775) pathway. However, no statistically significant enrichment of DEAS genes was found under the categories of cellular components, and KEGG pathway database (www.kegg.jp/kegg/kegg1.html) .
String database assigned 116 (75.8%) out of 153 DEAS genes identified in the posterior kidney of brown trout during PKD for the construction of a protein-protein interaction network (Fig. 5). There were 116 nodes and 302 edges in the protein-protein interaction network. The network had an average node degree of 5.21 and an average clustering coefficient of 0.326 (enriched P-value = 0.000741). The protein-protein interaction network demonstrated Ubiquitin-specific peptidase 25 (USP25) as the central node with 21 distinct interactions. Two local network clusters were predicted to be significantly enriched, (i) Cation transporting ATPase, C-terminus, and Sodium / potassium ATPase beta chain (CL:20005, FDR P-value = 0.0089), and (ii) mixed inclusion of Ion homeostasis, and EF-hand domain (CL:19893, FDR P-value = 0.0047) (Table 2 and Fig. 6 A-F). In addition, the protein-protein interaction network was significantly enriched for Salmonella infection pathway (dre05132, FDR P-value = 0.0357) under KEGG pathways  (Table 2 and Fig. 6 G-J). Abbreviations of proteins used in the protein-protein interaction network are listed in Supplementary File S4.
Validation of DEAS genes
The validation of a subset of the DEAS gene was done using reverse transcription PCR (RT-PCR). A3SS event was represented by protein kinase C binding protein 1 - like (prkcbp1l), and both A3SS long (230 bp) and A3SS short (152 bp) alternative spliced products were amplified (Fig. 7A). Bromodomain adjacent to zinc finger domain 2Ba (baz2ba) was validated as a representative of the A5SS event, and A5SS long (260 bp) and A5SS short (152 bp) alternative spliced products were amplified (Fig. 7B). SE event was represented by RAP1, GTP-GDP dissociation stimulator 1 (rap1gds1), and phosphoinositide-3-kinase adaptor protein 1 (pik3ap1) (Fig. 7C). The amplified PCR product size of the inclusion exon and skipped exon of rap1gds1 was 320 bp and 173 bp respectively. Similarly, pik3ap1 had a PCR product size of 655 bp (inclusion exon) and 229 bp (skipped exon). The alternative splicing event RI was represented by RAB guanine nucleotide exchange factor 1 (rabgef1) and both retention (246 bp) and non-retention (153 bp) isoforms were detected (Fig. 7D). The isoforms of myocyte enhancer factor 2d (mef2d) was validated as a representative of MXE alternative splice and both MXE exon 1 (241 bp) and MXE exon 2 (238 bp) were amplified (Fig. 7E). Furthermore, all the PCR products were sequenced, and the nucleotide sequences aligned exactly to their corresponding isoform sequences.
Fish are poikilothermic animals and their body temperature is constantly regulated by the temperature of the water. Furthermore, fish are constantly exposed to pathogens and other stressors in the aquatic environment. In order to cope with the unforeseeable environmental changes, evolutionarily fish might have developed a rapid regulatory molecular mechanism including post-transcriptional mechanism. Alternative splicing is a paramount post-transcriptional mechanism that regulates transcriptomics and proteomic diversity in all eukaryotic organisms . The frequency of alternative splicing is inversely proportional to the genome size of the fish and hence the fishes having smaller genomes have to depend more on the alternative splicing mechanism for the synthesis of a diverse range of proteins . The dynamics of the splicing pattern of an organism are highly influenced by environmental and physiological factors. Hence the splicing pattern and the functional role of spliced mRNA have always interested researchers. Moreover, RNA-seq based high throughput sequencing technologies have advanced the rapid identification and quantification of alternative splicing events . Investigation of alternative splicing during host-pathogen interaction can provide biological insights about the host response during disease pathogenesis. However, only a few investigations are done in fish to understand the alternative splicing mechanism during infectious diseases, particularly during bacterial [34, 45] and viral  infections. To date, there are no studies in fish to understand the alternative splicing pattern during parasitic infection.
The kidney of teleost fish has both excretory and immune functions. In addition, the anterior kidney serves as a haematopoietic organ and its functions are synonymous with the functions of the bone marrow of higher vertebrates. Apart from these functions, the kidney of teleost fish has an important physiological role in ion homeostasis and osmoregulation . Interestingly, T. bryosalmonae has been co-evolved with brown trout host and can undergo sporogenesis and persist for a long time in the kidney, and the parasitic spores are consistently shed via urine. In spite of the strong immune response, the T. bryosalmonae can proliferate in the kidney of the affected brown trout [20,21,22]. Hence kidney was chosen as the organ of interest for investigating the molecular response of brown trout during PKD.
The present study was aimed to investigate the alternative splicing mechanism in the posterior kidney of brown trout during PKD. The raw RNA-seq data generated from the PKD-affected brown trout from our previous experiment  and the recently published brown trout genomic sequence  were utilized in this study. In spite of the 19,722 DEAS genes identified during PKD pathogenesis in brown trout, only 153 were observed to be statistically significant. Moreover, we identified SE as the most common splicing pattern contributing 70.59% of the total DEAS event (Table 1 and Fig. 1). However, in the case of channel catfish, A3SS and A5SS were the most abundant DEAS event during E. ictaluri and F. columnare infections, respectively [34, 45].
Freshwater fish like brown trout dwells in a hyposmotic medium and the kidney of these animals have a major role in excreting osmotic water loads and re-absorption of filtered salts . The development of T. bryosalmonae in the kidney of PKD affected brown trout may affect the ion homeostasis of the host. In the present investigation, DEAS genes involved in the ion homeostasis, particularly in the sodium/potassium transporter activity were significantly enriched (Fig. 4). They are S-adenosylhomocysteine hydrolase-like protein 1 (ahcyl1), ATPase Na+/K+ transporting subunit alpha 3a (atp1a3a), ATPase Na+/K+ transporting subunit alpha 1a tandem duplicate 1 (atp1a1a.1), and ATPase Na+/K+ transporting subunit alpha 1a tandem duplicate 5 (atp1a1a.5).
AHCYL1 is a novel enzyme that belongs to the AHCY family of proteins responsible for the metabolism of S-adenosyl-L-homocysteine . Furthermore, AHCYL1 is involved in the inositol phospholipid signaling pathway and binds with inositol 1,4,5-trisphosphate receptor facilitating intracellular calcium (Ca2+) release leading to the activation of Ca2+ signalling . Various biological processes such as organogenesis, fertilization, and apoptosis are regulated by AHCYL1 and inositol 1,4,5-trisphosphate receptor-mediated Ca2+ signaling cascade. Experiments in zebrafish suggest that AHCYL1 has a significant function in the embryogenesis of zebrafish . Regulation of intracellular Ca2+ is important for the homeostasis of an organism and the differential regulation of AHCYL1 splice event in the present study may be a regulatory mechanism adapted by the brown host to attain homeostasis during PKD.
Sodium-Potassium ATPase is an enzyme protein present in the cell membrane of almost every cell of an animal. Moreover, this protein is highly conserved across fish and higher vertebrates . Across the cell membrane, the sodium-potassium ATPase protein facilitates the transportation of three sodium (Na+) ions out and two potassium (K+) ions into the cell with a single molecule of ATP hydrolysis . This Na+ / K+ transport has a crucial physiological function of ion homeostasis, transmembrane electrochemical Na+ gradient, nutrients transport across the cell, and cell signalling . Moreover, sodium-potassium ATPase is important for regulating the resting potential and volume of the cell. Sodium-Potassium ATPase protein is a heterodimer having alpha and beta subunits . In the present experiment, three alpha subunits (atp1a1a.1, atp1a1a.5, atp1a3a) were identified as DEAS genes during PKD pathogenesis in the kidney of brown trout. In zebrafish, a mutation in the atp1a1a.1 resulted in a slower functioning rate of the heart . Moreover, knockdown of atp1a3a in zebrafish led to dilation in the ventricle of the brain and abnormal motility . Hence sodium-Potassium ATPase is crucial for the normal physiology and biological functioning of all animals including fish.
In our previous study, transcriptome analysis in the posterior kidney of infected brown trout suggests that upregulated differentially expressed genes were enriched for the immune system process (58.54%) in the gene ontology analysis . However, in the present study, immune-related DEAS genes were not statistically enriched in the gene ontology analysis. Alternatively, sodium-potassium transporter and ion homeostasis were enriched in the gene ontology analysis of DEAS genes in the posterior kidney of brown trout during PKD. It is known that the kidney is responsible for maintaining relatively constant levels of key ions including sodium, potassium, and calcium. It is possible that the molecular function of sodium-potassium transporter activity and ion homeostasis of the kidney might have been affected in the PKD affected brown trout. Moreover, alternative splicing of genes related to sodium-potassium transporter activity and ion homeostasis may be a regulatory mechanism to maintain homeostasis  in the PKD affected brown trout. However, further studies are needed to test this hypothesis.
In the protein-protein interaction network analysis two local network clusters, namely cation transporting ATPase C-terminus and Sodium/potassium ATPase beta chain cluster (CL:20005), and mixed inclusion of Ion homeostasis and EF-hand domain cluster (CL:19893) were significantly enriched (Table 2, Figs. 5 and 6 A-F). This is similar to the results of functional enrichment analysis of DEAS genes where genes related to ion homeostasis were significantly enriched in the PKD affected brown trout (Fig. 4). Furthermore, proteins involved in the salmonella infection pathway such as filamin A alpha (flna), MYD88 innate immune signal transduction adaptor (myd88), mitogen-activated protein kinase 12a (mapk12a), and pleckstrin homology domain containing family M (plekhm2) were significantly enriched in the network analysis (Table 2 and Fig. 6 G-J). Moreover, it is important to highlight that the protein-protein interaction network is predicted based on the functional interaction of zebrafish proteins. Hence, this network analysis warrants further detailed investigation.
In conclusion, in this experiment, we determined the alternative splicing profile in the posterior kidney of brown trout during PKD. This is the first study in fish to explore the alternative splicing profile of the host during parasitic infection. Interestingly, genes involved in ion homeostasis were significantly enriched among the DEAS genes. Alternative splicing may have a vital role in maintaining the homeostasis of brown trout during the parasite development in the kidney. This study provides new biological insights into the alternative splicing pattern during PKD pathogenesis in brown trout. Of late, the importance and role of alternative splicing during host-pathogen interaction is being slowly recognized for developing therapeutic strategies. Hence, the insights gained from the present study can help to develop treatment measures for PKD.
Experiment and transcriptome data
This is a continuation of our previous work, where brown trout experimental infection and RNA-seq data generation are explained in detail . Briefly, specific pathogen-free brown trout were exposed with T. bryosalmonae spores. The posterior kidneys from the exposed and non-exposed control groups were collected at 2, 4, 6, 8, 10, 12, and 17 wpe. The posterior kidney samples were fixed in RNAlater (Sigma, Steinheim, Germany) and 10% neutral-buffered formalin (NBF) for further studies. RNA extraction was done using RNeasy Mini Kit (Qiagen, Hilden, Germany) and 12 cDNA libraries were prepared from the posterior kidney of T. bryosalmonae exposed brown trout (exposed group, n = 6) and non-exposed control brown trout (control group, n = 6) at 12 wpe using TruSeq RNA Library Prep Kit v2 (Illumina, San Diego, CA, USA). Subsequently, the sequencing of cDNA libraries (100-bp single reads) was performed in Illumina Hi-seq 2500 platform. The raw reads obtained were deposited in the NCBI database under the bioproject ID PRJNA542491. The animal experiment was conducted in accordance with the regulations and guidelines of the §26 Austrian Law for Animal Experiments, Tierversuchsgesetz 2012 with permit number BMWFW-68.205/0181-WF/V/3b/2017. All methods of this study are reported in accordance to the ARRIVE guidelines for animal research .
Mapping and alignment of RNA-seq reads
The NGS data was processed in the Galaxy server (https://usegalaxy.org/) . The raw reads were assessed for quality using the tool FastQC (version 0.11.8). Adapter removal and quality trimming were done using Trimmomatic . The trimmed reads were mapped to the brown trout genome (fSalTru1.1., ftp://ftp.ensembl.org/pub/release-100/fasta/salmo_trutta/dna/Salmo_trutta.fSalTru1.1.dna_sm.toplevel.fa.gz) obtained from Ensembl database using the splice aware mapping tool TopHat2 (Galaxy Version 2.1.1) . Following this Cufflinks 126.96.36.199 tool was used to assemble the mapped reads to the transcripts .
Alternative splicing analysis
The DEAS events between T. bryosalmonae exposed and control kidney were identified using rMATS turbo version 4.1.0 tool . BAM files obtained from TopHat2 were provided as input with default settings in rMATS. Significant DEAS events (delta PSI = 5%, FDR P-value < 0.05) were filtered using Maser package version 1.10.0 in R program. Chromosome-wise distribution of significant DEAS events in the brown trout genome during PKD was visualized using circlize package version 0.4.12 in R program. R programming tool version 4.1.0 was used for the data analysis.
Functional enrichment analysis
Functional enrichment analysis of DEAS events was done in g:Profiler web server . Initially the Ensembl gene identifiers of significant DEAS events of brown trout were converted to the corresponding orthologous Ensembl gene identifiers of the model organism zebrafish (Danio rerio) using g:Convert option. Subsequently, the converted gene identifiers were provided as input in g:GOSt tool for functional enrichment analysis (adjusted P-value < 0.05) for biological process, molecular function, and pathways.
Protein–protein interaction analysis of DEAS events
The protein–protein interaction network of DEAS events was constructed using STRING software version 11.0.b . The converted gene identifiers of DEAS events were given as queries against zebrafish in active interactive sources of textmining, experiments, databases, co-expression, neighbourhood, gene fusion and co-occurrence with an interaction confidence score of 0.15.
RNA extraction and reverse transcription PCR
RT-PCR was used to validate six representative DEAS genes (prkcbp1l, baz2ba, rap1gds1, pik3ap1, rabgef1, and mef2d). Briefly, total RNA was extracted from the posterior kidney of brown trout (exposed and control n = 6) with an on-column DNase I digestion step using RNeasy Mini Kit (Qiagen, Hilden, Germany). Quantification and quality check of the extracted RNA was done using NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Vienna, Austria). Subsequently, one μg of total RNA per sample was used as a template for cDNA synthesis using the iScript cDNA synthesis kit (Bio-Rad, Hercules, USA). For designing primers initially, the alternative spliced transcripts and the spliced sites were visualized in the brown trout genome using integrated genome browser (IGB) version 9.1.8 . Subsequently, the nucleic acid sequences were mined from IGB and used as input in the Primer3Plus tool to design specific primers . The details of the primers and their binding position are provided in Supplementary File S5 and Fig. 7 respectively. The PCR reaction contained 5.0 μL of 2x DreamTaq Green PCR master mix (Thermo Scientific), 0.5 μL of each primer (10 pmol/μL), 1.0 μL DEPC-treated sterile distilled water (Carl Roth, Karlsruhe, Germany), and 3.0 μL of 1:20-fold diluted cDNA. The PCR reaction was performed in C1000 Touch Thermal Cycler (Bio-Rad, München, Germany) under the conditions of initial denaturation at 95 °C for 5 min, 35 cycles of denaturation (95 °C for 30 s), annealing (60 °C for 30 s) and extension (72 °C for 30 s), and the final extension was done for 10 min at 72 °C. The PCR products were subjected to electrophorese using 2% agarose gel containing Gel Red Nucleic Acid Gel Stain (Thermo Scientific) and visualized using a Bio Imaging System (Bio-Rad, München, Germany). The PCR products were excised from the gel, purified, and sequenced from both 5′ and 3′ directions. Subsequently, the sequences were aligned to the corresponding DEAS transcript and the structures of the splice variants were validated.
Availability of data and materials
All raw RNA-seq reads are available in the NCBI Short Read Archive (SRA) database under the Bioproject ID PRJNA542491.
Alternative 3′ splice site
Alternative 5′ splice site
S-adenosylhomocysteine hydrolase-like protein 1
ATPase Na+/K+ transporting subunit alpha 1a tandem duplicate 1
ATPase Na+/K+ transporting subunit alpha 1a tandem duplicate 5
ATPase Na+/K+ transporting subunit alpha 3a
Bromodomain adjacent to zinc finger domain 2Ba
- Ca2+ :
Differently expressed alternatively spliced
Filamin A alpha
- K+ :
Integrated genome browser
- map k12a:
Mitogen-activated protein kinase 12a
Myocyte enhancer factor 2d
Mutually exclusive exons
MYD88 innate immune signal transduction adaptor
- Na+ :
Phosphoinositide-3-kinase adaptor protein 1
Proliferative kidney disease
Pleckstrin homology domain containing family M
Protein kinase C binding protein 1 – like
Premature termination codons
RAB guanine nucleotide exchange factor 1
RAP1 GTP-GDP dissociation stimulator 1
Reverse transcription PCR
Ubiquitin-specific peptidase 25
Weeks post exposure
Sudhagar A, Kumar G, El-Matbouli M. The malacosporean myxozoan parasite Tetracapsuloides bryosalmonae: a threat to wild salmonids. Pathogens. 2020;9(1):16.
Canning EU, Curry A, Feist SW, Longshaw M, Okamura B. Tetracapsula bryosalmonae n.sp. for PKX organism, the cause of PKD in salmonid fish. Bull Eur ass. Fish Pathol. 1999;19(5):203–6.
Morris DJ, Adams A. Transmission of Tetracapsuloides bryosalmonae (Myxozoa: Malacosporea), the causative organism of salmonid proliferative kidney disease, to the freshwater bryozoan Fredericella sultana. Parasitology. 2006;133:701–9.
Grabner DS, El-Matbouli M. Transmission of Tetracapsuloides bryosalmonae (Myxozoa: Malacosporea) to Fredericella sultana (Bryozoa: Phylactolaemata) by various fish species. Dis Aquat Org. 2008;79(2):133–9.
Bettge K, Wahli T, Segner H, Schmidt-Posthaus H. Proliferative kidney disease in rainbow trout: time- and temperature-related renal pathology and parasite distribution. Dis Aquat Org. 2009;83(1):67–76.
Waldner K, Bechter T, Auer S, Borgwardt F, El-Matbouli M, Unfer G. A brown trout (Salmo trutta) population faces devastating consequences due to proliferative kidney disease and temperature increase: a case study from Austria. Ecol Freshw Fish. 2020;29(3):465–76.
Sage JL. Economic contributions of the Yellowstone River to Park County, Montana. Insti Tourism Recreat Res Public. 2016;346:1–8.
Grabner DS, El-Matbouli M. Tetracapsuloides bryosalmonae (Myxozoa: Malacosporea) portal of entry into the fish host. Dis Aquat Org. 2010;90(3):197–206.
Kumar G, Abd-Elfattah A, Saleh M, El-Matbouli M. Fate of Tetracapsuloides bryosalmonae (Myxozoa) after infection of brown trout Salmo trutta and rainbow trout Oncorhynchus mykiss. Dis Aquat Org. 2013;107(1):9–18.
Soliman H, Kumar G, El-Matbouli M. Tetracapsuloides bryosalmonae persists in brown trout Salmo trutta for five years post exposure. Dis Aquat Org. 2018;127(2):151–6.
Okamura B, Hartikainen H, Schmidt-Posthaus H, Wahli T. Life cycle complexity, environmental change and the emerging status of salmonid proliferative kidney disease. Freshw Biol. 2011;56(4):735–53.
Bailey C, Segner H, Casanova-Nakayama A, Wahli T. Who needs the hotspot? The effect of temperature on the fish host immune response to Tetracapsuloides bryosalmonae the causative agent of proliferative kidney disease. Fish Shellfish Immunol. 2017;63:424–37.
Holland JW, Gould CRW, Jones CS, Noble LR, Secombes CJ. The expression of immune-regulatory genes in rainbow trout, Oncorhynchus mykiss, during a natural outbreak of proliferative kidney disease (PKD). Parasitology. 2003;126(7):S95–102.
Gorgoglione B, Wang T, Secombes CJ, Holland JW. Immune gene expression profiling of proliferative kidney disease in rainbow trout Oncorhynchus mykiss reveals a dominance of anti-inflammatory, antibody and Th cell-like activities. Vet Res. 2013;44:55.
Bailey C, Segner H, Wahli T. What goes around comes around: an investigation of resistance to proliferative kidney disease in rainbow trout Oncorhynchus mykiss (Walbaum) following experimental re-exposure. J Fish Dis. 2017;40(11):1599–612.
Abos B, Estensoro I, Perdiguero P, Faber M, Hu Y, Rosales PD, et al. Dysregulation of B cell activity during proliferative kidney disease in rainbow trout. Front Immunol. 2018;9:1203.
Kumar G, Abd-Elfattah A, El-Matbouli M. Identification of differentially expressed genes of brown trout (Salmo trutta) and rainbow trout (Oncorhynchus mykiss) in response to Tetracapsuloides bryosalmonae (Myxozoa). Parasitol Res. 2015;114(3):929–39.
Kotob MH, Kumar G, Saleh M, Gorgoglione B, Abdelzaher M. Differential modulation of host immune genes in the kidney and cranium of the rainbow trout (Oncorhynchus mykiss) in response to Tetracapsuloides bryosalmonae and Myxobolus cerebralis co-infections. Parasit Vectors. 2018;11:326.
Bailey C, Strepparava N, Wahli T, Segner H. Exploring the immune response, tolerance and resistance in proliferative kidney disease of salmonids. Dev Comp Immunol. 2019;90:165–75.
Kumar G, Abd-elfattah A, El-matbouli M. Differential modulation of host genes in the kidney of brown trout Salmo trutta during sporogenesis of Tetracapsuloides bryosalmonae (Myxozoa). Vet Res. 2014;45:101.
Sudhagar A, Ertl R, Kumar G, El MM. Transcriptome profiling of posterior kidney of brown trout, Salmo trutta, during proliferative kidney disease. Parasit Vectors. 2019;12:569.
Sudhagar A, El-matbouli M, Kumar G. Identification and expression profiling of toll-like receptors of brown trout (Salmo trutta) during proliferative kidney disease. Int J Mol Sci. 2020;21(11):3755.
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456:470–6.
Graveley BR. Alternative splicing: increasing diversity in the proteomic world. Trends Genet. 2001;17(2):100–7.
Lejeune F, Maquat LE. Mechanistic links between nonsense-mediated mRNA decay and pre-mRNA splicing in mammalian cells. Curr Opin Cell Biol. 2005;17(3):309–15.
Lewis BP, Green RE, Brenner SE. Evidence for the widespread coupling of alternative splicing and nonsense-mediated mRNA decay in humans. Proc Natl Acad Sci U S A. 2003;100(1):189–92.
Kelemen O, Convertini P, Zhang Z, Wen Y, Shen M, Falaleeva M, et al. Function of alternative splicing. Gene. 2013;514(1):1–30.
Chauhan K, Kalam H, Dutt R, Kumar D. RNA splicing: a new paradigm in host–pathogen interactions. J Mol Biol. 2019;431:1565–75.
Martinez NM, Lynch KW. Control of alternative splicing in immune responses: many regulators, many predictions, much still to learn. Immunol Rev. 2013;253(1):216–36.
Lynch KW. Consequences of regulated pre-mRNA splicing in the immune system. Nat Rev Immunol. 2004;4:931–40.
Chang MX, Zhang J. Alternative pre-mRNA splicing in mammals and teleost fish: a effective strategy for the regulation of immune responses against pathogen infection. Int J Mol Sci. 2017;18(7):1530.
Yang Y, Huang JM, Ju ZH, Li QL, Zhou L, Li RL, et al. Increased expression of a novel splice variant of the complement component 4 (C4A) gene in mastitis-infected dairy cattle. Genet Mol Res. 2012;11(3):2909–16.
Tan S, Wang W, Tian C, Niu D, Zhou T, Jin Y, et al. Heat stress induced alternative splicing in catfish as determined by transcriptome analysis. Comp Biochem Physiol Part D Genomics Proteomics. 2019;29:166–72.
Tan S, Wang W, Zhong X, Tian C, Niu D, Bao L, et al. Increased alternative splicing as a host response to Edwardsiella ictaluri infection in catfish. Mar Biotechnol. 2018;20:729–38.
Hansen T, Fjelldal PG, Lien S, Smith M, Corton C, Oliver K, et al. The genome sequence of the brown trout, Salmo trutta Linnaeus 1758. Wellcome Open Res. 2021;6:108.
Sudhagar A, Kumar G, El-Matbouli M. Transcriptome analysis based on RNA-Seq in understanding pathogenic mechanisms of diseases and the immune system of fish: a comprehensive review. Int J Mol Sci. 2018;19(1):245.
Shivam S, El-Matbouli M, Kumar G. Development of fish parasite vaccines in the OMICs era: Progress and opportunities. Vaccines. 2021;9(2):179.
Kumar G, Ertl R, Nilsen F, Bartholomew JL, El-Matbouli M. Data of de novo transcriptome assembly of the myxozoan parasite Tetracapsuloides bryosalmonae. Data Br. 2021;35:106831.
Faber M, Shaw S, Yoon S, de Paiva AE, Wang B, Qi Z, et al. Comparative transcriptomics and host-specific parasite gene expression profiles inform on drivers of proliferative kidney disease. Sci Rep. 2021;11:2149.
Kumar G, Ertl R, Bartholomew JL, El-Matbouli M. First transcriptome analysis of bryozoan Fredericella sultana, the primary host of myxozoan parasite Tetracapsuloides bryosalmonae. PeerJ. 2020;8:e9027.
Bailey C, von Siebenthal EW, Rehberger K, Segner H. Transcriptomic analysis of the impacts of ethinylestradiol (EE2) and its consequences for proliferative kidney disease outcome in rainbow trout (Oncorhynchus mykiss). Comp Biochem Physiol Part C Toxicol Pharmacol. 2019;222:31–48.
Kumar G, Ertl R, Bartholomew JL, El-Matbouli M. Transcriptome analysis elucidates the key responses of bryozoan Fredericella sultana during the development of Tetracapsuloides bryosalmonae (Myxozoa). Int J Mol Sci. 2020;21(16):5910.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Lu J, Peatman E, Wang W, Yang Q, Abernathy J, Wang S, et al. Alternative splicing in teleost fish genomes: same-species and cross-species analysis and comparisons. Mol Gen Genomics. 2010;283:531–9.
Tan S, Wang W, Tian C, Niu D, Zhou T, Yang Y, et al. Post-transcriptional regulation through alternative splicing after infection with Flavobacterium columnare in channel catfish (Ictalurus punctatus). Fish Shellfish Immunol. 2019;91:188–93.
Wan Q, Su J. Transcriptome analysis provides insights into the regulatory function of alternative splicing in antiviral immunity in grass carp (Ctenopharyngodon idella). Sci Rep. 2015;5:12946.
Takvam M, Wood CM, Kryvi H, Nilsen TO. Ion transporters and osmoregulation in the kidney of teleost fishes as a function of salinity. Front Physiol. 2021;12:664588.
Beyenbach KW. Kidneys sans glomeruli. Am J Physiol Physiol. 2004;286(5):F811–27.
Cooper BJ, Key B, Carter A, Angel NZ, Hart DNJ, Kato M. Suppression and overexpression of adenosylhomocysteine hydrolase-like protein 1 (AHCYL1) influences zebrafish embryo development: a possible role for AHCYL1 in inositol phospholipid signaling. J Biol Chem. 2006;281(32):22471–84.
Armesto P, Campinho MA, Rodríguez-Rúa A, Cousin X, Power DM, Manchado M, et al. Molecular characterization and transcriptional regulation of the Na+/K+ ATPase α subunit isoforms during development and salinity challenge in a teleost fish, the Senegalese sole (Solea senegalensis). Comp Biochem Physiol Part B Biochem Mol Biol. 2014;175:23–38.
Mobasheri A, Avila J, Cózar-Castellano I, Brownleader MD, Trevan M, Francis MJ, et al. Na+, K+-ATPase isozyme diversity; comparative biochemistry and physiological implications of novel functional interactions. Biosci Rep. 2000;20(2):51–91.
Lai F, Madan N, Ye Q, Duan Q, Li Z, Wang S, et al. Identification of a mutant α1 Na/K-ATPase that pumps but is defective in signal transduction. J Biol Chem. 2013;288(19):13295–304.
Blanco G, Koster JC, Mercer RW. The alpha subunit of the Na,K-ATPase specifically and stably associates into oligomers. Proc Natl Acad Sci U S A. 1994;91(18):8542–6.
Pott A, Bock S, Berger IM, Frese K, Dahme T, Keßler M, et al. Mutation of the Na+/K+ − ATPase Atp1a1a.1 causes QT interval prolongation and bradycardia in zebrafish. J Mol Cell Cardiol. 2018;120:42–52.
Doğanli C, Beck HC, Ribera AB, Oxvig C, Lykke-Hartmann K. α3Na+/K+-ATPase deficiency causes brain ventricle dilation and abrupt embryonic motility in zebrafish. J Biol Chem. 2013;288(13):8862–74.
Gamba G. Alternative splicing and diversity of renal transporters. Am J Physiol Physiol. 2001;281:F781–94.
Kilkenny C, Browne WJ, Cuthill IC, Emerson M, Altman DG. Improving bioscience research reporting: the ARRIVE guidelines for reporting animal research. PLoS Biol. 2010;8(6):e1000412.
Afgan E, Baker D, Batut B, van den Beek M, Bouvier D, Cech M, et al. The galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018;46(W1):W537–44.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
Shen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, et al. rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci U S A. 2014;111(51):E5593–601.
Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, et al. G:profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 2019;47(W1):W191–8.
Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49(D1):D605–12.
Freese NH, Norris DC, Loraine AE. Integrated genome browser: visual analytics platform for genomics. Bioinformatics. 2016;32:2089–95.
Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007;35:W71–4.
The first author acknowledges the Indian Council of Agriculture Research, New Delhi, India for providing PhD scholarship under Netaji Subhas—ICAR International Fellowship 2016–2017.
This study was funded by the Austrian Science Fund (FWF) Project No. P 30981-B32 to Gokhlesh Kumar.
Ethics approval and consent to participate
This study was approved by the institutional ethics committee of the University of Veterinary Medicine Vienna and the national authority, according to §26 of the Austrian Law for Animal Experiments, Tierversuchsgesetz 2012 under approval number BMWFW-68.205/0181-WF/V/3b/2017. All methods of this study are reported in accordance to the ARRIVE guidelines for animal research .
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Sudhagar, A., El-Matbouli, M. & Kumar, G. Genome-wide alternative splicing profile in the posterior kidney of brown trout (Salmo trutta) during proliferative kidney disease. BMC Genomics 23, 446 (2022). https://doi.org/10.1186/s12864-022-08685-4
- Post-transcriptional modification
- Tetracapsuloides bryosalmonae