A general test of association for quantitative traits in nuclear families

Causal interpretations of family GWAS in the presence of heterogeneous effects

Family-based genome-wide association studies (GWASs) are often claimed to provide an unbiased estimate of the average causal effects (or average treatment effects; ATEs) of alleles, on the basis of an analogy between the random transmission of alleles from parents to children and a randomized controlled trial. We show that this claim does not hold in general. Because Mendelian segregation only randomizes alleles among children of heterozygotes, the effects of alleles in the children of homozygotes are not observable. This feature will matter if an allele has different average effects in the children of homozygotes and heterozygotes, as can arise in the presence of gene-by-environment interactions, gene-by-gene interactions, or differences in linkage disequilibrium patterns. At a single locus, family-based GWAS can be thought of as providing an unbiased estimate of the average effect in the children of heterozygotes (i.e., a local average treatment effect; LATE). This interpretation does not extend to polygenic scores (PGSs), however, because different sets of SNPs are heterozygous in each family. Therefore, other than under specific conditions, the within-family regression slope of a PGS cannot be assumed to provide an unbiased estimate of the LATE for any subset or weighted average of families. In practice, the potential biases of a family-based GWAS are likely smaller than those that can arise from confounding in a standard, population-based GWAS, and so family studies remain important for the dissection of genetic contributions to phenotypic variation. Nonetheless, their causal interpretation is less straightforward than has been widely appreciated.

Interpreting population- and family-based genome-wide association studies in the presence of confounding

A central aim of genome-wide association studies (GWASs) is to estimate direct genetic effects: the causal effects on an individual's phenotype of the alleles that they carry. However, estimates of direct effects can be subject to genetic and environmental confounding and can also absorb the "indirect" genetic effects of relatives' genotypes. Recently, an important development in controlling for these confounds has been the use of within-family GWASs, which, because of the randomness of mendelian segregation within pedigrees, are often interpreted as producing unbiased estimates of direct effects. Here, we present a general theoretical analysis of the influence of confounding in standard population-based and within-family GWASs. We show that, contrary to common interpretation, family-based estimates of direct effects can be biased by genetic confounding. In humans, such biases will often be small per-locus, but can be compounded when effect-size estimates are used in polygenic scores (PGSs). We illustrate the influence of genetic confounding on population- and family-based estimates of direct effects using models of assortative mating, population stratification, and stabilizing selection on GWAS traits. We further show how family-based estimates of indirect genetic effects, based on comparisons of parentally transmitted and untransmitted alleles, can suffer substantial genetic confounding. We conclude that, while family-based studies have placed GWAS estimation on a more rigorous footing, they carry subtle issues of interpretation that arise from confounding.

Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer

INTRODUCTION

Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex.

METHODS

In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa).

RESULTS

Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer.

CONCLUSION

Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

Causal interpretations of family GWAS in the presence of heterogeneous effects

Family-based genome-wide association studies (GWAS) have emerged as a gold standard for assessing causal effects of alleles and polygenic scores. Notably, family studies are often claimed to provide an unbiased estimate of the average causal effect (or average treatment effect; ATE) of an allele, on the basis of an analogy between the random transmission of alleles from parents to children and a randomized controlled trial. Here, we show that this interpretation does not hold in general. Because Mendelian segregation only randomizes alleles among children of heterozygotes, the effects of alleles in the children of homozygotes are not observable. Consequently, if an allele has different average effects in the children of homozygotes and heterozygotes, as can arise in the presence of gene-by-environment interactions, gene-by-gene interactions, or differences in LD patterns, family studies provide a biased estimate of the average effect in the sample. At a single locus, family-based association studies can be thought of as providing an unbiased estimate of the average effect in the children of heterozygotes (i.e., a local average treatment effect; LATE). This interpretation does not extend to polygenic scores, however, because different sets of SNPs are heterozygous in each family. Therefore, other than under specific conditions, the within-family regression slope of a PGS cannot be assumed to provide an unbiased estimate for any subset or weighted average of families. Instead, family-based studies can be reinterpreted as enabling an unbiased estimate of the extent to which Mendelian segregation at loci in the PGS contributes to the population-level variance in the trait. Because this estimate does not include the between-family variance, however, this interpretation applies to only (roughly) half of the sample PGS variance. In practice, the potential biases of a family-based GWAS are likely smaller than those arising from confounding in a standard, population-based GWAS, and so family studies remain important for the dissection of genetic contributions to phenotypic variation. Nonetheless, the causal interpretation of family-based GWAS estimates is less straightforward than has been widely appreciated.

A Non-Coding Fc Gamma Receptor Cis-Regulatory Variant within the 1q23 Gene Cluster Is Associated with Plasmodium falciparum Infection in Children Residing in Burkina Faso

Antibodies play a crucial role in activating protective immunity against malaria by interacting with Fc-gamma receptors (FcγRs). Genetic variations in genes encoding FcγRs can affect immune cell responses to the parasite. In this study, our aim was to investigate whether non-coding variants that regulate FcγR expression could influence the prevalence of Plasmodium falciparum infection. Through bioinformatics approaches, we selected expression quantitative trait loci (eQTL) for FCGR2A, FCGR2B, FCGR2C, FCGR3A, and FCGR3B genes encoding FcγRs (FCGR), in whole blood. We prioritized two regulatory variants, rs2099684 and rs1771575, located in open genomic regions. These variants were identified using RegVar, ImmuNexUT, and transcription factor annotations specific to immune cells. In addition to these, we genotyped the coding variants FCGR2A/rs1801274 and FCGR2B/rs1050501 in 234 individuals from a malaria-endemic area in Burkina Faso. We conducted age and family-based analyses to evaluate associations with the prevalence of malarial infection in both children and adults. The analysis revealed that the regulatory rs1771575-CC genotype was predicted to influence FCGR2B/FCGR2C/FCGR3A transcripts in immune cells and was the sole variant associated with a higher prevalence of malarial infection in children. In conclusion, this study identifies the rs1771575 cis-regulatory variant affecting several FcγRs in myeloid and neutrophil cells and associates it with the inter-individual capacity of children living in Burkina Faso to control malarial infection.

A study of the genetic architecture of social responsiveness in families with parental schizophrenia or bipolar disorder and population-based controls

Twin-studies of social responsiveness have reported moderate to high heritabilities, but studies using parent-child data are lacking. Additionally, social impairments have been suggested as a vulnerability marker for schizophrenia and bipolar disorder, but the heritability of social responsiveness in this context is unknown. This study is part of the Danish High Risk and Resilience Study - VIA, comprising families with one parent with schizophrenia (n = 202) or bipolar disorder (n = 120) and population-based controls (PBC, n = 200). Social responsiveness was assessed with The Social Responsiveness Scale, Second Edition (SRS-2). Heritability was estimated from variance components, and a polygenic risk score (PRS) for autism spectrum disorder (ASD) was calculated to assess the genetic relationship between ASD and SRS-2. SRS-2 heritability was moderate to high and significantly different from zero in all groups when the children were rated by the primary caregiver. With teacher ratings, the heritability was lower and only significant in the full cohort and PBC. We found no significant association between SRS-2 and PRS for ASD. Our study confirms that social responsiveness is heritable, but that heritability estimates are affected by the child-respondent relation and familial risk of mental illness. This has implications for clinical practice and research using SRS-2 and provides insight on the familial transmission of mental illness.

