Genomic structure of growth hormone genes in chinook salmon (Oncorhynchus tshawytscha): presence of two functional genes, GH-I and GH-II, and a male-specific pseudogene, GH-psi

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

A nonfunctional copy of the salmonid sex-determining gene (sdY) is responsible for the “apparent” XY females in Chinook salmon, Oncorhynchus tshawytscha

Many salmonids have a male heterogametic (XX/XY) sex determination system, and they are supposed to have a conserved master sex-determining gene (sdY) that interacts at the protein level with Foxl2 leading to the blockage of the synergistic induction of Foxl2 and Nr5a1 of the cyp19a1a promoter. However, this hypothesis of a conserved master sex-determining role of sdY in salmonids is challenged by a few exceptions, one of them being the presence of naturally occurring “apparent” XY Chinook salmon, Oncorhynchus tshawytscha, females. Here, we show that some XY Chinook salmon females have a sdY gene (sdY-N183), with 1 missense mutation leading to a substitution of a conserved isoleucine to an asparagine (I183N). In contrast, Chinook salmon males have both a nonmutated sdY-I183 gene and the missense mutation sdY-N183 gene. The 3-dimensional model of SdY-I183N predicts that the I183N hydrophobic to hydrophilic amino acid change leads to a modification in the SdY β-sandwich structure. Using in vitro cell transfection assays, we found that SdY-I183N, like the wild-type SdY, is preferentially localized in the cytoplasm. However, compared to wild-type SdY, SdY-I183N is more prone to degradation, its nuclear translocation by Foxl2 is reduced, and SdY-I183N is unable to significantly repress the synergistic Foxl2/Nr5a1 induction of the cyp19a1a promoter. Altogether, our results suggest that the sdY-N183 gene of XY Chinook females is nonfunctional and that SdY-I183N is no longer able to promote testicular differentiation by impairing the synthesis of estrogens in the early differentiating gonads of wild Chinook salmon XY females.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Polymorphisms of the growth hormone gene and their association with growth traits and sex in Sarcocheilichthys sinensis

The growth hormone gene (gh) of Sarcocheilichthys sinensis was cloned and characterized in this study. The cDNA length of gh was 973 bp, containing a 5'-UTR of 15 bp, a 3'-UTR of 325 bp and an open reading frame of 633 bp. The genomic DNA of gh was 2135 bp in length containing five exons and four introns. The precursor peptide of gh contained 210 amino acids (aa), including a signal peptide of 22 aa (Met¹-Ala²²) and a mature region of 188 aa (Ser²³-Leu²¹⁰). The similarity and identity ranges of the gh precursor peptide with those of other cyprinids were 88.6%-99.0% and 84.8%-98.6%, respectively. The gh of S. sinensis expressed at the highest level in the pituitary, and its expression was also detected in muscle and brain. Six polymorphic sites were detected in intron 1 (g.51InDel, g.64InDel and g.242InDel), intron 2 (g.864T>C), intron 3 (g.1017InDel) and intron 4 (g.1541A>G). Among these sites, g.242InDel was significantly associated with condition factor, g.1541A>G was associated with all six growth traits, while g.864T>C was associated with sex. The data obtained herein provide useful information for further studies on the regulation mechanisms of growth and sexual growth differences in S. sinensis.

Discordance for genotypic sex in phenotypic female Atlantic salmon (Salmo salar) is related to a reduced sdY copy number

The master sex determinant in rainbow trout (Oncorhynchus mykiss), sexually dimorphic on the Y chromosome (sdY), is strongly but not perfectly associated with male phenotype in several other species from the family Salmonidae. Currently, the cause and implications of discordance for sdY-predicted genotypic sex and phenotypic sex in these species is unclear. Using an established multiplex PCR test for exons 2 and 3 of sdY, we demonstrated that sdY-predicted genotypic sex was discordant with histologically evidenced phenotypic sex in 4% of 176 Tasmanian Atlantic salmon. All discordant individuals were phenotypic females presenting a male genotype. Using real-time qPCR assays that we developed and validated for exons 2, 3 and 4 of sdY, all genotype-phenotype discordant females were confirmed to possess sdY, albeit at a reduced number of copies when compared to phenotypic males. The real-time qPCR assays also demonstrated reduced levels of sdY in 30% of phenotypic females that the established multiplex PCR-based test indicated to be devoid of sdY. These findings suggest sdY may be reduced in copy number or mosaicked in the genomic DNA of sdY-positive phenotypic female Atlantic salmon and highlight the importance of understanding the effects of reduced sdY copies on the development of phenotypic sex.

A newly developed genetic sex marker and its application to understanding chemically induced feminisation in roach (Rutilus rutilus)

Oestrogenic wastewater treatment works (WwTW) effluents discharged into UK rivers have been shown to affect sexual development, including inducing intersex, in wild roach (Rutilus rutilus). This can result in a reduced breeding capability with potential population level impacts. In the absence of a sex probe for roach it has not been possible to confirm whether intersex fish in the wild arise from genetic males or females, or whether sex reversal occurs in the wild, as this condition can be induced experimentally in controlled exposures to WwTW effluents and a steroidal oestrogen. Using restriction site-associated DNA sequencing (RAD-seq), we identified a candidate for a genetic sex marker and validated this marker as a sex probe through PCR analyses of samples from wild roach populations from nonpolluted rivers. We also applied the sex marker to samples from roach exposed experimentally to oestrogen and oestrogenic effluents to confirm suspected phenotypic sex reversal from males to females in some treatments, and also that sex-reversed males are able to breed as females. We then show, unequivocally, that intersex in wild roach populations results from feminisation of males, but find no strong evidence for complete sex reversal in wild roach at river sites contaminated with oestrogens. The discovered marker has utility for studies in roach on chemical effects, wild stock assessments, and reducing the number of fish used where only one sex is required for experimentation. Furthermore, we show that the marker can be applied nondestructively using a fin clip or skin swab, with animal welfare benefits.

