Gastrectomy (Gsx) is associated with altered emotional function and a predisposition to depression/anxiety disorders. Here we investigated the effects of Gsx on emotional reactivity in rats and explored the underlying neurobiological mechanisms. Gsx- and sham-operated rats were exposed to behavioural tests that explore anxiety- and depression-like behaviour (open field, black and white box, elevated plus maze, social interaction, forced swim) as well as memory (object recognition). The potential neurobiological mechanisms underlying these differences were explored by measuring (i) turnover of candidate neurotransmitter systems in the nucleus accumbens, (ii) hippocampal neurogenesis by BrdU labelling or by analysis of candidate genes involved in neuronal growth and (iii) changes in mRNA expression of candidate genes in dissected hippocampal and amygdala tissue. Data from individual behavioural tests as well as from multivariate analysis revealed differing emotional reactivity between Gsx- and sham-operated rats. Gsx rats showed reduced emotional reactivity in a new environment and decreased depression-like behaviour. Accumbal serotonin and dopamine turnover were both reduced in Gsx rats. Gsx also led to a memory deficit, although hippocampal neurogenesis was unaffected. Of the many candidate genes studied by real-time RT-PCR, we highlight a Gsx-associated decrease in expression of Egr-1, a transcription factor linked to neural plasticity and cognition, in the hippocampus and amygdala. Thus, Gsx induces an alteration of emotional reactivity and a memory/cognitive deficit that is associated with reduced turnover of serotonin and dopamine in the nucleus accumbens and decreased expression of Egr-1 in the hippocampus and amygdala.

Knowing the phase of marker genotype data can be useful in genome-wide association studies, because it makes it possible to use analysis frameworks that account for identity by descent or parent of origin of alleles and it can lead to a large increase in data quantities via genotype or sequence imputation. Long-range phasing and haplotype library imputation constitute a fast and accurate method to impute phase for SNP data.

Current genome-wide association studies (GWAS) are moving towards the use of large cohorts of primary cell lines to study a disease of interest and to assign biological relevance to the genetic signals identified. Here, we use a panel of human osteoblasts (HObs) to carry out a transcriptomic survey, similar to recent studies in lymphoblastoid cell lines (LCLs). The distinct nature of HObs and LCLs is reflected by the preferential grouping of cell type-specific genes within biologically and functionally relevant pathways unique to each tissue type. We performed cis-association analysis with SNP genotypes to identify genetic variations of transcript isoforms, and our analysis indicates that differential expression of transcript isoforms in HObs is also partly controlled by cis-regulatory genetic variants. These isoforms are regulated by genetic variants in both a tissue-specific and tissue-independent fashion, and these associations have been confirmed by RT-PCR validation. Our study suggests that multiple transcript isoforms are often present in both tissues and that genetic control may affect the relative expression of one isoform to another, rather than having an all-or-none effect. Examination of the top SNPs from a GWAS of bone mineral density show overlap with probeset associations observed in this study. The top hit corresponding to the FAM118A gene was tested for association studies in two additional clinical studies, revealing a novel transcript isoform variant. Our approach to examining transcriptome variation in multiple tissue types is useful for detecting the proportion of genetic variation common to different cell types and for the identification of cell-specific isoform variants that may be functionally relevant, an important follow-up step for GWAS.

Runs of homozygosity (ROH) are extended tracts of adjacent homozygous single nucleotide polymorphisms (SNPs) that are more common in unrelated individuals than previously thought. It has been proposed that estimating ROH on a genome-wide level, by making use of the genome-wide single nucleotide polymorphism (SNP) data, will enable to indentify recessive variants underlying complex traits. Here, we examined ROH larger than 1.5 Mb individually and in combination for association with survival in 5974 participants of the Rotterdam Study. In addition, we assessed the role of overall homozygosity, expressed as a percentage of the autosomal genome that is in ROH longer than 1.5 Mb, on survival during a mean follow-up period of 12 years. None of these measures of homozygosity was associated with survival to old age.

