Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

Estimating disease heritability from complex pedigrees allowing for ascertainment and covariates.

We propose TetraHer, a method for estimating the liability heritability of binary phenotypes. TetraHer has five key features. First, it can be applied to data from complex pedigrees that contain multiple types of relationships. Second, it can correct for ascertainment of cases or controls. Third, it can accommodate covariates. Fourth, it can model the contribution of common environment. Fifth, it produces a likelihood that can be used for significance testing. We first demonstrate the validity of TetraHer on simulated data. We then use TetraHer to estimate liability heritability for 229 codes from the tenth International Classification of Diseases (ICD-10). We identify 107 codes with significant heritability (p < 0.05/229), which can be used in future analyses for investigating the genetic architecture of human diseases.

Evaluating indirect genetic effects of siblings using singletons.

Estimating effects of parental and sibling genotypes (indirect genetic effects) can provide insight into how the family environment influences phenotypic variation. There is growing molecular genetic evidence for effects of parental phenotypes on their offspring (e.g. parental educational attainment), but the extent to which siblings affect each other is currently unclear. Here we used data from samples of unrelated individuals, without (singletons) and with biological full-siblings (non-singletons), to investigate and estimate sibling effects. Indirect genetic effects of siblings increase (or decrease) the covariance between genetic variation and a phenotype. It follows that differences in genetic association estimates between singletons and non-singletons could indicate indirect genetic effects of siblings if there is no heterogeneity in other sources of genetic association between singletons and non-singletons. We used UK Biobank data to estimate polygenic score (PGS) associations for height, BMI and educational attainment in self-reported singletons (N = 50,143) and non-singletons (N = 328,549). The educational attainment PGS association estimate was 12% larger (95% C.I. 3%, 21%) in the non-singleton sample than in the singleton sample, but the height and BMI PGS associations were consistent. Birth order data suggested that the difference in educational attainment PGS associations was driven by individuals with older siblings rather than firstborns. The relationship between number of siblings and educational attainment PGS associations was non-linear; PGS associations were 24% smaller in individuals with 6 or more siblings compared to the rest of the sample (95% C.I. 11%, 38%). We estimate that a 1 SD increase in sibling educational attainment PGS corresponds to a 0.025 year increase in the index individual's years in schooling (95% C.I. 0.013, 0.036). Our results suggest that older siblings may influence the educational attainment of younger siblings, adding to the growing evidence that effects of the environment on phenotypic variation partially reflect social effects of germline genetic variation in relatives.

Fetal alleles predisposing to metabolically favorable adiposity are associated with higher birth weight.

BACKGROUND: Higher birthweight is associated with higher adult body mass index (BMI). Alleles that predispose to greater adult adiposity might act in fetal life to increase fetal growth and birthweight. Whether there are fetal effects of recently identified adult metabolically favorable adiposity alleles on birthweight is unknown. AIM: We aimed to test the effect on birthweight of fetal genetic predisposition to higher metabolically favorable adult adiposity and compare that with the effect of fetal genetic predisposition to higher adult BMI. METHODS: We used published genome wide association study data (n = upto 406 063) to estimate fetal effects on birthweight (adjusting for maternal genotype) of alleles known to raise metabolically favorable adult adiposity or BMI. We combined summary data across single nucleotide polymorphisms (SNPs) with random effects meta-analyses. We performed weighted linear regression of SNP-birthweight effects against SNP-adult adiposity effects to test for a dose-dependent association. RESULTS: Fetal genetic predisposition to higher metabolically favorable adult adiposity and higher adult BMI were both associated with higher birthweight (3 g per effect allele (95% CI: 1-5) averaged over 14 SNPs; P = 0.002; 0.5 g per effect allele (95% CI: 0-1) averaged over 76 SNPs; P = 0.042, respectively). SNPs with greater effects on metabolically favorable adiposity tended to have greater effects on birthweight (R2 = 0.2912, P = 0.027). There was no dose-dependent association for BMI (R2 = -0.0019, P = 0.602). CONCLUSIONS: Fetal genetic predisposition to both higher adult metabolically favorable adiposity and BMI is associated with birthweight. Fetal effects of metabolically favorable adiposity-raising alleles on birthweight are modestly proportional to their effects on future adiposity, but those of BMI-raising alleles are not.

A Note on Modelling Bidirectional Feedback Loops in Mendelian Randomization Studies.