Loci associated with variation in gene expression and growth in juvenile salmon are influenced by the presence of a growth hormone transgene

BACKGROUND

Growth regulation is a complex process influenced by genetic and environmental factors. We examined differences between growth hormone (GH) transgenic (T) and non-transgenic (NT) coho salmon to elucidate whether the same loci were involved in controlling body size and gene expression phenotypes, and to assess whether physiological transformations occurring from GH transgenesis were under the influence of alternative pathways. The following genomic techniques were used to explore differences between size classes within and between transgenotypes (T vs. NT): RNA-Seq/Differentially Expressed Gene (DEG) analysis, quantitative PCR (qPCR) and OpenArray analysis, Genotyping-by-Sequencing, and Genome-Wide Association Study (GWAS).

RESULTS

DEGs identified in comparisons between the large and small tails of the size distributions of T and NT salmon (NT_Large, NT_Small, T_Large and T_Small) spanned a broad range of biological processes, indicating wide-spread influence of the transgene on gene expression. Overexpression of growth hormone led to differences in regulatory loci between transgenotypes and size classes. Expression levels were significantly greater in T fish at 16 of 31 loci and in NT fish for 10 loci. Eleven genes exhibited different mRNA levels when the interaction of size and transgenotype was considered (IGF1, IGFBP1, GH, C3-4, FAS, FAD6, GLUT1, G6PASE1, GOGAT, MID1IP1). In the GWAS, 649 unique SNPs were significantly associated with at least one study trait, with most SNPs associated with one of the following traits: C3_4, ELA1, GLK, IGF1, IGFBP1, IGFII, or LEPTIN. Only 1 phenotype-associated SNP was found in common between T and NT fish, and there were no SNPs in common between transgenotypes when size was considered.

CONCLUSIONS

Multiple regulatory loci affecting gene expression were shared between fast-growing and slow-growing fish within T or NT groups, but no such regulatory loci were found to be shared between NT and T groups. These data reveal how GH overexpression affects the regulatory responses of the genome resulting in differences in growth, physiological pathways, and gene expression in T fish compared with the wild type. Understanding the complexity of regulatory gene interactions to generate phenotypes has importance in multiple fields ranging from applications in selective breeding to quantifying influences on evolutionary processes.

Unstable Linkage of Molecular Markers with Sex Determination Gene in Pacific Salmon (Oncorhynchus spp.)

In the present study, we tested the congruence between the sdY sex-specific marker and other commonly used male markers, located on the Y-chromosome, with the sex phenotypes in 5 species of Pacific salmon in Asian waters, including Chinook, chum, sockeye, masu, and pink salmon. We found that the localization of the sex-specific marker of both males and females of these species is not consistent with the phenotypic sex. Also, no linkage was found between noncoding markers and the sdY gene in the same species samples. Possible genetic and physiological mechanisms underlying this discrepancy are discussed.

Comparison of the frequency of the growth hormone pseudogene between juvenile and adult female masu salmon Oncorhynchus masou

The proportions of individual masu salmon Oncorhynchus masou whose genotypic and phenotypic sex differed were compared among juvenile and adult fish in three rivers, and genotypically male but phenotypically female individuals were observed 6-16% more often among adults than among juveniles. This suggests that during the transition from juvenile to adult, survival rates of genotypically male but phenotypically female individuals are higher than those of normal females. In contrast, genotypically female but phenotypically male individuals were only found in the juvenile period, which suggests that they exhibit a decreased survival rate in comparison with normal males.

Comparative Analysis of the Shared Sex-Determination Region (SDR) among Salmonid Fishes

Salmonids present an excellent model for studying evolution of young sex-chromosomes. Within the genus, Oncorhynchus, at least six independent sex-chromosome pairs have evolved, many unique to individual species. This variation results from the movement of the sex-determining gene, sdY, throughout the salmonid genome. While sdY is known to define sexual differentiation in salmonids, the mechanism of its movement throughout the genome has remained elusive due to high frequencies of repetitive elements, rDNA sequences, and transposons surrounding the sex-determining regions (SDR). Despite these difficulties, bacterial artificial chromosome (BAC) library clones from both rainbow trout and Atlantic salmon containing the sdY region have been reported. Here, we report the sequences for these BACs as well as the extended sequence for the known SDR in Chinook gained through genome walking methods. Comparative analysis allowed us to study the overlapping SDRs from three unique salmonid Y chromosomes to define the specific content, size, and variation present between the species. We found approximately 4.1 kb of orthologous sequence common to all three species, which contains the genetic content necessary for masculinization. The regions contain transposable elements that may be responsible for the translocations of the SDR throughout salmonid genomes and we examine potential mechanistic roles of each one.

Development of a Chinook Salmon Sex Identification SNP Assay Based on the Growth Hormone Pseudogene

