Clock genes and their genomic distributions in three species of salmonid fishes: Associations with genes regulating sexual maturation and cell cycling

Rhythmic Clock Gene Expression in Atlantic Salmon Parr Brain

To better understand the complexity of clock genes in salmonids, a taxon with an additional whole genome duplication, an analysis was performed to identify and classify gene family members (clock, arntl, period, cryptochrome, nr1d, ror, and csnk1). The majority of clock genes, in zebrafish and Northern pike, appeared to be duplicated. In comparison to the 29 clock genes described in zebrafish, 48 clock genes were discovered in salmonid species. There was also evidence of species-specific reciprocal gene losses conserved to the Oncorhynchus sister clade. From the six period genes identified three were highly significantly rhythmic, and circadian in their expression patterns (per1a.1, per1a.2, per1b) and two was significantly rhythmically expressed (per2a, per2b). The transcriptomic study of juvenile Atlantic salmon (parr) brain tissues confirmed gene identification and revealed that there were 2,864 rhythmically expressed genes (p < 0.001), including 1,215 genes with a circadian expression pattern, of which 11 were clock genes. The majority of circadian expressed genes peaked 2 h before and after daylight. These findings provide a foundation for further research into the function of clock genes circadian rhythmicity and the role of an enriched number of clock genes relating to seasonal driven life history in salmonids.

Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing

BACKGROUND

Carp fish, rohu (Labeo rohita Ham.) is important freshwater aquaculture species of South-East Asia having seasonal reproductive rhythm. There is no holistic study at transcriptome level revealing key candidate genes involved in such circannual rhythm regulated by biological clock genes (BCGs). Seasonality manifestation has two contrasting phases of reproduction, i.e., post-spawning resting and initiation of gonadal activity appropriate for revealing the associated candidate genes. It can be deciphered by RNA sequencing of tissues involved in BPGL (Brain-Pituitary-Gonad-Liver) axis controlling seasonality. How far such BCGs of this fish are evolutionarily conserved across different phyla is unknown. Such study can be of further use to enhance fish productivity as seasonality restricts seed production beyond monsoon season.

RESULT

A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were generated using Illumina TruSeq. De-novo assembled BPGL tissues revealed 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs, 514 pathways, 5379 transcription factors, 187 mature miRNA which regulates candidate genes represented by 1576 differentially expressed transcripts are available in the form of web-genomic resources. Findings were validated by qPCR. This is the first report in carp fish having 32 BCGs, found widely conserved in fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is due to universal mechanism of rhythmicity in response to environment and earth rotation having adaptive and reproductive significance.

CONCLUSION

This study elucidates evolutionary conserved mechanism of photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish. Study reveals fish as a good model for research on biological clock besides its relevance in reproductive efficiency enhancement.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

DeepHBV: a deep learning model to predict hepatitis B virus (HBV) integration sites

Background

The hepatitis B virus (HBV) is one of the main causes of viral hepatitis and liver cancer. HBV integration is one of the key steps in the virus-promoted malignant transformation.

Results

An attention-based deep learning model, DeepHBV, was developed to predict HBV integration sites. By learning local genomic features automatically, DeepHBV was trained and tested using HBV integration site data from the dsVIS database. Initially, DeepHBV showed an AUROC of 0.6363 and an AUPR of 0.5471 for the dataset. The integration of genomic features of repeat peaks and TCGA Pan-Cancer peaks significantly improved model performance, with AUROCs of 0.8378 and 0.9430 and AUPRs of 0.7535 and 0.9310, respectively. The transcription factor binding sites (TFBS) were significantly enriched near the genomic positions that were considered. The binding sites of the AR-halfsite, Arnt, Atf1, bHLHE40, bHLHE41, BMAL1, CLOCK, c-Myc, COUP-TFII, E2A, EBF1, Erra, and Foxo3 were highlighted by DeepHBV in both the dsVIS and VISDB datasets, revealing a novel integration preference for HBV.

Conclusions

DeepHBV is a useful tool for predicting HBV integration sites, revealing novel insights into HBV integration-related carcinogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01869-8.

Size-selective mortality induces evolutionary changes in group risk-taking behaviour and the circadian system in a fish

