Exploring polygenic-environment and residual-environment interactions for depressive symptoms within the UK Biobank.

Substantial advances have been made in identifying genetic contributions to depression, but little is known about how the effect of genes can be modulated by the environment, creating a gene-environment interaction. Using multivariate reaction norm models (MRNMs) within the UK Biobank (N = 61294-91644), we investigate whether the polygenic and residual variance components of depressive symptoms are modulated by 17 a priori selected covariate traits-12 environmental variables and 5 biomarkers. MRNMs, a mixed-effects modelling approach, provide unbiased polygenic-covariate interaction estimates for a quantitative trait by controlling for outcome-covariate correlations and residual-covariate interactions. A continuous depressive symptom variable was the outcome in 17 MRNMs-one for each covariate trait. Each MRNM had a fixed-effects model (fixed effects included the covariate trait, demographic variables, and principal components) and a random effects model (where polygenic-covariate and residual-covariate interactions are modelled). Of the 17 selected covariates, 11 significantly modulate deviations in depressive symptoms through the modelled interactions, but no single interaction explains a large proportion of phenotypic variation. Results are dominated by residual-covariate interactions, suggesting that covariate traits (including neuroticism, childhood trauma, and BMI) typically interact with unmodelled variables, rather than a genome-wide polygenic component, to influence depressive symptoms. Only average sleep duration has a polygenic-covariate interaction explaining a demonstrably nonzero proportion of the variability in depressive symptoms. This effect is small, accounting for only 1.22% (95% confidence interval: [0.54, 1.89]) of variation. The presence of an interaction highlights a specific focus for intervention, but the negative results here indicate a limited contribution from polygenic-environment interactions.

Comparison of Helical Blade Systems for Osteoporotic Intertrochanteric Fractures Using Biomechanical Analysis and Clinical Assessments.

Background and Objectives: This study aimed to compare the biomechanical properties and outcomes of osteoporotic intertrochanteric fractures treated with two different helical blade systems, the trochanteric fixation nail-advanced (TFNA) and proximal femoral nail antirotation II (PFNA), to evaluate the efficacy and safety of the newly introduced TFNA system. Materials and Methods: A biomechanical comparison of the two helical blades was performed using uniaxial compression tests on polyurethane foam blocks of different densities. The peak resistance (PR) and accumulated resistance (AR) were measured during the 20 mm advancement through the test block. For clinical comparison, 63 osteoporotic intertrochanteric fractures treated with TFNA were identified and compared with the same number of fractures treated with PFNA using propensity score matching. Ambulatory status, medial migration, lateral sliding, fixation failure, and patient-reported outcomes were compared between the two groups over a minimum of 1 year's follow up. Results: The uniaxial compression test showed that a slightly, but significantly lower resistance was required to advance the TFNA through the test block compared with the PFNA (20 PCF, p = 0.017 and p = 0.026; 30 PCF, p = 0.007 and p = 0.001 for PR and AR, respectively). Clinically, the two groups showed no significant differences in post-operative ambulatory status and patient-reported outcomes. However, in TFNA groups, significantly more medial migration (TFNA, 0.75 mm; PFNA, 0.40 mm; p = 0.0028) and also, lateral sliding was noted (TFNA, 3.99 mm; PFNA, 1.80 mm; p = 0.004). Surgical failure occurred in four and two fractures treated with the TFNA and PFNA, respectively. Conclusions: The results of our study suggest that the newly introduced TFNA provides clinical outcomes comparable with those of the PFNA. However, inferior resistance to medial migration in the TFNA raises concerns regarding potential fixation failures.

Time taken to resume driving following hip arthroscopy.

BACKGROUND: Resuming driving is a common concern among patients undergoing hip arthroscopy. The present study aimed to assess whether patients who had undergone right hip arthroscopy presented with poorer driving performance than patients with normal hips and to analyze the time required to regain preoperative driving performance. METHODS: Forty-seven patients who had undergone right hip arthroscopy and consented to our test protocol were included in this study. Using an immersive driving simulator, the patients were tested for their brake reaction time (BRT), total brake time (TBT), and brake pedal depression (BPD) preoperatively and postoperatively. The first postoperative assessments were conducted when the patients could comfortably sit on the driving seat, and the follow-up assessments were conducted for 6 consecutive weeks at weekly intervals. The patients were divided into the following two groups based on the type of surgery that they underwent: the femoroacetabular impingement (FAI) surgery group and the simple hip arthroscopy (SA) group. Twenty healthy volunteers underwent driving assessments thrice at weekly intervals and constituted the control group. The braking parameters were compared between preoperative and postoperative measurements and among the FAI surgery, SA, and control groups. RESULTS: The preoperative braking parameters of the patients who underwent arthroscopy did not differ significantly from those of the controls (p = 0.373, 0.763, and 0.447 for the BRT, TBT, and BPD, respectively). All braking parameters returned to normal in 2 weeks in the FAI surgery group and in 1 week in the SA group. CONCLUSIONS: Our study suggests that the driving performance of patients who underwent right hip arthroscopy is comparable to that of individuals with normal hips and that the braking parameters may normalize to the preoperative state at 1 week after SA and 2 weeks after FAI surgery.

