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

Identification of Single-Nucleotide Polymorphisms in Differentially Expressed Genes Favoring Soybean Meal Tolerance in Higher-Growth Zebrafish (Danio rerio)

Genetic variability within the same fish species could confer soybean meal (SBM) tolerance in some individuals, thus favoring growth. This study investigates the single-nucleotide polymorphisms (SNPs) in differentially expressed genes (DEGs) favoring SBM tolerance in higher-growth zebrafish (Danio rerio). In a previous work, nineteen families of zebrafish were fed a fish meal diet (100FM control diet) or SBM-based diets supplemented with saponin (50SBM + 2SPN-experimental diet), from juvenile to adult stages. Individuals were selected from families with a genotype-by-environment interaction higher (170 ± 18 mg) or lower (76 ± 10 mg) weight gain on 50SBM + 2SPN in relation to 100FM. Intestinal transcriptomic analysis using RNA-seq revealed six hundred and sixty-five differentially expressed genes in higher-growth fish fed 50SBM + 2SPN diet. In this work, using these results, 47 SNPs in DEGs were selected. These SNPs were genotyped by Sequenom in 340 zebrafish that were fed with a 50SBM + 2SPN diet or with 100FM diet. Marker-trait analysis revealed 4 SNPs associated with growth in 3 immunity-related genes (aif1l, arid3c, and cst14b.2) in response to the 50SBM + 2SPN diet (p-value < 0.05). Two SNPs belonging to aif1l y arid3c produce a positive (+19 mg) and negative (-26 mg) effect on fish growth, respectively. These SNPs can be used as markers to improve the early selection of tolerant fish to SBM diet or other plant-based diets. These genes can be used as biomarkers to identify SNPs in commercial fish, thus contributing to the aquaculture sustainability.

Population correlates of rapid captive-induced maladaptation in a wild fish

Understanding the extent to which captivity generates maladaptation in wild species can inform species recovery programs and elucidate wild population responses to novel environmental change. Although rarely quantified, effective population size (N e ) and genetic diversity should influence the magnitude of plastic and genetic changes manifested in captivity that reduce wild fitness. Sexually dimorphic traits might also mediate consequences of captivity. To evaluate these relationships, we generated >600 full- and half-sibling families from nine wild brook trout populations, reared them for one generation under common, captive environmental conditions and contrasted several fitness-related traits in wild versus captive lines. We found substantial variation in lifetime success (lifetime survival and reproductive success) and life history traits among wild populations after just one captive generation (fourteen- and threefold ranges across populations, respectively). Populations with lower heterozygosity showed lower captive lifetime success, suggesting that captivity generates maladaptation within one generation. Greater male-biased mortality in captivity occurred in populations having disproportionately higher growth rates in males than females. Wild population N e and allelic diversity had little or no influence on captive trait expression and lifetime success. Our results have four conservation implications: (i) Trait values and lifetime success were highly variable across populations following one generation of captivity. (ii) Maladaptation induced by captive breeding might be particularly intense for the very populations practitioners are most interested in conserving, such as those with low heterozygosity. (iii) Maladaptive sex differences in captivity might be associated with population-dependent growth costs of reproduction. (iv) Heterozygosity can be a good indicator of short-term, intraspecific responses to novel environmental change.

Sex-Specific Co-expression Networks and Sex-Biased Gene Expression in the Salmonid Brook Charr Salvelinus fontinalis

Networks of co-expressed genes produce complex phenotypes associated with functional novelty. Sex differences in gene expression levels or in the structure of gene co-expression networks can cause sexual dimorphism and may resolve sexually antagonistic selection. Here we used RNA-sequencing in the salmonid Brook Charr Salvelinus fontinalis to characterize sex-specific co-expression networks in the liver of 47 female and 53 male offspring. In both networks, modules were characterized for functional enrichment, hub gene identification, and associations with 15 growth, reproduction, and stress-related phenotypes. Modules were then evaluated for preservation in the opposite sex, and in the congener Arctic Charr Salvelinus alpinus Overall, more transcripts were assigned to a module in the female network than in the male network, which coincided with higher inter-individual gene expression and phenotype variation in the females. Most modules were preserved between sexes and species, including those involved in conserved cellular processes (e.g., translation, immune pathways). However, two sex-specific male modules were identified, and these may contribute to sexual dimorphism. To compare with the network analysis, differentially expressed transcripts were identified between the sexes, revealing a total of 16% of expressed transcripts as sex-biased. For both sexes, there was no overrepresentation of sex-biased genes or sex-specific modules on the putative sex chromosome. Sex-biased transcripts were also not overrepresented in sex-specific modules, and in fact highly male-biased transcripts were enriched in preserved modules. Comparative network analysis and differential expression analyses identified different aspects of sex differences in gene expression, and both provided new insights on the genes underlying sexual dimorphism in the salmonid Brook Charr.