Intensive and trait-selective mortality of fish and wildlife can cause evolutionary changes in a range of life-history and behavioural traits. These changes might in turn alter the circadian system due to co-evolutionary mechanisms or correlated selection responses both at behavioural and molecular levels, with knock-on effects on daily physiological processes and behavioural outputs. We examined the evolutionary impact of size-selective mortality on group risk-taking behaviour and the circadian system in a model fish species. We exposed zebrafish Danio rerio to either large or small size-selective harvesting relative to a control over five generations, followed by eight generations during which harvesting was halted to remove maternal effects. Size-selective mortality affected fine-scale timing of behaviours. In particular, small size-selective mortality, typical of specialized fisheries and gape-limited predators targeting smaller size classes, increased group risk-taking behaviuor during feeding and after simulated predator attacks. Moreover, small size-selective mortality increased early peaks of daily activity as well as extended self-feeding daily activity to the photophase compared to controls. By contrast large size-selective mortality, typical of most wild capture fisheries, only showed an almost significant effect of decreasing group risk-taking behaviour during the habituation phase and no clear changes in fine-scale timing of daily behavioural rhythms compared to controls. We also found changes in the molecular circadian core clockwork in response to both size-selective mortality treatments. These changes disappeared in the clock output pathway because both size-selected lines showed similar transcription profiles. This switch downstream to the molecular circadian core clockwork also resulted in similar overall behavioural rhythms (diurnal swimming and self-feeding in the last hours of darkness) independent of the underlying molecular clock. To conclude, our experimental harvest left an asymmetrical evolutionary legacy in group risk-taking behaviour and in fine-scale daily behavioural rhythms. Yet, the overall timing of activity showed evolutionary resistance probably maintained by a molecular switch. Our experimental findings suggest that size-selective mortality can have consequences for behaviour and physiological processes.

The Molecular Evolution of Circadian Clock Genes in Spotted Gar (Lepisosteus oculatus)

Circadian rhythms are biological rhythms with a period of approximately 24 h. While canonical circadian clock genes and their regulatory mechanisms appear highly conserved, the evolution of clock gene families is still unclear due to several rounds of whole genome duplication in vertebrates. The spotted gar (Lepisosteus oculatus), as a non-teleost ray-finned fish, represents a fish lineage that diverged before the teleost genome duplication (TGD), providing an outgroup for exploring the evolutionary mechanisms of circadian clocks after whole-genome duplication. In this study, we interrogated the spotted gar draft genome sequences and found that spotted gar contains 26 circadian clock genes from 11 families. Phylogenetic analysis showed that 9 of these 11 spotted gar circadian clock gene families have the same number of genes as humans, while the members of the nfil3 and cry families are different between spotted gar and humans. Using phylogenetic and syntenic analyses, we found that nfil3-1 is conserved in vertebrates, while nfil3-2 and nfil3-3 are maintained in spotted gar, teleost fish, amphibians, and reptiles, but not in mammals. Following the two-round vertebrate genome duplication (VGD), spotted gar retained cry1a, cry1b, and cry2, and cry3 is retained in spotted gar, teleost fish, turtles, and birds, but not in mammals. We hypothesize that duplication of core clock genes, such as (nfil3 and cry), likely facilitated diversification of circadian regulatory mechanisms in teleost fish. We also found that the transcription factor binding element (Ahr::Arnt) is retained only in one of the per1 or per2 duplicated paralogs derived from the TGD in the teleost fish, implicating possible subfuctionalization cases. Together, these findings help decipher the repertoires of the spotted gar’s circadian system and shed light on how the vertebrate circadian clock systems have evolved.

Teleosts Genomics: Progress and Prospects in Disease Prevention and Control

Genome wide studies based on conventional molecular tools and upcoming omics technologies are beginning to gain functional applications in the control and prevention of diseases in teleosts fish. Herein, we provide insights into current progress and prospects in the use genomics studies for the control and prevention of fish diseases. Metagenomics has emerged to be an important tool used to identify emerging infectious diseases for the timely design of rational disease control strategies, determining microbial compositions in different aquatic environments used for fish farming and the use of host microbiota to monitor the health status of fish. Expounding the use of antimicrobial peptides (AMPs) as therapeutic agents against different pathogens as well as elucidating their role in tissue regeneration is another vital aspect of genomics studies that had taken precedent in recent years. In vaccine development, prospects made include the identification of highly immunogenic proteins for use in recombinant vaccine designs as well as identifying gene signatures that correlate with protective immunity for use as benchmarks in optimizing vaccine efficacy. Progress in quantitative trait loci (QTL) mapping is beginning to yield considerable success in identifying resistant traits against some of the highly infectious diseases that have previously ravaged the aquaculture industry. Altogether, the synopsis put forth shows that genomics studies are beginning to yield positive contribution in the prevention and control of fish diseases in aquaculture.