Detection of genomic regions underlying resistance to gastrointestinal parasites in Australian sheep.

BACKGROUND: This study aimed at identifying genomic regions that underlie genetic variation of worm egg count, as an indicator trait for parasite resistance in a large population of Australian sheep, which was genotyped with the high-density 600 K Ovine single nucleotide polymorphism array. This study included 7539 sheep from different locations across Australia that underwent a field challenge with mixed gastrointestinal parasite species. Faecal samples were collected and worm egg counts for three strongyle species, i.e. Teladorsagia circumcincta, Haemonchus contortus and Trichostrongylus colubriformis were determined. Data were analysed using genome-wide association studies (GWAS) and regional heritability mapping (RHM). RESULTS: Both RHM and GWAS detected a region on Ovis aries (OAR) chromosome 2 that was highly significantly associated with parasite resistance at a genome-wise false discovery rate of 5%. RHM revealed additional significant regions on OAR6, 18, and 24. Pathway analysis revealed 13 genes within these significant regions (SH3RF1, HERC2, MAP3K, CYFIP1, PTPN1, BIN1, HERC3, HERC5, HERC6, IBSP, SPP1, ISG20, and DET1), which have various roles in innate and acquired immune response mechanisms, as well as cytokine signalling. Other genes involved in haemostasis regulation and mucosal defence were also detected, which are important for protection of sheep against invading parasites. CONCLUSIONS: This study identified significant genomic regions on OAR2, 6, 18, and 24 that are associated with parasite resistance in sheep. RHM was more powerful in detecting regions that affect parasite resistance than GWAS. Our results support the hypothesis that parasite resistance is a complex trait and is determined by a large number of genes with small effects, rather than by a few major genes with large effects.

Using imputed whole-genome sequence data to improve the accuracy of genomic prediction for parasite resistance in Australian sheep.

BACKGROUND: This study aimed at (1) comparing the accuracies of genomic prediction for parasite resistance in sheep based on whole-genome sequence (WGS) data to those based on 50k and high-density (HD) single nucleotide polymorphism (SNP) panels; (2) investigating whether the use of variants within quantitative trait loci (QTL) regions that were selected from regional heritability mapping (RHM) in an independent dataset improved the accuracy more than variants selected from genome-wide association studies (GWAS); and (3) comparing the prediction accuracies between variants selected from WGS data to variants selected from the HD SNP panel. RESULTS: The accuracy of genomic prediction improved marginally from 0.16 ± 0.02 and 0.18 ± 0.01 when using all the variants from 50k and HD genotypes, respectively, to 0.19 ± 0.01 when using all the variants from WGS data. Fitting a GRM from the selected variants alongside a GRM from the 50k SNP genotypes improved the prediction accuracy substantially compared to fitting the 50k SNP genotypes alone. The gain in prediction accuracy was slightly more pronounced when variants were selected from WGS data compared to when variants were selected from the HD panel. When sequence variants that passed the GWAS [Formula: see text] threshold of 3 across the entire genome were selected, the prediction accuracy improved by 5% (up to 0.21 ± 0.01), whereas when selection was limited to sequence variants that passed the same GWAS [Formula: see text] threshold of 3 in regions identified by RHM, the accuracy improved by 9% (up to 0.25 ± 0.01). CONCLUSIONS: Our results show that through careful selection of sequence variants from the QTL regions, the accuracy of genomic prediction for parasite resistance in sheep can be improved. These findings have important implications for genomic prediction in sheep.

Accuracy of planar anteversion measurements using anteroposterior radiographs.

