Estimating cancer risk in carriers of Lynch syndrome variants in UK Biobank.

BackgroundLynch syndrome (LS) is an inherited cancer predisposition syndrome caused by genetic variants affecting DNA mismatch repair (MMR) genes MLH1, MSH2, MSH6 and PMS2 Cancer risk in LS is estimated from cohorts of individuals ascertained by individual or family history of cancer, which may upwardly bias estimates. METHODS: 830 carriers of pathogenic or likely pathogenic (path_MMR) MMR gene variants classified by InSiGHT were identified in 454 756 UK Biobank (UKB) participants using whole-exome sequence. Nelson-Aalen survival analysis was used to estimate cumulative incidence of colorectal, endometrial and breast cancer (BC). RESULTS: Cumulative incidence of colorectal and endometrial cancer (EC) by age 70 years was elevated in path_MMR carriers compared with non-carriers (colorectal: 11.8% (95% confidence interval (CI): 9.5% to 14.6%) vs 1.7% (95% CI: 1.6% to 1.7%), endometrial: 13.4% (95% CI: 10.2% to 17.6%) vs 1.0% (95% CI: 0.9% to 1.0%)), but the magnitude of this increase differed between genes. Cumulative BC incidence by age 70 years was not elevated in path_MMR carriers compared with non-carriers (8.9% (95% CI: 6.3% to 12.4%) vs 7.5% (95% CI: 7.4% to 7.6%)). Cumulative cancer incidence estimates in UKB were similar to estimates from the Prospective Lynch Syndrome Database for all genes and cancers, except there was no evidence for elevated EC risk in carriers of pathogenic PMS2 variants in UKB. CONCLUSION: These results support offering incidentally identified carriers of any path_MMR surveillance to manage colorectal cancer risk. Incidentally identified carriers of pathogenic variants in MLH1, MSH2 and MSH6 would also benefit from interventions to reduce EC risk. The results suggest that BC is not an LS-related cancer.

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep.

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.

Genome-wide methylation data improves dissection of the effect of smoking on body mass index.

Variation in obesity-related traits has a genetic basis with heritabilities between 40 and 70%. While the global obesity pandemic is usually associated with environmental changes related to lifestyle and socioeconomic changes, most genetic studies do not include all relevant environmental covariates, so the genetic contribution to variation in obesity-related traits cannot be accurately assessed. Some studies have described interactions between a few individual genes linked to obesity and environmental variables but there is no agreement on their total contribution to differences between individuals. Here we compared self-reported smoking data and a methylation-based proxy to explore the effect of smoking and genome-by-smoking interactions on obesity related traits from a genome-wide perspective to estimate the amount of variance they explain. Our results indicate that exploiting omic measures can improve models for complex traits such as obesity and can be used as a substitute for, or jointly with, environmental records to better understand causes of disease.

Linking protein to phenotype with Mendelian Randomization detects 38 proteins with causal roles in human diseases and traits.

To efficiently transform genetic associations into drug targets requires evidence that a particular gene, and its encoded protein, contribute causally to a disease. To achieve this, we employ a three-step proteome-by-phenome Mendelian Randomization (MR) approach. In step one, 154 protein quantitative trait loci (pQTLs) were identified and independently replicated. From these pQTLs, 64 replicated locally-acting variants were used as instrumental variables for proteome-by-phenome MR across 846 traits (step two). When its assumptions are met, proteome-by-phenome MR, is equivalent to simultaneously running many randomized controlled trials. Step 2 yielded 38 proteins that significantly predicted variation in traits and diseases in 509 instances. Step 3 revealed that amongst the 271 instances from GeneAtlas (UK Biobank), 77 showed little evidence of pleiotropy (HEIDI), and 92 evidence of colocalization (eCAVIAR). Results were wide ranging: including, for example, new evidence for a causal role of tyrosine-protein phosphatase non-receptor type substrate 1 (SHPS1; SIRPA) in schizophrenia, and a new finding that intestinal fatty acid binding protein (FABP2) abundance contributes to the pathogenesis of cardiovascular disease. We also demonstrated confirmatory evidence for the causal role of four further proteins (FGF5, IL6R, LPL, LTA) in cardiovascular disease risk.

Effect of genotyping strategies on the sustained benefit of single-step genomic BLUP over multiple generations.