Structural equation models (SEMs) involving feedback loops may offer advantages over standard instrumental variables estimators in terms of modelling causal effects in the presence of bidirectional relationships. In the following note, we show that in the case of a single "exposure" and "outcome" variable, modelling relationships using a SEM with a simple bidirectional linear feedback loop offers no advantage over traditional instrumental variables estimators in terms of consistency (i.e. both approaches yield consistent estimates of the causal effect, provided that causal estimates are obtained in both directions). In the case of finite samples, traditional IV estimators and SEM exhibited similar power across many of the conditions we examined, although which method performed best depended on the residual correlation between variables and the strength of the instruments. In particular, the power of SEM was insensitive to the residual correlation between variables, whereas the power of the Wald estimator/2SLS improved (deteriorated) relative to SEM as the magnitude of the residual correlation increased (decreased) assuming a positive causal effect of the exposure on the outcome. The power of SEM improved relative to the Wald estimator/2SLS as the instruments explained more residual variance in the "outcome" variable.

Using Genomic Structural Equation Modeling to Partition the Genetic Covariance Between Birthweight and Cardiometabolic Risk Factors into Maternal and Offspring Components in the Norwegian HUNT Study.

The Barker Hypothesis posits that adverse intrauterine environments result in fetal growth restriction and increased risk of cardiometabolic disease through developmental compensations. Here we introduce a new statistical model using the genomic SEM software that is capable of simultaneously partitioning the genetic covariation between birthweight and cardiometabolic traits into maternally mediated and offspring mediated contributions. We model the covariance between birthweight and later life outcomes, such as blood pressure, non-fasting glucose, blood lipids and body mass index in the Norwegian HUNT study, consisting of 15,261 mother-eldest offspring pairs with genetic and phenotypic data. Application of this model showed some evidence for maternally mediated effects of systolic blood pressure on offspring birthweight, and pleiotropy between birthweight and non-fasting glucose mediated through the offspring genome. This underscores the importance of genetic links between birthweight and cardiometabolic phenotypes and offer alternative explanations to environmentally based hypotheses for the phenotypic correlation between these variables.

Estimating indirect parental genetic effects on offspring phenotypes using virtual parental genotypes derived from sibling and half sibling pairs.

Indirect parental genetic effects may be defined as the influence of parental genotypes on offspring phenotypes over and above that which results from the transmission of genes from parents to their children. However, given the relative paucity of large-scale family-based cohorts around the world, it is difficult to demonstrate parental genetic effects on human traits, particularly at individual loci. In this manuscript, we illustrate how parental genetic effects on offspring phenotypes, including late onset conditions, can be estimated at individual loci in principle using large-scale genome-wide association study (GWAS) data, even in the absence of parental genotypes. Our strategy involves creating "virtual" mothers and fathers by estimating the genotypic dosages of parental genotypes using physically genotyped data from relative pairs. We then utilize the expected dosages of the parents, and the actual genotypes of the offspring relative pairs, to perform conditional genetic association analyses to obtain asymptotically unbiased estimates of maternal, paternal and offspring genetic effects. We apply our approach to 19066 sibling pairs from the UK Biobank and show that a polygenic score consisting of imputed parental educational attainment SNP dosages is strongly related to offspring educational attainment even after correcting for offspring genotype at the same loci. We develop a freely available web application that quantifies the power of our approach using closed form asymptotic solutions. We implement our methods in a user-friendly software package IMPISH (IMputing Parental genotypes In Siblings and Half Siblings) which allows users to quickly and efficiently impute parental genotypes across the genome in large genome-wide datasets, and then use these estimated dosages in downstream linear mixed model association analyses. We conclude that imputing parental genotypes from relative pairs may provide a useful adjunct to existing large-scale genetic studies of parents and their offspring.

Mendelian randomization analysis of factors related to ovulation and reproductive function and endometrial cancer risk.