BACKGROUND: Several methods using simple anteroposterior (AP) radiographs have been suggested for the measurement of anteversion of the cup component after total hip arthroplasty. Herein, we compared six widely used anteversion measurement methods using two different types of AP radiograph, the conventional pelvis AP and hip-centered AP radiographs, to identify the measurement method and the type of radiograph that would provide the highest accuracy and reliability. METHODS: We developed two custom-made bi-planar anteversion measurement models for the validation test. The models were designed for pelvis AP and hip-centered AP radiographs, respectively. The radiographs were acquired using the inclination angles of both models, changing from 10° to 70° at 10° increments. For each inclination angle, anteversion was changed from 0° to 30° at 5° increments. The measurements were obtained independently by two orthopedic surgeons blinded from each other's measurements, using the methods of 1) Pradhan et al., 2) Lewinnek et al., 3) Widmer et al., 4) Liaw et al., 5) Hassan et al., and 6) Ackland et al. The measurements were repeated after 2 months. The accuracy, compared with that of the reference angle, and intra-observer and inter-observer reliabilities of each method were calculated. RESULTS: The highest accuracy was found when the method of Liaw et al. was used with hip-centered AP radiographs, which showed a difference of 1.37° ± 1.73 from the reference angle. Moreover, regardless of the type of radiograph, the methods by Pradhan et al., Lewinnek et al., and Liaw et al. showed excellent correlations with the reference anteversion. However, substantial differences were found when the methods by Widmer et al., Hassan et al., and Ackland et al. were used, regardless of the type of radiograph used. When anteversion was measured in an inclination between 30° and 50°, the method of Pradhan et al., when used with pelvis AP radiographs, showed the highest accuracy (1.23° ± 0.92°). We also found no significant difference in anteversions between the measurements made on pelvic and hip-centered AP radiographs. Both interobserver and intraobserver reliabilities were high for all the measurements tested. CONCLUSIONS: The methods by Pradhan et al., Liaw et al., and Lewinnek et al. may provide relatively accurate anteversion measurements with high reliability, regardless of the type of radiograph.

Effect of selection and selective genotyping for creation of reference on bias and accuracy of genomic prediction.

Reference populations for genomic selection usually involve selected individuals, which may result in biased prediction of estimated genomic breeding values (GEBV). In a simulation study, bias and accuracy of GEBV were explored for various genetic models with individuals selectively genotyped in a typical nucleus breeding program. We compared the performance of three existing methods, that is, Best Linear Unbiased Prediction of breeding values using pedigree-based relationships (PBLUP), genomic relationships for genotyped animals only (GBLUP) and a Single-Step approach (SSGBLUP) using both. For a scenario with no-selection and random mating (RR), prediction was unbiased. However, lower accuracy and bias were observed for scenarios with selection and random mating (SR) or selection and positive assortative mating (SA). As expected, bias disappeared when all individuals were genotyped and used in GBLUP. SSGBLUP showed higher accuracy compared to GBLUP, and bias of prediction was negligible with SR. However, PBLUP and SSGBLUP still showed bias in SA due to high inbreeding. SSGBLUP and PBLUP were unbiased provided that inbreeding was accounted for in the relationship matrices. Selective genotyping based on extreme phenotypic contrasts increased the prediction accuracy, but prediction was biased when using GBLUP. SSGBLUP could correct the biasedness while gaining higher accuracy than GBLUP. In a typical animal breeding program, where it is too expensive to genotype all animals, it would be appropriate to genotype phenotypically contrasting selection candidates and use a Single-Step approach to obtain accurate and unbiased prediction of GEBV.

Dominance genetic variation contributes little to the missing heritability for human complex traits.

For human complex traits, non-additive genetic variation has been invoked to explain "missing heritability," but its discovery is often neglected in genome-wide association studies. Here we propose a method of using SNP data to partition and estimate the proportion of phenotypic variance attributed to additive and dominance genetic variation at all SNPs (hSNP(2) and δSNP(2)) in unrelated individuals based on an orthogonal model where the estimate of hSNP(2) is independent of that of δSNP(2). With this method, we analyzed 79 quantitative traits in 6,715 unrelated European Americans. The estimate of δSNP(2) averaged across all the 79 quantitative traits was 0.03, approximately a fifth of that for additive variation (average hSNP(2) = 0.15). There were a few traits that showed substantial estimates of δSNP(2), none of which were replicated in a larger sample of 11,965 individuals. We further performed genome-wide association analyses of the 79 quantitative traits and detected SNPs with genome-wide significant dominance effects only at the ABO locus for factor VIII and von Willebrand factor. All these results suggest that dominance variation at common SNPs explains only a small fraction of phenotypic variation for human complex traits and contributes little to the missing narrow-sense heritability problem.