The third generation selective estrogen receptor modulators lasofoxifene (las) and bazedoxifene (bza) are indicated for treatment of postmenopausal osteoporosis. 17β-Estradiol (E2) and the second generation SERM raloxifene (ral) have major effects on the immune system, particularly on B cells. Treatment with E2 or ral inhibits B lymphopoiesis and treatment with E2, but not ral, stimulates antibody production. The effects of las and bza on the immune system have not been studied. Therefore, the aim of this study was to investigate their role in B cell development, maturation, and function. C57BL/6 mice were sham-operated or ovariectomized (ovx) and treated with vehicle, E2, ral, las, or bza. All substances increased total bone mineral density in ovx mice, as measured by peripheral quantitative computed tomography. In uterus, bza alone lacked agonistic effect in ovx mice and even acted as an antagonist in sham mice. As expected, E2 decreased B cell numbers at all developmental stages from pre-BI cells (in bone marrow) to transitional 1 (T1) B cells (in spleen) and increased marginal zone (MZ) B cells as determined by flow cytometry. However, treatment with las or bza only decreased the last stages of bone marrow B cell development and splenic T1 B cells, but had no effect MZ B cells. E2 increased antibody-producing cells quantified by ELISPOT, but las or bza did not. In conclusion, las and bza differ from E2 by retaining normal number of cells at most B cell stages during B lymphopoiesis and maturation and by not increasing antibody-producing cells.

Hydroxysteroid (17β)-dehydrogenase type 1 (HSD17B1) catalyzes the conversion of low active 17-ketosteroids, androstenedione (A-dione) and estrone (E1) to highly active 17-hydroxysteroids, testosterone (T) and E2, respectively. In this study, the importance of HSD17B1 in ovarian estrogen production was determined using Hsd17b1 knockout (HSD17B1KO) mice. In these mice, the ovarian HSD17B enzyme activity was markedly reduced, indicating a central role of HSD17B1 in ovarian physiology. The lack of Hsd17b activity resulted in increased ovarian E1:E2 and A-dione:T ratios, but we also observed reduced progesterone concentration in HSD17B1KO ovaries. Accordingly with the altered steroid production, altered expression of Star, Cyp11a1, Lhcgr, Hsd17b7, and especially Cyp17a1 was observed. The ovaries of HSD17B1KO mice presented with all stages of folliculogenesis, while the corpus luteum structure was less defined and number reduced. Surprisingly, bundles of large granular cells of unknown origin appeared in the stroma of the KO ovaries. The HSD17B1KO mice presented with severe subfertility and failed to initiate pseudopregnancy. However, the HSD17B1KO females presented with normal estrous cycle defined by vaginal smears and normal puberty appearance. This study indicates that HSD17B1 is a key enzyme in ovarian steroidogenesis and has a novel function in initiation and stabilization of pregnancy.

The Generation Scotland Scottish Family Health Study (GS:SFHS) includes 23,960 participants from across Scotland with records for many health-related traits and environmental covariates. Genotypes at ~700 K SNPs are currently available for 10,000 participants. The cohort was designed as a resource for genetic and health related research and the study of complex traits. In this study we developed a suite of analyses to disentangle the genomic differentiation within GS:SFHS individuals to describe and optimise the sample and methods for future analyses.

Variants associated with blood lipid levels may be population-specific. To identify low-frequency variants associated with this phenotype, population-specific reference panels may be used. Here we impute nine large Dutch biobanks (~35,000 samples) with the population-specific reference panel created by the Genome of The Netherlands Project and perform association testing with blood lipid levels. We report the discovery of five novel associations at four loci (P value <6.61 × 10(-4)), including a rare missense variant in ABCA6 (rs77542162, p.Cys1359Arg, frequency 0.034), which is predicted to be deleterious. The frequency of this ABCA6 variant is 3.65-fold increased in the Dutch and its effect (βLDL-C=0.135, βTC=0.140) is estimated to be very similar to those observed for single variants in well-known lipid genes, such as LDLR.

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.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: β = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Treatment with anti-inflammatory glucocorticoids is associated with osteoporosis. Many of the treated patients are postmenopausal women, who even without treatment have an increased risk of osteoporosis. Lymphocytes have been shown to play a role in postmenopausal and arthritis-induced osteoporosis, and they are targeted by glucocorticoids. The aim of this study was to investigate the mechanisms behind effects of glucocorticoids on bone during health and menopause, focusing on lymphocytes. Female C57BL/6 or SCID mice were therefore sham-operated or ovariectomized and 2 weeks later treatment with dexamethasone (dex), the nonsteroidal anti-inflammatory drug carprofen, or vehicle was started and continued for 2.5 weeks. At the termination of experiments, femurs were phenotyped using peripheral quantitative computed tomography and high-resolution micro-computed tomography, and markers of bone turnover were analyzed in serum. T and B lymphocyte populations in bone marrow and spleen were analyzed by flow cytometry. Dex-treated C57BL/6 mice had increased trabecular bone mineral density, but lower cortical content and thickness compared with vehicle-treated mice. The dex-treated mice also had lower levels of bone turnover markers and markedly decreased numbers of spleen T and B lymphocytes. In contrast, these effects could not be repeated when mice were treated with the nonsteroidal anti-inflammatory drug carprofen. In addition, dex did not increase trabecular bone in ovariectomized SCID mice lacking functional T and B lymphocytes. In contrast to most literature, the results from this study indicate that treatment with dex increased trabecular bone density, which may indicate that this effect is associated with corticosteroid-induced alterations of the lymphocyte populations.