BACKGROUND: Observational epidemiological studies suggest a link between several factors related to ovulation and reproductive function and endometrial cancer (EC) risk; however, it is not clear whether these relationships are causal, and whether the risk factors act independently of each other. The aim of this study was to investigate putative causal relationships between the number of live births, age at last live birth, and years ovulating and EC risk.  METHODS: We conducted a series of observational analyses to investigate various risk factors and EC risk in the UK Biobank (UKBB). Additionally, multivariate analysis was performed to elucidate the relationship between the number of live births, age at last live birth, and years ovulating and other related factors such as age at natural menopause, age at menarche, and body mass index (BMI). Secondly, we used Mendelian randomization (MR) to assess if these observed relationships were causal. Genome-wide significant single nucleotide polymorphisms (SNPs) were extracted from previous studies of woman's number of live births, age at menopause and menarche, and BMI. We conducted a genome-wide association analysis using the UKBB to identify SNPs associated with years ovulating, years using the contraceptive pill, and age at last live birth. RESULTS: We found evidence for a causal effect of the number of live births (inverse variance weighted (IVW) odds ratio (OR): 0.537, p = 0.006), the number of years ovulating (IVW OR: 1.051, p = 0.014), in addition to the known risk factors BMI, age at menarche, and age at menopause on EC risk in the univariate MR analyses. Due to the close relationships between these factors, we followed up with multivariable MR (MVMR) analysis. Results from the MVMR analysis showed that number of live births had a causal effect on EC risk (OR: 0.783, p = 0.036) independent of BMI, age at menarche and age at menopause. CONCLUSIONS: MVMR analysis showed that the number of live births causally reduced the risk of EC.

Exploring the causal effect of maternal pregnancy adiposity on offspring adiposity: Mendelian randomisation using polygenic risk scores.

BACKGROUND: Greater maternal adiposity before or during pregnancy is associated with greater offspring adiposity throughout childhood, but the extent to which this is due to causal intrauterine or periconceptional mechanisms remains unclear. Here, we use Mendelian randomisation (MR) with polygenic risk scores (PRS) to investigate whether associations between maternal pre-/early pregnancy body mass index (BMI) and offspring adiposity from birth to adolescence are causal. METHODS: We undertook confounder adjusted multivariable (MV) regression and MR using mother-offspring pairs from two UK cohorts: Avon Longitudinal Study of Parents and Children (ALSPAC) and Born in Bradford (BiB). In ALSPAC and BiB, the outcomes were birthweight (BW; N = 9339) and BMI at age 1 and 4 years (N = 8659 to 7575). In ALSPAC only we investigated BMI at 10 and 15 years (N = 4476 to 4112) and dual-energy X-ray absorptiometry (DXA) determined fat mass index (FMI) from age 10-18 years (N = 2659 to 3855). We compared MR results from several PRS, calculated from maternal non-transmitted alleles at between 29 and 80,939 single nucleotide polymorphisms (SNPs). RESULTS: MV and MR consistently showed a positive association between maternal BMI and BW, supporting a moderate causal effect. For adiposity at most older ages, although MV estimates indicated a strong positive association, MR estimates did not support a causal effect. For the PRS with few SNPs, MR estimates were statistically consistent with the null, but had wide confidence intervals so were often also statistically consistent with the MV estimates. In contrast, the largest PRS yielded MR estimates with narrower confidence intervals, providing strong evidence that the true causal effect on adolescent adiposity is smaller than the MV estimates (Pdifference = 0.001 for 15-year BMI). This suggests that the MV estimates are affected by residual confounding, therefore do not provide an accurate indication of the causal effect size. CONCLUSIONS: Our results suggest that higher maternal pre-/early-pregnancy BMI is not a key driver of higher adiposity in the next generation. Thus, they support interventions that target the whole population for reducing overweight and obesity, rather than a specific focus on women of reproductive age.

Higher maternal adiposity reduces offspring birthweight if associated with a metabolically favourable profile.