Simultaneous discovery, estimation and prediction analysis of complex traits using a bayesian mixture model.

Gene discovery, estimation of heritability captured by SNP arrays, inference on genetic architecture and prediction analyses of complex traits are usually performed using different statistical models and methods, leading to inefficiency and loss of power. Here we use a Bayesian mixture model that simultaneously allows variant discovery, estimation of genetic variance explained by all variants and prediction of unobserved phenotypes in new samples. We apply the method to simulated data of quantitative traits and Welcome Trust Case Control Consortium (WTCCC) data on disease and show that it provides accurate estimates of SNP-based heritability, produces unbiased estimators of risk in new samples, and that it can estimate genetic architecture by partitioning variation across hundreds to thousands of SNPs. We estimated that, depending on the trait, 2,633 to 9,411 SNPs explain all of the SNP-based heritability in the WTCCC diseases. The majority of those SNPs (>96%) had small effects, confirming a substantial polygenic component to common diseases. The proportion of the SNP-based variance explained by large effects (each SNP explaining 1% of the variance) varied markedly between diseases, ranging from almost zero for bipolar disorder to 72% for type 1 diabetes. Prediction analyses demonstrate that for diseases with major loci, such as type 1 diabetes and rheumatoid arthritis, Bayesian methods outperform profile scoring or mixed model approaches.

Performance of risk prediction for inflammatory bowel disease based on genotyping platform and genomic risk score method.

BACKGROUND: Predicting risk of disease from genotypes is being increasingly proposed for a variety of diagnostic and prognostic purposes. Genome-wide association studies (GWAS) have identified a large number of genome-wide significant susceptibility loci for Crohn's disease (CD) and ulcerative colitis (UC), two subtypes of inflammatory bowel disease (IBD). Recent studies have demonstrated that including only loci that are significantly associated with disease in the prediction model has low predictive power and that power can substantially be improved using a polygenic approach. METHODS: We performed a comprehensive analysis of risk prediction models using large case-control cohorts genotyped for 909,763 GWAS SNPs or 123,437 SNPs on the custom designed Immunochip using four prediction methods (polygenic score, best linear genomic prediction, elastic-net regularization and a Bayesian mixture model). We used the area under the curve (AUC) to assess prediction performance for discovery populations with different sample sizes and number of SNPs within cross-validation. RESULTS: On average, the Bayesian mixture approach had the best prediction performance. Using cross-validation we found little differences in prediction performance between GWAS and Immunochip, despite the GWAS array providing a 10 times larger effective genome-wide coverage. The prediction performance using Immunochip is largely due to the power of the initial GWAS for its marker selection and its low cost that enabled larger sample sizes. The predictive ability of the genomic risk score based on Immunochip was replicated in external data, with AUC of 0.75 for CD and 0.70 for UC. CD patients with higher risk scores demonstrated clinical characteristics typically associated with a more severe disease course including ileal location and earlier age at diagnosis. CONCLUSIONS: Our analyses demonstrate that the power of genomic risk prediction for IBD is mainly due to strongly associated SNPs with considerable effect sizes. Additional SNPs that are only tagged by high-density GWAS arrays and low or rare-variants over-represented in the high-density region on the Immunochip contribute little to prediction accuracy. Although a quantitative assessment of IBD risk for an individual is not currently possible, we show sufficient power of genomic risk scores to stratify IBD risk among individuals at diagnosis.

Effect of Coronary Artery Calcification Score by Lifestyle and Correlation With Coronary Artery Stenosis by Multidetector Computed Tomography.

