Finding the adaptive needles in a population-structured haystack: a case study in a New Zealand mollusc

Genetic adaptation to future environmental conditions is crucial to help species persist as the climate changes. Genome scans are powerful tools to understand adaptive landscapes, enabling us to correlate genetic diversity with environmental gradients while disentangling neutral from adaptive variation. However, low gene flow can lead to both local adaptation and highly structured populations, and is a major confounding factor for genome scans, resulting in an inflated number of candidate loci.

Here, we compared candidate locus detection in a marine mollusc (Onithochiton neglectus), taking advantage of a natural geographical contrast in the levels of genetic structure between its populations. O. neglectus is endemic to New Zealand and distributed throughout an environmental gradient from the subtropical north to the subantarctic south. Due to a brooding developmental mode, populations tend to be locally isolated. However, adult hitchhiking on rafting kelp increases connectivity among southern populations.

We applied two genome scans for outliers (Bayescan and PCAdapt) and two genotype–environment association (GEA) tests (BayeScEnv and RDA). To limit issues with false positives, we combined results using the geometric mean of q‐values and performed association tests with random environmental variables. This novel approach is a compromise between stringent and relaxed approaches widely used before, and allowed us to classify candidate loci as low confidence or high confidence.

Genome scans for outliers detected a large number of significant outliers in strong and moderately structured populations. No high‐confidence GEA loci were detected in the context of strong population structure. However, 86 high‐confidence loci were associated predominantly with latitudinally varying abiotic factors in the less structured southern populations. This suggests that the degree of connectivity driven by kelp rafting over the southern scale may be insufficient to counteract local adaptation in this species.

Our study supports the expectation that genome scans may be prone to errors in highly structured populations. Nonetheless, it also empirically demonstrates that careful statistical controls enable the identification of candidate loci that invite more detailed investigations. Ultimately, genome scans are valuable tools to help guide further research aiming to determine the potential of non‐model species to adapt to future environments.

Genetic Regulation of Biomarkers as Stress Proxies in Dairy Cows

Stress in livestock reduces productivity and is a welfare concern. At a physiological level, stress is associated with the activation of inflammatory responses and increased levels of harmful reactive oxygen species. Biomarkers that are indicative of stress could facilitate the identification of more stress-resilient animals. We examined twenty-one metabolic, immune response, and liver function biomarkers that have been associated with stress in 416 Italian Simmental and 436 Italian Holstein cows which were genotyped for 150K SNPs. Single-SNP and haplotype-based genome-wide association studies were carried out to assess whether the variation in the levels in these biomarkers is under genetic control and to identify the genomic loci involved. Significant associations were found for the plasma levels of ceruloplasmin (Bos taurus chromosome 1-BTA1), paraoxonase (BTA4) and γ-glutamyl transferase (BTA17) in the individual breed analysis that coincided with the position of the genes coding for these proteins, suggesting that their expression is under cis-regulation. A meta-analysis of both breeds identified additional significant associations with paraoxonase on BTA 16 and 26. Finding genetic associations with variations in the levels of these biomarkers suggests that the selection for high or low levels of expression could be achieved rapidly. Whether the level of expression of the biomarkers correlates with the response to stressful situations has yet to be determined.

Trait differentiation and adaptation of plants along elevation gradients

Studies of genetic adaptation in plant populations along elevation gradients in mountains have a long history, but there has until now been neither a synthesis of how frequently plant populations exhibit adaptation to elevation nor an evaluation of how consistent underlying trait differences across species are. We reviewed studies of adaptation along elevation gradients (i) from a meta-analysis of phenotypic differentiation of three traits (height, biomass and phenology) from plants growing in 70 common garden experiments; (ii) by testing elevation adaptation using three fitness proxies (survival, reproductive output and biomass) from 14 reciprocal transplant experiments; (iii) by qualitatively assessing information at the molecular level, from 10 genomewide surveys and candidate gene approaches. We found that plants originating from high elevations were generally shorter and produced less biomass, but phenology did not vary consistently. We found significant evidence for elevation adaptation in terms of survival and biomass, but not for reproductive output. Variation in phenotypic and fitness responses to elevation across species was not related to life history traits or to environmental conditions. Molecular studies, which have focussed mainly on loci related to plant physiology and phenology, also provide evidence for adaptation along elevation gradients. Together, these studies indicate that genetically based trait differentiation and adaptation to elevation are widespread in plants. We conclude that a better understanding of the mechanisms underlying adaptation, not only to elevation but also to environmental change, will require more studies combining the ecological and molecular approaches.

Nonsyndromic cleft palate: An association study at GWAS candidate loci in a multiethnic sample

BACKGROUND

Nonsyndromic cleft palate only (nsCPO) is a common and multifactorial form of orofacial clefting. In contrast to successes achieved for the other common form of orofacial clefting, that is, nonsyndromic cleft lip with/without cleft palate (nsCL/P), genome wide association studies (GWAS) of nsCPO have identified only one genome wide significant locus. Aim of the present study was to investigate whether common variants contribute to nsCPO and, if so, to identify novel risk loci.

METHODS

We genotyped 33 SNPs at 27 candidate loci from 2 previously published nsCPO GWAS in an independent multiethnic sample. It included: (i) a family-based sample of European ancestry (n = 212); and (ii) two case/control samples of Central European (n = 94/339) and Arabian ancestry (n = 38/231), respectively. A separate association analysis was performed for each genotyped dataset, and meta-analyses were performed.

RESULTS

After association analysis and meta-analyses, none of the 33 SNPs showed genome-wide significance. Two variants showed nominally significant association in the imputed GWAS dataset and exhibited a further decrease in p-value in a European and an overall meta-analysis including imputed GWAS data, respectively (rs395572: P_MetaEU  = 3.16 × 10^-4 ; rs6809420: P_MetaAll  = 2.80 × 10^-4 ).

CONCLUSION

Our findings suggest that there is a limited contribution of common variants to nsCPO. However, the individual effect sizes might be too small for detection of further associations in the present sample sizes. Rare variants may play a more substantial role in nsCPO than in nsCL/P, for which GWAS of smaller sample sizes have identified genome-wide significant loci. Whole-exome/genome sequencing studies of nsCPO are now warranted.