Genetic variation within a stick-insect species associated with community-level traits.

Phenotypic variation within species can affect the ecological dynamics of populations and communities. Characterizing the genetic variation underlying such effects can help parse the roles of genetic evolution and plasticity in "eco-evolutionary dynamics" and inform how genetic variation may shape patterns of evolution. Here, we employ genome-wide association (GWA) methods in Timema cristinae stick insects and their co-occurring arthropod communities to identify genetic variation associated with community-level traits. Previous studies have shown that maladaptation (i.e., imperfect crypsis) of T. cristinae can reduce the abundance and species richness of other arthropods due to an increase in bird predation. Whether genetic variation that is independent of crypsis has similar effects is unknown and was tested here using genome-wide genotyping-by-sequencing data of stick insects, arthropod community information, and GWA mapping with Bayesian sparse linear mixed models. We find associations between genetic variation in stick insects and arthropod community traits. However, these associations disappear when host-plant traits are accounted for. We thus use path analysis to disentangle interrelationships among stick-insect genetic variation, host-plant traits, and community traits. This revealed that host-plant size has large effects on arthropod communities, while genetic variation in stick insects has a smaller, but still significant effect. Our findings demonstrate that (1) genetic variation in a species can be associated with community-level traits but that (2) interrelationships among multiple factors may need to be analyzed to disentangle whether such associations represent causal relationships. This work helps to build a framework for genomic studies of eco-evolutionary dynamics.

A pituitary gene network linking vgll3 to regulators of sexual maturation in male Atlantic salmon.

Age at maturity is a key life history trait and a significant contributor to life history strategy variation. The maturation process is influenced by genetic and environmental factors, but specific causes of variation in maturation timing remain elusive. In many species, the increase in the regulatory gonadotropin-releasing hormone 1 (GnRH1) marks the onset of puberty. Atlantic salmon, however, lacks the gnrh1 gene, suggesting gnrh3 and/or other regulatory factors are involved in the maturation process. Earlier research in Atlantic salmon has found a strong association between alternative alleles of vgll3 and maturation timing. Recently we reported strong induction of gonadotropin genes (fshb and lhb) in the pituitary of Atlantic salmon homozygous for the early maturation allele (E) of vgll3. The induction of gonadotropins was accompanied by increased expression of their direct upstream regulators, c-jun and sf1 (nr5a1b) but the regulatory connection between vgll3 and these regulators has never been investigated in any organism. In this study, we investigated the potential regulatory connection between vgll3 genotypes and these regulators through a stepwise approach of identifying a gene regulatory network (GRN) containing c-jun and sf1, and transcription factor motif enrichment analysis. We found a GRN containing c-jun with predicted upstream regulators, e2f1, egr1, foxj1 and klf4, to be differentially expressed in the pituitary. Finally, we suggest a vgll3 and Hippo pathway -dependent model for transcriptional regulation of c-jun and sf1 in the pituitary, which may have broader implications across vertebrates.

Life-history genotype explains variation in migration activity in Atlantic salmon (Salmo salar).

One of the most well-known life-history continuums is the fast-slow axis, where 'fast' individuals mature earlier than 'slow' individuals. 'Fast' individuals are predicted to be more active than 'slow' individuals because high activity is required to maintain a fast life-history strategy. Recent meta-analyses revealed mixed evidence for such integration. Here, we test whether known life-history genotypes differ in activity expression by using Atlantic salmon (Salmo salar) as a model. In salmon, variation in Vgll3, a transcription cofactor, explains approximately 40% of variation in maturation timing. We predicted that the allele related to early maturation (vgll3*E) would be associated with higher activity. We used an automated surveillance system to follow approximately 1900 juveniles including both migrants and non-migrants (i.e. smolt and parr fish, respectively) in semi-natural conditions over 31 days (approx. 580 000 activity measurements). In migrants, but not in non-migrants, vgll3 explained variation in activity according to our prediction in a sex-dependent manner. Specifically, in females the vgll3*E allele was related to increasing activity, whereas in males the vgll3*L allele (later maturation allele) was related to increasing activity. These sex-dependent effects might be a mechanism maintaining within-population genetic life-history variation.

Refining the genomic location of single nucleotide polymorphism variation affecting Atlantic salmon maturation timing at a key large-effect locus.