Genetic Ancestry Estimates within Dutch Family Units and Across Genotyping Arrays: Insights from Empirical Analysis Using Two Estimation Methods

Accurate inference of genetic ancestry is crucial for population-based association studies, accounting for population heterogeneity and structure. This study analyzes genome-wide SNP data from the Netherlands Twin Register to compare genetic ancestry estimates. The focus is on the comparison of ancestry estimates between family members and individuals genotyped on multiple arrays (Affymetrix 6.0, Affymetrix Axiom, and Illumina GSA). Two conventional methods, principal component analysis and ADMIXTURE, were implemented to estimate ancestry, each serving its specific purpose, rather than for direct comparison. The results reveal that as the degree of genetic relatedness decreases, the Euclidean distances of genetic ancestry estimates between family members significantly increase (empirical p < 0.001), regardless of the estimation method and genotyping array. Ancestry estimates among individuals genotyped on multiple arrays also show statistically significant differences (empirical p < 0.001). Additionally, this study investigates the relationship between the ancestry estimates of non-identical twin offspring with ancestrally diverse parents and those with ancestrally similar parents. The results indicate a statistically significant weak correlation between the variation in ancestry estimates among offspring and differences in ancestry estimates among parents (Spearman's rho: 0.07, p = 0.005). This study highlights the utility of current methods in inferring genetic ancestry, emphasizing the importance of reference population composition in determining ancestry estimates.

Leveraging family data to design Mendelian randomization that is provably robust to population stratification

Mendelian randomization (MR) has emerged as a powerful approach to leverage genetic instruments to infer causality between pairs of traits in observational studies. However, the results of such studies are susceptible to biases owing to weak instruments, as well as the confounding effects of population stratification and horizontal pleiotropy. Here, we show that family data can be leveraged to design MR tests that are provably robust to confounding from population stratification, assortative mating, and dynastic effects. We show in simulations that our approach, MR-Twin, is robust to confounding from population stratification and is not affected by weak instrument bias, whereas standard MR methods yield inflated false positive rates. We then conduct an exploratory analysis of MR-Twin and other MR methods applied to 121 trait pairs in the UK Biobank data set. Our results suggest that confounding from population stratification can lead to false positives for existing MR methods, whereas MR-Twin is immune to this type of confounding, and that MR-Twin can help assess whether traditional approaches may be inflated owing to confounding from population stratification.

Genetics of destemming in pepper: A step towards mechanical harvesting

Introduction: The majority of peppers in the US for fresh market and processing are handpicked, and harvesting can account for 20-50% of production costs. Innovation in mechanical harvesting would increase availability; lower the costs of local, healthy vegetable products; and perhaps improve food safety and expand markets. Most processed peppers require removal of pedicels (stem and calyx) from the fruit, but lack of an efficient mechanical process for this operation has hindered adoption of mechanical harvest. In this paper, we present characterization and advancements in breeding green chile peppers for mechanical harvesting. Specifically, we describe inheritance and expression of an easy-destemming trait derived from the landrace UCD-14 that facilitates machine harvest of green chiles. Methods: A torque gauge was used for measuring bending forces similar to those of a harvester and applied to two biparental populations segregating for destemming force and rate. Genotyping by sequencing was used to generate genetic maps for quantitative trait locus (QTL) analyses. Results: A major destemming QTL was found on chromosome 10 across populations and environments. Eight additional population and/or environment-specific QTL were also identified. Chromosome 10 QTL markers were used to help introgress the destemming trait into jalapeño-type peppers. Low destemming force lines combined with improvements in transplant production enabled mechanical harvest of destemmed fruit at a rate of 41% versus 2% with a commercial jalapeńo hybrid. Staining for the presence of lignin at the pedicel/fruit boundary indicated the presence of an abscission zone and homologs of genes known to affect organ abscission were found under several QTL, suggesting that the easy-destemming trait may be due to the presence and activation of a pedicel/fruit abscission zone. Conclusion: Presented here are tools to measure the easy-destemming trait, its physiological basis, possible molecular pathways, and expression of the trait in various genetic backgrounds. Mechanical harvest of destemmed mature green chile fruits was achieved by combining easy-destemming with transplant management.

Interpreting population and family-based genome-wide association studies in the presence of confounding

A central aim of genome-wide association studies (GWASs) is to estimate direct genetic effects: the causal effects on an individual's phenotype of the alleles that they carry. However, estimates of direct effects can be subject to genetic and environmental confounding, and can also absorb the 'indirect' genetic effects of relatives' genotypes. Recently, an important development in controlling for these confounds has been the use of within-family GWASs, which, because of the randomness of Mendelian segregation within pedigrees, are often interpreted as producing unbiased estimates of direct effects. Here, we present a general theoretical analysis of the influence of confounding in standard population-based and within-family GWASs. We show that, contrary to common interpretation, family-based estimates of direct effects can be biased by genetic confounding. In humans, such biases will often be small per-locus, but can be compounded when effect size estimates are used in polygenic scores. We illustrate the influence of genetic confounding on population- and family-based estimates of direct effects using models of assortative mating, population stratification, and stabilizing selection on GWAS traits. We further show how family-based estimates of indirect genetic effects, based on comparisons of parentally transmitted and untransmitted alleles, can suffer substantial genetic confounding. In addition to known biases that can arise in family-based GWASs when interactions between family members are ignored, we show that biases can also arise from gene-by-environment (G×E) interactions when parental genotypes are not distributed identically across interacting environmental and genetic backgrounds. We conclude that, while family-based studies have placed GWAS estimation on a more rigorous footing, they carry subtle issues of interpretation that arise from confounding and interactions.

Leveraging family data to design Mendelian Randomization that is provably robust to population stratification

Mendelian Randomization (MR) has emerged as a powerful approach to leverage genetic instruments to infer causality between pairs of traits in observational studies. However, the results of such studies are susceptible to biases due to weak instruments as well as the confounding effects of population stratification and horizontal pleiotropy. Here, we show that family data can be leveraged to design MR tests that are provably robust to confounding from population stratification, assortative mating, and dynastic effects. We demonstrate in simulations that our approach, MR-Twin, is robust to confounding from population stratification and is not affected by weak instrument bias, while standard MR methods yield inflated false positive rates. We applied MR-Twin to 121 trait pairs in the UK Biobank dataset and found that MR-Twin identifies likely causal trait pairs and does not identify trait pairs that are unlikely to be causal. Our results suggest that confounding from population stratification can lead to false positives for existing MR methods, while MR-Twin is immune to this type of confounding.

A family-based study of genetic and epigenetic effects across multiple neurocognitive, motor, social-cognitive and social-behavioral functions

Many psychiatric and neurodevelopmental disorders are known to be heritable, but studies trying to elucidate the genetic architecture of such traits often lag behind studies of somatic traits and diseases. The reasons as to why relatively few genome-wide significant associations have been reported for such traits have to do with the sample sizes needed for the detection of small effects, the difficulty in defining and characterizing the phenotypes, partially due to overlaps in affected underlying domains (which is especially true for cognitive phenotypes), and the complex genetic architectures of the phenotypes, which are not wholly captured in traditional case-control GWAS designs. We aimed to tackle the last two issues by performing GWASs of eight quantitative neurocognitive, motor, social-cognitive and social-behavioral traits, which may be considered endophenotypes for a variety of psychiatric and neurodevelopmental conditions, and for which we employed models capturing both general genetic association and parent-of-origin effects, in a family-based sample comprising 402 children and their parents (mostly family trios). We identified 48 genome-wide significant associations across several traits, of which 3 also survived our strict study-wide quality criteria. We additionally performed a functional annotation of implicated genes, as most of the 48 associations were with variants within protein-coding genes. In total, our study highlighted associations with five genes (TGM3, CACNB4, ANKS1B, CSMD1 and SYNE1) associated with measures of working memory, processing speed and social behavior. Our results thus identify novel associations, including previously unreported parent-of-origin associations with relevant genes, and our top results illustrate new potential gene → endophenotype → disorder pathways.

Dissecting Complex Traits Using Omics Data: A Review on the Linear Mixed Models and Their Application in GWAS

Genome-wide association study (GWAS) is the most popular approach to dissecting complex traits in plants, humans, and animals. Numerous methods and tools have been proposed to discover the causal variants for GWAS data analysis. Among them, linear mixed models (LMMs) are widely used statistical methods for regulating confounding factors, including population structure, resulting in increased computational proficiency and statistical power in GWAS studies. Recently more attention has been paid to pleiotropy, multi-trait, gene-gene interaction, gene-environment interaction, and multi-locus methods with the growing availability of large-scale GWAS data and relevant phenotype samples. In this review, we have demonstrated all possible LMMs-based methods available in the literature for GWAS. We briefly discuss the different LMM methods, software packages, and available open-source applications in GWAS. Then, we include the advantages and weaknesses of the LMMs in GWAS. Finally, we discuss the future perspective and conclusion. The present review paper would be helpful to the researchers for selecting appropriate LMM models and methods quickly for GWAS data analysis and would benefit the scientific society.

Estimating effects of parents' cognitive and non-cognitive skills on offspring education using polygenic scores

Understanding how parents' cognitive and non-cognitive skills influence offspring education is essential for educational, family and economic policy. We use genetics (GWAS-by-subtraction) to assess a latent, broad non-cognitive skills dimension. To index parental effects controlling for genetic transmission, we estimate indirect parental genetic effects of polygenic scores on childhood and adulthood educational outcomes, using siblings (N = 47,459), adoptees (N = 6407), and parent-offspring trios (N = 2534) in three UK and Dutch cohorts. We find that parental cognitive and non-cognitive skills affect offspring education through their environment: on average across cohorts and designs, indirect genetic effects explain 36-40% of population polygenic score associations. However, indirect genetic effects are lower for achievement in the Dutch cohort, and for the adoption design. We identify potential causes of higher sibling- and trio-based estimates: prenatal indirect genetic effects, population stratification, and assortative mating. Our phenotype-agnostic, genetically sensitive approach has established overall environmental effects of parents' skills, facilitating future mechanistic work.

Genetic analyses of circulating PUFA-derived mediators identifies heritable dihydroxyeicosatrienoic acid species

Estimates of heritability are the first step in identifying a trait with substantial variation due to genetic factors. Large-scale genetic analyses can identify the DNA variants that influence the levels of circulating lipid species and the statistical technique Mendelian randomisation can use these DNA variants to address potential causality of these lipids in disease. We estimated the heritability of plasma eicosanoids, octadecanoids and docosanoids to identify those lipid species with substantial heritability. We analysed plasma lipid mediators in 31 White British families (196 participants) ascertained for high blood pressure and deeply clinically and biochemically phenotyped over a 25-year period. We found that the dihydroxyeicosatrienoic acid (DHET) species, 11,12-DHET and 14,15-DHET, products of arachidonic acid metabolism by cytochrome P450 (CYP) monooxygenase and soluble epoxide hydrolase (sEH), exhibited substantial heritability (h2 = 33%-37%; Padj<0.05). Identification of these two heritable bioactive lipid species allows for future large-scale, targeted, lipidomics-genomics analyses to address causality in cardiovascular and other diseases.

Within-sibship genome-wide association analyses decrease bias in estimates of direct genetic effects

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies

Genome-wide association studies (GWAS) are effectively applied to detect the marker trait associations (MTAs) using whole genome-wide variants for complex quantitative traits in different crop species. GWAS has been applied in wheat for different quality, biotic and abiotic stresses, and agronomic and yield-related traits. Predictions for marker-trait associations are controlled with the development of better statistical models taking population structure and familial relatedness into account. In this review, we have provided a detailed overview of the importance of association mapping, population design, high-throughput genotyping and phenotyping platforms, advancements in statistical models and multiple threshold comparisons, and recent GWA studies conducted in wheat. The information about MTAs utilized for gene characterization and adopted in breeding programs is also provided. In the literature that we surveyed, as many as 86,122 wheat lines have been studied under various GWA studies reporting 46,940 loci. However, further utilization of these is largely limited. The future breakthroughs in area of genomic selection, multi-omics-based approaches, machine, and deep learning models in wheat breeding after exploring the complex genetic structure with the GWAS are also discussed. This is a most comprehensive study of a large number of reports on wheat GWAS and gives a comparison and timeline of technological developments in this area. This will be useful to new researchers or groups who wish to invest in GWAS.

Assortative mating and within-spouse pair comparisons

Spousal comparisons have been proposed as a design that can both reduce confounding and estimate effects of the shared adulthood environment. However, assortative mating, the process by which individuals select phenotypically (dis)similar mates, could distort associations when comparing spouses. We evaluated the use of spousal comparisons, as in the within-spouse pair (WSP) model, for aetiological research such as genetic association studies. We demonstrated that the WSP model can reduce confounding but may be susceptible to collider bias arising from conditioning on assorted spouse pairs. Analyses using UK Biobank spouse pairs found that WSP genetic association estimates were smaller than estimates from random pairs for height, educational attainment, and BMI variants. Within-sibling pair estimates, robust to demographic and parental effects, were also smaller than random pair estimates for height and educational attainment, but not for BMI. WSP models, like other within-family models, may reduce confounding from demographic factors in genetic association estimates, and so could be useful for triangulating evidence across study designs to assess the robustness of findings. However, WSP estimates should be interpreted with caution due to potential collider bias.

Heritability and family-based GWAS analyses of the N-acyl ethanolamine and ceramide plasma lipidome

Signalling lipids of the N-acyl ethanolamine (NAE) and ceramide (CER) classes have emerged as potential biomarkers of cardiovascular disease (CVD). We sought to establish the heritability of plasma NAEs (including the endocannabinoid anandamide) and CERs, to identify common DNA variants influencing the circulating concentrations of the heritable lipids, and assess causality of these lipids in CVD using 2-sample Mendelian randomization (2SMR). Nine NAEs and 16 CERs were analyzed in plasma samples from 999 members of 196 British Caucasian families, using targeted ultra-performance liquid chromatography with tandem mass spectrometry. All lipids were significantly heritable (h2 = 36-62%). A missense variant (rs324420) in the gene encoding the enzyme fatty acid amide hydrolase (FAAH), which degrades NAEs, associated at genome-wide association study (GWAS) significance (P < 5 × 10-8) with four NAEs (DHEA, PEA, LEA and VEA). For CERs, rs680379 in the SPTLC3 gene, which encodes a subunit of the rate-limiting enzyme in CER biosynthesis, associated with a range of species (e.g. CER[N(24)S(19)]; P = 4.82 × 10-27). We observed three novel associations between SNPs at the CD83, SGPP1 and DEGS1 loci, and plasma CER traits (P < 5 × 10-8). 2SMR in the CARDIoGRAMplusC4D cohorts (60 801 cases; 123 504 controls) and in the DIAGRAM cohort (26 488 cases; 83 964 controls), using the genetic instruments from our family-based GWAS, did not reveal association between genetically determined differences in CER levels and CVD or diabetes. Two of the novel GWAS loci, SGPP1 and DEGS1, suggested a casual association between CERs and a range of haematological phenotypes, through 2SMR in the UK Biobank, INTERVAL and UKBiLEVE cohorts (n = 110 000-350 000).

Variants in regulatory elements of PDE4D associate with major mental illness in the Finnish population

We have previously reported a replicable association between variants at the PDE4D gene and familial schizophrenia in a Finnish cohort. In order to identify the potential functional mutations underlying these previous findings, we sequenced 1.5 Mb of the PDE4D genomic locus in 20 families (consisting of 96 individuals and 79 independent chromosomes), followed by two stages of genotyping across 6668 individuals from multiple Finnish cohorts for major mental illnesses. We identified 4570 SNPs across the PDE4D gene, with 380 associated to schizophrenia (p ≤ 0.05). Importantly, two of these variants, rs35278 and rs165940, are located at transcription factor-binding sites, and displayed replicable association in the two-stage enlargement of the familial schizophrenia cohort (combined statistics for rs35278 p = 0.0012; OR = 1.18, 95% CI: 1.06-1.32; and rs165940 p = 0.0016; OR = 1.27, 95% CI: 1.13-1.41). Further analysis using additional cohorts and endophenotypes revealed that rs165940 principally associates within the psychosis (p = 0.025, OR = 1.18, 95% CI: 1.07-1.30) and cognitive domains of major mental illnesses (g-score p = 0.044, β = -0.033). Specifically, the cognitive domains represented verbal learning and memory (p = 0.0091, β = -0.044) and verbal working memory (p = 0.0062, β = -0.036). Moreover, expression data from the GTEx database demonstrated that rs165940 significantly correlates with the mRNA expression levels of PDE4D in the cerebellum (p-value = 0.04; m-value = 0.9), demonstrating a potential functional consequence for this variant. Thus, rs165940 represents the most likely functional variant for major mental illness at the PDE4D locus in the Finnish population, increasing risk broadly to psychotic disorders.

Innovative approach to identify multigenomic and environmental interactions associated with birth defects in family-based hybrid designs

Genes, including those with transgenerational effects, work in concert with behavioral, environmental, and social factors via complex biological networks to determine human health. Understanding complex relationships between causal factors underlying human health is an essential step towards deciphering biological mechanisms. We propose a new analytical framework to investigate the interactions between maternal and offspring genetic variants or their surrogate single nucleotide polymorphisms (SNPs) and environmental factors using family-based hybrid study design. The proposed approach can analyze diverse genetic and environmental factors and accommodate samples from a variety of family units, including case/control-parental triads, and case/control-parental dyads, while minimizing potential bias introduced by population admixture. Comprehensive simulations demonstrated that our innovative approach outperformed the log-linear approach, the best available method for case-control family data. The proposed approach had greater statistical power and was capable to unbiasedly estimate the maternal and child genetic effects and the effects of environmental factors, while controlling the Type I error rate against population stratification. Using our newly developed approach, we analyzed the associations between maternal and fetal SNPs and obstructive and conotruncal heart defects, with adjustment for demographic and lifestyle factors and dietary supplements. Fourteen and 11 fetal SNPs were associated with obstructive and conotruncal heart defects, respectively. Twenty-seven and 17 maternal SNPs were associated with obstructive and conotruncal heart defects, respectively. In addition, maternal body mass index was a significant risk factor for obstructive defects. The proposed approach is a powerful tool for interrogating the etiological mechanism underlying complex traits.

Notes on Three Decades of Methodology Workshops

Since 1987, a group of behavior geneticists have been teaching an annual methodology workshop on how to use state-of-the-art methods to analyze genetically informative data. In the early years, the focus was on analyzing twin and family data, using information of their known genetic relatedness to infer the role of genetic and environmental factors on phenotypic variation. With the rapid evolution of genotyping and sequencing technology and availability of measured genetic data, new methods to detect genetic variants associated with human traits were developed and became the focus of workshop teaching in alternate years. Over the years, many of the methodological advances in the field of statistical genetics have been direct outgrowths of the workshop, as evidence by the software and methodological publications authored by workshop faculty. We provide data and demographics of workshop attendees and evaluate the impact of the methodology workshops on scientific output in the field by evaluating the number of papers applying specific statistical genetic methodologies authored by individuals who have attended workshops.

Disease Expression and Familial Transmission of Fuchs Endothelial Corneal Dystrophy With and Without CTG18.1 Expansion

Purpose

To characterize inheritance, penetrance, and trinucleotide repeat expansion stability in Fuchs endothelial corneal dystrophy (FECD).

Methods

One thousand unrelated and related subjects with and without FECD were prospectively recruited. CTG18.1 repeat length (CTG18.1L) was determined via short tandem repeat assay and Southern blotting of leukocyte DNA. Multivariable logistic regression and generalized estimating equation models were employed.

Results

There were 546 unrelated FECD cases (67.6% female; 70 ± 10 years) and 235 controls (63.8% female; 73 ± 8 years; all ≥ 50 years). CTG18.1 expansion (CTG18.1exp+) was observed in 424 (77.7%) cases and 18 (7.7%) controls (P = 2.48 × 10-44). CTG18.1 expansion was associated with FECD severity (P = 5.62 × 10-7). The family arm of the study included 331 members from 112 FECD-affected families; 87 families were CTG18.1exp+. Autosomal dominant inheritance with variable expression of FECD was observed, regardless of expansion status. FECD penetrance of CTG18.1 expansion increased with age, ranging from 44.4% in the youngest (19-46 years) to 86.2% in the oldest (64-91 years) age quartiles. Among 62 parent-offspring transmissions of CTG18.1exp+, 48 (77.4%) had a change in CTG18.1L ≤ 10 repeats, and eight (12.9%) were ≥50 repeats, including five large expansions (∼1000-2000 repeats) that contracted. Among 44 offspring who did not inherit the CTG18.1exp+ allele, eight (18.2%) exhibited FECD.

Conclusions

CTG18.1 expansion was highly associated with FECD but demonstrated incomplete penetrance. CTG18.1L instability occurred in a minority of parent-offspring transmissions, with large expansions exhibiting contraction. The observation of FECD without CTG18.1 expansion among family members in CTG18.1exp+ families highlights the complexity of the relationship between the FECD phenotype and CTG18.1 expansion.

LGR4 Gene Polymorphisms Are Associated With Bone and Obesity Phenotypes in Chinese Female Nuclear Families

OBJECTIVE

The current study was conducted to determine whether peak bone mineral density (BMD) and obesity phenotypes are associated with certain LGR4 gene polymorphisms found in Chinese nuclear families with female children.

METHODS

A total of 22 single nucleotide polymorphisms (SNPs) located in and around the LGR4 gene were identified in 1,300 subjects who were members of 390 Chinese nuclear families with female children. Then, BMD readings of the femoral neck, total hip, and lumbar spine as well as measurements of the total lean mass (TLM), total fat mass (TFM), and trunk fat mass were obtained via dual-energy X-ray absorptiometry. The quantitative transmission disequilibrium test was used to analyze the associations between specific SNPs and LGR4 haplotypes and peak BMD as well as between LGR4 haplotypes and TLM, percent lean mass, TFM, percent fat mass, trunk fat mass, and body mass index (BMI).

RESULTS

Here, rs7936621 was significantly associated with the BMD values for the total hip and lumbar spine, while rs10835171 and rs6484295 were associated with the trunk fat mass and BMI, respectively. Regarding the haplotypes, we found significant associations between GAA in block 2 and trunk fat mass and BMI, between AGCGT in block 3 and total hip BMD, between TGCTCC in block 5 and femoral neck BMD, and between TACTTC in block 5 and both lumbar spine and femoral neck BMD (all P-values < 0.05).

CONCLUSION

Genetic variations of the LGR4 gene are related to peak BMD, BMI, and trunk fat mass.

The rs1634330 Polymorphisms in the SOST Gene Are Associated with Body Composition in Chinese Nuclear Families with Male Offspring

OBJECTIVE

The purpose of this study was to explore the effect of the SOST gene polymorphisms on body composition in Chinese nuclear families with male offspring.

METHODS

1,016 individuals were recruited from 335 Chinese nuclear families with male offspring. The nuclear families consist of at least one male offspring aged 18 to 44. We genotyped the 10 tagged single-nucleotide polymorphisms (SNPs) in SOST gene (rs7220711, rs865429, rs851057, rs1708635, rs2023794, rs1234612, rs74252774, rs1634330, rs851058, and rs1513670) in all the above people. We used dual-energy X-ray absorptiometry to measure the composition of the human body. The quantitative transmission disequilibrium test (QTDT) was used to analyze the associations of the SNPs with the body composition.

RESULTS

QTDT analysis showed that rs1634330 was significantly associated with trunk LM (P < 0.05). However, haplotypes were not found to be significantly associated with the body composition in the within-family association. The 1000 permutations were consistent with these within-family association results.

CONCLUSIONS

Our results showed that the genetic variation in the SOST gene may contribute to variations in the body composition of Chinese male offspring.

The Boulder Workshop Question Box

The International Statistical Genetics Workshop (commonly referred to as the "Boulder Workshop") has been held annually in Boulder, Colorado almost every year since 1990. A staple feature of each workshop has been the presence of a "question box" (either a physical box or an online virtual one) where workshop participants are given the opportunity of asking questions to the faculty. In this manuscript, we have compiled a list of ten "classic" questions that have appeared in one form or another across multiple workshops and our attempts at answering them.

Statistically efficient association analysis of quantitative traits with haplotypes and untyped SNPs in family studies

BACKGROUND

Associations between haplotypes and quantitative traits provide valuable information about the genetic basis of complex human diseases. Haplotypes also provide an effective way to deal with untyped SNPs. Two major challenges arise in haplotype-based association analysis of family data. First, haplotypes may not be inferred with certainty from genotype data. Second, the trait values within a family tend to be correlated because of common genetic and environmental factors.

RESULTS

To address these challenges, we present an efficient likelihood-based approach to analyzing associations of quantitative traits with haplotypes or untyped SNPs. This approach properly accounts for within-family trait correlations and can handle general pedigrees with arbitrary patterns of missing genotypes. We characterize the genetic effects on the quantitative trait by a linear regression model with random effects and develop efficient likelihood-based inference procedures. Extensive simulation studies are conducted to examine the performance of the proposed methods. An application to family data from the Childhood Asthma Management Program Ancillary Genetic Study is provided. A computer program is freely available.

CONCLUSIONS

Results from extensive simulation studies show that the proposed methods for testing the haplotype effects on quantitative traits have correct type I error rates and are more powerful than some existing methods.

Heritability of haemodynamics in the ascending aorta

Blood flow in the vasculature can be characterised by dimensionless numbers commonly used to define the level of instabilities in the flow, for example the Reynolds number, Re. Haemodynamics play a key role in cardiovascular disease (CVD) progression. Genetic studies have identified mechanosensitive genes with causal roles in CVD. Given that CVD is highly heritable and abnormal blood flow may increase risk, we investigated the heritability of fluid metrics in the ascending aorta calculated using patient-specific data from cardiac magnetic resonance (CMR) imaging. 341 participants from 108 British Caucasian families were phenotyped by CMR and genotyped for 557,124 SNPs. Flow metrics were derived from the CMR images to provide some local information about blood flow in the ascending aorta, based on maximum values at systole at a single location, denoted max, and a 'peak mean' value averaged over the area of the cross section, denoted pm. Heritability was estimated using pedigree-based (QTDT) and SNP-based (GCTA-GREML) methods. Estimates of Reynolds number based on spatially averaged local flow during systole showed substantial heritability ([Formula: see text], [Formula: see text]), while the estimated heritability for Reynolds number calculated using the absolute local maximum velocity was not statistically significant (12-13%; [Formula: see text]). Heritability estimates of the geometric quantities alone; e.g. aortic diameter ([Formula: see text], [Formula: see text]), were also substantially heritable, as described previously. These findings indicate the potential for the discovery of genetic factors influencing haemodynamic traits in large-scale genotyped and phenotyped cohorts where local spatial averaging is used, rather than instantaneous values. Future Mendelian randomisation studies of aortic haemodynamic estimates, which are swift to derive in a clinical setting, will allow for the investigation of causality of abnormal blood flow in CVD.

Copy number variants in siblings of Mexican origin concordant for schizophrenia or bipolar disorder

Schizophrenia (SCZ) and bipolar disorder (BD) cause similar symptomatology. A correlation between these disorders has been found. We aimed to explore shared CNVs between SCZ and BD, in 35 sibpairs diagnosed with SCZ and 21 sibpairs diagnosed with BD. CNV calling was performed using data derived of Psycharray, by PennCNV. We did not find any shared CNVs between individuals diagnosed with BD and SCZ, neither with psychotic symptoms in individuals with BD. Nevertheless, we found a significant higher CNV burden in early-onset SCZ. This is one of the first's studies analyzing shared CNVs between SCZ and BD in Mexican population.

Detecting X-linked common and rare variant effects in family-based sequencing studies

The breakthroughs in next generation sequencing have allowed us to access data consisting of both common and rare variants, and in particular to investigate the impact of rare genetic variation on complex diseases. Although rare genetic variants are thought to be important components in explaining genetic mechanisms of many diseases, discovering these variants remains challenging, and most studies are restricted to population-based designs. Further, despite the shift in the field of genome-wide association studies (GWAS) towards studying rare variants due to the "missing heritability" phenomenon, little is known about rare X-linked variants associated with complex diseases. For instance, there is evidence that X-linked genes are highly involved in brain development and cognition when compared with autosomal genes; however, like most GWAS for other complex traits, previous GWAS for mental diseases have provided poor resources to deal with identification of rare variant associations on X-chromosome. In this paper, we address the two issues described above by proposing a method that can be used to test X-linked variants using sequencing data on families. Our method is much more general than existing methods, as it can be applied to detect both common and rare variants, and is applicable to autosomes as well. Our simulation study shows that the method is efficient, and exhibits good operational characteristics. An application to the University of Miami Study on Genetics of Autism and Related Disorders also yielded encouraging results.

Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses

Estimates from Mendelian randomization studies of unrelated individuals can be biased due to uncontrolled confounding from familial effects. Here we describe methods for within-family Mendelian randomization analyses and use simulation studies to show that family-based analyses can reduce such biases. We illustrate empirically how familial effects can affect estimates using data from 61,008 siblings from the Nord-Trøndelag Health Study and UK Biobank and replicated our findings using 222,368 siblings from 23andMe. Both Mendelian randomization estimates using unrelated individuals and within family methods reproduced established effects of lower BMI reducing risk of diabetes and high blood pressure. However, while Mendelian randomization estimates from samples of unrelated individuals suggested that taller height and lower BMI increase educational attainment, these effects were strongly attenuated in within-family Mendelian randomization analyses. Our findings indicate the necessity of controlling for population structure and familial effects in Mendelian randomization studies.

Quantitative genome-wide association analyses of receptive language in the Danish High Risk and Resilience Study

BACKGROUND

One of the most basic human traits is language. Linguistic ability, and disability, have been shown to have a strong genetic component in family and twin studies, but molecular genetic studies of language phenotypes are scarce, relative to studies of other cognitive traits and neurodevelopmental phenotypes. Moreover, most genetic studies examining such phenotypes do not incorporate parent-of-origin effects, which could account for some of the heritability of the investigated trait. We performed a genome-wide association study of receptive language, examining both child genetic effects and parent-of-origin effects.

RESULTS

Using a family-based cohort with 400 children with receptive language scores, we found a genome-wide significant paternal parent-of-origin effect with a SNP, rs11787922, on chromosome 9q21.31, whereby the T allele reduced the mean receptive language score by ~ 23, constituting a reduction of more than 1.5 times the population SD (P = 1.04 × 10-8). We further confirmed that this association was not driven by broader neurodevelopmental diagnoses in the child or a family history of psychiatric diagnoses by incorporating covariates for the above and repeating the analysis.

CONCLUSIONS

Our study reports a genome-wide significant association for receptive language skills; to our knowledge, this is the first documented genome-wide significant association for this phenotype. Furthermore, our study illustrates the importance of considering parent-of-origin effects in association studies, particularly in the case of cognitive or neurodevelopmental traits, in which parental genetic data are not always incorporated.

Design of experiments for fine-mapping quantitative trait loci in livestock populations

BACKGROUND

Single nucleotide polymorphisms (SNPs) which capture a significant impact on a trait can be identified with genome-wide association studies. High linkage disequilibrium (LD) among SNPs makes it difficult to identify causative variants correctly. Thus, often target regions instead of single SNPs are reported. Sample size has not only a crucial impact on the precision of parameter estimates, it also ensures that a desired level of statistical power can be reached. We study the design of experiments for fine-mapping of signals of a quantitative trait locus in such a target region.

METHODS

A multi-locus model allows to identify causative variants simultaneously, to state their positions more precisely and to account for existing dependencies. Based on the commonly applied SNP-BLUP approach, we determine the z-score statistic for locally testing non-zero SNP effects and investigate its distribution under the alternative hypothesis. This quantity employs the theoretical instead of observed dependence between SNPs; it can be set up as a function of paternal and maternal LD for any given population structure.

RESULTS

We simulated multiple paternal half-sib families and considered a target region of 1 Mbp. A bimodal distribution of estimated sample size was observed, particularly if more than two causative variants were assumed. The median of estimates constituted the final proposal of optimal sample size; it was consistently less than sample size estimated from single-SNP investigation which was used as a baseline approach. The second mode pointed to inflated sample sizes and could be explained by blocks of varying linkage phases leading to negative correlations between SNPs. Optimal sample size increased almost linearly with number of signals to be identified but depended much stronger on the assumption on heritability. For instance, three times as many samples were required if heritability was 0.1 compared to 0.3. An R package is provided that comprises all required tools.

CONCLUSIONS

Our approach incorporates information about the population structure into the design of experiments. Compared to a conventional method, this leads to a reduced estimate of sample size enabling the resource-saving design of future experiments for fine-mapping of candidate variants.

Within family Mendelian randomization studies

Mendelian randomization (MR) is increasingly used to make causal inferences in a wide range of fields, from drug development to etiologic studies. Causal inference in MR is possible because of the process of genetic inheritance from parents to offspring. Specifically, at gamete formation and conception, meiosis ensures random allocation to the offspring of one allele from each parent at each locus, and these are unrelated to most of the other inherited genetic variants. To date, most MR studies have used data from unrelated individuals. These studies assume that genotypes are independent of the environment across a sample of unrelated individuals, conditional on covariates. Here we describe potential sources of bias, such as transmission ratio distortion, selection bias, population stratification, dynastic effects and assortative mating that can induce spurious or biased SNP-phenotype associations. We explain how studies of related individuals such as sibling pairs or parent-offspring trios can be used to overcome some of these sources of bias, to provide potentially more reliable evidence regarding causal processes. The increasing availability of data from related individuals in large cohort studies presents an opportunity to both overcome some of these biases and also to evaluate familial environmental effects.

WGS-based telomere length analysis in Dutch family trios implicates stronger maternal inheritance and a role for RRM1 gene

Telomere length (TL) regulation is an important factor in ageing, reproduction and cancer development. Genetic, hereditary and environmental factors regulating TL are currently widely investigated, however, their relative contribution to TL variability is still understudied. We have used whole genome sequencing data of 250 family trios from the Genome of the Netherlands project to perform computational measurement of TL and a series of regression and genome-wide association analyses to reveal TL inheritance patterns and associated genetic factors. Our results confirm that TL is a largely heritable trait, primarily with mother's, and, to a lesser extent, with father's TL having the strongest influence on the offspring. In this cohort, mother's, but not father's age at conception was positively linked to offspring TL. Age-related TL attrition of 40 bp/year had relatively small influence on TL variability. Finally, we have identified TL-associated variations in ribonuclease reductase catalytic subunit M1 (RRM1 gene), which is known to regulate telomere maintenance in yeast. We also highlight the importance of multivariate approach and the limitations of existing tools for the analysis of TL as a polygenic heritable quantitative trait.

Loss of ZnT8 function protects against diabetes by enhanced insulin secretion

A rare loss-of-function allele p.Arg138* in SLC30A8 encoding the zinc transporter 8 (ZnT8), which is enriched in Western Finland, protects against type 2 diabetes (T2D). We recruited relatives of the identified carriers and showed that protection was associated with better insulin secretion due to enhanced glucose responsiveness and proinsulin conversion, particularly when compared with individuals matched for the genotype of a common T2D-risk allele in SLC30A8, p.Arg325. In genome-edited human induced pluripotent stem cell (iPSC)-derived β-like cells, we establish that the p.Arg138* allele results in reduced SLC30A8 expression due to haploinsufficiency. In human β cells, loss of SLC30A8 leads to increased glucose responsiveness and reduced KATP channel function similar to isolated islets from carriers of the T2D-protective allele p.Trp325. These data position ZnT8 as an appealing target for treatment aimed at maintaining insulin secretion capacity in T2D.

Marked variation in heritability estimates of left ventricular mass depending on modality of measurement

Left ventricular (LV) hypertrophy is a strong risk factor for heart failure and cardiovascular death. ECG measures of LV mass are estimated as heritable in twin and family-based analyses and heritability estimates of LV mass measured by echocardiography are lower. We hypothesised that CMR-derived measurements, being more precise than echocardiographic measurements, would advance our understanding of heritable LV traits. We phenotyped 116 British families (427 individuals) by CMR and ECG, and undertook heritability analyses using variance-components (QTDT) and GWAS SNP-based (GCTA-GREML) methods. ECG-based traits such as LV mass and Sokolow-Lyon duration showed substantial estimates of heritability (60%), whereas CMR-derived LV mass was only modestly heritable (20%). However, the ECG LV mass was positively correlated with the lateral diameter of the chest (rho = 0.67), and adjustment for this attenuated the heritability estimate (42%). Finally, CMR-derived right ventricular mass showed considerable heritability (44%). Heritability estimates of LV phenotypes show substantial variation depending on the modality of measurement, being greater when measured by ECG than CMR. This may reflect the differences between electrophysiological as opposed to anatomical hypertrophy. However, ECG LV hypertrophy traits are likely to be influenced by genetic association with anthropometric measures, inflating their overall measured heritability.

GWAS Analysis Indicated Importance of NF-κB Signaling Pathway in Host Resistance Against Motile Aeromonas Septicemia Disease in Catfish

Motile Aeromonas septicemia (MAS) disease caused by a bacterial pathogen, Aeromonas hydrophila, is an emerging but severe disease of catfish. Genetic enhancement of disease resistance is considered to be effective to control the disease. To provide an insight into the genomic basis of MAS disease resistance, in this study, we conducted a genome-wide association study (GWAS) to identify quantitative trait loci (QTL). A total of 1820 interspecific backcross catfish of 7 families were challenged with A. hydrophila, and 382 phenotypic extremes were selected for genotyping with the catfish 690 K SNP arrays. Three QTL on linkage group (LG) 2, 26 and 29 were identified to be significantly associated with MAS resistance. Within these regions, a total of 24 genes had known functions in immunity, 10 of which were involved in NF-κB signaling pathway, suggesting the importance of NF-κB signaling pathway in MAS resistance. In addition, three suggestively significant QTL were identified on LG 11, 17, and 20. The limited numbers of QTL involved in MAS resistance suggests that marker-assisted selection may be a viable approach for catfish breeding.

High-Density Linkage Map and QTLs for Growth in Snapper (Chrysophrys auratus)

Characterizing the genetic variation underlying phenotypic traits is a central objective in biological research. This research has been hampered in the past by the limited genomic resources available for most non-model species. However, recent advances in sequencing technologies and related genotyping methods are rapidly changing this. Here we report the use of genome-wide SNP data from the ecologically and commercially important marine fish species Chrysophrys auratus (snapper) to 1) construct the first linkage map for this species, 2) scan for growth QTL, and 3) search for putative candidate genes in the surrounding QTL regions. The newly constructed linkage map contained ∼11K SNP markers and is one of the densest maps to date in the fish family Sparidae. Comparisons with genome scaffolds of the recently assembled snapper genome indicated that marker placement was mostly consistent between the scaffolds and linkage map (R = 0.7), but that at fine scales (< 5 cM) some precision limitations occurred. Of the 24 linkage groups, which likely reflect the 24 chromosomes of this species, three were found to contain QTL with genome-wide significance for growth-related traits. A scan of 13 candidate growth genes located the growth hormone, myogenin, and parvalbumin genes within 5.3, 9.6, and 25.0 cM of these QTL, respectively. The linkage map and QTL found in this study will advance the investigation of genome structure and aquaculture breeding efforts in this and related species.

