Changes in gene expression as biochemical adaptations to environmental change: a tribute to Peter Hochachka

Adaptive Differences in Gene Expression in Farm-Impacted Seedbeds of the Native Blue Mussel Mytilus chilensis

The study of adaptive population differences is relevant for evolutionary biology, as it evidences the power of selective local forces relative to gene flow in maintaining adaptive phenotypes and their underlying genetic determinants. However, human-mediated hybridization through habitat translocations, a common and recurrent aquaculture practice where hybrids could eventually replace local genotypes, risk populations' ability to cope with perturbations. The endemic marine mussel Mytilus chilensis supports a booming farming industry in the inner sea of Chiloé Island, southern Chile, which entirely relies on artificially collected seeds from natural beds that are translocated to ecologically different fattening centers. A matter of concern is how farm-impacted seedbeds will potentially cope with environmental shifts and anthropogenic perturbations. This study provides the first de novo transcriptome of M. chilensis; assembled from tissue samples of mantles and gills of individuals collected in ecologically different farm-impacted seedbeds, Cochamó (41°S) and Yaldad (43°S). Both locations and tissue samples differentially expressed transcripts (DETs) in candidate adaptive genes controlling multiple fitness traits, involved with metabolism, genetic and environmental information processing, and cellular processes. From 189,743 consensus contigs assembled: 1,716 (Bonferroni p value ≤ 0.05) were DETs detected in different tissues of samples from different locations, 210 of them (fold change ≥ | 100|) in the same tissue of samples from a different location, and 665 (fold change ≥ | 4|) regardless of the tissue in samples from a different location. Site-specific DETs in Cochamó (169) and Yaldad (150) in candidate genes controlling tolerance to temperature and salinity shifts, and biomineralization exhibit a high number of nucleotide genetic variants with regular occurrence (frequency > 99%). This novel M. chilensis transcriptome should help assessing and monitoring the impact of translocations in wild and farm-impacted mussel beds in Chiloé Island. At the same time, it would help designing effective managing practices for conservation, and translocation traceability.

Patterns of alternative splicing in response to cold acclimation in fish

Phenotypic plasticity is an important aspect of an organism's response to environmental change that often requires the modulation of gene expression. These changes in gene expression can be quantitative, as a result of increases or decreases in the amounts of specific transcripts, or qualitative, as a result of the expression of alternative transcripts from the same gene (e.g. via alternative splicing of pre-mRNAs). Although the role of quantitative changes in gene expression in phenotypic plasticity is well known, relatively few studies have examined the role of qualitative changes. Here, we use skeletal muscle RNA-seq data from Atlantic killifish (Fundulus heteroclitus), threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio) to investigate the extent of qualitative changes in gene expression in response to cold acclimation. Fewer genes demonstrated alternative splicing than differential expression as a result of cold acclimation; however, differences in splicing were detected for 426 to 866 genes depending on species, indicating that large numbers of qualitative changes in gene expression are associated with cold acclimation. Many of these alternatively spliced genes were also differentially expressed, and there was functional enrichment for involvement in muscle contraction among the genes demonstrating qualitative changes in response to cold acclimation. Additionally, there was a common group of 29 genes with cold-acclimation-mediated changes in splicing in all three species, suggesting that there may be a set of genes with expression patterns that respond qualitatively to prolonged exposure to cold temperatures across fishes.

Differential survivorship of congeneric ornamental fishes under forecasted climate changes are related to anaerobic potential

Two Amazonian closely related tetras – cardinal Paracheirodon axelrodi and green neon P. simulans – were artificially acclimatized to environmental chambers mimicking future climate change scenarios (mild, moderate and extreme), using a microcosm facility. P. simulans survived (100%) to all scenarios after 30 days exposure, while P. axelrodi presented decreasing survival percentages according to environmental severity. These differences may be the reflection of distinct natural acclimatization to microhabitats between the species, which differ in thermal conditions. Survival responses might be related to differences in relative gene expression of lactate dehydrogenase (Ldh), suggesting that P. axelrodi anaerobic potential is lower or non-existent compared to P. simulans, not tolerating long-term thermal challenges. Accordingly, increases in temperature and in CO₂ levels caused increases in energy demand and resulted in activation of the anaerobic pathway, as demonstrated by the higher enzyme levels measured in head and tail portions of both species. Sustained anaerobic glycolysis is possible when fish live in challenging environments (low oxygen or high temperature). Our results clearly show that P. simulans has a larger scope for survival to higher energy demands due to its increased anaerobic potential compared to P. axelrodi.