Efforts to understand the genetic underpinnings of phenotypic variation are becoming more and more frequent in molecular ecology. Such efforts often lead to the identification of candidate regions showing signals of association and/or selection. These regions may contain multiple genes and therefore validation of which genes are actually responsible for the signal is required. In Atlantic salmon (Salmo salar), a large-effect locus for maturation timing, an ecologically important trait, occurs in a genomic region including two genes, vgll3 and akap11, but data for clearly determining which of the genes (or both) contribute to the association have been lacking. Here, we take advantage of natural recombination events detected between the two candidate genes in a salmon broodstock to reduce linkage disequilibrium at the locus, thus enabling delineation of the influence of variation at these two genes on early maturation. By rearing 5,895 males to maturation age, of which 81% had recombinant vgll3/akap11 allelic combinations, we found that vgll3 single nucleotide polymorphism (SNP) variation was strongly associated with early maturation, whereas there was little or no association between akap11 SNP variation and early maturation. These findings provide strong evidence supporting vgll3 as the primary candidate gene in the chromosome 25 locus for influencing early maturation. This will help guide future research for understanding the genetic processes controlling early maturation. This also exemplifies the utility of natural recombinants to more precisely map causal variation underlying ecologically important phenotypic diversity.

Dissecting the loci underlying maturation timing in Atlantic salmon using haplotype and multi-SNP based association methods.

Characterizing the role of different mutational effect sizes in the evolution of fitness-related traits has been a major goal in evolutionary biology for a century. Such characterization in a diversity of systems, both model and non-model, will help to understand the genetic processes underlying fitness variation. However, well-characterized genetic architectures of such traits in wild populations remain uncommon. In this study, we used haplotype-based and multi-SNP Bayesian association methods with sequencing data for 313 individuals from wild populations to test the mutational composition of known candidate regions for sea age at maturation in Atlantic salmon (Salmo salar). We detected an association at five loci out of 116 candidates previously identified in an aquaculture strain with maturation timing in wild Atlantic salmon. We found that at four of these five loci, variation explained by the locus was predominantly driven by a single SNP suggesting the genetic architecture of this trait includes multiple loci with simple, non-clustered alleles and a locus with potentially more complex alleles. This highlights the diversity of genetic architectures that can exist for fitness-related traits. Furthermore, this study provides a useful multi-SNP framework for future work using sequencing data to characterize genetic variation underlying phenotypes in wild populations.

Strong regulatory effects of vgll3 genotype on reproductive axis gene expression in juvenile male Atlantic salmon.

Age at maturity is a major contributor to the diversity of life history strategies in organisms. The process of maturation is influenced by both genetics and the environment, and includes changes in levels of sex hormones and behavior, but the specific factors leading to variation in maturation timing are poorly understood. gnrh1 regulates the transcription of gonadotropin genes at pubertal onset in many species, but this gene is lacking in certain teleost species including Atlantic salmon (Salmo salar), which raises the possibility of the involvement of other important regulatory factors during this process. Earlier research has reported a strong association of alternative alleles of the vgll3 gene with maturation timing in Atlantic salmon, suggesting it as a potential candidate regulating reproductive axis genes. Here, we investigated the expression of reproductive axis genes in one-year-old Atlantic salmon males with immature gonads and different vgll3 genotypes during the spawning period. We detected strong vgll3 genotype-dependent differential expression of reproductive axis genes (such as fshb, lhb, amh and igf3) tested in the pituitary, and testis. In addition, we observed differential expression of jun (ap1) and nr5a1b (sf1), potential upstream regulators of gonadotropins in the pituitary, as well as axin2, id3, insl3, itch, ptgs2a and ptger4b, the downstream targets of amh and igf3 in the testis. Hereby, we provide evidence of strong vgll3 genotype-dependent transcriptional regulation of reproductive axis genes prior to sexual maturation and suggest alternative models for distinct actions of vgll3 genotypes on the related molecular processes.

Beyond large-effect loci: large-scale GWAS reveals a mixed large-effect and polygenic architecture for age at maturity of Atlantic salmon.