Genotypic sex identification assays can provide valuable information about fish populations when phenotypic sex determination is difficult. Here we describe the development of a TaqMan® assay (Ots_SexID) designed to identify the genotypic sex by targeting a region previously examined in the growth hormone pseudogene for winter-run Chinook salmon (Oncorhynchus tshawytscha) collected from the Sacramento River and spawned at the Livingston Stone National Fish Hatchery. Accuracy of the marker was assessed by comparing genotypic sex assignments for Chinook salmon spawned at Livingston Stone National Fish hatchery in 2012 (n = 84) with phenotypic sex recorded during spawning. Genotypic sex was observed to be concordant with phenotypic sex identified using Ots_SexID in 83/84 individuals, suggesting that the assay could be used to predict phenotypic sex with ∼︀99% accuracy. To evaluate the utility of the TaqMan assay in other parts of the species’ range, we examined collections from 29 other populations ranging from Alaska to California. Genotypic sex assignments based on the assay were generally concordant with observed phenotypes, but there were some strong exceptions. These results suggest that the new assay will be very useful for Sacramento River winter-run Chinook salmon, but also highlight the importance of thoroughly testing any genotypic sex identification assay before application in a population of interest.

Reduced oceanic growth of growth hormone pseudogene-positive female masu salmon Oncorhynchus masou

This study compared the growth rates of female masu salmon Oncorhynchus masou, who possessed a male-specific gene marker, the growth hormone pseudogene (GHp), and normal females, as estimated from their scale growth. There was a difference between the growth rates of GHp-positive females and those of normal females of the same age during the ocean period, although their growth rates during the river period were similar. These results suggest that GHp-positive salmonid females exhibit male-like characteristics such as reduced feeding activity during the ocean period, which depresses their growth.

High-throughput sequencing and pathway analysis reveal alteration of the pituitary transcriptome by 17α-ethynylestradiol (EE2) in female coho salmon, Oncorhynchus kisutch

Considerable research has been done on the effects of endocrine disrupting chemicals (EDCs) on reproduction and gene expression in the brain, liver and gonads of teleost fish, but information on impacts to the pituitary gland are still limited despite its central role in regulating reproduction. The aim of this study was to further our understanding of the potential effects of natural and synthetic estrogens on the brain-pituitary-gonad axis in fish by determining the effects of 17α-ethynylestradiol (EE2) on the pituitary transcriptome. We exposed sub-adult coho salmon (Oncorhynchus kisutch) to 0 or 12 ng EE2/L for up to 6 weeks and effects on the pituitary transcriptome of females were assessed using high-throughput Illumina(®) sequencing, RNA-Seq and pathway analysis. After 1 or 6 weeks, 218 and 670 contiguous sequences (contigs) respectively, were differentially expressed in pituitaries of EE2-exposed fish relative to control. Two of the most highly up- and down-regulated contigs were luteinizing hormone β subunit (241-fold and 395-fold at 1 and 6 weeks, respectively) and follicle-stimulating hormone β subunit (-3.4-fold at 6 weeks). Additional contigs related to gonadotropin synthesis and release were differentially expressed in EE2-exposed fish relative to controls. These included contigs involved in gonadotropin releasing hormone (GNRH) and transforming growth factor-β signaling. There was an over-representation of significantly affected contigs in 33 and 18 canonical pathways at 1 and 6 weeks, respectively, including circadian rhythm signaling, calcium signaling, peroxisome proliferator-activated receptor (PPAR) signaling, PPARα/retinoid x receptor α activation, and netrin signaling. Network analysis identified potential interactions between genes involved in circadian rhythm and GNRH signaling, suggesting possible effects of EE2 on timing of reproductive events.

Characterization of the OmyY1 Region on the Rainbow Trout Y Chromosome

We characterized the male-specific region on the Y chromosome of rainbow trout, which contains both sdY (the sex-determining gene) and the male-specific genetic marker, OmyY1. Several clones containing the OmyY1 marker were screened from a BAC library from a YY clonal line and found to be part of an 800 kb BAC contig. Using fluorescence in situ hybridization (FISH), these clones were localized to the end of the short arm of the Y chromosome in rainbow trout, with an additional signal on the end of the X chromosome in many cells. We sequenced a minimum tiling path of these clones using Illumina and 454 pyrosequencing. The region is rich in transposons and rDNA, but also appears to contain several single-copy protein-coding genes. Most of these genes are also found on the X chromosome; and in several cases sex-specific SNPs in these genes were identified between the male (YY) and female (XX) homozygous clonal lines. Additional genes were identified by hybridization of the BACs to the cGRASP salmonid 4x44K oligo microarray. By BLASTn evaluations using hypothetical transcripts of OmyY1-linked candidate genes as query against several EST databases, we conclude at least 12 of these candidate genes are likely functional, and expressed.

Coding Gene SNP Mapping Reveals QTL Linked to Growth and Stress Response in Brook Charr (Salvelinus fontinalis)