Heritability estimation and differential analysis of count data with generalized linear mixed models in genomic sequencing studies

Motivation

Genomic sequencing studies, including RNA sequencing and bisulfite sequencing studies, are becoming increasingly common and increasingly large. Large genomic sequencing studies open doors for accurate molecular trait heritability estimation and powerful differential analysis. Heritability estimation and differential analysis in sequencing studies requires the development of statistical methods that can properly account for the count nature of the sequencing data and that are computationally efficient for large datasets.

Results

Here, we develop such a method, PQLseq (Penalized Quasi-Likelihood for sequencing count data), to enable effective and efficient heritability estimation and differential analysis using the generalized linear mixed model framework. With extensive simulations and comparisons to previous methods, we show that PQLseq is the only method currently available that can produce unbiased heritability estimates for sequencing count data. In addition, we show that PQLseq is well suited for differential analysis in large sequencing studies, providing calibrated type I error control and more power compared to the standard linear mixed model methods. Finally, we apply PQLseq to perform gene expression heritability estimation and differential expression analysis in a large RNA sequencing study in the Hutterites.

Availability and implementation

PQLseq is implemented as an R package with source code freely available at www.xzlab.org/software.html and https://cran.r-project.org/web/packages/PQLseq/index.html.

Supplementary information

Supplementary data are available at Bioinformatics online.