Several infrequent genetic polymorphisms in the SERPINA1 gene are known to substantially reduce concentration of alpha1-antitrypsin (AAT) in the blood. Since low AAT serum levels fail to protect pulmonary tissue from enzymatic degradation, these polymorphisms also increase the risk for early onset chronic obstructive pulmonary disease (COPD). The role of more common SERPINA1 single nucleotide polymorphisms (SNPs) in respiratory health remains poorly understood. We present here an agnostic investigation of genetic determinants of circulating AAT levels in a general population sample by performing a genome-wide association study (GWAS) in 1392 individuals of the SAPALDIA cohort. Five common SNPs, defined by showing minor allele frequencies (MAFs) >5%, reached genome-wide significance, all located in the SERPINA gene cluster at 14q32.13. The top-ranking genotyped SNP rs4905179 was associated with an estimated effect of β = -0.068 g/L per minor allele (P = 1.20*10(-12)). But denser SERPINA1 locus genotyping in 5569 participants with subsequent stepwise conditional analysis, as well as exon-sequencing in a subsample (N = 410), suggested that AAT serum level is causally determined at this locus by rare (MAF<1%) and low-frequent (MAF 1-5%) variants only, in particular by the well-documented protein inhibitor S and Z (PI S, PI Z) variants. Replication of the association of rs4905179 with AAT serum levels in the Copenhagen City Heart Study (N = 8273) was successful (P<0.0001), as was the replication of its synthetic nature (the effect disappeared after adjusting for PI S and Z, P = 0.57). Extending the analysis to lung function revealed a more complex situation. Only in individuals with severely compromised pulmonary health (N = 397), associations of common SNPs at this locus with lung function were driven by rarer PI S or Z variants. Overall, our meta-analysis of lung function in ever-smokers does not support a functional role of common SNPs in the SERPINA gene cluster in the general population.

Polycystic ovary syndrome (PCOS) is a reproductive and metabolic disorder associated with obesity and insulin resistance that often precedes the development of type-2 diabetes. Rats continuously exposed to dihydrotestosterone from prepuberty display typical reproductive and metabolic PCOS characteristics including anovulation, polycystic ovaries, insulin resistance, and obesity. Our aim was to investigate if resveratrol improves reproductive and metabolic functions in PCOS rats. The effect was compared to exercise. Control and PCOS rats were treated with vehicle or resveratrol (400 mg · kg(-1) · day(-1)) for 5-6 weeks. Another group of PCOS rats received vehicle treatment and exercised for 5-6 weeks. Insulin sensitivity was determined by euglycemic-hyperinsulinemic clamp. The glucose infusion rate was lower in the PCOS-vehicle group compared to control-vehicle rats (P < 0.05). Exercise increased insulin sensitivity compared with PCOS-vehicle rats (P < 0.05), but resveratrol did not. Resveratrol treatment and exercise resulted in smaller adipocytes, upregulated estrogen-related receptor α gene expression in subcutaneous fat, and improved estrus cyclicity in the previously acyclic PCOS rats. Although resveratrol had positive effects on adiposity and cyclicity in a similar manner to exercise, resveratrol does not seem to be a good candidate for treating insulin resistance associated with PCOS because no improvement in insulin sensitivity was observed in PCOS rats on normal chow.