Transcriptome response of Atlantic salmon (Salmo salar) to competition with ecologically similar non-native species

Non-native species may be introduced either intentionally or unintentionally, and their impact can range from benign to highly disruptive. Non-native salmonids were introduced into Lake Ontario, Canada, to provide recreational fishing opportunities; however, the establishment of those species has been proposed as a significant barrier to the reintroduction of native Atlantic salmon (Salmo salar) due to intense interspecific competition. In this study, we compared population differences of Atlantic salmon in transcriptome response to interspecific competition. We reared Atlantic salmon from two populations (LaHave River and Sebago Lake) with fish of each of three non-native salmonids (Chinook salmon Oncorhynchus tshawytscha, rainbow trout O. mykiss, and brown trout S. trutta) in artificial streams. We used RNA-seq to assess transcriptome differences between the Atlantic salmon populations and the responses of these populations to the interspecific competition treatments after 10 months of competition in the stream tanks. We found that population differences in gene expression were generally greater than the effects of interspecific competition. Interestingly, we found that the two Atlantic salmon populations exhibited similar responses to interspecific competition based on functional gene ontologies, but the specific genes within those ontologies were different. Our transcriptome analyses suggest that the most stressful competitor (as measured by the highest number of differentially expressed genes) differs between the two study populations. Our transcriptome characterization highlights the importance of source population selection for conservation applications, as organisms with different evolutionary histories can possess different transcriptional responses to the same biotic stressors. The results also indicate that generalized predictions of the response of native species to interactions with introduced species may not be appropriate without incorporating potential population-specific response to introduced species.

Transcriptional mechanisms underlying life-history responses to climate change in the three-spined stickleback

Phenotypic plasticity, the ability of an organism to express different phenotypes depending on the environment, provides an important mechanism by which an animal population can persist under rapid climate change. We experimentally tested both life‐history and transcriptional responses of an ecological model species, the three‐spined stickleback, to warm acclimation at the southern edge of its European range. We explored cross‐environment genetic correlations of key life‐history traits in male sticklebacks exposed to long‐term temperature changes to examine whether the plasticity pattern was variable among genotypes by using a character‐state approach. We also studied gene expression plasticity by analysing both whole‐transcriptome and candidate gene expression in brain and liver. Male sticklebacks that developed under warmer conditions during winter were smaller in size and invested less in nuptial coloration at the beginning of the breeding season, showing similar responses across different genotypes. The lack of genetic variation in life‐history responses may limit any future evolution of the thermal reaction norm in the study population. After long‐term exposure to increased winter temperatures, genes responsible for several metabolic and oxidation–reduction processes were upregulated, and some hormone genes involved in growth and reproduction were downregulated in the brain. In the liver, there was no significantly represented gene ontology by the differentially expressed genes. Since a higher temperature leads to a higher resting metabolic rate, living in warmer environments may incur higher energetic costs for ectotherms to maintain cellular homoeostasis, resulting in negative consequences for life‐history traits. The expression of genes related to metabolism, cellular homoeostasis and regulatory signalling may underlie temperature‐induced changes in life history.

Genetic architecture of gene transcription in two Atlantic salmon (Salmo salar) populations

Gene expression regulation has an important role in short-term acclimation and long-term adaptation to changing environments. However, the genetic architecture of gene expression has received much less attention than that of traditional phenotypic traits. In this study, we used a 5 × 5 full-factorial breeding design within each of two Atlantic salmon (Salmo salar) populations to characterize the genetic architecture of gene transcription. The two populations (LaHave and Sebago) are being used for reintroduction efforts into Lake Ontario, Canada. We used high-throughput quantitative real-time PCR to measure gene transcription levels for 22 genes in muscle tissue of Atlantic salmon fry. We tested for population differences in gene transcription and partitioned the transcription variance into additive genetic, non-additive genetic and maternal effects within each population. Interestingly, average additive genetic effects for gene transcription were smaller than those reported for traditional phenotypic traits in salmonids, suggesting that the evolutionary potential of gene transcription is lower than that of traditional traits. Contrary to expectations for early life stage traits, maternal effects were small. In general, the LaHave population had higher additive genetic effects for gene transcription than the Sebago population had, indicating that the LaHave fish have a higher adaptive potential to respond to the novel selection pressures associated with reintroduction into a novel environment. This study highlights not only the profound variation in gene transcription possible among salmonid populations but also the among-population variation in the underlying genetic architecture of such traits.