Family-based association tests for rare variants with censored traits

We propose a set of family-based burden and kernel tests for censored traits (FamBAC and FamKAC). Here, censored traits refer to time-to-event outcomes, for instance, age-at-onset of a disease. To model censored traits in family-based designs, we used the frailty model, which incorporated not only fixed genetic effects of rare variants in a region of interest but also random polygenic effects shared within families. We first partitioned genotype scores of rare variants into orthogonal between- and within-family components, and then derived their corresponding efficient score statistics from the frailty model. Finally, FamBAC and FamKAC were constructed by aggregating the weighted efficient scores of the within-family components across rare variants and subjects. FamBAC collapsed rare variants within subject first to form a burden test that followed a chi-squared distribution; whereas FamKAC was a variant component test following a mixture of chi-squared distributions. For FamKAC, p-values can be computed by permutation tests or for computational efficiency by approximation methods. Through simulation studies, we showed that type I error was correctly controlled by FamBAC for various variant weighting schemes (0.0371 to 0.0527). However, FamKAC type I error rates based on approximation methods were deflated (max 0.0376) but improved by permutation tests. Our simulations also demonstrated that burden test FamBAC had higher power than kernel test FamKAC when high proportion (e.g. ≥ 80%) of causal variants had effects in the same direction. In contrast, when the effects of causal variants on the censored trait were in mixed directions, FamKAC outperformed FamBAC and had comparable or higher power than an existing method, RVFam. Our proposed framework has the flexibility to accommodate general nuclear families, and can be used to analyze sequence data for censored traits such as age-at-onset of a complex disease of interest.