Cold temperature represses daily rhythms in the liver transcriptome of a stenothermal teleost under decreasing day length

The climate-change-driven increase in temperature is occurring rapidly and decreasing the predictability of seasonal rhythms at high latitudes. It is therefore urgent to understand how a change in the relationship between photoperiod and temperature can affect ectotherms in these environments. We tested whether temperature affects daily rhythms of transcription in a cold-adapted salmonid using high-throughput RNA sequencing. Arctic char (Salvelinus alpinus) from a subarctic population were reared at a high and a low temperature (15 and 8°C) for 1 month under natural, decreasing day length during late summer. Liver transcriptomes were compared between samples collected in the middle and towards the end of the light period and in the middle of the dark period. Daily variation in transcription was lower in fish from the low temperature compared with strong daily variation in warm-acclimated fish, suggesting that cold temperatures dampen the cycling of transcriptional rhythms under a simultaneously decreasing day length. Different circadian clock genes had divergent expression patterns, responding either by decreased expression or by increased rhythmicity at 15°C compared with 8°C. The results point out mechanisms that can affect the ability of fish to adapt to increasing temperatures caused by climate change.

Expressional changes of AMH signaling system in the quail testis induced by photoperiod

Gonadal sex differentiation proceeds by the interplay of various genes including the transcription factors and secretory factors in a complex network. The sex-differentiating genes are expressed not only during early sex differentiation but also throughout the gonadal development and even in the adult gonads. In addition, the evidence that they actually function in the adult gonads have been accumulated from the studies using the conditional knockout mice. However, many previous studies were focused on one single gene though those genes function in a network. In this study, the expressions of various sex-differentiating genes were analyzed simultaneously in the adult testis of the Japanese quail (Coturnix japonica), whose testicular functions are dramatically changed by altering the photoperiod, to elucidate the roles of them in the adult gonad. Anti-Müllerian hormone (AMH) was significantly upregulated in the regressed testis induced by the short-day condition. The expressions of the transcription factors that promoteAMHexpression in mammals (SF1,SOX9,WT1andGATA4) were also increased in the regressed testis. Moreover, AMH receptor (AMHR2) showed similar expression pattern to its ligand. We also analyzed the expressions of other transforming growth factor beta (TGFB) superfamily members and their receptors. The expressions of the ligands and receptors of TGFB family, and follistatin and betaglycan in addition to inhibin subunits were increased in the regressed testis. These results suggest that AMH is involved in the adult testicular functions of the Japanese quail together with other TGFB superfamily members.

Quantifying the intraindividual variation of antimüllerian hormone in the ovarian cycle

OBJECTIVE

To quantify intraindividual variability of antimüllerian hormone (AMH) as analytical and biological coefficients of variation and assess the effects of variation on clinical classification.

DESIGN

Retrospective cohort study.

SETTING

Not applicable.

PATIENT(S)

Thirty-eight women referred by general practitioners.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

Total intraindividual variability (CV_W), analytical (CV_A) and biological variability (CV_I) for each woman and for AMH ranges: low (<5 pmol/L), reduced (5-10), moderate (>10-30) and high (>30 pmol/L), with calculation of proportion of women crossing clinical cutoffs and expected variability around each cutoff.

RESULT(S)

Cycling women (n = 38) contributed 238 blood samples (average 6 samples each). The average total intraindividual AMH variability was 20% (range: 2.1% to 73%). Biological variation was 19% (range: 0 to 71%) and at least twice the analytical variation of 6.9% (range: 4.5% to 16%). Reclassification rates were highest in women with low (33%) or reduced AMH (67%) levels. Expected variations around the 5, 10, and 30 pmol/L cutoffs were 3-7, 7-13, and 20-40 pmol/L, respectively. In a woman with mean AMH in the 10-30 pmol/L range, the span of results that could occur was 7-40 pmol/L.

CONCLUSION(S)