Sex Chromosome Evolution, Heterochiasmy, and Physiological QTL in the Salmonid Brook Charr Salvelinus fontinalis

Whole-genome duplication (WGD) can have large impacts on genome evolution, and much remains unknown about these impacts. This includes the mechanisms of coping with a duplicated sex determination system and whether this has an impact on increasing the diversity of sex determination mechanisms. Other impacts include sexual conflict, where alleles having different optimums in each sex can result in sequestration of genes into nonrecombining sex chromosomes. Sex chromosome development itself may involve sex-specific recombination rate (i.e., heterochiasmy), which is also poorly understood. The family Salmonidae is a model system for these phenomena, having undergone autotetraploidization and subsequent rediploidization in most of the genome at the base of the lineage. The salmonid master sex determining gene is known, and many species have nonhomologous sex chromosomes, putatively due to transposition of this gene. In this study, we identify the sex chromosome of Brook Charr Salvelinus fontinalis and compare sex chromosome identities across the lineage (eight species and four genera). Although nonhomology is frequent, homologous sex chromosomes and other consistencies are present in distantly related species, indicating probable convergence on specific sex and neo-sex chromosomes. We also characterize strong heterochiasmy with 2.7-fold more crossovers in maternal than paternal haplotypes with paternal crossovers biased to chromosome ends. When considering only rediploidized chromosomes, the overall heterochiasmy trend remains, although with only 1.9-fold more recombination in the female than the male. Y chromosome crossovers are restricted to a single end of the chromosome, and this chromosome contains a large interspecific inversion, although its status between males and females remains unknown. Finally, we identify quantitative trait loci (QTL) for 21 unique growth, reproductive, and stress-related phenotypes to improve knowledge of the genetic architecture of these traits important to aquaculture and evolution.

Salmonid Chromosome Evolution as Revealed by a Novel Method for Comparing RADseq Linkage Maps

Whole genome duplication (WGD) can provide material for evolutionary innovation. Family Salmonidae is ideal for studying the effects of WGD as the ancestral salmonid underwent WGD relatively recently, ∼65 Ma, then rediploidized and diversified. Extensive synteny between homologous chromosome arms occurs in extant salmonids, but each species has both conserved and unique chromosome arm fusions and fissions. Assembly of large, outbred eukaryotic genomes can be difficult, but structural rearrangements within such taxa can be investigated using linkage maps. RAD sequencing provides unprecedented ability to generate high-density linkage maps for nonmodel species, but can result in low numbers of homologous markers between species due to phylogenetic distance or differences in library preparation. Here, we generate a high-density linkage map (3,826 markers) for the Salvelinus genera (Brook Charr S. fontinalis), and then identify corresponding chromosome arms among the other available salmonid high-density linkage maps, including six species of Oncorhynchus, and one species for each of Salmo, Coregonus, and the nonduplicated sister group for the salmonids, Northern Pike Esox lucius for identifying post-duplicated homeologs. To facilitate this process, we developed MapComp to identify identical and proximate (i.e. nearby) markers between linkage maps using a reference genome of a related species as an intermediate, increasing the number of comparable markers between linkage maps by 5-fold. This enabled a characterization of the most likely history of retained chromosomal rearrangements post-WGD, and several conserved chromosomal inversions. Analyses of RADseq-based linkage maps from other taxa will also benefit from MapComp, available at: https://github.com/enormandeau/mapcomp/

Population size, habitat fragmentation, and the nature of adaptive variation in a stream fish

Whether and how habitat fragmentation and population size jointly affect adaptive genetic variation and adaptive population differentiation are largely unexplored. Owing to pronounced genetic drift, small, fragmented populations are thought to exhibit reduced adaptive genetic variation relative to large populations. Yet fragmentation is known to increase variability within and among habitats as population size decreases. Such variability might instead favour the maintenance of adaptive polymorphisms and/or generate more variability in adaptive differentiation at smaller population size. We investigated these alternative hypotheses by analysing coding-gene, single-nucleotide polymorphisms associated with different biological functions in fragmented brook trout populations of variable sizes. Putative adaptive differentiation was greater between small and large populations or among small populations than among large populations. These trends were stronger for genetic population size measures than demographic ones and were present despite pronounced drift in small populations. Our results suggest that fragmentation affects natural selection and that the changes elicited in the adaptive genetic composition and differentiation of fragmented populations vary with population size. By generating more variable evolutionary responses, the alteration of selective pressures during habitat fragmentation may affect future population persistence independently of, and perhaps long before, the effects of demographic and genetic stochasticity are manifest.