Reduced signal for polygenic adaptation of height in UK Biobank

Several recent papers have reported strong signals of selection on European polygenic height scores. These analyses used height effect estimates from the GIANT consortium and replication studies. Here, we describe a new analysis based on the the UK Biobank (UKB), a large, independent dataset. We find that the signals of selection using UKB effect estimates are strongly attenuated or absent. We also provide evidence that previous analyses were confounded by population stratification. Therefore, the conclusion of strong polygenic adaptation now lacks support. Moreover, these discrepancies highlight (1) that methods for correcting for population stratification in GWAS may not always be sufficient for polygenic trait analyses, and (2) that claims of differences in polygenic scores between populations should be treated with caution until these issues are better understood.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Metabolic GWAS-based dissection of genetic bases underlying the diversity of plant metabolism

Plants have served as sources providing humans with metabolites for food and nutrition, biomaterials for living, and treatment for pain and disease. Plants produce a huge array of metabolites, with an immense diversity at both the population and individual levels. Dissection of the genetic bases for metabolic diversity has attracted increasing research attention. The concept of genome-wide association study (GWAS) was extended to studies on the diversity of plant metabolome that benefitted from the development of mass-spectrometry-based analytical systems and genome sequencing technologies. Metabolic genome-wide association study (mGWAS) is one of the most powerful tools for global identification of genetic determinants for diversity of plant metabolism. Recently, mGWAS has been performed for various species with continuous improvements, providing deeper insights into the genetic bases of metabolic diversity. In this review, we discuss fully the achievements to date and remaining challenges that are associated with both mGWAS and mGWAS-based multi-dimensional analysis. We begin with a summary of GWAS and its development based on statistical methods and populations. As variation in targeted traits is essential for GWAS, we review metabolic diversity and its rise at both the population and individual levels. Subsequently, the application of mGWAS for plants and its corresponding achievements are fully discussed. We address the current knowledge on mGWAS-based multi-dimensional analysis and emerging insights into the diversity of metabolism.