Growth performance and reduced stress response are traits of major interest in fish production. Growth and stress-related quantitative trait loci (QTL) have been already identified in several salmonid species, but little effort has been devoted to charrs (genus Salvelinus). Moreover, most QTL studies to date focused on one or very few traits, and little investigation has been devoted to QTL identification for gene expression. Here, our objective was to identify QTL for 27 phenotypes related to growth and stress responses in brook charr (Salvelinus fontinalis), which is one of the most economically important freshwater aquaculture species in Canada. Phenotypes included 12 growth parameters, six blood and plasma variables, three hepatic variables, and one plasma hormone level as well as the relative expression measurements of five genes of interest linked to growth regulation. QTL analysis relied on a linkage map recently built from S. fontinalis consisting of both single-nucleotide polymorphism (SNP, n = 266) and microsatellite (n =81) markers in an F₂ interstrain hybrid population (n = 171). We identified 63 growth-related QTL and four stress-related QTL across 18 of the 40 linkage groups of the brook charr linkage map. Percent variance explained, confidence interval, and allelic QTL effects also were investigated to provide insight into the genetic architecture of growth- and stress-related QTL. QTL related to growth performance and stress response that were identified could be classified into two groups: (1) a group composed of the numerous, small-effect QTL associated with some traits related to growth (i.e., weight) that may be under the control of a large number of genes or pleiotropic genes, and (2) a group of less numerous QTL associated with growth (i.e., gene expression) and with stress-related QTL that display a larger effect, suggesting that these QTL are under the control of a limited number of genes of major effect. This study represents a first step toward the identification of genes potentially linked to phenotypic variation of growth and stress response in brook charr. The ultimate goal is to provide new tools for developing Molecular Assisted Selection for this species.

Variable incidences and morphological characteristics of female masu salmon Oncorhynchus masou with growth hormone pseudogene

The genetic sex of mature masu salmon Oncorhynchus masou (212 males and 243 females) collected from three rivers located in the mid-western part of Hokkaido was determined using the male-specific genetic marker growth hormone pseudogene (GHp). A total of 72 phenotypic females were found to have GHp, whereas none of the phenotypic males lacked GHp. The occurrence of females with incongruence between genotypic and phenotypic sex varied with the river and reached the highest rate of 67%. In contrast, none of the phenotypic males were found to have a female genotype. Comparison of adult body size showed that females with GHp were significantly smaller than those without GHp. Moreover, comparison of secondary characteristics revealed that both upper jaw length and head length were greater in females with GHp than in those without GHp. These results suggest that genetic changes (e.g. transposition and mutation) may also affect morphological characteristics.

Disruption of the salmon reproductive endocrine axis through prolonged nutritional stress: changes in circulating hormone levels and transcripts for ovarian genes involved in steroidogenesis and apoptosis

Mechanisms regulating the normal progression of ovarian follicular growth versus onset of atresia in fishes are poorly understood. To gain a better understanding of these processes, we exposed immature female coho salmon (Oncorhynchus kisutch) to prolonged fasting to induce follicular atresia and monitored body growth, development of the ovarian follicles, changes in reproductive hormones, and transcripts for ovarian genes. Prolonged fasting reduced body and ovary weight and increased the appearance of atretic follicles relative to normally fed controls. Endocrine analyses showed that fasting reduced plasma insulin-like growth factor 1 (IGF1), estradiol-17β (E2), and pituitary, but not plasma, levels of follicle-stimulating hormone (FSH). Transcripts for ovarian fsh receptor (fshr) and steroidogenesis-related genes, such as steroidogenic acute regulatory protein (star), 3β-hydroxysteroid dehydrogenase (hsd3b), and P450 aromatase (cyp19a1a) were significantly lower in fasted fish. Ovarian expression of apoptosis-related genes, such as Fas-associated death domain (fadd), caspase 8 (casp8), caspase 3 (casp3), and caspase 9 (casp9) were significantly elevated in fasted fish compared to fed fish, indicating that apoptosis is involved in the process of atresia in this species. Interestingly, some genes such as fadd, casp8, casp3, and hsd3b, were differentially expressed prior to increases in the number of atretic follicles and reductions in hormone levels induced by fasting, and may therefore have potential as early indicators of atresia. Together these results suggest that prolonged nutritional stress may disrupt the reproductive system and induce follicular atresia in part via reductions in ovarian IGF and FSH signaling, and downstream effects on steroidogenesis-related genes and E2 production.

Screening and identification of male-specific DNA fragments in common carps Cyprinus carpio using suppression subtractive hybridization

In this study, a sex subtractive genomic DNA library was constructed using suppression subtractive hybridization (SSH) between male and female Cyprinus carpio. Twenty-two clones with distinguishable hybridization signals were selected and sequenced. The specific primers were designed based on the sequence data. Those primers were then used to amplify the sex-specific fragments from the genomic DNA of male and female carp. The amplified fragments from two clones showed specificity to males but not to females, which were named as Ccmf2 [387 base pairs (bp)] and Ccmf3 (183 bp), respectively. The sex-specific pattern was analysed in a total of 40 individuals from three other different C. carpio. stocks and grass carp Ctenopharyngodon idella using Ccmf2 and Ccmf3 as dot-blotting probes. The results revealed that the molecular diversity exists on the Y chromosome of C. carpio. No hybridization signals, however, were detected from individuals of C. idella, suggesting that the two sequences are specific to C. carpio. No significant homologous sequences of Ccmf2 and Ccmf3 were found in GenBank. Therefore, it was interpreted that the results as that Ccmf2 and Ccmf3 are two novel male-specific sequences; and both fragments could be used as markers to rapidly and accurately identify the genetic sex of part of C. carpio. This may provide a very efficient selective tool for practically breeding monosex female populations in aquacultural production.