Role of genomics and transcriptomics in selection of reintroduction source populations

The use and importance of reintroduction as a conservation tool to return a species to its historical range from which it has been extirpated will increase as climate change and human development accelerate habitat loss and population extinctions. Although the number of reintroduction attempts has increased rapidly over the past 2 decades, the success rate is generally low. As a result of population differences in fitness-related traits and divergent responses to environmental stresses, population performance upon reintroduction is highly variable, and it is generally agreed that selecting an appropriate source population is a critical component of a successful reintroduction. Conservation genomics is an emerging field that addresses long-standing challenges in conservation, and the potential for using novel molecular genetic approaches to inform and improve conservation efforts is high. Because the successful establishment and persistence of reintroduced populations is highly dependent on the functional genetic variation and environmental stress tolerance of the source population, we propose the application of conservation genomics and transcriptomics to guide reintroduction practices. Specifically, we propose using genome-wide functional loci to estimate genetic variation of source populations. This estimate can then be used to predict the potential for adaptation. We also propose using transcriptional profiling to measure the expression response of fitness-related genes to environmental stresses as a proxy for acclimation (tolerance) capacity. Appropriate application of conservation genomics and transcriptomics has the potential to dramatically enhance reintroduction success in a time of rapidly declining biodiversity and accelerating environmental change.

Transcriptional response of stress-regulated genes to industrial effluent exposure in the cockle Cerastoderma glaucum

This study assessed the responses of molecular biomarkers and heavy metal levels in Cerastoderma glaucum exposed for 1 week to two industrial effluents (1%) discharged into the Tunisian coastal area, F1 and F2, produced by different units of production of a phosphate treatment plant. A significant uptake of metals (Cd, Cu, Zn, and Ni) was observed in exposed cockles compared to controls, with an uptake higher for F1 than for F2. A decrease in LT50 (stress on stress test) was also observed after an exposure to the effluent F1. Treatments resulted in different patterns of messenger RNA (mRNA) expression of the different genes tested in this report. Gene transcription monitoring performed on seven genes potentially involved in the tolerance to metal exposure showed that for both exposures, mechanisms are rapidly and synchronically settled down to prevent damage to cellular components, by (1) handling and exporting out metal ions through the up-regulation of ATP-binding cassette xenobiotic transporter (ABCB1) and metallothionein (MT), (2) increasing the mRNA expression of antioxidant enzymes (catalase (CAT), superoxide dismutases, CuZnSOD and MnSOD), (3) protecting and/or repairing proteins through the expression of heat shock protein 70 (HSP70) mRNAs, and (4) increasing ATP production (through the up-regulation of cytochrome c oxidase 1 (CO1)) to provide energy for cells to tolerate stress exposure. The tools developed may be useful both for future control strategies and for the use of the cockle C. glaucum as a sentinel species.

Transcriptional response of stress-regulated genes to cadmium exposure in the cockle Cerastoderma glaucum from the gulf of Gabès area (Tunisia)

This study investigates cadmium effects on key messenger RNA (mRNA) expression (MT, MnSOD, CuZnSOD, CAT, ABCB1, HSP70, and CO1) by qPCR in the cockle Cerastoderma glaucum after chronic exposure to two high but environmentally relevant concentrations of CdCl2 (50 μg/L and 5 mg/L) for 12 h to 18 days. Cd accumulation measured in cockles' tissues is significantly higher in both treatment conditions compared to controls and in a dose-dependent manner. Stress on stress tests performed at different times of the experiment clearly demonstrated that exposure to both concentrations of Cd significantly affects cockle survival time in air. Important changes in gene transcription were also highlighted. In particular, MT, HSP70, CAT, and CuZnSOD seem to be relevant biomarkers of Cd exposure because (1) their mRNA levels increase upon exposure and (2) they are highly correlated to Cd accumulation in tissues. Results may be useful for control strategies and for the use of cockles as sentinel organisms.