The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels

Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Familial resemblances in human plasma metabolites are attributable to both genetic and common environmental effects

Metabolites are of great importance for understanding the pathogenesis of several diseases. Understanding the genetic contribution to metabolite concentrations may provide insights into mechanisms of complex diseases. Several studies have investigated heritability of metabolites but none investigated potential influences of genetic and environmental factors on the relationship between metabolites and cardiometabolic (CM) risk factors. Thus, we tested the hypothesis that both genetic and common environmental effects contribute to the variance of plasma metabolite concentrations and that shared genetic and environmental effects explain their phenotypic correlations with CM risk factors. To test this hypothesis, variance component method and bivariate genetic analysis were performed in a family-based sample of 48 French Canadians from 16 families. Familial resemblances were computed for all 147 detected metabolites and 9 (acetylornithine, acylcarnitine C9, arginine, phosphatidylcholine acyl-alkyl C36:4, serotonin, lysophosphatidylcholine acyl C20:4, citrulline, asymmetric dimethylarginine, phosphatidylcholine acyl-alkyl C36:5) showed a significant familial effect (55.7%, 18.7%, and 37.0% for maximal heritability, genetic heritability, and common environmental effect, respectively). Citrulline, phosphatidylcholine acyl-alkyl C36:4, phosphatidylcholine acyl-alkyl C36:5, and serotonin had significant phenotypic correlations with CM risk factors. Citrulline had a positive genetic correlation with apolipoprotein B100, while phosphatidylcholine acyl-alkyl C36:5 had a positive environmental correlation with total cholesterol. In conclusion, familial resemblances in metabolite concentrations were mainly attributable to common environmental effect when considering metabolites with a significant familial effect. Common genetic and environmental factors may also influence the relationship between metabolites and CM risk factors.

