PolarMorphism enables discovery of shared genetic variants across multiple traits from GWAS summary statistics

Pleiotropic SNPs are associated with multiple traits. Such SNPs can help pinpoint biological processes with an effect on multiple traits or point to a shared etiology between traits. We present PolarMorphism, a new method for the identification of pleiotropic SNPs from GWAS summary statistics. PolarMorphism can be readily applied to more than two traits or whole trait domains. PolarMorphism makes use of the fact that trait-specific SNP effect sizes can be seen as Cartesian coordinates and can thus be converted to polar coordinates r (distance from the origin) and theta (angle with the Cartesian x-axis). r describes the overall effect of a SNP, while theta describes the extent to which a SNP is shared. r and theta are used to determine the significance of SNP sharedness, resulting in a p-value per SNP that can be used for further analysis. We apply PolarMorphism to a large collection of publicly available GWAS summary statistics enabling the construction of a pleiotropy network that shows the extent to which traits share SNPs. This network shows how PolarMorphism can be used to gain insight into relationships between traits and trait domains. Furthermore, pathway analysis of the newly discovered pleiotropic SNPs demonstrates that analysis of more than two traits simultaneously yields more biologically relevant results than the combined results of pairwise analysis of the same traits. Finally, we show that PolarMorphism is more efficient and more powerful than previously published methods.

Evaluation of methods incorporating biological function and GWAS summary statistics to accelerate discovery

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 18 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

Genome-wide superior alleles, haplotypes and candidate genes associated with tolerance on sodic-dispersive soils in wheat (Triticum aestivum L.)

Key message The pleiotropic SNPs/haplotypes, overlapping genes (metal ion binding, photosynthesis), and homozygous/biallelic SNPs and transcription factors (HTH myb-type and BHLH) hold great potential for improving wheat yield potential on sodic-dispersive soils. Abstract Sodic-dispersive soils have multiple subsoil constraints including poor soil structure, alkaline pH and subsoil toxic elemental ion concentration, affecting growth and development in wheat. Tolerance is required at all developmental stages to enhance wheat yield potential on such soils. An in-depth investigation of genome-wide associations was conducted using a field phenotypic data of 206 diverse Focused Identification of Germplasm Strategy (FIGS) wheat lines for two consecutive years from different sodic and non-sodic plots and the exome targeted genotyping by sequencing (tGBS) assay. A total of 39 quantitative trait SNPs (QTSs), including 18 haplotypes were identified on chromosome 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 5A, 5D, 6B, 7A, 7B, 7D for yield and yield-components tolerance. Among these, three QTSs had common associations for multiple traits, indicating pleiotropism and four QTSs had close associations for multiple traits, within 32.38 Mb. The overlapping metal ion binding (Mn, Ca, Zn and Al) and photosynthesis genes and transcription factors (PHD-, Dof-, HTH myb-, BHLH-, PDZ_6-domain) identified are known to be highly regulated during germination, maximum stem elongation, anthesis, and grain development stages. The homozygous/biallelic SNPs having allele frequency above 30% were identified for yield and crop establishment/plants m−2. These SNPs correspond to HTH myb-type and BHLH transcription factors, brassinosteroid signalling pathway, kinase activity, ATP and chitin binding activity. These resources are valuable in haplotype-based breeding and genome editing to improve yield potential on sodic-dispersive soils.

Identification of Grain Size-Related QTLs in Korean japonica Rice Using Genome Resequencing and High-Throughput Image Analysis