Development of intestinal ion-transporting mechanisms during smoltification and seawater acclimation in Atlantic salmon Salmo salar

This study investigated the expression of ion transporters involved in intestinal fluid absorption and presents evidence for developmental changes in abundance and tissue distribution of these transporters during smoltification and seawater (SW) acclimation of Atlantic salmon Salmo salar. Emphasis was placed on Na(+) , K(+) -ATPase (NKA) and Na(+) , K(+) , Cl(-) co-transporter (NKCC) isoforms, at both transcriptional and protein levels, together with transcription of chloride channel genes. The nka α1c was the dominant isoform at the transcript level in both proximal and distal intestines; also, it was the most abundant isoform expressed in the basolateral membrane of enterocytes in the proximal intestine. This isoform was also abundantly expressed in the distal intestine in the lower part of the mucosal folds. The protein expression of intestinal Nkaα1c increased during smoltification. Immunostaining was localized to the basal membrane of the enterocytes in freshwater (FW) fish, and re-distributed to a lateral position after SW entry. Two other Nka isoforms, α1a and α1b, were expressed in the intestine but were not regulated to the same extent during smoltification and subsequent SW transfer. Their localization in the intestinal wall indicates a house-keeping function in excitatory tissues. The absorptive form of the NKCC-like isoform (sub-apically located NKCC2 and/or Na(+) , Cl(-) co-transporter) increased during smoltification and further after SW transfer. The cellular distribution changed from a diffuse expression in the sub-apical regions during smoltification to clustering of the transporters closer to the apical membrane after entry to SW. Furthermore, transcript abundance indicates that the mechanisms necessary for exit of chloride ions across the basolateral membrane and into the lateral intercellular space are present in the form of one or more of three different chloride channels: cystic fibrosis transmembrane conductance regulator I and II and chloride channel 3.

Effect of heat stress and recovery on viability, oxidative damage, and heat shock protein expression in hepatic cells of grass carp (Ctenopharyngodon idellus)

In this study, we investigated the effects of hyperthermia and recovery on cell viability, lactate dehydrogenase (LDH) activity, superoxide dismutase (SOD) activity, malondialdehyde (MDA), total antioxidant capacity (T-AOC), and heat shock protein (HSP60, 70, and 90) mRNA expression in the hepatic cells of the grass carp, Ctenopharyngodon idellus. Triplicate groups of cultured cells were exposed to 30, 32, or 34 °C for 0.5 h and then immediately incubated at 27 °C in 5 % CO2 for 6, 12, 24, or 48 h. Hyperthermia stress greatly reduced cell viability and increased LDH release. Cell damage declined after recovery. Hyperthermia stress increased the lipid peroxide levels and reduced the antioxidant capacity (e.g., reduced SOD and T-AOC) of the cells. However, oxidative damage declined as the recovery period increased, and the levels of MDA, SOD, and T-AOC were restored. After cells were exposed to 32 °C, the expression of HSP60 after recovery for 1, 2, and 4 h (P < 0.05), the expression of HSP70 after recovery for 0.5 and 1 h (P < 0.01), and the expression of HSP90 throughout recovery were significantly higher (P < 0.01) than the prestress levels. During the recovery period, the variations in HSP gene expression reflected the transition period from a state of cellular growth to one of the cellular repairs. In conclusion, hyperthermia depresses cell viability, induces oxidative damage, and increases HSP expression, which plays an important role during hyperthermic stress in grass carp hepatic cells.

Tilapia (Oreochromis mossambicus) brain cells respond to hyperosmotic challenge by inducing myo-inositol biosynthesis