AIMS/HYPOTHESIS: Higher maternal BMI during pregnancy is associated with higher offspring birthweight, but it is not known whether this is solely the result of adverse metabolic consequences of higher maternal adiposity, such as maternal insulin resistance and fetal exposure to higher glucose levels, or whether there is any effect of raised adiposity through non-metabolic (e.g. mechanical) factors. We aimed to use genetic variants known to predispose to higher adiposity, coupled with a favourable metabolic profile, in a Mendelian randomisation (MR) study comparing the effect of maternal 'metabolically favourable adiposity' on offspring birthweight with the effect of maternal general adiposity (as indexed by BMI). METHODS: To test the causal effects of maternal metabolically favourable adiposity or general adiposity on offspring birthweight, we performed two-sample MR. We used variants identified in large, published genetic-association studies as being associated with either higher adiposity and a favourable metabolic profile, or higher BMI (n = 442,278 and n = 322,154 for metabolically favourable adiposity and BMI, respectively). We then extracted data on the metabolically favourable adiposity and BMI variants from a large, published genetic-association study of maternal genotype and offspring birthweight controlling for fetal genetic effects (n = 406,063 with maternal and/or fetal genotype effect estimates). We used several sensitivity analyses to test the reliability of the results. As secondary analyses, we used data from four cohorts (total n = 9323 mother-child pairs) to test the effects of maternal metabolically favourable adiposity or BMI on maternal gestational glucose, anthropometric components of birthweight and cord-blood biomarkers. RESULTS: Higher maternal adiposity with a favourable metabolic profile was associated with lower offspring birthweight (-94 [95% CI -150, -38] g per 1 SD [6.5%] higher maternal metabolically favourable adiposity, p = 0.001). By contrast, higher maternal BMI was associated with higher offspring birthweight (35 [95% CI 16, 53] g per 1 SD [4 kg/m2] higher maternal BMI, p = 0.0002). Sensitivity analyses were broadly consistent with the main results. There was evidence of outlier SNPs for both exposures; their removal slightly strengthened the metabolically favourable adiposity estimate and made no difference to the BMI estimate. Our secondary analyses found evidence to suggest that a higher maternal metabolically favourable adiposity decreases pregnancy fasting glucose levels while a higher maternal BMI increases them. The effects on neonatal anthropometric traits were consistent with the overall effect on birthweight but the smaller sample sizes for these analyses meant that the effects were imprecisely estimated. We also found evidence to suggest that higher maternal metabolically favourable adiposity decreases cord-blood leptin while higher maternal BMI increases it. CONCLUSIONS/INTERPRETATION: Our results show that higher adiposity in mothers does not necessarily lead to higher offspring birthweight. Higher maternal adiposity can lead to lower offspring birthweight if accompanied by a favourable metabolic profile. DATA AVAILABILITY: The data for the genome-wide association studies (GWAS) of BMI are available at https://portals.broadinstitute.org/collaboration/giant/index.php/GIANT_consortium_data_files . The data for the GWAS of body fat percentage are available at https://walker05.u.hpc.mssm.edu .

Within family Mendelian randomization studies.

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

Genome-Wide Association Study Identifies Two Novel Loci Associated with Female Stress and Urgency Urinary Incontinence.

PURPOSE: Genome-wide association studies have not identified replicable genetic risk loci for stress or urgency urinary incontinence. MATERIALS AND METHODS: We carried out a discovery stage, case control, genome-wide association study in 3 independent discovery cohorts of European women (8,979) for stress incontinence, urgency incontinence, and any incontinence phenotypes. We conducted replication in 6 additional studies of European ancestry (4,069). We collected bladder biopsies from women with incontinence (50) to further investigate bladder expression of implicated genes and pathways and used symptom questionnaires for phenotyping. We conducted meta-analyses using inverse variance fixed effects models and whole transcriptome analyses using Affymetrix® arrays with replication with TaqMan® polymerase chain reaction. RESULTS: In the discovery stage, we identified 16 single nucleotide polymorphisms genotyped or imputed at 5 loci that reached genome-wide significance (p <5×10-8). In replication, rs138724718 on chromosome 2 near the macrophage receptor with collagenous structure (MARCO) gene (replication p=0.003) was associated with stress incontinence. In addition, rs34998271 on chromosome 6 near the endothelin 1 (EDN1) gene (replication p=0.0008) was associated with urgency incontinence. In combined meta-analyses of discovery and replication cohorts, associations with genome-wide significance for these 2 single nucleotide polymorphisms were confirmed. Transcriptomics analyses showed differential expression of 7 of 19 genes in the endothelin pathway between stress and urgency incontinence (p <0.0001). CONCLUSIONS: We uncovered 2 new risk loci near the genes endothelin 1 (EDN1), associated with urgency incontinence, and macrophage receptor with collagenous structure (MARCO), associated with stress incontinence. These loci are biologically plausible given their roles in smooth muscle contraction and innate host defense, respectively.

The Boulder Workshop Question Box.

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

Modeling Parent-Specific Genetic Nurture in Families with Missing Parental Genotypes: Application to Birthweight and BMI.