OBJECTIVE: This study examined the effect of coronary artery calcification score by lifestyle and correlation with coronary artery stenosis in persons who underwent coronary artery computed tomography (CT) angiography among health examinees for heart diseases in Korea. METHODS AND MATERIALS: The study included 506 subjects (256 men and 250 women) who underwent coronary artery CT angiography among health examines for heart diseases at the Incheon Branch of the Korea Association of Health Promotion between January 2, 2014, and December 31, 2014. The demographical variables of the subjects were determined by frequency analysis, and the difference by sex was compared and analyzed using χ independence test. Independent 2-sample t test was performed to determine any difference in main factors by coronary artery calcification. RESULTS: According to the results, 175 (34.6%) had calcification, men showed statistically higher scores than women, and calcification seemed higher in those who were older, taller, heavier, and thicker in waist. Regarding blood pressure, calcification was shown if contraction phase and relaxation blood pressure was higher, blood sugar before meal was higher, and neutral fat was higher. By lifestyle, calcification seemed to be higher in those with more alcohol drinking per week, long past smoking years, and higher smoking amount per day in the past and present. In addition, coronary artery stenosis rate showed statistical correlation with calcification from the left anterior descending artery, right coronary artery, left circumflex artery, and left main coronary artery in sequence. CONCLUSIONS: In conclusion, coronary artery calcification score CT is deemed to be a suitable method for the estimation of coronary artery stenosis with short examination time, low radiation exposure, and noninvasive method.

Statistical power to detect genetic (co)variance of complex traits using SNP data in unrelated samples.

We have recently developed analysis methods (GREML) to estimate the genetic variance of a complex trait/disease and the genetic correlation between two complex traits/diseases using genome-wide single nucleotide polymorphism (SNP) data in unrelated individuals. Here we use analytical derivations and simulations to quantify the sampling variance of the estimate of the proportion of phenotypic variance captured by all SNPs for quantitative traits and case-control studies. We also derive the approximate sampling variance of the estimate of a genetic correlation in a bivariate analysis, when two complex traits are either measured on the same or different individuals. We show that the sampling variance is inversely proportional to the number of pairwise contrasts in the analysis and to the variance in SNP-derived genetic relationships. For bivariate analysis, the sampling variance of the genetic correlation additionally depends on the harmonic mean of the proportion of variance explained by the SNPs for the two traits and the genetic correlation between the traits, and depends on the phenotypic correlation when the traits are measured on the same individuals. We provide an online tool for calculating the power of detecting genetic (co)variation using genome-wide SNP data. The new theory and online tool will be helpful to plan experimental designs to estimate the missing heritability that has not yet been fully revealed through genome-wide association studies, and to estimate the genetic overlap between complex traits (diseases) in particular when the traits (diseases) are not measured on the same samples.

Estimation and partitioning of (co)heritability of inflammatory bowel disease from GWAS and immunochip data.

As custom arrays are cheaper than generic GWAS arrays, larger sample size is achievable for gene discovery. Custom arrays can tag more variants through denser genotyping of SNPs at associated loci, but at the cost of losing genome-wide coverage. Balancing this trade-off is important for maximizing experimental designs. We quantified both the gain in captured SNP-heritability at known candidate regions and the loss due to imperfect genome-wide coverage for inflammatory bowel disease using immunochip (iChip) and imputed GWAS data on 61,251 and 38.550 samples, respectively. For Crohn's disease (CD), the iChip and GWAS data explained 19 and 26% of variation in liability, respectively, and SNPs in the densely genotyped iChip regions explained 13% of the SNP-heritability for both the iChip and GWAS data. For ulcerative colitis (UC), the iChip and GWAS data explained 15 and 19% of variation in liability, respectively, and the dense iChip regions explained 10 and 9% of the SNP-heritability in the iChip and the GWAS data. From bivariate analyses, estimates of the genetic correlation in risk between CD and UC were 0.75 (SE 0.017) and 0.62 (SE 0.042) for the iChip and GWAS data, respectively. We also quantified the SNP-heritability of genomic regions that did or did not contain the previous 163 GWAS hits for CD and UC, and SNP-heritability of the overlapping loci between the densely genotyped iChip regions and the 163 GWAS hits. For both diseases, over different genomic partitioning, the densely genotyped regions on the iChip tagged at least as much variation in liability as in the corresponding regions in the GWAS data, however a certain amount of tagged SNP-heritability in the GWAS data was lost using the iChip due to the low coverage at unselected regions. These results imply that custom arrays with a GWAS backbone will facilitate more gene discovery, both at associated and novel loci.

Genome-wide association study reveals greater polygenic loading for schizophrenia in cases with a family history of illness.