Growth-related quantitative trait loci in domestic and wild rainbow trout (Oncorhynchus mykiss)

Background

Somatic growth is a complex process that involves the action and interaction of genes and environment. A number of quantitative trait loci (QTL) previously identified for body weight and condition factor in rainbow trout (Oncorhynchus mykiss), and two other salmonid species, were used to further investigate the genetic architecture of growth-influencing genes in this species. Relationships among previously mapped candidate genes for growth and their co-localization to identified QTL regions are reported. Furthermore, using a comparative genomic analysis of syntenic rainbow trout linkage group clusters to their homologous regions within model teleost species such as zebrafish, stickleback and medaka, inferences were made regarding additional possible candidate genes underlying identified QTL regions.

Results

Body weight (BW) QTL were detected on the majority of rainbow trout linkage groups across 10 parents from 3 strains. However, only 10 linkage groups (i.e., RT-3, -6, -8, -9, -10, -12, -13, -22, -24, -27) possessed QTL regions with chromosome-wide or genome-wide effects across multiple parents. Fewer QTL for condition factor (K) were identified and only six instances of co-localization across families were detected (i.e. RT-9, -15, -16, -23, -27, -31 and RT-2/9 homeologs). Of note, both BW and K QTL co-localize on RT-9 and RT-27. The incidence of epistatic interaction across genomic regions within different female backgrounds was also examined, and although evidence for interaction effects within certain QTL regions were evident, these interactions were few in number and statistically weak. Of interest, however, was the fact that these predominantly occurred within K QTL regions. Currently mapped growth candidate genes are largely congruent with the identified QTL regions. More QTL were detected in male, compared to female parents, with the greatest number evident in an F₁ male parent derived from an intercross between domesticated and wild strain of rainbow trout which differed strongly in growth rate.

Conclusions

Strain background influences the degree to which QTL effects are evident for growth-related genes. The process of domestication (which primarily selects faster growing fish) may largely reduce the genetic influences on growth-specific phenotypic variation. Although heritabilities have been reported to be relatively high for both BW and K growth traits, the genetic architecture of K phenotypic variation appears less defined (i.e., fewer major contributing QTL regions were identified compared with BW QTL regions).

Gene structure and functional characterization of growth hormone in dogfish, Squalus acanthias

Dogfish (Squalus acanthias) growth hormone (GH) was identified by cDNA cloning and protein purification from the pituitary gland. Dogfish GH cDNA encoded a prehormone of 210 amino acids (aa). Sequence analysis of purified GH revealed that the prehormone is composed of a signal peptide of 27 aa and a mature protein of 183 aa. Dogfish GH showed 94% sequence identity with blue shark GH, and also showed 37-66%, 26%, and 48-67% sequence identity with GH from osteichtyes, an agnathan, and tetrapods. The site of production was identified through immunocytochemistry to be cells of the proximal pars distalis of the pituitary gland. Dogfish GH stimulates both insulin-like growth factor-I and II mRNA levels in dogfish liver in vitro. The dogfish GH gene consisted of five exons and four introns, the same as in lamprey, teleosts such as cypriniforms and siluriforms, and tetrapods. The 5'-flanking region within 1082 bp of the transcription start site contained consensus sequences for the TATA box, Pit-1/GHF-1, CRE, TRE, and ERE. These results show that the endocrine mechanism for growth stimulation by the GH-IGF axis was established at an early stage of vertebrate evolution, and that the 5-exon-type gene organization might reflect the structure of the ancestral gene for the GH gene family.

The growth hormone-encoding gene isolated and characterized from Labeo rohita Hamilton is expressed in CHO cells under the control of constitutive promoters in 'autotransgene' constructs

The growth hormone (GH) gene along with its regulatory sequences has been isolated from the blood and pituitary gland of Labeo rohita. This GH gene is approximately 2.8 kb long and consists of five exons and four introns of varying sizes with AG/TA in its exon-intron junctions. The promoter has a single cyclic AMP response unit (CRE) element, TATA, CAT and several Pit 1 binding sequences. The 1169-bp gene transcript starts 54 bp upstream of the ATG initiation codon and has two polyadenylation signals, ATTAAA, after the TAG stop codon. The mature mRNA has the poly (A) tail inserted 16 bp downstream of the second polyadenylation signal. Four chimeric 'autotransgenes' were constructed having either histone 3 or beta-actin promoter and cDNA or the total GH gene. The functionality of the individual components of the autotransgene was determined in the Chinese hamster ovary (CHO) cells by transfection experiments. Based on the results, the transcription of the GH gene is initiated at the transcription start signal of the respective promoters and terminates at the 3' regulatory sequence of the GH gene. Expression of GH in CHO cells shows that the fish promoters are active, the splicing signal is recognized, and the mRNA produced is stable and translated. The GH protein produced is effectively translocated and secreted into the medium. These results indicate the usefulness of CHO cells in determining the property of individual components of autotransgenes constructed from L. rohita and overall functional commonality between fish and mammal.

Unilateral ovariectomy increases egg size and reduces follicular atresia in the semelparous coho salmon, Oncorhynchus kisutch