Disaggregation and estimation of genetic effects from offspring and parents has long been of interest to statistical geneticists. Recently, technical and methodological advances have made the genome-wide and loci-specific estimation of direct offspring and parental genetic nurture effects more possible. However, unbiased estimation using these methods requires datasets where both parents and at least one child have been genotyped, which are relatively scarce. Our group has recently developed a method and accompanying software (IMPISH; Hwang et al. in PLoS Genet 16:e1009154, 2020) which is able to impute missing parental genotypes from observed data on sibships and estimate their effects on an offspring phenotype conditional on the effects of genetic transmission. However, this method is unable to disentangle maternal and paternal effects, which may differ in magnitude and direction. Here, we introduce an extension to the original IMPISH routine which takes advantage of all available nuclear families to impute parent-specific missing genotypes and obtain asymptotically unbiased estimates of genetic effects on offspring phenotypes. We apply this this method to data from related individuals in the UK Biobank, showing concordance with previous estimates of maternal genetic effects on offspring birthweight. We also conduct the first GWAS jointly estimating offspring-, maternal-, and paternal-specific genetic effects on body-mass index.

Direct and Indirect Effects of Maternal, Paternal, and Offspring Genotypes: Trio-GCTA.

Indirect genetic effects from relatives may result in misleading quantifications of heritability, but can also be of interest in their own right. In this paper we propose Trio-GCTA, a model for separating direct and indirect genetic effects when genome-wide single nucleotide polymorphism data have been collected from parent-offspring trios. The model is applicable to phenotypes obtained from any of the family members. We discuss appropriate parameter interpretations and apply the method to three exemplar phenotypes: offspring birth weight, maternal relationship satisfaction, and paternal body-mass index, using real data from the Norwegian Mother, Father and Child Cohort Study (MoBa).

The Augmented Classical Twin Design: Incorporating Genome-Wide Identity by Descent Sharing Into Twin Studies in Order to Model Violations of the Equal Environments Assumption.

The Classical Twin Method (CTM) compares the similarity of monozygotic (MZ) twins with that of dizygotic (DZ) twins to make inferences about the relative importance of genes and environment in the etiology of individual differences. The design has been applied to thousands of traits across the biomedical, behavioral and social sciences and is arguably the most widely used natural experiment known to science. The fundamental assumption of the CTM is that trait relevant environmental covariation within MZ pairs is the same as that found within DZ pairs, so that zygosity differences in within-pair variance must be due to genetic factors uncontaminated by the environment. This equal environments assumption (EEA) has been, and still is hotly contested, and has been mentioned as a possible contributing factor to the missing heritability conundrum. In this manuscript, we introduce a new model for testing the EEA, which we call the Augmented Classical Twin Design which uses identity by descent (IBD) sharing between DZ twin pairs to estimate separate environmental variance components for MZ and DZ twin pairs, and provides a test of whether these are equal. We show through simulation that given large samples of DZ twin pairs, the model provides unbiased estimates of variance components and valid tests of the EEA under strong assumptions (e.g. no epistatic variance, IBD sharing in DZ twins estimated accurately etc.) which may not hold in reality. Sample sizes in excess of 50,000 DZ twin pairs with genome-wide genetic data are likely to be required in order to detect substantial violations of the EEA with moderate power. Consequently, we recommend that the Augmented Classical Twin Design only be applied to datasets with very large numbers of DZ twin pairs (> 50,000 DZ twin pairs), and given the strong assumptions relating to the absence of epistatic variance, appropriate caution be exercised regarding interpretation of the results.

Plasma Cortisol, Aldosterone, and Ascorbic Acid Concentrations in Patients with Septic Shock Do Not Predict Treatment Effect of Hydrocortisone on Mortality. A Nested Cohort Study.