Population size is weakly related to quantitative genetic variation and trait differentiation in a stream fish

How population size influences quantitative genetic variation and differentiation among natural, fragmented populations remains unresolved. Small, isolated populations might occupy poor quality habitats and lose genetic variation more rapidly due to genetic drift than large populations. Genetic drift might furthermore overcome selection as population size decreases. Collectively, this might result in directional changes in additive genetic variation (VA ) and trait differentiation (QST ) from small to large population size. Alternatively, small populations might exhibit larger variation in VA and QST if habitat fragmentation increases variability in habitat types. We explored these alternatives by investigating VA and QST using nine fragmented populations of brook trout varying 50-fold in census size N (179-8416) and 10-fold in effective number of breeders, Nb (18-135). Across 15 traits, no evidence was found for consistent differences in VA and QST with population size and almost no evidence for increased variability of VA or QST estimates at small population size. This suggests that (i) small populations of some species may retain adaptive potential according to commonly adopted quantitative genetic measures and (ii) populations of varying sizes experience a variety of environmental conditions in nature, however extremely large studies are likely required before any firm conclusions can be made.

Multidisciplinary population monitoring when demographic data are sparse: a case study of remote trout populations

The potential of genetic, genomic, and phenotypic metrics for monitoring population trends may be especially high in isolated regions, where traditional demographic monitoring is logistically difficult and only sporadic sampling is possible. This potential, however, is relatively underexplored empirically. Over eleven years, we assessed several such metrics along with traditional ecological knowledge and catch data in a socioeconomically important trout species occupying a large, remote lake. The data revealed largely stable characteristics in two populations over 2–3 generations, but possible contemporary changes in a third population. These potential shifts were suggested by reduced catch rates, reduced body size, and changes in selection implied at one gene-associated single nucleotide polymorphism. A demographic decline in this population, however, was ambiguously supported, based on the apparent lack of temporal change in effective population size, and corresponding traditional knowledge suggesting little change in catch. We illustrate how the pluralistic approach employed has practicality for setting future monitoring efforts of these populations, by guiding monitoring priorities according to the relative merits of different metrics and availability of resources. Our study also considers some advantages and disadvantages to adopting a pluralistic approach to population monitoring where demographic data are not easily obtained.

Stocking impacts the expression of candidate genes and physiological condition in introgressed brook charr (Salvelinus fontinalis) populations

Translocation of plants and animal populations between environments is one of the major forms of anthropogenic perturbation experienced by pristine populations, and consequently, human-mediated hybridization by stocking practices between wild and exogenous conspecifics is of increasing concern. In this study, we compared the expression of seven candidate genes involved in multifactorial traits and regulatory pathways for growth as a function of level of introgressive hybridization between wild and domestic brook charr to test the null hypothesis of no effect of introgression on wild fish. Our analyses revealed that the expression of two of the genes tested, cytochrome c oxidase VIIa and the growth hormone receptor isoform I, was positively correlated with the level of introgression. We also observed a positive relationship between the extent of introgression and physiological status quantified by the Fulton's condition index. The expression of other genes was influenced by other variables, including year of sampling (reflecting different thermal conditions), sampling method and lake of origin. This is the first demonstration in nature that introgression from stocked populations has an impact on the expression of genes playing a role in important biological functions that may be related with fitness in wild introgressed populations.

Localization of association signal from risk and protective variants in sequencing studies

Aggregating information across multiple variants in a gene or region can improve power for rare variant association testing. Power is maximized when the aggregation region contains many causal variants and few neutral variants. In this paper, we present a method for the localization of the association signal in a region using a sliding-window based approach to rare variant association testing in a region. We first introduce a novel method for analysis of rare variants, the Difference in Minor Allele Frequency test (DMAF), which allows combined analysis of common and rare variants, and makes no assumptions about the direction of effects. In whole-region analyses of simulated data with risk and protective variants, DMAF and other methods which pool data across individuals were found to outperform methods which pool data across variants. We then implement a sliding-window version of DMAF, using a step-down permutation approach to control type I error with the testing of multiple windows. In simulations, the sliding-window DMAF improved power to detect a causal sub-region, compared to applying DMAF to the whole region. Sliding-window DMAF was also effective in localizing the causal sub-region. We also applied the DMAF sliding-window approach to test for an association between response to the drug gemcitabine and variants in the gene FKBP5 sequenced in 91 lymphoblastoid cell lines derived from white non-Hispanic individuals. The application of the sliding-window test procedure detected an association in a sub-region spanning an exon and two introns, when rare and common variants were analyzed together.