Unilateral ovariectomy (ULO, removal of one ovary) is a powerful technique for studying aspects of reproductive physiology, including follicular recruitment and growth. To examine effects of ULO for the first time in a semelparous species, coho salmon (Oncorhynchus kisutch) were unilaterally ovariectomized during mid-vitellogenesis approximately 3 months before spawning. At termination of the study (79 days post-surgery), single ovaries of ULO fish were gravimetrically equivalent to paired ovaries of sham surgery, control fish. There was no evidence of recruitment of new vitellogenic follicles. Instead, the dramatic increase in ovary mass was attributable to hypertrophy of existing vitellogenic follicles (33% increase in volume) and increased fecundity achieved through a greater than two-fold reduction in follicular atresia. The composition of whole ovaries on a dry weight basis from ULO fish was greater in protein, but lower in lipid than that of control fish. Expressing the data on a per follicle basis, however, showed that follicles of ULO fish contained more protein, ash, water, and lipid. The results indicate that ULO of coho salmon induces compensatory hypertrophy of existing vitellogenic follicles, while maximizing fecundity through reduction of atresia. Thus, 3 months before spawning, coho salmon exhibit the ability to adjust final egg size and number when faced with significant depletion of ovarian follicles. This in vivo system provides a platform for further study of physiological mechanisms regulating follicular growth and atresia, and the trade-off between egg size and egg number.

Isolation, characterization and comparison of Atlantic and Chinook salmon growth hormone 1 and 2

Background

Growth hormone (GH) is an important regulator of skeletal growth, as well as other adapted processes in salmonids. The GH gene (gh) in salmonids is represented by duplicated, non-allelic isoforms designated as gh1 and gh2. We have isolated and characterized gh-containing bacterial artificial chromosomes (BACs) of both Atlantic and Chinook salmon (Salmo salar and Oncorhynchus tshawytscha) in order to further elucidate our understanding of the conservation and regulation of these loci.

Results

BACs containing gh1 and gh2 from both Atlantic and Chinook salmon were assembled, annotated, and compared to each other in their coding, intronic, regulatory, and flanking regions. These BACs also contain the genes for skeletal muscle sodium channel oriented in the same direction. The sequences of the genes for interferon alpha-1, myosin alkali light chain and microtubule associated protein Tau were also identified, and found in opposite orientations relative to gh1 and gh2. Viability of each of these genes was examined by PCR. We show that transposon insertions have occurred differently in the promoters of gh, within and between each species. Other differences within the promoters and intronic and 3'-flanking regions of the four gh genes provide evidence that they have distinct regulatory modes and possibly act to function differently and/or during different times of salmonid development.

Conclusion

A core proximal promoter for transcription of both gh1 and gh2 is conserved between the two species of salmon. Nevertheless, transposon integration and regulatory element differences do exist between the promoters of gh1 and gh2. Additionally, organization of transposon families into the BACs containing gh1 and for the BACs containing gh2, are very similar within orthologous regions, but much less clear conservation is apparent in comparisons between the gh1- and gh2-containing paralogous BACs for the two fish species. This is consistent with the hypothesis that a burst of transposition activity occurred during the speciation events which led to Atlantic and Pacific salmon. The Chinook and other Oncorhynchus GH1s are strikingly different in comparison to the other GHs and this change is not apparent in the surrounding non-coding sequences.

Y-specific sequences and polymorphisms in rainbow trout and Chinook salmon

Improved methods for genetically sexing salmonids and for characterization of Y-chromosome homologies between species can contribute to understanding the evolution of sex chromosomes and sex-determining mechanisms. In this study we have explored 12.5 kb of Y-chromosome-specific sequence flanking the previously described OtY2 locus in Chinook salmon (Oncorhynchus tshawytscha) and 21 kb of homologous rainbow trout (Oncorhynchus mykiss) Y-chromosome-specific sequence. This is the first confirmed Y-specific sequence for rainbow trout. New Y-specific markers are described for Chinook salmon (OtY3) and rainbow trout (OmyY1), which are readily detected by PCR assays and are advantageous because they also produce autosomal control amplification products. Additionally, AFLP analysis of Chinook salmon yielded another potential Y-chromosome marker. These descriptions will facilitate genotypic sexing and should be useful for population studies of Y-chromosome polymorphisms and for future studies to characterize what appears to be a common sex-determining mechanism between these species.

Characterization of a chromosomal rearrangement responsible for producing "apparent" XY-female fall-run Chinook salmon in California

Fluorescence in situ hybridization (FISH) was used to identify the X and Y chromosomes of offspring produced by normal and "apparent" XY-female fall-run Chinook salmon (Oncorhynchus tshawytscha) from California. FISH experiments were performed using probes to 2 sex-linked loci, growth hormone pseudogene (GH-Psi), and OtY1, as well as a probe to a sex-linked microsatellite (Omy7INRA). Comparison of FISH staining patterns between the offspring produced by normal and apparent XY-females revealed that the apparent XY-female examined transmitted a "Y-like" chromosome with an attenuated OtY1 and GH-Psi signal to half of its offspring. Segregation analysis of microsatellites derived from rainbow trout (Oncorhynchus mykiss) with respect to phenotypic sex was carried out for 2 normal and 2 apparent XY-female Chinook salmon families. Inheritance patterns of Omy7INRA were consistent with this locus being closely linked to GH-Psi in males and in apparent XY-females carrying the Y-like chromosome. Another microsatellite locus (Omm1077) was closely linked to the primary sex-determining locus (SEX) in males but not to GH-Psi/OtY1 in apparent XY-females. The FISH analyses suggest that apparent XY-female fall-run Chinook salmon in California are not the product of a Y chromosome to autosome translocation. Despite the combined FISH and inheritance analyses, we were unable to differentiate between 2 alternative explanations for apparent XY-females, namely, recombination of markers between the sex chromosomes, or a Y chromosome with a dysfunctional or missing sex-determining region.