Total variation in AMH was 20%, and the majority of this was biological. Changes in AMH resulted in reclassification in 29% of women and occurred most frequently in those with low and reduced AMH. In cycling women, the variability in AMH should be considered by clinicians, especially if a result is close to a clinical cutoff.

Integrated gene mapping and synteny studies give insights into the evolution of a sex proto-chromosome in Solea senegalensis

The evolution of genes related to sex and reproduction in fish shows high plasticity and, to date, the sex determination system has only been identified in a few species. Solea senegalensis has 42 chromosomes and an XX/XY chromosome system for sex determination, while related species show the ZZ/ZW system. Next-generation sequencing (NGS), multi-color fluorescence in situ hybridization (mFISH) techniques, and bioinformatics analysis have been carried out, with the objective of revealing new information about sex determination and reproduction in S. senegalensis. To that end, several bacterial artificial chromosome (BAC) clones that contain candidate genes involved in such processes (dmrt1, dmrt2, dmrt3, dmrt4, sox3, sox6, sox8, sox9, lh, cyp19a1a, amh, vasa, aqp3, and nanos3) were analyzed and compared with the same region in other related species. Synteny studies showed that the co-localization of dmrt1-dmrt2-drmt3 in the largest metacentric chromosome of S. senegalensis is coincident with that found in the Z chromosome of Cynoglossus semilaevis, which would potentially make this a sex proto-chromosome. Phylogenetic studies show the close proximity of S. senegalensis to Oryzias latipes, a species with an XX/XY system and a sex master gene. Comparative mapping provides evidence of the preferential association of these candidate genes in particular chromosome pairs. By using the NGS and mFISH techniques, it has been possible to obtain an integrated genetic map, which shows that 15 out of 21 chromosome pairs of S. senegalensis have at least one BAC clone. This result is important for distinguishing those chromosome pairs of S. senegalensis that are similar in shape and size. The mFISH analysis shows the following co-localizations in the same chromosomes: dmrt1-dmrt2-dmrt3, dmrt4-sox9-thrb, aqp3-sox8, cyp19a1a-fshb, igsf9b-sox3, and lysg-sox6.

The role of Amh signaling in teleost fish--Multiple functions not restricted to the gonads

This review summarizes the important role of Anti-Müllerian hormone (Amh) during gonad development in fishes. This Tgfβ-domain bearing hormone was named after one of its known functions, the induction of the regression of Müllerian ducts in male mammalian embryos. Later in development it is involved in male and female gonad differentiation and extragonadal expression has been reported in mammals as well. Teleosts lack Müllerian ducts, but they have amh orthologous genes. amh expression is reported from 21 fish species and possible regulatory interactions with further factors like sex steroids and gonadotropic hormones are discussed. The gonadotropin Fsh inhibits amh expression in all fish species studied. Sex steroids show no consistent influence on amh expression. Amh is produced in male Sertoli cells and female granulosa cells and inhibits germ cell proliferation and differentiation as well as steroidogenesis in both sexes. Therefore, Amh might be a central player in gonad development and a target of gonadotropic Fsh. Furthermore, there is evidence that an Amh-type II receptor is involved in germ cell regulation. Amh and its corresponding type II receptor are also present in brain and pituitary, at least in some teleosts, indicating additional roles of Amh effects in the brain-pituitary-gonadal axis. Unraveling Amh signaling is important in stem cell research and for reproduction as well as for aquaculture and in environmental science.

Genome-wide association study (GWAS) for growth rate and age at sexual maturation in Atlantic salmon (Salmo salar)

Early sexual maturation is considered a serious drawback for Atlantic salmon aquaculture as it retards growth, increases production times and affects flesh quality. Although both growth and sexual maturation are thought to be complex processes controlled by several genetic and environmental factors, selection for these traits has been continuously accomplished since the beginning of Atlantic salmon selective breeding programs. In this genome-wide association study (GWAS) we used a 6.5K single-nucleotide polymorphism (SNP) array to genotype ∼480 individuals from the Cermaq Canada broodstock program and search for SNPs associated with growth and age at sexual maturation. Using a mixed model approach we identified markers showing a significant association with growth, grilsing (early sexual maturation) and late sexual maturation. The most significant associations were found for grilsing, with markers located in Ssa10, Ssa02, Ssa13, Ssa25 and Ssa12, and for late maturation with markers located in Ssa28, Ssa01 and Ssa21. A lower level of association was detected with growth on Ssa13. Candidate genes, which were linked to these genetic markers, were identified and some of them show a direct relationship with developmental processes, especially for those in association with sexual maturation. However, the relatively low power to detect genetic markers associated with growth (days to 5 kg) in this GWAS indicates the need to use a higher density SNP array in order to overcome the low levels of linkage disequilibrium observed in Atlantic salmon before the information can be incorporated into a selective breeding program.