This study aimed to determine the regulation of the de novo myo-inositol biosynthetic (MIB) pathway in Mozambique tilapia (Oreochromis mossambicus) brain following acute (25 ppt) and chronic (30, 60 and 90 ppt) salinity acclimations. The MIB pathway plays an important role in accumulating the compatible osmolyte, myo-inositol, in cells in response to hyperosmotic challenge and consists of two enzymes, myo-inositol phosphate synthase and inositol monophosphatase. In tilapia brain, MIB enzyme transcriptional regulation was found to robustly increase in a time (acute acclimation) or dose (chronic acclimation) dependent manner. Blood plasma osmolality and Na(+) and Cl(-) concentrations were also measured and significantly increased in response to both acute and chronic salinity challenges. Interestingly, highly significant positive correlations were found between MIB enzyme mRNA and blood plasma osmolality in both acute and chronic salinity acclimations. Additionally, a mass spectrometry assay was established and used to quantify total myo-inositol concentration in tilapia brain, which closely mirrored the hyperosmotic MIB pathway induction. Thus, myo-inositol is a major compatible osmolyte that is accumulated in brain cells when exposed to acute and chronic hyperosmotic challenge. These data show that the MIB pathway is highly induced in response to environmental salinity challenge in tilapia brain and that this induction is likely prompted by increases in blood plasma osmolality. Because the MIB pathway uses glucose-6-phosphate as a substrate and large amounts of myo-inositol are being synthesized, our data also illustrate that the MIB pathway likely contributes to the high energetic demand posed by salinity challenge.

Role of hypoxia-inducible factor α in response to hypoxia and heat shock in the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas inhabits the intertidal zone and shows tolerance to stress conditions such as hypoxia and heat shock. Although some information is available about the genes expressed in response to hypoxia, little is known about the molecular mechanism of the regulation of their expression in mollusks, including the Pacific oyster. Hypoxia-inducible factor 1α (HIF-1α) is a master regulator of hypoxia-responsive transcription. In this study, we cloned HIF-α from the oyster and investigated its response to unique stress conditions, including air exposure, for the first time in mollusks. The cDNA of oyster Hif-α is 3,182 bp long, of which 2,094 bp encodes a protein of 698 amino acid residues. Northern and Western blot analysis showed that expression of oyster HIF-α mRNA and protein were induced by air exposure, and that expression was induced periodically during air exposure. In addition, induction of Hif-α mRNA increased by a maximum 8.0-fold by heat shock. Under heat shock at 35°C (lethal temperature for the oyster), however, it was induced later than at 30°C. After recovery from hypoxia and/or heat shock, Hif-α mRNA also upregulated. These data suggest that the oyster has a strategy to induce Hif-α mRNA in order to survive hypoxia and heat shock, and that HIF signaling is necessary for recovery from stress.

Novel isoforms of heat shock transcription factor 1 are induced by hypoxia in the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas inhabits the intertidal zone and shows tolerance to various stress conditions such as hypoxia and heat shock. However, little is known about the cellular mechanism of responses to these stresses. Heat shock transcription factor 1 (HSF1) regulates the transcription of several genes, including heat shock proteins (HSPs). In this study, we cloned HSF1 from the oyster and investigated its response to air-exposure. The cDNA of oyster Hsf1 contains 2,931 bp, of which 1,389 bp encode a protein of 463 amino acid residues. Moreover, we found that the oyster has seven novel alternatively spliced isoforms, Hsf1b-h, consisting of insertion A with 48 bp, insertion B with 42 bp and/or insertion C with 42 bp. We determined the sequences of oyster genomic DNA containing Hsf1 insertions A, B and C. The results indicated that eight isoforms of Hsf1 are generated from a single Hsf1 gene by alternative splicing without frameshifting. Real-time PCR analysis showed that Hsf1a is expressed constitutively, and the expression of Hsf1b-h and Hsp70 mRNA is induced by air exposure and/or hypoxia. In addition, we found that 11 putative hypoxia response elements, which are hypoxia-inducible factor 1 (HIF-1) binding sequences, are located in the Hsf1 promoter region. These data suggest that the oyster has an HIF-HSF pathway in which HSPs are induced in an HIF-dependent manner, and that it also has a novel mechanism of alternative splicing of Hsf1 in response to hypoxia.

Population-specific gene expression responses to hybridization between farm and wild Atlantic salmon