Genome-wide association studies (GWAS) have provided valuable insights into the genetic basis of complex traits. However, they have explained relatively little trait heritability. Recently, we proposed a new analytical approach called regional heritability mapping (RHM) that captures more of the missing genetic variation. This method is applicable both to related and unrelated populations. Here, we demonstrate the power of RHM in comparison with single-SNP GWAS and gene-based association approaches under a wide range of scenarios with variable numbers of quantitative trait loci (QTL) with common and rare causal variants in a narrow genomic region. Simulations based on real genotype data were performed to assess power to capture QTL variance, and we demonstrate that RHM has greater power to detect rare variants and/or multiple alleles in a region than other approaches. In addition, we show that RHM can capture more accurately the QTL variance, when it is caused by multiple independent effects and/or rare variants. We applied RHM to analyze three biometrical eye traits for which single-SNP GWAS have been published or performed to evaluate the effectiveness of this method in real data analysis and detected some additional loci which were not detected by other GWAS methods. RHM has the potential to explain some of missing heritability by capturing variance caused by QTL with low MAF and multiple independent QTL in a region, not captured by other GWAS methods. RHM analyses can be implemented using the software REACTA (http://www.epcc.ed.ac.uk/projects-portfolio/reacta).

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged > 55 years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey–Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey–Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han–Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p < 5 × 10− 8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p = 4.6 × 10− 8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98–1.14; p = 0.17), displaying high degree of heterogeneity (I2 = 57%; Qhet p = 0.0006). Under Han–Eskin alternative random effects model the summary effect was significant (p = 0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size > 1.25) may still be consistent with an effect size < 1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Screening gene function in vivo is a powerful approach to discover novel drug targets. We present high-throughput screening (HTS) data for 3 762 distinct global gene knockout (KO) mouse lines with viable adult homozygous mice generated using either gene-trap or homologous recombination technologies. Bone mass was determined from DEXA scans of male and female mice at 14 weeks of age and by microCT analyses of bones from male mice at 16 weeks of age. Wild-type (WT) cagemates/littermates were examined for each gene KO. Lethality was observed in an additional 850 KO lines. Since primary HTS are susceptible to false positive findings, additional cohorts of mice from KO lines with intriguing HTS bone data were examined. Aging, ovariectomy, histomorphometry and bone strength studies were performed and possible non-skeletal phenotypes were explored. Together, these screens identified multiple genes affecting bone mass: 23 previously reported genes (Calcr, Cebpb, Crtap, Dcstamp, Dkk1, Duoxa2, Enpp1, Fgf23, Kiss1/Kiss1r, Kl (Klotho), Lrp5, Mstn, Neo1, Npr2, Ostm1, Postn, Sfrp4, Slc30a5, Slc39a13, Sost, Sumf1, Src, Wnt10b), five novel genes extensively characterized (Cldn18, Fam20c, Lrrk1, Sgpl1, Wnt16), five novel genes with preliminary characterization (Agpat2, Rassf5, Slc10a7, Slc26a7, Slc30a10) and three novel undisclosed genes coding for potential osteoporosis drug targets.

Extraversion is a relatively stable and heritable personality trait associated with numerous psychosocial, lifestyle and health outcomes. Despite its substantial heritability, no genetic variants have been detected in previous genome-wide association (GWA) studies, which may be due to relatively small sample sizes of those studies. Here, we report on a large meta-analysis of GWA studies for extraversion in 63,030 subjects in 29 cohorts. Extraversion item data from multiple personality inventories were harmonized across inventories and cohorts. No genome-wide significant associations were found at the single nucleotide polymorphism (SNP) level but there was one significant hit at the gene level for a long non-coding RNA site (LOC101928162). Genome-wide complex trait analysis in two large cohorts showed that the additive variance explained by common SNPs was not significantly different from zero, but polygenic risk scores, weighted using linkage information, significantly predicted extraversion scores in an independent cohort. These results show that extraversion is a highly polygenic personality trait, with an architecture possibly different from other complex human traits, including other personality traits. Future studies are required to further determine which genetic variants, by what modes of gene action, constitute the heritable nature of extraversion.

The bone-sparing effect of estrogens is mediated primarily via estrogen receptor (ER)α, which stimulates gene transcription through activation function (AF)-1 and AF-2. The role of ERαAF-1 for the estradiol (E2) effects is tissue specific. The selective ER modulators (SERMs) raloxifene (Ral), lasofoxifene (Las), and bazedoxifene (Bza) can be used to treat postmenopausal osteoporosis. They all reduce the risk for vertebral fractures, whereas Las and partly Bza, but not Ral, reduce the risk for nonvertebral fractures. Here, we have compared the tissue specificity of Ral, Las, and Bza and evaluated the role of ERαAF-1 for the effects of these SERMs, with an emphasis on bone parameters. We treated ovariectomized (OVX) wild-type (WT) mice and OVX mice lacking ERαAF-1 (ERαAF-1(0)) with E2, Ral, Las, or Bza. All three SERMs increased trabecular bone mass in the axial skeleton. In the appendicular skeleton, only Las increased the trabecular bone volume/tissue volume and trabecular number, whereas both Ral and Las increased the cortical bone thickness and strength. However, Ral also increased cortical porosity. The three SERMs had only a minor effect on uterine weight. Notably, all evaluated effects of these SERMs were absent in ovx ERαAF-1(0) mice. In conclusion, all SERMs had similar effects on axial bone mass. However, the SERMs had slightly different effects on the appendicular skeleton since only Las increased the trabecular bone mass and only Ral increased the cortical porosity. Importantly, all SERM effects require a functional ERαAF-1 in female mice. These results could lead to development of more specific treatments for osteoporosis.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.