Divergence in expression of candidate genes for the smoltification process between juvenile resident rainbow and anadromous steelhead trout

Rainbow and steelhead trout (Oncorhynchus mykiss), among other salmonid fishes, exhibit tremendous life history diversity, foremost of which is variation in migratory propensity. While some individuals possess the ability to undertake an anadromous marine migration, others remain resident in freshwater throughout their life cycle. Those that will migrate undergo tremendous physiological, morphological, and behavioral transformations in a process called smoltification which transitions freshwater-adapted parr to marine-adapted smolts. While the behavior, ecology, and physiology of smoltification are well described, our understanding of the proximate genetic mechanisms that trigger the process are not well known. Quantitative genetic analyses have identified several genomic regions associated with smoltification and migration-related traits within this species. Here we investigate the divergence in gene expression of 18 functional and positional candidate genes for the smoltification process in the brain, gill, and liver tissues of migratory smolts, resident parr, and precocious mature male trout at the developmental stage of out-migration. Our analysis reveals several genes differentially expressed between life history classes and validates the candidate nature of several genes in the parr-smolt transformation including Clock1α, FSHβ, GR, GH2, GHR1, GHR2, NDK7, p53, SC6a7, Taldo1, THRα, THRβ, and Vdac2.

Molecular markers for variation in spawning date in a hatchery population of rainbow trout (Oncorhynchus mykiss)

We examined the distribution of alleles at 63 microsatellite loci distributed across 29 linkage groups in broodstock females from a commercial population of rainbow trout spawning on different dates throughout the season (August to January). A total of 368 females, 184 and 117 females from each of the tail-ends of the spawning distribution and a subsample of 67 females spawning in the middle, were used to detect marker-trait associations. Twenty-one loci in a subset of genomic regions (RT-5, 7, 8, 10, 12, 14, 15, 22, 23, 24, 25, 29, 30, and 31) were significantly associated with variation in spawning date. Many of these markers localize to regions with known spawning date quantitative trait loci based on previous studies. An individual assignment analysis was used to test how well the molecular data could be used to assign individuals to their correct spawning group, and markers were given a ranking reflecting their contribution to the accuracy of assignment. The top 15 ranked markers were successful at assigning the majority of females to the correct spawning group based on genotype with an average accuracy of 76 %. The most likely genes that could contribute to these differences in spawning date are discussed. Together, these data indicate that the loci could be incorporated into a selection index with phenotype data to increase the accuracy of selection for spawning date.

Mapping and Expression of Candidate Genes for Development Rate in Rainbow Trout (Oncorhynchus mykiss)

Development rate has important implications for individual fitness and physiology. In salmonid fishes, development rate correlates with many traits later in life, including life-history diversity, growth, and age and size at sexual maturation. In rainbow trout (Oncorhynchus mykiss), a quantitative trait locus for embryonic development rate has been detected on chromosome 5 across populations. However, few candidate genes have been identified within this region. In this study, we use gene mapping, gene expression, and quantitative genetic methods to further identify the genetic basis of embryonic developmental rate in O. mykiss Among the genes located in the region of the major development rate quantitative trait locus (GHR1, Clock1a, Myd118-1, and their paralogs), all were expressed early in embryonic development (fertilization through hatch), but none were differentially expressed between individuals with the fast- or slow-developing alleles for a major embryonic development rate quantitative trait locus. In a follow-up study of migratory and resident rainbow trout from natural populations in Alaska, we found significant additive variation in development rate and, moreover, found associations between development rate and allelic variation in all 3 candidate genes within the quantitative trait locus for embryonic development. The mapping of these genes to this region and associations in multiple populations provide positional candidates for further study of their roles in growth, development, and life-history diversity in this model salmonid.