Genome-wide association studies (GWAS) of schizophrenia have yielded more than 100 common susceptibility variants, and strongly support a substantial polygenic contribution of a large number of small allelic effects. It has been hypothesized that familial schizophrenia is largely a consequence of inherited rather than environmental factors. We investigated the extent to which familiality of schizophrenia is associated with enrichment for common risk variants detectable in a large GWAS. We analyzed single nucleotide polymorphism (SNP) data for cases reporting a family history of psychotic illness (N = 978), cases reporting no such family history (N = 4,503), and unscreened controls (N = 8,285) from the Psychiatric Genomics Consortium (PGC1) study of schizophrenia. We used a multinomial logistic regression approach with model-fitting to detect allelic effects specific to either family history subgroup. We also considered a polygenic model, in which we tested whether family history positive subjects carried more schizophrenia risk alleles than family history negative subjects, on average. Several individual SNPs attained suggestive but not genome-wide significant association with either family history subgroup. Comparison of genome-wide polygenic risk scores based on GWAS summary statistics indicated a significant enrichment for SNP effects among family history positive compared to family history negative cases (Nagelkerke's R(2 ) = 0.0021; P = 0.00331; P-value threshold <0.4). Estimates of variability in disease liability attributable to the aggregate effect of genome-wide SNPs were significantly greater for family history positive compared to family history negative cases (0.32 and 0.22, respectively; P = 0.031). We found suggestive evidence of allelic effects detectable in large GWAS of schizophrenia that might be specific to particular family history subgroups. However, consideration of a polygenic risk score indicated a significant enrichment among family history positive cases for common allelic effects. Familial illness might, therefore, represent a more heritable form of schizophrenia, as suggested by previous epidemiological studies.

Runs of homozygosity implicate autozygosity as a schizophrenia risk factor.

Autozygosity occurs when two chromosomal segments that are identical from a common ancestor are inherited from each parent. This occurs at high rates in the offspring of mates who are closely related (inbreeding), but also occurs at lower levels among the offspring of distantly related mates. Here, we use runs of homozygosity in genome-wide SNP data to estimate the proportion of the autosome that exists in autozygous tracts in 9,388 cases with schizophrenia and 12,456 controls. We estimate that the odds of schizophrenia increase by ~17% for every 1% increase in genome-wide autozygosity. This association is not due to one or a few regions, but results from many autozygous segments spread throughout the genome, and is consistent with a role for multiple recessive or partially recessive alleles in the etiology of schizophrenia. Such a bias towards recessivity suggests that alleles that increase the risk of schizophrenia have been selected against over evolutionary time.

Estimating missing heritability for disease from genome-wide association studies.

Genome-wide association studies are designed to discover SNPs that are associated with a complex trait. Employing strict significance thresholds when testing individual SNPs avoids false positives at the expense of increasing false negatives. Recently, we developed a method for quantitative traits that estimates the variation accounted for when fitting all SNPs simultaneously. Here we develop this method further for case-control studies. We use a linear mixed model for analysis of binary traits and transform the estimates to a liability scale by adjusting both for scale and for ascertainment of the case samples. We show by theory and simulation that the method is unbiased. We apply the method to data from the Wellcome Trust Case Control Consortium and show that a substantial proportion of variation in liability for Crohn disease, bipolar disorder, and type I diabetes is tagged by common SNPs.

Research review: Polygenic methods and their application to psychiatric traits.

BACKGROUND: Despite evidence from twin and family studies for an important contribution of genetic factors to both childhood and adult onset psychiatric disorders, identifying robustly associated specific DNA variants has proved challenging. In the pregenomics era the genetic architecture (number, frequency and effect size of risk variants) of complex genetic disorders was unknown. Empirical evidence for the genetic architecture of psychiatric disorders is emerging from the genetic studies of the last 5 years. METHODS AND SCOPE: We review the methods investigating the polygenic nature of complex disorders. We provide mini-guides to genomic profile (or polygenic) risk scoring and to estimation of variance (or heritability) from common SNPs; a glossary of key terms is also provided. We review results of applications of the methods to psychiatric disorders and related traits and consider how these methods inform on missing heritability, hidden heritability and still-missing heritability. FINDINGS: Genome-wide genotyping and sequencing studies are providing evidence that psychiatric disorders are truly polygenic, that is they have a genetic architecture of many genetic variants, including risk variants that are both common and rare in the population. Sample sizes published to date are mostly underpowered to detect effect sizes of the magnitude presented by nature, and these effect sizes may be constrained by the biological validity of the diagnostic constructs. CONCLUSIONS: Increasing the sample size for genome wide association studies of psychiatric disorders will lead to the identification of more associated genetic variants, as already found for schizophrenia. These loci provide the starting point of functional analyses that might eventually lead to new prevention and treatment options and to improved biological validity of diagnostic constructs. Polygenic analyses will contribute further to our understanding of complex genetic traits as sample sizes increase and as sample resources become richer in phenotypic descriptors, both in terms of clinical symptoms and of nongenetic risk factors.

Alcohol consumption and the risk of all-cause and cause-specific mortality-a linear and nonlinear Mendelian randomization study.

BACKGROUND: Many observational studies support light-to-moderate alcohol intake as potentially protective against premature death. We used a genetic approach to evaluate the linear and nonlinear relationships between alcohol consumption and mortality from different underlying causes. METHODS: We used data from 278 093 white-British UK Biobank participants, aged 37-73 years at recruitment and with data on alcohol intake, genetic variants, and mortality. Habitual alcohol consumption was instrumented by 94 variants. Linear Mendelian randomization (MR) analyses were conducted using five complementary approaches, and nonlinear MR analyses by the doubly-ranked method. RESULTS: There were 20 834 deaths during the follow-up (median 12.6 years). In conventional analysis, the association between alcohol consumption and mortality outcomes was 'J-shaped'. In contrast, MR analyses supported a positive linear association with premature mortality, with no evidence for curvature (Pnonlinearity ≥ 0.21 for all outcomes). The odds ratio [OR] for each standard unit increase in alcohol intake was 1.27 (95% confidence interval [CI] 1.16-1.39) for all-cause mortality, 1.30 (95% CI 1.10-1.53) for cardiovascular disease, 1.20 (95% CI 1.08-1.33) for cancer, and 2.06 (95% CI 1.36-3.12) for digestive disease mortality. These results were consistent across pleiotropy-robust methods. There was no clear evidence for an association between alcohol consumption and mortality from respiratory diseases or COVID-19 (1.32, 95% CI 0.96-1.83 and 1.46, 95% CI 0.99-2.16, respectively; Pnonlinearity ≥ 0.21). CONCLUSION: Higher levels of genetically predicted alcohol consumption had a strong linear association with an increased risk of premature mortality with no evidence for any protective benefit at modest intake levels.

Genetic and brain similarity independently predict childhood anthropometrics and neighborhood socioeconomic conditions.

Linking the developing brain with individual differences in clinical and demographic traits is challenging due to the substantial interindividual heterogeneity of brain anatomy and organization. Here we employ an integrative approach that parses individual differences in both cortical thickness and common genetic variants, and assess their effects on a wide set of childhood traits. The approach uses a linear mixed model framework to obtain the unique effects of each type of similarity, as well as their covariance. We employ this approach in a sample of 7760 unrelated children in the ABCD cohort baseline sample (mean age 9.9, 46.8% female). In general, associations between cortical thickness similarity and traits were limited to anthropometrics such as height, weight, and birth weight, as well as a marker of neighborhood socioeconomic conditions. Common genetic variants explained significant proportions of variance across nearly all included outcomes, although estimates were somewhat lower than previous reports. No significant covariance of the effects of genetic and cortical thickness similarity was found. The present findings highlight the connection between anthropometrics as well as neighborhood socioeconomic conditions and the developing brain, which appear to be independent from individual differences in common genetic variants in this population-based sample.

A better coefficient of determination for genetic profile analysis.

Genome-wide association studies have facilitated the construction of risk predictors for disease from multiple Single Nucleotide Polymorphism markers. The ability of such "genetic profiles" to predict outcome is usually quantified in an independent data set. Coefficients of determination (R(2) ) have been a useful measure to quantify the goodness-of-fit of the genetic profile. Various pseudo-R(2) measures for binary responses have been proposed. However, there is no standard or consensus measure because the concept of residual variance is not easily defined on the observed probability scale. Unlike other nongenetic predictors such as environmental exposure, there is prior information on genetic predictors because for most traits there are estimates of the proportion of variation in risk in the population due to all genetic factors, the heritability. It is this useful ability to benchmark that makes the choice of a measure of goodness-of-fit in genetic profiling different from that of nongenetic predictors. In this study, we use a liability threshold model to establish the relationship between the observed probability scale and underlying liability scale in measuring R(2) for binary responses. We show that currently used R(2) measures are difficult to interpret, biased by ascertainment, and not comparable to heritability. We suggest a novel and globally standard measure of R(2) that is interpretable on the liability scale. Furthermore, even when using ascertained case-control studies that are typical in human disease studies, we can obtain an R(2) measure on the liability scale that can be compared directly to heritability.