BACKGROUND: Single-step genomic best linear unbiased prediction (ssGBLUP) allows the inclusion of information from genotyped and ungenotyped individuals in a single analysis. This avoids the need to genotype all candidates with the potential benefit of reducing overall costs. The aim of this study was to assess the effect of genotyping strategies, the proportion of genotyped candidates and the genotyping criterion to rank candidates to be genotyped, when using ssGBLUP evaluation. A simulation study was carried out assuming selection over several discrete generations where a proportion of the candidates were genotyped and evaluation was done using ssGBLUP. The scenarios compared were: (i) three genotyping strategies defined by their protocol for choosing candidates to be genotyped (RANDOM: candidates were chosen at random; TOP: candidates with the best genotyping criterion were genotyped; and EXTREME: candidates with the best and worse criterion were genotyped); (ii) eight proportions of genotyped candidates (p); and (iii) two genotyping criteria to rank candidates to be genotyped (candidates' own phenotype or estimated breeding values). The criteria of the comparison were the cumulated gain and reliability of the genomic estimated breeding values (GEBV). RESULTS: The genotyping strategy with the greatest cumulated gain was TOP followed by RANDOM, with EXTREME behaving as RANDOM at low p and as TOP with high p. However, the reliability of GEBV was higher with RANDOM than with TOP. This disparity between the trend of the gain and the reliability is due to the TOP scheme genotyping the candidates with the greater chances of being selected. The extra gain obtained with TOP increases when the accuracy of the selection criterion to rank candidates to be genotyped increases. CONCLUSIONS: The best strategy to maximise genetic gain when only a proportion of the candidates are to be genotyped is TOP, since it prioritises the genotyping of candidates which are more likely to be selected. However, the strategy with the greatest GEBV reliability does not achieve the largest gain, thus reliability cannot be considered as an absolute and sufficient criterion for determining the scheme which maximises genetic gain.

The genetic landscape of Scotland and the Isles.

Britain and Ireland are known to show population genetic structure; however, large swathes of Scotland, in particular, have yet to be described. Delineating the structure and ancestry of these populations will allow variant discovery efforts to focus efficiently on areas not represented in existing cohorts. Thus, we assembled genotype data for 2,554 individuals from across the entire archipelago with geographically restricted ancestry, and performed population structure analyses and comparisons to ancient DNA. Extensive geographic structuring is revealed, from broad scales such as a NE to SW divide in mainland Scotland, through to the finest scale observed to date: across 3 km in the Northern Isles. Many genetic boundaries are consistent with Dark Age kingdoms of Gaels, Picts, Britons, and Norse. Populations in the Hebrides, the Highlands, Argyll, Donegal, and the Isle of Man show characteristics of isolation. We document a pole of Norwegian ancestry in the north of the archipelago (reaching 23 to 28% in Shetland) which complements previously described poles of Germanic ancestry in the east, and "Celtic" to the west. This modern genetic structure suggests a northwestern British or Irish source population for the ancient Gaels that contributed to the founding of Iceland. As rarer variants, often with larger effect sizes, become the focus of complex trait genetics, more diverse rural cohorts may be required to optimize discoveries in British and Irish populations and their considerable global diaspora.

New App-Based Dietary and Lifestyle Intervention on Weight Loss and Cardiovascular Health.

Consumer digital technology is rapidly evolving, allowing users to manage their health in a simple, non-invasive manner. However, there are few studies revealing if using digital technology as part of an intervention really has an impact in consumer health compared with traditional strategies. The objective of the current study is to compare two groups (MTB; n = 18, 46.1 ± 10.4 years and MTBAPP; n = 19, 45.3 ± 6.40 years) of overweight, prehypertensive individuals in losing weight and lowering their blood pressure. Both were provided with nutritionist-guided recommendations, a wearable tracking device and a dietary supplement that has previously been proven to help lose body weight and lower blood pressure. In addition, one of the groups (MTBAPP) used a mobile app specifically designed for the intervention. Blood pressure, body composition, triglyceride level, peak expiratory flow, forced expiratory volume in the first second and maximum oxygen volume were measured at different time points. In addition, participants were monitored with an activity bracelet throughout the intervention. As a result, both groups significantly lost body weight, while the group using the app additionally improved blood pressure levels and lowered fat mass. Furthermore, the app users significantly increased the number of daily steps and decreased sedentary time. In conclusion, the addition of a mobile app with daily reminders to follow healthy lifestyle recommendations increased physical activity and overall improved blood pressure and fat mass levels when compared with a group performing the same intervention but in absence of the mobile application.

Genomic regions influencing intramuscular fat in divergently selected rabbit lines.

Intramuscular fat (IMF) is one of the main meat quality traits for breeding programmes in livestock species. The main objective of this study was to identify genomic regions associated with IMF content comparing two rabbit populations divergently selected for this trait, and to generate a list of putative candidate genes. Animals were genotyped using the Affymetrix Axiom OrcunSNP Array (200k). After quality control, the data involved 477 animals and 93 540 SNPs. Two methods were used in this research: single marker regressions with the data adjusted by genomic relatedness, and a Bayesian multiple marker regression. Associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU8 and OCU13. The highest value for the percentage of the genomic variance explained by a genomic region was found in two consecutive genomic windows on OCU8 (7.34%). Genes in the associated regions of OCU1 and OCU8 presented biological functions related to the control of adipose cell function, lipid binding, transportation and localisation (APOLD1, PLBD1, PDE6H, GPRC5D and GPRC5A) and lipid metabolic processes (MTMR2). The EWSR1 gene, underlying the OCU13 region, is linked to the development of brown adipocytes. The findings suggest that there is a large component of polygenic effect behind the differences in IMF content in these two lines, as the variance explained by most of the windows was low. The genomic regions of OCU1, OCU8 and OCU13 revealed novel candidate genes. Further studies would be needed to validate the associations and explore their possible application in selection programmes.

Correction: Pedigree- and SNP-Associated Genetics and Recent Environment are the Major Contributors to Anthropometric and Cardiometabolic Trait Variation.

[This corrects the article DOI: 10.1371/journal.pgen.1005804.].

Genomic analysis of family data reveals additional genetic effects on intelligence and personality.

Pedigree-based analyses of intelligence have reported that genetic differences account for 50-80% of the phenotypic variation. For personality traits these effects are smaller, with 34-48% of the variance being explained by genetic differences. However, molecular genetic studies using unrelated individuals typically report a heritability estimate of around 30% for intelligence and between 0 and 15% for personality variables. Pedigree-based estimates and molecular genetic estimates may differ because current genotyping platforms are poor at tagging causal variants, variants with low minor allele frequency, copy number variants, and structural variants. Using ~20,000 individuals in the Generation Scotland family cohort genotyped for ~700,000 single-nucleotide polymorphisms (SNPs), we exploit the high levels of linkage disequilibrium (LD) found in members of the same family to quantify the total effect of genetic variants that are not tagged in GWAS of unrelated individuals. In our models, genetic variants in low LD with genotyped SNPs explain over half of the genetic variance in intelligence, education, and neuroticism. By capturing these additional genetic effects our models closely approximate the heritability estimates from twin studies for intelligence and education, but not for neuroticism and extraversion. We then replicated our finding using imputed molecular genetic data from unrelated individuals to show that ~50% of differences in intelligence, and ~40% of the differences in education, can be explained by genetic effects when a larger number of rare SNPs are included. From an evolutionary genetic perspective, a substantial contribution of rare genetic variants to individual differences in intelligence, and education is consistent with mutation-selection balance.

Pedigree- and SNP-Associated Genetics and Recent Environment are the Major Contributors to Anthropometric and Cardiometabolic Trait Variation.

Genome-wide association studies have successfully identified thousands of loci for a range of human complex traits and diseases. The proportion of phenotypic variance explained by significant associations is, however, limited. Given the same dense SNP panels, mixed model analyses capture a greater proportion of phenotypic variance than single SNP analyses but the total is generally still less than the genetic variance estimated from pedigree studies. Combining information from pedigree relationships and SNPs, we examined 16 complex anthropometric and cardiometabolic traits in a Scottish family-based cohort comprising up to 20,000 individuals genotyped for ~520,000 common autosomal SNPs. The inclusion of related individuals provides the opportunity to also estimate the genetic variance associated with pedigree as well as the effects of common family environment. Trait variation was partitioned into SNP-associated and pedigree-associated genetic variation, shared nuclear family environment, shared couple (partner) environment and shared full-sibling environment. Results demonstrate that trait heritabilities vary widely but, on average across traits, SNP-associated and pedigree-associated genetic effects each explain around half the genetic variance. For most traits the recently-shared environment of couples is also significant, accounting for ~11% of the phenotypic variance on average. On the other hand, the environment shared largely in the past by members of a nuclear family or by full-siblings, has a more limited impact. Our findings point to appropriate models to use in future studies as pedigree-associated genetic effects and couple environmental effects have seldom been taken into account in genotype-based analyses. Appropriate description of the trait variation could help understand causes of intra-individual variation and in the detection of contributing loci and environmental factors.

A meta-analysis of 120 246 individuals identifies 18 new loci for fibrinogen concentration.

Genome-wide association studies have previously identified 23 genetic loci associated with circulating fibrinogen concentration. These studies used HapMap imputation and did not examine the X-chromosome. 1000 Genomes imputation provides better coverage of uncommon variants, and includes indels. We conducted a genome-wide association analysis of 34 studies imputed to the 1000 Genomes Project reference panel and including ∼120 000 participants of European ancestry (95 806 participants with data on the X-chromosome). Approximately 10.7 million single-nucleotide polymorphisms and 1.2 million indels were examined. We identified 41 genome-wide significant fibrinogen loci; of which, 18 were newly identified. There were no genome-wide significant signals on the X-chromosome. The lead variants of five significant loci were indels. We further identified six additional independent signals, including three rare variants, at two previously characterized loci: FGB and IRF1. Together the 41 loci explain 3% of the variance in plasma fibrinogen concentration.

Identification of a novel locus on chromosome 2q13, which predisposes to clinical vertebral fractures independently of bone density.

OBJECTIVES: To identify genetic determinants of susceptibility to clinical vertebral fractures, which is an important complication of osteoporosis. METHODS: Here we conduct a genome-wide association study in 1553 postmenopausal women with clinical vertebral fractures and 4340 controls, with a two-stage replication involving 1028 cases and 3762 controls. Potentially causal variants were identified using expression quantitative trait loci (eQTL) data from transiliac bone biopsies and bioinformatic studies. RESULTS: A locus tagged by rs10190845 was identified on chromosome 2q13, which was significantly associated with clinical vertebral fracture (P=1.04×10-9) with a large effect size (OR 1.74, 95% CI 1.06 to 2.6). Bioinformatic analysis of this locus identified several potentially functional SNPs that are associated with expression of the positional candidate genes TTL (tubulin tyrosine ligase) and SLC20A1 (solute carrier family 20 member 1). Three other suggestive loci were identified on chromosomes 1p31, 11q12 and 15q11. All these loci were novel and had not previously been associated with bone mineral density or clinical fractures. CONCLUSION: We have identified a novel genetic variant that is associated with clinical vertebral fractures by mechanisms that are independent of BMD. Further studies are now in progress to validate this association and evaluate the underlying mechanism.

Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis.

BACKGROUND: Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes. METHODS & FINDINGS: Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants. CONCLUSIONS: As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Genomic prediction of complex human traits: relatedness, trait architecture and predictive meta-models.

We explore the prediction of individuals' phenotypes for complex traits using genomic data. We compare several widely used prediction models, including Ridge Regression, LASSO and Elastic Nets estimated from cohort data, and polygenic risk scores constructed using published summary statistics from genome-wide association meta-analyses (GWAMA). We evaluate the interplay between relatedness, trait architecture and optimal marker density, by predicting height, body mass index (BMI) and high-density lipoprotein level (HDL) in two data cohorts, originating from Croatia and Scotland. We empirically demonstrate that dense models are better when all genetic effects are small (height and BMI) and target individuals are related to the training samples, while sparse models predict better in unrelated individuals and when some effects have moderate size (HDL). For HDL sparse models achieved good across-cohort prediction, performing similarly to the GWAMA risk score and to models trained within the same cohort, which indicates that, for predicting traits with moderately sized effects, large sample sizes and familial structure become less important, though still potentially useful. Finally, we propose a novel ensemble of whole-genome predictors with GWAMA risk scores and demonstrate that the resulting meta-model achieves higher prediction accuracy than either model on its own. We conclude that although current genomic predictors are not accurate enough for diagnostic purposes, performance can be improved without requiring access to large-scale individual-level data. Our methodologically simple meta-model is a means of performing predictive meta-analysis for optimizing genomic predictions and can be easily extended to incorporate multiple population-level summary statistics or other domain knowledge.

Meta-analysis of 49†549 individuals imputed with the 1000 Genomes Project reveals an exonic damaging variant in ANGPTL4 determining fasting TG levels.

BACKGROUND: So far, more than 170 loci have been associated with circulating lipid levels through genome-wide association studies (GWAS). These associations are largely driven by common variants, their function is often not known, and many are likely to be markers for the causal variants. In this study we aimed to identify more new rare and low-frequency functional variants associated with circulating lipid levels. METHODS: We used the 1000 Genomes Project as a reference panel for the imputations of GWAS data from ∼60 000 individuals in the discovery stage and ∼90 000 samples in the replication stage. RESULTS: Our study resulted in the identification of five new associations with circulating lipid levels at four loci. All four loci are within genes that can be linked biologically to lipid metabolism. One of the variants, rs116843064, is a damaging missense variant within the ANGPTL4 gene. CONCLUSIONS: This study illustrates that GWAS with high-scale imputation may still help us unravel the biological mechanism behind circulating lipid levels.

A general approach for haplotype phasing across the full spectrum of relatedness.

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally 'unrelated' individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

The effect of increasing lifespan and recycling rate on carbon storage in wood products from theoretical model to application for the European wood sector.

The use of wood products is often promoted as a climate change mitigation option to reduce atmospheric carbon dioxide concentrations. In previous literature, we identified longevity and recycling rate as two determining factors that influence the carbon stock in wood products, but no studies have predicted the effect of improved wood use on carbon storage over time. In this study, we aimed at evaluating changes in the lifespan and the recycling rate as two options for enhancing carbon stock in wood products for different time horizons. We first explored the behaviour over time of both factors in a theoretical simulation, and then calculated their effect for the European wood sector of the future. The theoretical simulation shows that the carbon stock in wood products increases linearly when increasing the average lifespan of wood products and exponentially when improving the recycling rate. The emissions savings under the current use of wood products in Europe in 2030 were estimated at 57.65 Mt carbon dioxide (CO2) per year. This amount could be increased 5 Mt CO2 if average lifespan increased 19.54 % or if recycling rate increased 20.92 % in 2017. However, the combination of both strategies could increase the emissions saving almost 5 Mt CO2 more by 2030. Incrementing recycling rate of paper and paperboard is the best short-term strategy (2030) to reduce emissions, but elongating average lifespan of wood-based panels is a better strategy for longer term periods (2046).

Modelling carbon stocks and fluxes in the wood product sector: a comparative review.

In addition to forest ecosystems, wood products are carbon pools that can be strategically managed to mitigate climate change. Wood product models (WPMs) simulating the carbon balance of wood production, use and end of life can complement forest growth models to evaluate the mitigation potential of the forest sector as a whole. WPMs can be used to compare scenarios of product use and explore mitigation strategies. A considerable number of WPMs have been developed in the last three decades, but there is no review available analysing their functionality and performance. This study analyses and compares 41 WPMs. One surprising initial result was that we discovered the erroneous implementation of a few concepts and assumptions in some of the models. We further described and compared the models using six model characteristics (bucking allocation, industrial processes, carbon pools, product removal, recycling and substitution effects) and three model-use characteristics (system boundaries, model initialization and evaluation of results). Using a set of indicators based on the model characteristics, we classified models using a hierarchical clustering technique and differentiated them according to their increasing degrees of complexity and varying levels of user support. For purposes of simulating carbon stock in wood products, models with a simple structure may be sufficient, but to compare climate change mitigation options, complex models are needed. The number of models has increased substantially over the last ten years, introducing more diversity and accuracy. Calculation of substitution effects and recycling has also become more prominent. However, the lack of data is still an important constraint for a more realistic estimation of carbon stocks and fluxes. Therefore, if the sector wants to demonstrate the environmental quality of its products, it should make it a priority to provide reliable life cycle inventory data, particularly regarding aspects of time and location.

Meta-analysis of genome-wide association studies identifies six new Loci for serum calcium concentrations.

Calcium is vital to the normal functioning of multiple organ systems and its serum concentration is tightly regulated. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤ 21,679 additional individuals. Seven loci (six new regions) in association with serum calcium were identified and replicated. Rs1570669 near CYP24A1 (P = 9.1E-12), rs10491003 upstream of GATA3 (P = 4.8E-09) and rs7481584 in CARS (P = 1.2E-10) implicate regions involved in Mendelian calcemic disorders: Rs1550532 in DGKD (P = 8.2E-11), also associated with bone density, and rs7336933 near DGKH/KIAA0564 (P = 9.1E-10) are near genes that encode distinct isoforms of diacylglycerol kinase. Rs780094 is in GCKR. We characterized the expression of these genes in gut, kidney, and bone, and demonstrate modulation of gene expression in bone in response to dietary calcium in mice. Our results shed new light on the genetics of calcium homeostasis.