Investigation of parent-of-origin effects induced by fenofibrate treatment on triglycerides levels

Background

Genome-wide association studies performed on triglycerides (TGs) have not accounted for epigenetic mechanisms that may partially explain trait heritability.

Results

Parent-of-origin (POO) effect association analyses using an agnostic approach or a candidate approach were performed for pretreatment TG levels, posttreatment TG levels, and pre- and posttreatment TG-level differences in the real GAW20 family data set. We detected 22 genetic variants with suggestive POO effects with at least 1 phenotype (P ≤ 10^− 5). We evaluated the association of these 22 significant genetic variants showing POO effects with close DNA methylation probes associated with TGs. A total of 18 DNA methylation probes located in the vicinity of the 22 SNPs were associated with at least 1 phenotype and 6 SNP-probe pairs were associated with DNA methylation probes at the nominal level of P < 0.05, among which 1 pair presented evidence of POO effect. Our analyses identified a paternal effect of SNP rs301621 on the difference between pre- and posttreatment TG levels (P = 1.2 × 10^− 5). This same SNP showed evidence for a maternal effect on methylation levels of a nearby probe (cg10206250; P = 0.01). Using a causal inference test we established that the observed POO effect of rs301621 was not mediated by DNA methylation at cg10206250.

Conclusions