Grain size is a key factor influencing the grain yield in rice. To identify the as-yet-unknown genes regulating grain size in Korean japonica rice, we developed a recombinant inbred line population (n = 162) from a cross between Odae (large-grain) and Joun (small-grain), and measured six traits including the thousand-grain weights of unhulled and hulled seeds, grain area, grain length, grain width and grain length-to-width ratio using high-throughput image analysis at the F8 and F9 generations. A genetic map was constructed using 248 kompetitive allele-specific PCR (KASP) markers that were polymorphic between the parental genotypes, and 29 QTLs affecting the six traits were identified, of which 15 were stable in both F8 and F9 generations. Notably, three QTL clusters affecting multiple traits were detected on chromosomes 6, 7 and 11. We analyzed whole-genome resequencing data of Odae and Joun, and selected candidate genes for the stable QTLs in the identified clusters that have high- or moderate-impact variations between Odae and Joun and encode proteins the families of which have been reported to be related to grain size regulation. These results will facilitate the identification of genes underlying the QTLs and promote molecular breeding of high-yielding Korean japonica rice varieties.

Dissecting supergene: evidence for recombination suppression among multiple functional loci within inversions.

Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between set of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining an landmark-free phenotyping algorithm with multivariate genome wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We showed that the supergene controlling the striking wing-pattern polymorphism displayed by this species contains many independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are associated with colour variations only in morphs controlled by inversions, indicating that they were recruited after the formation of these inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations.

Characterizing the genetic basis of trait evolution in the Mexican cavefish

Evolution in response to a change in ecology often coincides with various morphological, physiological, and behavioral traits. For most organisms little is known about the genetic and functional relationship between evolutionarily derived traits, representing a critical gap in our understanding of adaptation The Mexican tetra, Astyanax mexicanus, consists of largely independent populations of fish that inhabit at least 30 caves in Northeast Mexico, and a surface fish population, that inhabits the rivers of Mexico and Southern Texas. The recent application of molecular genetic approaches combined with behavioral phenotyping have established A. mexicanus as a model for studying the evolution of complex traits. Cave populations of A. mexicanus are interfertile with surface populations and have evolved numerous traits including eye degeneration, insomnia, albinism and enhanced mechanosensory function. The infertility of different populations from the same species provides a unique opportunity to define the genetic relationship between evolved traits and assess the co-evolution of behavioral and morphological traits with one another. To define the relationships between morphological and behavioral traits, we developed a pipeline to test individual fish for multiple traits. This pipeline confirmed differences in locomotor activity, prey capture, and startle reflex between surface and cavefish populations. To measure the relationship between traits, individual F2 hybrid fish were characterized for locomotor behavior, prey-capture behavior, startle reflex and morphological attributes. Analysis revealed an association between body length and slower escape reflex, suggesting a trade-off between increased size and predator avoidance in cavefish. Overall, there were few associations between individual behavioral traits, or behavioral and morphological traits, suggesting independent genetic changes underlie the evolution of behavioral and morphological traits. Taken together, this approach provides a novel system to identify genes that underlie naturally occurring genetic variation in morphological and behavioral traits.

Plasticity across levels: relating epigenomic, transcriptomic, and phenotypic responses to osmotic stress in a halotolerant microalga

Phenotypic plasticity, the ability of a given genotype to produce alternative phenotypes in response to its environment of development, is an important mechanism for coping with variable environments. While the mechanisms underlying phenotypic plasticity are diverse, their relative contributions need to be investigated quantitatively to better understand the evolvability of plasticity across biological levels. This requires relating plastic responses of the epigenome, transcriptome, and organismal phenotype, and how they vary with the genotype. Here we carried out this approach for responses to osmotic stress in Dunaliella salina, a green microalga that is a model organism for salinity tolerance. We compared two strains that show markedly different demographic responses to osmotic stress, and showed that these phenotypic responses involve strain- and environment-specific variation in gene expression levels, but a relative low - but significant - effect of strain x environment interaction. We also found an important genotype effect on the genome-wide methylation pattern, but little contribution from environmental conditions to the latter. However, we did detect a significant marginal effect of epigenetic variation on gene expression, beyond the influence of genetic differences on epigenetic state, and we showed that hypomethylated regions are correlated with higher gene expression. Our results indicate that epigenetic mechanisms are either not involved in the rapid plastic response to environmental change in this species, or involve only few changes in trans that are sufficient to trigger concerted changes in the expression of many genes, and phenotypic responses by multiple traits.