Rationale: Whether biomarkers can identify subgroups of patients with septic shock with differential treatment responses to hydrocortisone is unknown.Objectives: To determine if there is heterogeneity in effect for hydrocortisone on mortality, shock resolution, and other clinical outcomes based on baseline cortisol, aldosterone, and ascorbic acid concentrations.Methods: From May 2014 to April 2017, we obtained serum samples from 529 patients with septic shock from 22 ICUs in Australia and New Zealand.Measurements and Main Results: There were no significant interactions between the association with 90-day mortality and treatment with either hydrocortisone or placebo for total cortisol (odds ratio [OR], 1.09; 95% confidence interval [CI], 1.02-1.16 vs. OR, 1.07; 95% CI, 1.00-1.13; P = 0.70), free cortisol (OR, 1.20; 95% CI, 1.04-1.38 vs. OR, 1.16; 95% CI, 1.02-1.32; P = 0.75), aldosterone (OR, 1.01; 95% CI, 0.97-1.05 vs. OR, 1.01; 95% CI, 0.98-1.04; P = 0.99), or ascorbic acid (OR, 1.11; 95% CI, 0.89-1.39 vs. OR, 1.05; 95% CI, 0.91-1.22; P = 0.70), respectively. Similar results were observed for the association with shock resolution. Elevated free cortisol was significantly associated with 90-day mortality (OR, 1.13; 95% CI, 1.00-1.27; P = 0.04), but total cortisol, aldosterone, and ascorbic acid were not.Conclusions: In patients with septic shock, there was no heterogeneity in effect of adjunctive hydrocortisone on mortality, shock resolution, or other clinical outcomes based on cortisol, aldosterone, and ascorbic acid concentrations. Plasma aldosterone and ascorbic acid concentrations are not associated with outcome.

Development of a polygenic risk score to improve screening for fracture risk: A genetic risk prediction study.

BACKGROUND: Since screening programs identify only a small proportion of the population as eligible for an intervention, genomic prediction of heritable risk factors could decrease the number needing to be screened by removing individuals at low genetic risk. We therefore tested whether a polygenic risk score for heel quantitative ultrasound speed of sound (SOS)-a heritable risk factor for osteoporotic fracture-can identify low-risk individuals who can safely be excluded from a fracture risk screening program. METHODS AND FINDINGS: A polygenic risk score for SOS was trained and selected in 2 separate subsets of UK Biobank (comprising 341,449 and 5,335 individuals). The top-performing prediction model was termed "gSOS", and its utility in fracture risk screening was tested in 5 validation cohorts using the National Osteoporosis Guideline Group clinical guidelines (N = 10,522 eligible participants). All individuals were genome-wide genotyped and had measured fracture risk factors. Across the 5 cohorts, the average age ranged from 57 to 75 years, and 54% of studied individuals were women. The main outcomes were the sensitivity and specificity to correctly identify individuals requiring treatment with and without genetic prescreening. The reference standard was a bone mineral density (BMD)-based Fracture Risk Assessment Tool (FRAX) score. The secondary outcomes were the proportions of the screened population requiring clinical-risk-factor-based FRAX (CRF-FRAX) screening and BMD-based FRAX (BMD-FRAX) screening. gSOS was strongly correlated with measured SOS (r2 = 23.2%, 95% CI 22.7% to 23.7%). Without genetic prescreening, guideline recommendations achieved a sensitivity and specificity for correct treatment assignment of 99.6% and 97.1%, respectively, in the validation cohorts. However, 81% of the population required CRF-FRAX tests, and 37% required BMD-FRAX tests to achieve this accuracy. Using gSOS in prescreening and limiting further assessment to those with a low gSOS resulted in small changes to the sensitivity and specificity (93.4% and 98.5%, respectively), but the proportions of individuals requiring CRF-FRAX tests and BMD-FRAX tests were reduced by 37% and 41%, respectively. Study limitations include a reliance on cohorts of predominantly European ethnicity and use of a proxy of fracture risk. CONCLUSIONS: Our results suggest that the use of a polygenic risk score in fracture risk screening could decrease the number of individuals requiring screening tests, including BMD measurement, while maintaining a high sensitivity and specificity to identify individuals who should be recommended an intervention.

Genome-wide survey of parent-of-origin effects on DNA methylation identifies candidate imprinted loci in humans.

Genomic imprinting is an epigenetic mechanism leading to parent-of-origin silencing of alleles. So far, the precise number of imprinted regions in humans is uncertain. In this study, we leveraged genome-wide DNA methylation in whole blood measured longitudinally at three time points (birth, childhood and adolescence) and genome-wide association studies (GWAS) data in 740 mother-child duos from the Avon Longitudinal Study of parents and children to identify candidate imprinted loci. We reasoned that cis-meQTLs at genomic regions that were imprinted would show strong evidence of parent-of-origin associations with DNA methylation, enabling the detection of imprinted regions. Using this approach, we identified genome-wide significant cis-meQTLs that exhibited parent-of-origin effects (POEs) at 82 loci, 34 novel and 48 regions previously implicated in imprinting (3.7-10

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero.

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (ß = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.