We performed POO effect association analyses of SNPs with TGs, as well as association analyses of SNPs with DNA methylation probes. These analyses, which were followed by a causal inference test, established that the paternal effect at the SNP rs301621 is induced by treatment and is not mediated by methylation level at cg10206250.

Application of novel and existing methods to identify genes with evidence of epigenetic association: results from GAW20

BACKGROUND

The rise in popularity and accessibility of DNA methylation data to evaluate epigenetic associations with disease has led to numerous methodological questions. As part of GAW20, our working group of 8 research groups focused on gene searching methods.

RESULTS

Although the methods were varied, we identified 3 main themes within our group. First, many groups tackled the question of how best to use pedigree information in downstream analyses, finding that (a) the use of kinship matrices is common practice, (b) ascertainment corrections may be necessary, and (c) pedigree information may be useful for identifying parent-of-origin effects. Second, many groups also considered multimarker versus single-marker tests. Multimarker tests had modestly improved power versus single-marker methods on simulated data, and on real data identified additional associations that were not identified with single-marker methods, including identification of a gene with a strong biological interpretation. Finally, some of the groups explored methods to combine single-nucleotide polymorphism (SNP) and DNA methylation into a single association analysis.

CONCLUSIONS

A causal inference method showed promise at discovering new mechanisms of SNP activity; gene-based methods of summarizing SNP and DNA methylation data also showed promise. Even though numerous questions still remain in the analysis of DNA methylation data, our discussions at GAW20 suggest some emerging best practices.