Because of intrinsic differences in their genetic architectures, wild populations invaded by domesticated individuals could experience population-specific consequences following introgression by genetic material of domesticated origin. Expression levels of 16 000 transcripts were quantified by microarrays in liver tissue from farm, wild, and farm-wild backcross (i.e. F1 farm-wild hybrid × wild; total n = 50) Atlantic salmon (Salmo salar) raised under common environmental conditions. The wild populations and farm strain originated from three North American rivers in eastern Canada (Stewiacke, Tusket, and Saint John rivers, respectively). Analysis of variance revealed 177 transcripts with different expression levels among the five strains compared. Five times more of these transcripts were differentiated between farmed parents and Tusket backcrosses (n = 53) than between Stewiacke backcrosses and their farmed parents (n = 11). Altered biological processes in backcrosses also differed between populations both in number and in the type of processes impacted (metabolism vs immunity). Over-dominant gene expression regulation in backcrosses varied considerably between populations (23% in Stewiacke vs 44% in Tusket). Hence, the consequences of introgression of farm genetic material on gene expression depended on population-specific genetic architectures. These results support the need to evaluate impacts of farm-wild genetic interactions at the population scale.

Animal performance and stress: responses and tolerance limits at different levels of biological organisation

Recent advances in molecular biology and the use of DNA microarrays for gene expression profiling are providing new insights into the animal stress response, particularly the effects of stress on gene regulation. However, interpretation of the complex transcriptional changes that occur during stress still poses many challenges because the relationship between changes at the transcriptional level and other levels of biological organisation is not well understood. To confront these challenges, a conceptual model linking physiological and transcriptional responses to stress would be helpful. Here, we provide the basis for one such model by synthesising data from organismal, endocrine, cellular, molecular, and genomic studies. We show using available examples from ectothermic vertebrates that reduced oxygen levels and oxidative stress are common to many stress conditions and that the responses to different types of stress, such as environmental, handling and confinement stress, often converge at the challenge of dealing with oxygen imbalance and oxidative stress. As a result, a common set of stress responses exists that is largely independent of the type of stressor applied. These common responses include the repair of DNA and protein damage, cell cycle arrest or apoptosis, changes in cellular metabolism that reflect the transition from a state of cellular growth to one of cellular repair, the release of stress hormones, changes in mitochondrial densities and properties, changes in oxygen transport capacities and changes in cardio-respiratory function. Changes at the transcriptional level recapitulate these common responses, with many stress-responsive genes functioning in cell cycle control, regulation of transcription, protein turnover, metabolism, and cellular repair. These common transcriptional responses to stress appear coordinated by only a limited number of stress-inducible and redox-sensitive transcription factors and signal transduction pathways, such as the immediate early genes c-fos and c-jun, the transcription factors NFkappaB and HIF-1alpha, and the JNK and p38 kinase signalling pathways. As an example of environmental stress responses, we present temperature response curves at organismal, cellular and molecular levels. Acclimation and physiological adjustments that can shift the threshold temperatures for the onset of these responses are discussed and include, for example, adjustments of the oxygen delivery system, the heat shock response, cellular repair system, and transcriptome. Ultimately, however, an organism's ability to cope with environmental change is largely determined by its ability to maintain aerobic scope and to prevent loss in performance. These systemic constraints can determine an organism's long-term survival well before cellular and molecular functions are disturbed. The conceptual model we propose here discusses some of the crosslinks between responses at different levels of biological organisation and the central role of oxygen balance and oxidative stress in eliciting these responses with the aim to help the interpretation of environmental genomic data in the context of organismal function and performance.

Expression and water calcium dependence of calcium transporter isoforms in zebrafish gill mitochondrion-rich cells

Background

Freshwater fish absorb Ca²⁺ predominantly from ambient water, and more than 97% of Ca²⁺ uptake is achieved by active transport through gill mitochondrion-rich (MR) cells. In the current model for Ca²⁺ uptake in gill MR cells, Ca²⁺ passively enters the cytosol via the epithelium Ca²⁺ channel (ECaC), and then is extruded into the plasma through the basolateral Na⁺/Ca²⁺ exchanger (NCX) and plasma membrane Ca²⁺-ATPase (PMCA). However, no convincing molecular or cellular evidence has been available to support the role of specific PMCA and/or NCX isoforms in this model. Zebrafish (Danio rerio) is a good model for analyzing isoforms of a gene because of the plentiful genomic databases and expression sequence tag (EST) data.

Results