BACKGROUND: Understanding genetic architecture is essential for determining how traits will change in response to evolutionary processes such as selection, genetic drift and/or gene flow. In Atlantic salmon, age at maturity is an important life history trait that affects factors such as survival, reproductive success, and growth. Furthermore, age at maturity can seriously impact aquaculture production. Therefore, characterizing the genetic architecture that underlies variation in age at maturity is of key interest. RESULTS: Here, we refine our understanding of the genetic architecture for age at maturity of male Atlantic salmon using a genome-wide association study of 11,166 males from a single aquaculture strain, using imputed genotypes at 512,397 single nucleotide polymorphisms (SNPs). All individuals were genotyped with a 50K SNP array and imputed to higher density using parents genotyped with a 930K SNP array and pedigree information. We found significant association signals on 28 of 29 chromosomes (P-values: 8.7 × 10-133-9.8 × 10-8), including two very strong signals spanning the six6 and vgll3 gene regions on chromosomes 9 and 25, respectively. Furthermore, we identified 116 independent signals that tagged 120 candidate genes with varying effect sizes. Five of the candidate genes found here were previously associated with age at maturity in other vertebrates, including humans. DISCUSSION: These results reveal a mixed architecture of large-effect loci and a polygenic component that consists of multiple smaller-effect loci, suggesting a more complex genetic architecture of Atlantic salmon age at maturity than previously thought. This more complex architecture will have implications for selection on this key trait in aquaculture and for management of wild salmon populations.

Ancient chromosomal rearrangement associated with local adaptation of a postglacially colonized population of Atlantic Cod in the northwest Atlantic.

Intraspecific diversity is central to the management and conservation of exploited species, yet knowledge of how this diversity is distributed and maintained in the genome of many marine species is lacking. Recent advances in genomic analyses allow for genome-wide surveys of intraspecific diversity and offer new opportunities for exploring genomic patterns of divergence. Here, we analysed genome-wide polymorphisms to measure genetic differentiation between an offshore migratory and a nonmigratory population and to define conservation units of Atlantic Cod (Gadus morhua) in coastal Labrador. A total of 141 individuals, collected from offshore sites and from a coastal site within Gilbert Bay, Labrador, were genotyped using an ~11k single nucleotide polymorphism array. Analyses of population structure revealed strong genetic differentiation between migratory offshore cod and nonmigratory Gilbert Bay cod. Genetic differentiation was elevated for loci within a chromosomal rearrangement found on linkage group 1 (LG1) that coincides with a previously found double inversion associated with migratory and nonmigratory ecotype divergence of cod in the northeast Atlantic. This inverted region includes several genes potentially associated with adaptation to differences in salinity and temperature, as well as influencing migratory behaviour. Our work provides evidence that a chromosomal rearrangement on LG1 is associated with parallel patterns of divergence between migratory and nonmigratory ecotypes on both sides of the Atlantic Ocean.

Evolution repeats itself in replicate long-term studies in the wild.

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

Modular chromosome rearrangements reveal parallel and nonparallel adaptation in a marine fish.

Genomic architecture and standing variation can play a key role in ecological adaptation and contribute to the predictability of evolution. In Atlantic cod (Gadus morhua), four large chromosomal rearrangements have been associated with ecological gradients and migratory behavior in regional analyses. However, the degree of parallelism, the extent of independent inheritance, and functional distinctiveness of these rearrangements remain poorly understood. Here, we use a 12K single nucleotide polymorphism (SNP) array to demonstrate extensive individual variation in rearrangement genotype within populations across the species range, suggesting that local adaptation to fine-scale ecological variation is enabled by rearrangements with independent inheritance. Our results demonstrate significant association of rearrangements with migration phenotype and environmental gradients across the species range. Individual rearrangements exhibit functional modularity, but also contain loci showing multiple environmental associations. Clustering in genetic distance trees and reduced differentiation within rearrangements across the species range are consistent with shared variation as a source of contemporary adaptive diversity in Atlantic cod. Conversely, we also find that haplotypes in the LG12 and LG1 rearranged region have diverged across the Atlantic, despite consistent environmental associations. Exchange of these structurally variable genomic regions, as well as local selective pressures, has likely facilitated individual diversity within Atlantic cod stocks. Our results highlight the importance of genomic architecture and standing variation in enabling fine-scale adaptation in marine species.

A migration-associated supergene reveals loss of biocomplexity in Atlantic cod.

Chromosome structural variation may underpin ecologically important intraspecific diversity by reducing recombination within supergenes containing linked, coadapted alleles. Here, we confirm that an ancient chromosomal rearrangement is strongly associated with migratory phenotype and individual genetic structure in Atlantic cod (Gadus morhua) across the Northwest Atlantic. We reconstruct trends in effective population size over the last century and reveal declines in effective population size matching onset of industrialized harvest (after 1950). We find different demographic trajectories between individuals homozygous for the chromosomal rearrangement relative to heterozygous or homozygous individuals for the noninverted haplotype, suggesting different selective histories across the past 150 years. These results illustrate how chromosomal structural diversity can mediate fine-scale genetic, phenotypic, and demographic variation in a highly connected marine species and show how overfishing may have led to loss of biocomplexity within Northern cod stock.