Identification of the sex chromosome pair in chum salmon (Oncorhynchus keta) and pink salmon (Oncorhynchus gorbuscha)

Fluorescence in situ hybridization (FISH) using a probe to the male-specific GH-Y (growth hormone pseudogene) was used to identify the Y chromosome in the karyotypes of chum salmon (Oncorhynchus keta) and pink salmon (Oncorhynchus gorbuscha). The sex chromosome pair is a small acrocentric chromosome pair in chum salmon and the smallest metacentric chromosome pair in pink salmon. Both of these chromosome pairs are morphologically different from the sex chromosome pairs in chinook salmon (Oncorhynchus tshawytscha) and coho salmon (Oncorhynchus kisutch). The 5S rRNA genes are on multiple chromosome pairs including the sex chromosome pair in chum salmon, but at the centromeres of two autosomal metacentric pairs in pink salmon. The sex chromosome pairs and the chromosomal locations of the 5S rDNA appear to be different in all five of the North American Pacific salmon species and rainbow trout. The implications of these results for evolution of sex chromosomes in salmonids are discussed.

Assessment of environmental estrogens and the intersex/sex reversal capacity for chinook salmon (Oncorhynchus tshawytscha) in primary and final municipal wastewater effluents

A trickling filter/solid contact (TF/SC) biological secondary treatment plant with chlorine disinfection serving a suburban population of 740,000 was assessed for environmental estrogens. Weekly grab samples were taken at established sampling points and analyzed for various pertinent environmental estrogens including industrial chemicals, and natural and synthetic steroidal estrogens. Additionally, human estrogen receptor (hER) activity and capacity to elicit intersex/sex reversal for the wastewater was monitored using a recombinant yeast assay and whole fish exposures, respectively. hER activity levels varied from 76 to 106 ng/L E2 equivalents in the primary effluent, and were reduced by 25% by biological treatment. For the primary and final effluent no evidence of sex reversal or intersex was apparent in any of the treatment groups (1%, 3%, 10%, or 100%) based on genetic sex determinations and histological examination of the gonads in alevin from 28 d exposed chinook salmon (Oncorhynchus tshawytscha) eggs.

Transgene constructs in coho salmon (Oncorhynchus kisutch) are repeated in a head-to-tail fashion and can be integrated adjacent to horizontally-transmitted parasite DNA

Currently, little information is available regarding the molecular organization of integrated transgenes in genetically-engineered fish. We performed a detailed structural analysis of an inserted transgene in one strain (M77) of transgenic coho salmon (Oncorhynchus kisutch) containing a salmon growth hormone gene construct (OnMTGH1). Microinjected DNA was found to have inserted into a single site in the coho salmon genome, and was organized with four complete internal copies and two partial terminal copies of the OnMTGH1 construct. All construct copies were organized in a direct-tandem (head-to-tail) repeat fashion in strain M77 and five additional strains (one also possessed a second recombinant junction fragment). For strain M77, the junctions between the transgene insert and the insertion point within the wild-type genome were cloned from strain-specific cosmid libraries and sequenced, revealing that the transgene insertion was accompanied by a deletion of 587 bp of wild-type DNA as well as a small insertion (19 bp) of unknown DNA upstream and a 14 bp direct- tandem duplication of sequence downstream. Upstream and downstream wild-type DNA sequence contained several repetitive sequence elements based on Southern blot analysis and homology to repetitive sequences in GenBank. In the downstream flank, a pseudogene sequence was also identified which has high homology to the CA membrane protein gene from Schistosoma japonicum, a parasite closely related to Sanguinicola sp. parasites which infect salmonids. Whether the presence of an inserted transgene and the presence of potentially horizontally-transmitted DNA are indicative of a genomic region with a predisposition for insertion of foreign DNA requires further study. The information derived from this transgene structure provides information useful for comparison to other transgenic organisms and for determination of the mechanism of transgene integration in lower vertebrates.

Genomic structure of the sea lamprey growth hormone-encoding gene

Growth hormone (GH) belongs to a family of pituitary hormones together with prolactin and somatolactin. In our previous study, GH and its cDNA were identified in the pituitary gland of the sea lamprey, Petromyzon marinus, an extant representative of the most ancient class of vertebrates, and isolated GH stimulated expression of insulin-like growth factor in the liver. The evidence suggests that GH is the ancestral hormone in the molecular evolution of the GH/PRL/SL family and that the endocrine mechanism for growth stimulation was established at an early stage in the evolution of vertebrates. To further understand the molecular evolution of the GH/PRL/SL gene family, we report the genomic structure of sea lamprey GH including its 5'-flanking region, being cloned by PCR using specific primers prepared from its cDNA. The sea lamprey GH gene consists of 13,604 bp, making it the largest of all the GH genes. The 5'-flanking region within 697 bp contains consensus sequences for a TATA box, two Pit-1/GHF-1, three TRE, and a CRE. The sea lamprey GH gene consists of five exons and four introns, the same as in mammals, birds, and teleosts such as cypriniforms and siluriforms with the exception of some teleosts such as salmoniforms, percififorms, and tetradontiforms, in which there is an additional intron in the 5th exon. The 5-exon-type gene organization might reflect the structure of the ancestral gene for the GH/PRL/SL gene family.