Using a strategy of BLAST from the zebrafish genome database (Sanger Institute), 6 isoforms of PMCAs (PMCA1a, PMCA1b, PMCA2, PMCA3a, PMCA3b, and PMCA4) and 7 isoforms of NCXs (NCX1a, NCX1b, NCX2a, NCX2b, NCX3, NCX4a, and NCX4b) were identified. In the reverse-transcriptase polymerase chain reaction (RT-PCR) analysis, 5 PMCAs and 2 NCXs were ubiquitously expressed in various tissues including gills. Triple fluorescence in situ hybridization and immunocytochemistry showed the colocalization of zecac, zpmca2, and zncx1b mRNAs in a portion of gill MR cells (using Na⁺-K⁺-ATPase as the marker), implying a subset of ionocytes specifically responsible for the transepithelial Ca²⁺ uptake in zebrafish gills. The gene expressions in gills of high- or low-Ca²⁺-acclimated zebrafish by quantitative real-time PCR analysis showed that zecac was the only gene regulated in response to environmental Ca²⁺ levels, while zpmcas and zncxs remained steady.

Conclusion

The present study provides molecular evidence for the specific isoforms of Ca²⁺ transporters, zECaC, zPMCA2, and zNCX1b, supporting the current Ca²⁺ uptake model, in which ECaC may play a role as the major regulatory target for this mechanism during environmental challenge.

Molecular cloning and expression of the duplicated thyroid hormone responsive spot 14 (THRSP) genes in ducks

Thyroid hormone responsive Spot 14 (THRSP) is suggested as a transcription factor involved in the regulation of adipogenic enzymes by 3 thyroid response elements in the promoter region. In the chicken genome, THRSP gene was identified to duplicate into 2 paralogs, THRSPalpha and THRSPbeta. In the current study, cDNA sequences of the duplicated duck THRSP genes were cloned by real-time PCR and rapid amplification of cDNA ends. Duck THRSPalpha and THRSPbeta were predicted to encode peptides with 133 amino acids, which had 74 and 68% sequence identity at cDNA level, 78 and 74% identity at amino acid level to the chicken counterparts, respectively. A high percentage (73.1%) of G and C nucleotides were found in the 3' untranslated region of duck THRSPbeta cDNA. Although a low similarity of peptide composition was shared between ducks and mammals, and a moderate similarity was shared between ducks and chickens, many predicted properties of THRSP, including the pI, subcellular localization and functional domains seemed to be highly conserved. The present study demonstrated that the duck THRSP gene duplicates into the 2 paralogs as in chickens. Phylogenetic analysis indicated that the duplication for THRSP paralogs appeared to have taken place preceding the chicken-duck split, and the diverging rate between THRSP paralogs seemed faster in the chicken genome than that in the duck genome. Expression analysis by real-time quantitative PCR showed that THRSP paralogs in ducks were more actively transcribed in fat tissues (i.e., s.c. fat and abdominal fat) than in liver, and the mRNA concentrations of THRSPbeta were higher than that of THRSPalpha in liver and s.c. fat.

Alternate mRNA processing of the hepatocyte nuclear factor genes and its role in monogenic diabetes

Variation in mRNA processing has the capacity to exert fine control over gene expression in most cell types. The hepatic nuclear factor genes, like approximately 74% of the genome, produce multiple transcripts. Hepatic nuclear factor isoforms exhibit both spatial and temporal variation in expression. In this review, the known isoforms of the hepatocyte nuclear factor-1α, hepatocyte nuclear factor-1β and hepatocyte nuclear factor-4α genes are described and their properties are compared. Finally, data are discussed regarding the influence of hepatocyte nuclear factor-1α alternate mRNA processing on the clinical phenotype of maturity-onset diabetes of the young.

Isomers of the TCF1 gene encoding hepatocyte nuclear factor-1 alpha show differential expression in the pancreas and define the relationship between mutation position and clinical phenotype in monogenic diabetes