Rapid parallel evolution of standing variation in a single, complex, genomic region is associated with life history in steelhead/rainbow trout

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. Such changes may be population-specific or, alternatively, may involve parallel evolution of the same genomic region in multiple populations, if that region contains genes or co-adapted gene complexes affecting the selected trait(s). Both quantitative and population genetic approaches have identified associations between specific genomic regions and the anadromous (steelhead) and resident (rainbow trout) life-history strategies of Oncorhynchus mykiss. Here, we use genotype data from 95 single nucleotide polymorphisms and show that the distribution of variation in a large region of one chromosome, Omy5, is strongly associated with life-history differentiation in multiple above-barrier populations of rainbow trout and their anadromous steelhead ancestors. The associated loci are in strong linkage disequilibrium, suggesting the presence of a chromosomal inversion or other rearrangement limiting recombination. These results provide the first evidence of a common genomic basis for life-history variation in O. mykiss in a geographically diverse set of populations and extend our knowledge of the heritable basis of rapid adaptation of complex traits in novel habitats.

Genotypic sex determination in teleosts: insights from the testis-determining amhy gene

The master sex-determining genes identified so far in fishes are clearly not conserved, as evidenced by several unrelated genes reported to play critical roles in sex determination. In this study, we reviewed the molecular process of sex determination in the Patagonian pejerrey Odontesthes hatcheri, an emerging model due to the recent discovery that a Y-chromosome linked, duplicated copy of the anti-Müllerian hormone gene, amhy plays a pivotal role in sex determination. A comparative analysis with other newly found sex-determining genes of teleost fish, DMY/dmrt1bY, sdY, amhr2, and gsdf(Y) is performed and alternative ideas are proposed to explain the mechanism involved in the rise of various types of non-homologous sex-determining genes.

A Y-linked anti-Müllerian hormone duplication takes over a critical role in sex determination

Gonadal sex determination in vertebrates generally follows a sequence of genetically programmed events. In what is seemingly becoming a pattern, all confirmed or current candidate "master" sex-determining genes reported in this group, e.g., SRY in eutherian mammals, DMY/dmrt1bY in medaka, DM-W in the African clawed frog, and DMRT1 in chicken encode transcription factors. In contrast, here we show that a male-specific, duplicated copy of the anti-Müllerian hormone (amh) is implicated in testicular development of the teleost fish Patagonian pejerrey (Odontesthes hatcheri). The gene, termed amhy because it is found in a single metacentric/submetacentric chromosome of XY individuals, is expressed much earlier than the autosomal amh (6 d after fertilization vs. 12 wk after fertilization) and is localized to presumptive Sertoli cells of XY males during testicular differentiation. Moreover, amhy knockdown in XY embryos resulted in the up-regulation of foxl2 and cyp19a1a mRNAs and the development of ovaries. These results are evidence of a functional amh duplication in vertebrates and suggest that amhy may be the master sex-determining gene in this species. If confirmed, this would be a unique instance of a hormone-related gene, a member of the TGF-β superfamily, in such a role.

Genetic architecture of body weight, condition factor and age of sexual maturation in Icelandic Arctic charr (Salvelinus alpinus)

The high commercial value from the aquaculture of salmonid fishes has prompted many studies into the genetic architecture of complex traits and the need to identify genomic regions that have repeatable associations with trait variation both within and among species. We searched for quantitative trait loci (QTL) for body weight (BW), condition factor (CF) and age of sexual maturation (MAT) in families of Arctic charr (Salvelinus alpinus) from an Icelandic breeding program. QTL with genome-wide significance were detected for each trait on multiple Arctic charr (AC) linkage groups (BW: AC-4, AC-20; CF: AC-7, AC-20, AC-23, AC-36; MAT: AC-13/34, AC-39). In addition to the genome-wide significant QTL for both BW and CF on AC-20, linkage groups AC-4, AC-7, AC-8, and AC-16 contain QTL for both BW and CF with chromosome-wide significance. These regions had effects (albeit weaker) on MAT with the exception of the region on AC-8. Comparisons with a North American cultured strain of Arctic charr, as well as North American populations of Atlantic salmon (Salmo salar), and rainbow trout (Oncorhynchus mykiss), reveal some conservation in QTL location and structure, particularly with respect to the joint associations of QTL influencing BW and CF. The detection of some differences in genetic architecture between the two aquaculture strains of Arctic charr may be reflective of the differential evolutionary histories experienced by these fishes, and illustrates the importance of including different strains to investigate genetic variation in a species where the intent is to use that variation in selective breeding programs.