The dawn and evolution of hormones in the adenohypophysis

The adenohypophysial hormones have been believed to have evolved from several ancestral genes by duplication followed by evolutionary divergence. To understand the origin and evolution of the endocrine systems in vertebrates, we have characterized adenohypophysial hormones in an agnathan, the sea lamprey Petromyzon marinus. In gnathostomes, adrenocorticotropin (ACTH) and melanotropin (MSH) together with beta-endorphins (beta-END) are encoded in a single gene, designated as proopiomelanocortin (POMC), however in sea lamprey, ACTH and MSH are encoded in two distinct genes, proopoicortin (POC) gene and proopiomelanotropin (POM) gene, respectively. The POC and POM genes are expressed specifically in the rostral pars distalis (RPD) and the pars intermedia (PI), respectively. Consequently, the final products from both tissues are the same in all vertebrates, i.e., ACTH from the PD and MSH from the PI. The POMC gene might have been established in the early stages of invertebrate evolution by internal gene duplication of the MSH domains. The ancestral gene might be then inherited in lobe-finned fish and tetrapods, while internal duplication and deletion of MSH domains as well as duplication of whole POMC gene took place in lamprey and gnathostome fish. Sea lamprey growth hormone (GH) is expressed in the cells of the dorsal half of the proximal pars distalis (PPD) and stimulates the expression of an insulin-like growth factor (IGF) gene in the liver as in other vertebrates. Its gene consists of 5 exons and 4 introns spanning 13.6 kb, which is the largest gene among known GH genes. GH appears to be the only member of the GH family in the sea lamprey, which suggests that GH is the ancestral hormone of the GH family that originated first in the molecular evolution of the GH family in vertebrates and later, probably during the early evolution of gnathostomes. The other member of the gene family, PRL and SL, appeared by gene duplication. A beta-chain cDNA belonging to the gonadotropin (GTH) and thyrotropin (TSH) family was cloned. It is expressed in cells of the ventral half of PPD. Since the expression of this gene is stimulated by lamprey gonadotropin-releasing hormone, it was assigned to be a GTHbeta. This GTHbeta is far removed from beta-subunits of LH, FSH, and TSH in an unrooted tree derived from phylogenetic analysis, and takes a position as an out group, suggesting that lampreys have a single GTH gene, which duplicated after the agnathans and prior to the evolution of gnathostomes to give rise to LH and FSH.

Cloning, expression, and biological activity of growth hormone in bullfrog (Rana catesbeiana)

The cDNA (GenBank, AY251538) encoding bullfrog growth hormone (fGH) was cloned by RT-PCR from the total RNA of pituitary glands. Its sequence encoded a putative polypeptide of 215 amino acids, including a signal peptide of 25 amino acids with no change to those of other previous reported bullfrog GHs. The fGH precursor shares 98.1, 96.3, and 95.3% homologies to those of other bullfrog GHs (AAB24792, AAB19428, CAA31038) in amino-acid sequence and its nucleotide sequences of the coding region shares 99.1 and 98.5% homologies to those of previous bullfrog GH genes (S52027 and X12520). The fGH cDNA was also efficiently expressed in Escherichia coli carrying a plasmid pGfGH in which the cDNA was under the control of GST promoter of pGEX1-lambdaT. The expressed fusion protein GST-fGH is comprised about 29.3% of the total cellular protein in such bacteria. The purified GST-fGH cannot only showed a obvious dose-response curve when it reacted with the hepatic membrane receptor proteins from bullfrog, but also significantly increased the body weight and length of bullfrog after twice injection and such effects lasted two or three weeks after the last injection with purified GST-fGH.

Previtellogenic oocyte growth in salmon: relationships among body growth, plasma insulin-like growth factor-1, estradiol-17beta, follicle-stimulating hormone and expression of ovarian genes for insulin-like growth factors, steroidogenic-acute regulatory protein and receptors for gonadotropins, growth hormone, and somatolactin

Body growth during critical periods is known to be an important factor in determining the age of maturity and fecundity in fish. However, the endocrine mechanisms controlling oogenesis in fish and the effects of growth on this process are poorly understood. In this study interactions between the growth and reproductive systems were examined by monitoring changes in various components of the FSH-ovary axis, plasma insulin-like growth factor 1 (Igf1), and ovarian gene expression in relation to body and previtellogenic oocyte growth in coho salmon. Samples were collected from females during two hypothesized critical periods when growth influences maturation in this species. Body growth during the fall-spring months was strongly related to the degree of oocyte development, with larger fish possessing more advanced oocytes than smaller, slower growing fish. The accumulation of cortical alveoli in the oocytes was associated with increases in plasma and pituitary FSH, plasma estradiol-17beta, and ovarian steroidogenic acute regulatory protein (star) gene expression, whereas ovarian transcripts for growth hormone receptor and somatolactin receptor decreased. As oocytes accumulated lipid droplets, a general increase occurred in plasma Igf1 and components of the FSH-ovary axis, including plasma FSH, estradiol-17beta, and ovarian mRNAs for gonadotropin receptors, star, igf1, and igf2. A consistent positive relationship between plasma Igf1, estradiol-17beta, and pituitary FSH during growth in the spring suggests that these factors are important links in the mechanism by which body growth influences the rate of oocyte development.