The generation of multiple transcripts by mRNA processing has the potential to moderate differences in gene expression both between tissues and at different stages of development. Where gene function is compromised by mutation, the presence of multiple isoforms may influence the resulting phenotype. Heterozygous mutations in the transcription factor hepatocyte nuclear factor-1 alpha (HNF1A or TCF1 gene) result in early-onset diabetes as a result of pancreatic beta-cell dysfunction. We investigated the expression of the three alternatively processed isoforms of the HNF1A gene and their impact on the phenotype associated with mutations. Real-time PCR demonstrated variation in tissue expression of HNF1A isomers: HNF1A(A), with the lowest transactivation activity compared with the truncated isoforms HNF1A(B) and HNF1A(C), is the major isomer in liver (54%) and kidney (67%) but not in adult pancreas (24%) and islets (26%). However, in fetal pancreas HNF1A(A) is the major transcript (84%), which supports developmental regulation of isomer expression. We examined whether the isomers affected by the mutation altered the diabetes phenotype in 564 subjects with 123 mutations in HNF1A. Mutations that affected only isomer HNF1A(A) (exons 8-10) were diagnosed later (25.5 years) than mutations affecting all three isomers (exons 1-6) (18.0 years) (P=0.006). This first genotype/phenotype relationship described for patients with HNF1A mutations is explained by isomer structure and not by either mutation type or functional domain. We conclude that all three isomers may be critical for beta-cell function and could play a role in both the developing and mature beta cell.

Tissue-specific distribution of an alternatively spliced COL4A5 isoform and non-random X chromosome inactivation reflect phenotypic variation in heterozygous X-linked Alport syndrome

A novel type of hereditary transmission of COL4A5 in a Japanese family with X-linked Alport syndrome was detected through analysis of cDNA sequences and an X-chromosome inactivation assay. A female patient with moderately altered renal function, who was diagnosed with Alport syndrome by renal biopsy, and her mother, who was undergoing maintenance haemodialysis, showed similar tissue-specific expression of a truncated isoform of COL4A5, which was generated by alternative splicing without a splice-site mutation. Expression of the truncated isoform occurred in the renal specimen derived from the patient, but not in specimens from controls. Genomic analysis revealed two point mutations (c.4821 + 121, T>C; c.4822-151_150, ins T) in intron 49 of COL4A5 from the patient. The peripheral blood mononuclear cells of the patient and her mother showed non-random lyonization. While the females showed only renal impairment, an affected male in the same family suffered from severe renal insufficiency, visual defect and hearing disturbances. Hence, we suggest that this type of heredity COL4A5 presents with phenotypic variation in female heterozygous X-linked Alport syndrome patients.

Gene expression after freshwater transfer in gills and opercular epithelia of killifish: insight into divergent mechanisms of ion transport

We have explored the molecular basis for differences in physiological function between the gills and opercular epithelium of the euryhaline killifish Fundulus heteroclitus. These tissues are functionally similar in seawater, but in freshwater the gills actively absorb Na+ but not Cl-, whereas the opercular epithelium actively absorbs Cl- but not Na+. These differences in freshwater physiology are likely due to differences in absolute levels of gene expression (measured using real-time PCR), as several proteins important for Na+ transport, namely Na+,H+-exchanger 2 (NHE2), carbonic anhydrase 2 (CA2), Na+,HCO3- cotransporter 1, and V-type H+-ATPase, were expressed at 3- to over 30-fold higher absolute levels in the gills. In gills, transfer from 10% seawater to freshwater increased the activity of Na+,K+-ATPase by twofold (from 12 h to 7 days), increased the expression of NHE2 (at 12 h) and CA2 (from 12 h to 7 days), and decreased the expression of NHE3 (from 12 h to 3 days). In opercular epithelium, NHE2 was not expressed; furthermore, Na+,K+-ATPase activity was unchanged after transfer to freshwater, CA2 mRNA levels decreased, and NHE3 levels increased. Consistent with their functional similarities in seawater, killifish gills and opercular epithelium expressed Na+,K+-ATPase alpha 1a, Na+,K+,2Cl- cotransporter 1 (NKCC1), cystic fibrosis transmembrane conductance regulator (CFTR) Cl- channel and the signalling protein 14-3-3a at similar absolute levels. Furthermore, NKCC1 and CFTR were suppressed equally in each tissue after freshwater transfer, and 14-3-3a mRNA increased in both. These results provide insight into the mechanisms of ion transport by killifish gills and opercular epithelia, and demonstrate a potential molecular basis for the differences in physiological function between these two organs.