Leaf Trait Networks Based on Global Data: Representing Variation and Adaptation in Plants

The interdependence of multiple traits allows plants to perform multiple functions. Acquiring an accurate representation of the interdependence of plant traits could advance our understanding of the adaptative strategies of plants. However, few studies focus on complex relationships among multiple traits. Here, we proposed use of leaf trait networks (LTNs) to capture the complex relationships among traits, allowing us to visualize all relationships and quantify how they differ through network parameters. We established LTNs using six leaf economic traits. It showed that significant differences in LTNs of different life forms and growth forms. The trait relationships of broad-leaved trees were tighter than conifers; thus, broad-leaved trees could be more efficient than conifers. The trait relationships of shrubs were tighter than trees because shrubs require multiple traits to co-operate efficiently to perform multiple functions for thriving in limited resources. Furthermore, leaf nitrogen concentration and life span had the highest centrality in LTNs; consequently, the environmental selection of these two traits might impact the whole phenotype. In conclusion, LTNs are useful tools for identifying key traits and quantifying the interdependence of multiple traits.

Genome-Wide Superior Alleles, Haplotypes and Candidate Genes Associated With Tolerance on Sodic-Dispersive Soils in Wheat (Triticum Aestivum L.)

Sodic-dispersive soils have multiple subsoil constraints including poor soil structure, alkaline pH and subsoil toxic elemental ion concentration, affecting growth and development in wheat. Tolerance is required at all developmental stages to enhance wheat yield potential on such soils. An in-depth investigation of genome-wide associations was conducted using a field phenotypic data of 206 diverse Focused Identification of Germplasm Strategy (FIGS) wheat lines for two consecutive years from different sodic and non-sodic plots and the exome targeted genotyping by sequencing (tGBS) assay. A total of 39 quantitative trait SNPs (QTSs), including 18 haplotypes were identified on chromosome 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 5A, 5D, 6B, 7A, 7B, 7D for yield and yield-components tolerance. Among these, three QTSs had common associations for multiple traits, indicating pleiotropism and four QTSs had close associations for multiple traits, within 32.38 Mb. The overlapping metal ion binding (Mn, Ca, Zn and Al) and photosynthesis genes, and transcription factors (PHD-, Dof-, HTH myb-, BHLH-, PDZ_6-domain) identified are known to be highly regulated during germination, maximum stem elongation, anthesis, and grain development stages. The homozygous/biallelic SNPs having allele frequency above 30% were identified for yield and crop establishment/plants m-2. These SNPs correspond to HTH myb-type and BHLH transcription factors, brassinosteroid signaling pathway, kinase activity, ATP and chitin binding activity. These resources are valuable in haplotype-based breeding and genome editing to improve yield potential on sodic-dispersive soils.

Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with a lifetime risk of one in 350 people and an unmet need for disease-modifying therapies. We conducted a cross-ancestry genome-wide association study (GWAS) including 29,612 patients with ALS and 122,656 controls, which identified 15 risk loci. When combined with 8,953 individuals with whole-genome sequencing (6,538 patients, 2,415 controls) and a large cortex-derived expression quantitative trait locus (eQTL) dataset (MetaBrain), analyses revealed locus-specific genetic architectures in which we prioritized genes either through rare variants, short tandem repeats or regulatory effects. ALS-associated risk loci were shared with multiple traits within the neurodegenerative spectrum but with distinct enrichment patterns across brain regions and cell types. Of the environmental and lifestyle risk factors obtained from the literature, Mendelian randomization analyses indicated a causal role for high cholesterol levels. The combination of all ALS-associated signals reveals a role for perturbations in vesicle-mediated transport and autophagy and provides evidence for cell-autonomous disease initiation in glutamatergic neurons.