Epigenetic and genetic components of height regulation

Mapping Epigenetic Gene Variant Dynamics: Comparative Analysis of Frequency, Functional Impact and Trait Associations in African and European Populations

Epigenetic modifications influence gene expression levels, impact organismal traits, and play a role in the development of diseases. Therefore, variants in genes involved in epigenetic processes are likely to be important in disease susceptibility, and the frequency of variants may vary between populations with African and European ancestries. Here, we analyse an integrated dataset to define the frequencies, associated traits, and functional impact of epigenetic gene variants among individuals of African and European ancestry represented in the UK Biobank. We find that the frequencies of 88.4% of epigenetic gene variants significantly differ between these groups. Furthermore, we find that the variants are associated with many traits and diseases, and some of these associations may be population-specific owing to allele frequency differences. Additionally, we observe that variants associated with traits are significantly enriched for quantitative trait loci that affect DNA methylation, chromatin accessibility, and gene expression. We find that methylation quantitative trait loci account for 71.2% of the variants influencing gene expression. Moreover, variants linked to biomarker traits exhibit high correlation. We therefore conclude that epigenetic gene variants associated with traits tend to differ in their allele frequencies among African and European populations and are enriched for QTLs.

Characterising the genetic architecture of changes in adiposity during adulthood using electronic health records

Obesity is a heritable disease, characterised by excess adiposity that is measured by body mass index (BMI). While over 1,000 genetic loci are associated with BMI, less is known about the genetic contribution to adiposity trajectories over adulthood. We derive adiposity-change phenotypes from 24.5 million primary-care health records in over 740,000 individuals in the UK Biobank, Million Veteran Program USA, and Estonian Biobank, to discover and validate the genetic architecture of adiposity trajectories. Using multiple BMI measurements over time increases power to identify genetic factors affecting baseline BMI by 14%. In the largest reported genome-wide study of adiposity-change in adulthood, we identify novel associations with BMI-change at six independent loci, including rs429358 (APOE missense variant). The SNP-based heritability of BMI-change (1.98%) is 9-fold lower than that of BMI. The modest genetic correlation between BMI-change and BMI (45.2%) indicates that genetic studies of longitudinal trajectories could uncover novel biology of quantitative traits in adulthood.

Genome-wide analyses identify 21 infertility loci and over 400 reproductive hormone loci across the allele frequency spectrum

Genome-wide association studies (GWASs) may help inform treatments for infertility, whose causes remain unknown in many cases. Here we present GWAS meta-analyses across six cohorts for male and female infertility in up to 41,200 cases and 687,005 controls. We identified 21 genetic risk loci for infertility (P≤5E-08), of which 12 have not been reported for any reproductive condition. We found positive genetic correlations between endometriosis and all-cause female infertility (rg=0.585, P=8.98E-14), and between polycystic ovary syndrome and anovulatory infertility (rg=0.403, P=2.16E-03). The evolutionary persistence of female infertility-risk alleles in EBAG9 may be explained by recent directional selection. We additionally identified up to 269 genetic loci associated with follicle-stimulating hormone (FSH), luteinising hormone, oestradiol, and testosterone through sex-specific GWAS meta-analyses (N=6,095-246,862). While hormone-associated variants near FSHB and ARL14EP colocalised with signals for anovulatory infertility, we found no rg between female infertility and reproductive hormones (P>0.05). Exome sequencing analyses in the UK Biobank (N=197,340) revealed that women carrying testosterone-lowering rare variants in GPC2 were at higher risk of infertility (OR=2.63, P=1.25E-03). Taken together, our results suggest that while individual genes associated with hormone regulation may be relevant for fertility, there is limited genetic evidence for correlation between reproductive hormones and infertility at the population level. We provide the first comprehensive view of the genetic architecture of infertility across multiple diagnostic criteria in men and women, and characterise its relationship to other health conditions.

A population-based survey of FBN1 variants in Iceland reveals underdiagnosis of Marfan syndrome

Marfan syndrome (MFS) is an autosomal dominant condition characterized by aortic aneurysm, skeletal abnormalities, and lens dislocation, and is caused by variants in the FBN1 gene. To explore causes of MFS and the prevalence of the disease in Iceland we collected information from all living individuals with a clinical diagnosis of MFS in Iceland (n = 32) and performed whole-genome sequencing of those who did not have a confirmed genetic diagnosis (27/32). Moreover, to assess a potential underdiagnosis of MFS in Iceland we attempted a genotype-based approach to identify individuals with MFS. We interrogated deCODE genetics' database of 35,712 whole-genome sequenced individuals to search for rare sequence variants in FBN1. Overall, we identified 15 pathogenic or likely pathogenic variants in FBN1 in 44 individuals, only 22 of whom were previously diagnosed with MFS. The most common of these variants, NM_000138.4:c.8038 C > T p.(Arg2680Cys), is present in a multi-generational pedigree, and was found to stem from a single forefather born around 1840. The p.(Arg2680Cys) variant associates with a form of MFS that seems to have an enrichment of abdominal aortic aneurysm, suggesting that this may be a particularly common feature of p.(Arg2680Cys)-associated MFS. Based on these combined genetic and clinical data, we show that MFS prevalence in Iceland could be as high as 1/6,600 in Iceland, compared to 1/10,000 based on clinical diagnosis alone, which indicates underdiagnosis of this actionable genetic disorder.

The complex role of MEG3: An emerging long non-coding RNA in breast cancer

MEG3, a significant long non-coding RNA (lncRNA), substantially functions in diverse biological processes, particularly breast cancer (BC) development. Within the imprinting DLK-MEG3 region on human chromosomal region 14q32.3, MEG3 spans 35 kb and encompasses ten exons. It exerts regulatory effects through intricate interactions with miRNAs, proteins, and epigenetic modifications. MEG3's multifaceted function in BC is evident in gene expression modulation, osteogenic tissue differentiation, and involvement in bone-related conditions. Its role as a tumor suppressor is highlighted by its influence on miR-182 and miRNA-29 expression in BC. Additionally, MEG3 is implicated in acute myocardial infarction and endothelial cell function, emphasising cell-specific regulatory mechanisms. MEG3's impact on gene activity encompasses transcriptional and post-translational adjustments, including DNA methylation, histone modifications, and interactions with transcription factors. MEG3 dysregulation is linked to unfavourable outcomes and drug resistance. Notably, higher MEG3 expression is associated with enhanced survival in BC patients. Overcoming challenges such as unravelling context-specific interactions, understanding epigenetic control, and translating findings into clinical applications is imperative. Prospective endeavours involve elucidating underlying mechanisms, exploring epigenetic alterations, and advancing MEG3-based diagnostic and therapeutic approaches. A comprehensive investigation into broader signaling networks and rigorous clinical trials are pivotal. Rigorous validation through functional and molecular analyses will shed light on MEG3's intricate contribution to BC progression.

Height-Related Polygenic Variants Are Associated with Metabolic Syndrome Risk and Interact with Energy Intake and a Rice-Main Diet to Influence Height in KoGES

Adult height is inversely related to metabolic syndrome (MetS) risk, but its genetic impacts have not been revealed. The present study aimed to examine the hypothesis that adult height-related genetic variants interact with lifestyle to influence adult height and are associated with MetS risk in adults aged >40 in Korea during 2010–2014. Participants were divided into short stature (SS; control) and tall stature (TS; case) by the 85th percentile of adult height. The genetic variants linked to adult height were screened from a genome-wide association study in a city hospital-based cohort (n = 58,701) and confirmed in Ansan/Ansung plus rural cohorts (n = 13,783) among the Korean Genome and Epidemiology Study. Genetic variants that interacted with each other were identified using the generalized multifactor dimensionality reduction (GMDR) analysis. The interaction between the polygenic risk score (PRS) of the selected genetic variants and lifestyles was examined. Adult height was inversely associated with MetS, cardiovascular diseases, and liver function. The PRS, including zinc finger and BTB domain containing 38 (ZBTB38)_rs6762722, polyadenylate-binding protein-interacting protein-2B (PAIP2B)_rs13034890, carboxypeptidase Z (CPZ)_rs3756173, and latent-transforming growth factor beta-binding protein-1 (LTBP1)_rs4630744, was positively associated with height by 1.29 times and inversely with MetS by 0.894 times after adjusting for covariates. In expression quantitative trait loci, the gene expression of growth/differentiation factor-5 (GDF5)_rs224331, non-SMC condensin I complex subunit G (NCAPG)_rs2074974, ligand-dependent nuclear receptor corepressor like (LCORL)_rs7700107, and insulin-like growth factor-1 receptor (IGF1R)_rs2871865 was inversely linked to their risk allele in the tibial nerve and brain. The gene expression of PAIP2B_rs13034890 and a disintegrin and metalloproteinase with thrombospondin motifs-like-3 (ADAMTSL3)_rs13034890 was positively related to it. The PRS was inversely associated with MetS, hyperglycemia, HbA1c, and white blood cell counts. The wild type of GDF5_rs224331 (Ala276) lowered binding energy with rugosin A, D, and E (one of the hydrolyzable tannins) but not the mutated one (276Ser) in the in-silico analysis. The PRS interacted with energy intake and rice-main diet; PRS impact was higher in the high energy intake and the low rice-main diet. In conclusion, the PRS for adult height interacted with energy intake and diet patterns to modulate height and was linked to height and MetS by modulating their expression in the tibial nerve and brain.

Differences in wasting assessment through Middle-Upper Arm Circumference (MUAC) adjusted by sex, age and geographic origin for children aged 6-59 months: New reference based on anthropometric surveys from 22 low-and-middle-income countries

OBJECTIVES

The Middle Upper Arm Circumference (MUAC) bracelet is a widely used instrument in public health assessments and humanitarian assistance projects. The WHO guidelines present a universal cut-off point of 115 mm to determine whether a child has severe acute malnutrition. The objective of this study is to analyze the existing differences in the MUAC for boys and girls aged between 6 and 59 months, from 22 countries distributed in three different continents, in contrast to the use of this single cut-off point. In addition, the creation of MUAC growth charts is presented for reference use.

MATERIALS AND METHODS

This study was carried out with a database developed by Action Against Hunger, composed, after the data pre-processing phase, of 97 921 individuals without anthropometric failure from African, Asian, and American continents. MUAC measurements were compared between countries, dividing by sex and age groups. A k-means method was used to create country clusters to allow comparisons and the variability was resumed using a Principal Component Analysis. For each cluster, growth curves were created and smoothed using the LOESS method.

RESULTS

Our research has revealed the existence of differences in the MUAC between countries in both, males and females, although with different trends. The evidence was confirmed with the creation of two clusters using the k-means method, which, when graphically represented by the Principal Component Analysis, showed that the MUAC was clearly different. There were also differences between males and females within each cluster, where growth curves did not overlap in any age group.

CONCLUSIONS

All statistical analysis indicate that there are differences in the MUAC values for children without anthropometric failure between countries, but also between sexes. With this research, a new reference is proposed that consider the existing variability between human populations to improve the precision in the determination of severe acute malnutrition in children.

Genetic Influences on Fetal Alcohol Spectrum Disorder

Fetal alcohol spectrum disorder (FASD) encompasses the range of deleterious outcomes of prenatal alcohol exposure (PAE) in the affected offspring, including developmental delay, intellectual disability, attention deficits, and conduct disorders. Several factors contribute to the risk for and severity of FASD, including the timing, dose, and duration of PAE and maternal factors such as age and nutrition. Although poorly understood, genetic factors also contribute to the expression of FASD, with studies in both humans and animal models revealing genetic influences on susceptibility. In this article, we review the literature related to the genetics of FASD in humans, including twin studies, candidate gene studies in different populations, and genetic testing identifying copy number variants. Overall, these studies suggest different genetic factors, both in the mother and in the offspring, influence the phenotypic outcomes of PAE. While further work is needed, understanding how genetic factors influence FASD will provide insight into the mechanisms contributing to alcohol teratogenicity and FASD risk and ultimately may lead to means for early detection and intervention.

The genetic architecture of changes in adiposity during adulthood

Obesity is a heritable disease, characterised by excess adiposity that is measured by body mass index (BMI). While over 1,000 genetic loci are associated with BMI, less is known about the genetic contribution to adiposity trajectories over adulthood. We derive adiposity-change phenotypes from 1.5 million primary-care health records in over 177,000 individuals in UK Biobank to study the genetic architecture of weight-change. Using multiple BMI measurements over time increases power to identify genetic factors affecting baseline BMI. In the largest reported genome-wide study of adiposity-change in adulthood, we identify novel associations with BMI-change at six independent loci, including rs429358 (a missense variant in APOE). The SNP-based heritability of BMI-change (1.98%) is 9-fold lower than that of BMI, and higher in women than in men. The modest genetic correlation between BMI-change and BMI (45.2%) indicates that genetic studies of longitudinal trajectories could uncover novel biology driving quantitative trait values in adulthood.

Meta-analysis fine-mapping is often miscalibrated at single-variant resolution

Meta-analysis is pervasively used to combine multiple genome-wide association studies (GWASs). Fine-mapping of meta-analysis studies is typically performed as in a single-cohort study. Here, we first demonstrate that heterogeneity (e.g., of sample size, phenotyping, imputation) hurts calibration of meta-analysis fine-mapping. We propose a summary statistics-based quality-control (QC) method, suspicious loci analysis of meta-analysis summary statistics (SLALOM), that identifies suspicious loci for meta-analysis fine-mapping by detecting outliers in association statistics. We validate SLALOM in simulations and the GWAS Catalog. Applying SLALOM to 14 meta-analyses from the Global Biobank Meta-analysis Initiative (GBMI), we find that 67% of loci show suspicious patterns that call into question fine-mapping accuracy. These predicted suspicious loci are significantly depleted for having nonsynonymous variants as lead variant (2.7×; Fisher's exact p = 7.3 × 10-4). We find limited evidence of fine-mapping improvement in the GBMI meta-analyses compared with individual biobanks. We urge extreme caution when interpreting fine-mapping results from meta-analysis of heterogeneous cohorts.

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

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

Parent-of-Origin inference for biobanks

Identical genetic variations can have different phenotypic effects depending on their parent of origin. Yet, studies focusing on parent-of-origin effects have been limited in terms of sample size due to the lack of parental genomes or known genealogies. We propose a probabilistic approach to infer the parent-of-origin of individual alleles that does not require parental genomes nor prior knowledge of genealogy. Our model uses Identity-By-Descent sharing with second- and third-degree relatives to assign alleles to parental groups and leverages chromosome X data in males to distinguish maternal from paternal groups. We combine this with robust haplotype inference and haploid imputation to infer the parent-of-origin for 26,393 UK Biobank individuals. We screen 99 phenotypes for parent-of-origin effects and replicate the discoveries of 6 GWAS studies, confirming signals on body mass index, type 2 diabetes, standing height and multiple blood biomarkers, including the known maternal effect at the MEG3/DLK1 locus on platelet phenotypes. We also report a novel maternal effect at the TERT gene on telomere length, thereby providing new insights on the heritability of this phenotype. All our summary statistics are publicly available to help the community to better characterize the molecular mechanisms leading to parent-of-origin effects and their implications for human health.

Detection of Parent-of-Origin Effects for the Variants Associated With Behavioral Disinhibition in the MCTFR Data

Behavioral disinhibition is one of the important characteristics of many mental diseases. It has been reported in literature that serious behavioral disinhibition will affect people's health and greatly reduce people's quality of life. Meanwhile, behavioral disinhibition can easily lead to illegal drug abuse and violent crimes, etc., which will bring great harm to the society. At present, large-scale genome-wide association analysis has identified many loci associated with behavioral disinhibition. However, these studies have not incorporated the parent-of-origin effects (POE) into analysis, which may ignore or underestimate the genetic effects of loci on behavioral disinhibition. Therefore, in this article, we analyzed the five phenotypes related to behavioral disinhibition in the Minnesota Center for Twin and Family Research data (nicotine, alcohol consumption, alcohol dependence, illicit drugs, and non-substance use related behavioral disinhibition), to further explore the POE of variants on behavioral disinhibition. We applied a linear mixed model to test for the POE at a genome-wide scale on five transformed phenotypes, and found nine SNPs with statistically significant POE at the significance level of 5 × 10-8. Among them, SNPs rs4141854, rs9394515, and rs4711553 have been reported to be associated with two neurological disorders (restless legs syndrome and Tourette's syndrome) which are related to behavioral disinhibition; SNPs rs12960235 and rs715351 have been found to be associated with head and neck squamous cell carcinoma, skin cancer and type I diabetes, while both SNPs have not been identified to be related to behavioral disinhibition in literature; SNPs rs704833, rs6837925, rs1863548, and rs11067062 are novel loci identified in this article, and their function annotations have not been reported in literature. Follow-up study in molecular genetics is needed to verify whether they are surely related to behavioral disinhibition.

Novel missense ACAN gene variants linked to familial osteochondritis dissecans cluster in the C-terminal globular domain of aggrecan

The cartilage aggrecan proteoglycan is crucial for both skeletal growth and articular cartilage function. A number of aggrecan (ACAN) gene variants have been linked to skeletal disorders, ranging from short stature to severe chondrodyplasias. Osteochondritis dissecans is a disorder where articular cartilage and subchondral bone fragments come loose from the articular surface. We previously reported a missense ACAN variant linked to familial osteochondritis dissecans, with short stature and early onset osteoarthritis, and now describe three novel ACAN gene variants from additional families with this disorder. Like the previously described variant, these are autosomal dominant missense variants, resulting in single amino acid residue substitutions in the C-type lectin repeat of the aggrecan G3 domain. Functional studies showed that neither recombinant variant proteins, nor full-length variant aggrecan proteoglycan from heterozygous patient cartilage, were secreted to the same level as wild-type aggrecan. The variant proteins also showed decreased binding to known cartilage extracellular matrix ligands. Mapping these and other ACAN variants linked to hereditary skeletal disorders showed a clustering of osteochondritis dissecans-linked variants to the G3 domain. Taken together, this supports a link between missense ACAN variants affecting the aggrecan G3 domain and hereditary osteochondritis dissecans.

Genome-wide survey of parent-of-origin-specific associations across clinical traits derived from electronic health records

Parent-of-origin (PoO) effects refer to the differential phenotypic impacts of genetic variants dependent on their parental inheritance due to imprinting. While PoO effects can influence complex traits, they may be poorly captured by models that do not differentiate the parental origin of the variant. The aim of this study was to conduct a genome-wide screen for PoO effects on a broad range of clinical traits derived from electronic health records (EHR) in the DiscovEHR study enriched with familial relationships. Using pairwise kinship estimates from genetic data and demographic data, we identified 22,051 offspring among 134,049 individuals in the DiscovEHR study. PoO of ~9 million variants was assigned in the offspring by comparing offspring and parental genotypes and haplotypes. We then performed genome-wide PoO association analyses across 154 quantitative and 611 binary traits extracted from EHR. Of the 732 significant PoO associations identified (p < 5 × 10^-8), we attempted to replicate 274 PoO associations in the UK Biobank study with 5,015 offspring and replicated 9 PoO associations (p < 0.05). In summary, our study implements a bioinformatic and statistical approach to examine PoO effects genome-wide in a large population study enriched with familial relationships and systematically characterizes PoO effects on hundreds of clinical traits derived from EHR. Our results suggest that, while the statistical power to detect PoO effects remains modest yet, accurately modeling PoO effects has the potential to find new associations that may have been missed by the standard additive model, further enhancing the mechanistic understanding of genetic influence on complex traits.

Long-read sequencing of 3,622 Icelanders provides insight into the role of structural variants in human diseases and other traits

Long-read sequencing (LRS) promises to improve the characterization of structural variants (SVs). We generated LRS data from 3,622 Icelanders and identified a median of 22,636 SVs per individual (a median of 13,353 insertions and 9,474 deletions). We discovered a set of 133,886 reliably genotyped SV alleles and imputed them into 166,281 individuals to explore their effects on diseases and other traits. We discovered an association of a rare deletion in PCSK9 with lower low-density lipoprotein (LDL) cholesterol levels, compared to the population average. We also discovered an association of a multiallelic SV in ACAN with height; we found 11 alleles that differed in the number of a 57-bp-motif repeat and observed a linear relationship between the number of repeats carried and height. These results show that SVs can be accurately characterized at the population scale using LRS data in a genome-wide non-targeted approach and demonstrate how SVs impact phenotypes.

Impact of early life famine exposure on adulthood anthropometry among survivors of the 1983-1985 Ethiopian Great famine: a historical cohort study

BACKGROUND

Nutritional insult in early life brings adaptive changes in body structure and functioning that could remain throughout the affected individual's life course. The long term impact of early life famine exposure on adulthood anthropometric measurements has been recorded in previous studies. However, the results were contradictory. Hence, we extend this study to examine the impact of famine exposure during early life on adulthood's anthropometry among survivors of the 1983-85 Ethiopian great famine.

METHODS

A total of 1384 adult men and women survived from 1983 to 85 Ethiopian great famine were included in the study. Famine exposure status was classified into five groups: early life-exposed, prenatal-exposed, postnatal-exposed, adolescence-exposed, and non-exposed based on self-reported age and birthdate of the participants. Prenatal, post-natal, and adolescence exposed groups were considered as early life exposed. Following a standard procedure, anthropometric measurements were taken. A linear regression analysis was used to analyze the impact of famine exposure on adult anthropometric measurements adjusted for all possible covariates. The effect of famine exposure on overweight, general obesity, and abdominal obesity was examined using multinomial and binary logistic regression analysis.

RESULT

Compared to non-exposed groups, adult height was lower by 1.83 cm (β = - 1.83; 95% CI: - 3.05, - 0.58), 1.35 cm (β = - 1.35; 95% CI: - 2.56, - 0.14) and 2.07 cm (β = - 2.07 cm; 95% CI: - 3.31, - 0.80) among early life, prenatal and post-natal exposed groups, respectively. Likewise, famine exposure during early life (β = 0.02; 95% CI: 0.01, 0.03), prenatal (β = 0.03; 95% CI: 0.02, 0.03) and post-natal life (β = 0.02; 95% CI: 0.02, 0.03) was positively associated with increased waist to height ratio. However, none of the above exposures resulted in a significant association with body mass index (P > 0. 05). Additionally, exposure to famine during early stage of life was not associated with increased risk of overweight, general obesity and abdominal obesity in adults.

CONCLUSION

Decreased adult height and increased waist-to-height ratio were associated with early life exposure to famine, particularly prenatal and post-natal exposure. These results therefore underscore the significance of avoiding undernutrition in early life, which tends to be important for achieving once potential adult height and to minimize the increased risk of anthropometric markers of abdominal obesity such as waist to height ratio in later life.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease

Autoimmune thyroid disease is the most common autoimmune disease and is highly heritable¹. Here, by using a genome-wide association study of 30,234 cases and 725,172 controls from Iceland and the UK Biobank, we find 99 sequence variants at 93 loci, of which 84 variants are previously unreported^2-7. A low-frequency (1.36%) intronic variant in FLT3 (rs76428106-C) has the largest effect on risk of autoimmune thyroid disease (odds ratio (OR) = 1.46, P = 2.37 × 10^-24). rs76428106-C is also associated with systemic lupus erythematosus (OR = 1.90, P = 6.46 × 10^-4), rheumatoid factor and/or anti-CCP-positive rheumatoid arthritis (OR = 1.41, P = 4.31 × 10^-4) and coeliac disease (OR = 1.62, P = 1.20 × 10^-4). FLT3 encodes fms-related tyrosine kinase 3, a receptor that regulates haematopoietic progenitor and dendritic cells. RNA sequencing revealed that rs76428106-C generates a cryptic splice site, which introduces a stop codon in 30% of transcripts that are predicted to encode a truncated protein, which lacks its tyrosine kinase domains. Each copy of rs76428106-C doubles the plasma levels of the FTL3 ligand. Activating somatic mutations in FLT3 are associated with acute myeloid leukaemia⁸ with a poor prognosis and rs76428106-C also predisposes individuals to acute myeloid leukaemia (OR = 1.90, P = 5.40 × 10^-3). Thus, a predicted loss-of-function germline mutation in FLT3 causes a reduction in full-length FLT3, with a compensatory increase in the levels of its ligand and an increased disease risk, similar to that of a gain-of-function mutation.

The Challenge of Defining and Investigating the Causes of Idiopathic Short Stature and Finding an Effective Therapy

Idiopathic short stature (ISS) comprises a wide range of conditions associated with short stature that elude the conventional diagnostic work-up and are often caused by still largely unknown genetic variants. In the last decade, the improvement of diagnostic techniques has led to the discovery of causal mutations in genes involved in the function of the growth hormone (GH)/insulin-like growth factor-I (IGF-I) axis as well as in growth plate physiology. However, many cases of ISS remain idiopathic. In the future, the more frequent identification of the underlying causes will allow a better stratification of subjects and offer a tailored management. GH therapy has been proposed and approved in some countries for the treatment of children with ISS. To improve the efficacy of GH therapy, trials with GH combined with GnRH agonists, aromatase inhibitors, and even IGF-I have been conducted. This review aims to revise the current definition of ISS and discuss the management of children with ISS on the basis of the most recent evidence.

Searching for parent-of-origin effects on cardiometabolic traits in imprinted genomic regions

Cardiometabolic traits pose a major global public health burden. Large-scale genome-wide association studies (GWAS) have identified multiple loci accounting for up to 30% of the genetic variance in complex traits such as cardiometabolic traits. However, the contribution of parent-of-origin effects (POEs) to complex traits has been largely ignored in GWAS. Family-based studies enable the assessment of POEs in genetic association analyses. We investigated POEs on a range of complex traits in 3 family-based studies. The discovery phase was carried out in large pedigrees from the Kibbutzim Family Study (n = 901 individuals) and in 872 parent-offspring trios from the Jerusalem Perinatal Study. Focusing on imprinted genomic regions, we examined parent-specific associations with 12 complex traits (i.e., body-size, blood pressure, lipids), mostly cardiometabolic risk traits. Forty five of the 11,967 SNPs initially found to have POE were evaluated for replication (p value < 1 × 10-4) in Framingham Heart Study families (max n = 8000 individuals). Three common variants yielded evidence of POE in the meta-analysis. Two variants, located on chr6 in the HLA region, showed a paternal effect on height (rs1042136: βpaternal = -0.023, p value = 1.5 × 10-8 and rs1431403: βpaternal = -0.011, p value = 5.4 × 10-6). The corresponding maternally-derived effects were statistically nonsignificant. The variant rs9332053, located on chr13 in RCBTB2 gene, demonstrated a maternal effect on hip circumference (βmaternal = -4.24, p value = 9.6 × 10-6; βpaternal = 1.29, p value = 0.23). These findings provide evidence for the utility of incorporating POEs into association studies of cardiometabolic traits, especially anthropometric traits. The study highlights the benefits of using family-based data for deciphering the genetic architecture of complex traits.

Eighty-eight variants highlight the role of T cell regulation and airway remodeling in asthma pathogenesis

Asthma is one of the most common chronic diseases affecting both children and adults. We report a genome-wide association meta-analysis of 69,189 cases and 702,199 controls from Iceland and UK biobank. We find 88 asthma risk variants at 56 loci, 19 previously unreported, and evaluate their effect on other asthma and allergic phenotypes. Of special interest are two low frequency variants associated with protection against asthma; a missense variant in TNFRSF8 and 3‘ UTR variant in TGFBR1. Functional studies show that the TNFRSF8 variant reduces TNFRSF8 expression both on cell surface and in soluble form, acting as loss of function. eQTL analysis suggests that the TGFBR1 variant acts through gain of function and together with an intronic variant in a downstream gene, SMAD3, points to defective TGFβR1 signaling as one of the biological perturbations increasing asthma risk. Our results increase the number of asthma variants and implicate genes with known role in T cell regulation, inflammation and airway remodeling in asthma pathogenesis.

Asthma is a common allergic airway disease with significant inter-individual heterogeneity. Here, Olafsdottir et al. report a genome-wide meta-analysis of two large population-based cohorts to identify sequence variants that associate with asthma risk and perform follow-up functional analyses on a protective loss-of-function variant in TNFRSF8.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions

SUMMARY

We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci.

AVAILABILITY AND IMPLEMENTATION

ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Parental height modifies the association between linear growth and neurodevelopment in infancy

AIM

To estimate the extent to which maternal and paternal height modify the association between length-for-age Z-score (LAZ) and neurodevelopmental outcomes assessed by the Bayley Scales of Infant and Toddler Development, 3rd Edition (Bayley-III).

METHODS

Baseline data from a clinical trial in 600 Nepalese infants aged six to 11 months with LAZ less than -1 were utilised. The primary exposure was the LAZ score, interaction variables were maternal and paternal height, and the outcomes were Bayley-III cognitive, language and motor scaled scores. Linear regression and generalised additive model (GAM) were used to identify potential interactions.

RESULTS

Linear regression analysis stratified by parental height categories showed that association between unit increase in LAZ and cognitive scaled score differed across maternal (normal height: ß 1.16, 95% CI; 0.75, 1.57 and short height: ß 0.67, 95% CI; 0.28, 1.05) and paternal (normal height: ß 1.32, 95% CI; 0.91, 1.72 and short height: ß 0.61, 95% CI; 0.03, 1.18) height categories. Maternal height also modified the association between LAZ and fine motor scaled score.

CONCLUSION

The association between LAZ and neurodevelopmental outcomes was attenuated when maternal and paternal height was taken into account. Parental stature should be considered when using LAZ as a proxy for neurodevelopment among infants.

Complex Phenotypes: Mechanisms Underlying Variation in Human Stature

PURPOSE OF REVIEW

The goal of the review is to provide a comprehensive overview of the current understanding of the mechanisms underlying variation in human stature.

RECENT FINDINGS

Human height is an anthropometric trait that varies considerably within human populations as well as across the globe. Historically, much research focus was placed on understanding the biology of growth plate chondrocytes and how modifications to core chondrocyte proliferation and differentiation pathways potentially shaped height attainment in normal as well as pathological contexts. Recently, much progress has been made to improve our understanding regarding the mechanisms underlying the normal and pathological range of height variation within as well as between human populations, and today, it is understood to reflect complex interactions among a myriad of genetic, environmental, and evolutionary factors. Indeed, recent improvements in genetics (e.g., GWAS) and breakthroughs in functional genomics (e.g., whole exome sequencing, DNA methylation analysis, ATAC-sequencing, and CRISPR) have shed light on previously unknown pathways/mechanisms governing pathological and common height variation. Additionally, the use of an evolutionary perspective has also revealed important mechanisms that have shaped height variation across the planet. This review provides an overview of the current knowledge of the biological mechanisms underlying height variation by highlighting new research findings on skeletal growth control with an emphasis on previously unknown pathways/mechanisms influencing pathological and common height variation. In this context, this review also discusses how evolutionary forces likely shaped the genomic architecture of height across the globe.

Characterizing rare and low-frequency height-associated variants in the Japanese population

Human height is a representative phenotype to elucidate genetic architecture. However, the majority of large studies have been performed in European population. To investigate the rare and low-frequency variants associated with height, we construct a reference panel (N = 3,541) for genotype imputation by integrating the whole-genome sequence data from 1,037 Japanese with that of the 1000 Genomes Project, and perform a genome-wide association study in 191,787 Japanese. We report 573 height-associated variants, including 22 rare and 42 low-frequency variants. These 64 variants explain 1.7% of the phenotypic variance. Furthermore, a gene-based analysis identifies two genes with multiple height-increasing rare and low-frequency nonsynonymous variants (SLC27A3 and CYP26B1; P_SKAT-O < 2.5 × 10⁻⁶). Our analysis shows a general tendency of the effect sizes of rare variants towards increasing height, which is contrary to findings among Europeans, suggesting that height-associated rare variants are under different selection pressure in Japanese and European populations.

Thousands of genetic loci are known to associate with human height, but these are mainly based on studies in European ancestry populations. Here, Akiyama et al. construct a genotype reference panel for the Japanese population followed by GWAS and report 573 height associated variants in 191,787 Japanese.

Profiling the genome-wide landscape of tandem repeat expansions

Tandem repeat (TR) expansions have been implicated in dozens of genetic diseases, including Huntington's Disease, Fragile X Syndrome, and hereditary ataxias. Furthermore, TRs have recently been implicated in a range of complex traits, including gene expression and cancer risk. While the human genome harbors hundreds of thousands of TRs, analysis of TR expansions has been mainly limited to known pathogenic loci. A major challenge is that expanded repeats are beyond the read length of most next-generation sequencing (NGS) datasets and are not profiled by existing genome-wide tools. We present GangSTR, a novel algorithm for genome-wide genotyping of both short and expanded TRs. GangSTR extracts information from paired-end reads into a unified model to estimate maximum likelihood TR lengths. We validate GangSTR on real and simulated data and show that GangSTR outperforms alternative methods in both accuracy and speed. We apply GangSTR to a deeply sequenced trio to profile the landscape of TR expansions in a healthy family and validate novel expansions using orthogonal technologies. Our analysis reveals that healthy individuals harbor dozens of long TR alleles not captured by current genome-wide methods. GangSTR will likely enable discovery of novel disease-associated variants not currently accessible from NGS.

Idiopathic short stature and growth hormone sensitivity in prepubertal children

OBJECTIVE

We compared growth hormone sensitivity to an insulin-like growth factor I (IGF-I) generation test in children with idiopathic short stature (ISS) and of normal stature (NS) across the birthweight range.

METHODS

Forty-six prepubertal children (~7.1 years) born at term were studied: ISS (n = 23; 74% boys) and NS (n = 23; 57% boys). Children underwent a modified IGF-I generation test with recombinant human growth hormone (rhGH; 0.05 mg/kg/d) over four consecutive days. Hormonal concentrations were measured at baseline and day 5.

RESULTS

Children with idiopathic short stature were 1.90 SDS lighter (P < 0.0001) but had 4.5% more body fat (P = 0.0007) than NS children. Overall, decreasing birthweight SDS across the normal range (-1.9 to +1.5 SDS) was associated with lower percentage IGF-I response to rhGH stimulation in univariable (r = 0.45; P = 0.002) and multivariable models (β = 24.6; P = 0.006). Plasma IGF-I concentrations rose in both groups with rhGH stimulation (P < 0.0001). GHBP levels (P = 0.002) were suppressed in ISS children (-19%; P = 0.029) but increased among NS children (+18%; P = 0.028), with contrasting responses also observed for leptin and IGFBP-1. Further, the increase in insulin concentrations in response to rhGH stimulation was ~3-fold greater in NS children (142% vs 50%; P = 0.006).

CONCLUSIONS

A progressive decrease in birthweight SDS was associated with a reduction in GH sensitivity in both NS and ISS children. Thus, the lower IGF-I response to rhGH stimulation in association with decreasing birthweight indicates that the ISS children at the lower end of the birthweight spectrum may have partial GH resistance, which may contribute to their poorer growth.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10-42, β = -0.090) and confers risk of hip fracture (P = 1.0 × 10-8, OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

Sequence variation at ANAPC1 accounts for 24% of the variability in corneal endothelial cell density

The corneal endothelium is vital for transparency and proper hydration of the cornea. Here, we conduct a genome-wide association study of corneal endothelial cell density (cells/mm2), coefficient of cell size variation (CV), percentage of hexagonal cells (HEX) and central corneal thickness (CCT) in 6,125 Icelanders and find associations at 10 loci, including 7 novel. We assess the effects of these variants on various ocular biomechanics such as corneal hysteresis (CH), as well as eye diseases such as glaucoma and corneal dystrophies. Most notably, an intergenic variant close to ANAPC1 (rs78658973[A], frequency = 28.3%) strongly associates with decreased cell density and accounts for 24% of the population variance in cell density (β = -0.77 SD, P = 1.8 × 10-314) and associates with increased CH (β = 0.19 SD, P = 2.6 × 10-19) without affecting risk of corneal diseases and glaucoma. Our findings indicate that despite correlations between cell density and eye diseases, low cell density does not increase the risk of disease.

Characterizing mutagenic effects of recombination through a sequence-level genetic map

Genetic diversity arises from recombination and de novo mutation (DNM). Using a combination of microarray genotype and whole-genome sequence data on parent-child pairs, we identified 4,531,535 crossover recombinations and 200,435 DNMs. The resulting genetic map has a resolution of 682 base pairs. Crossovers exhibit a mutagenic effect, with overrepresentation of DNMs within 1 kilobase of crossovers in males and females. In females, a higher mutation rate is observed up to 40 kilobases from crossovers, particularly for complex crossovers, which increase with maternal age. We identified 35 loci associated with the recombination rate or the location of crossovers, demonstrating extensive genetic control of meiotic recombination, and our results highlight genes linked to the formation of the synaptonemal complex as determinants of crossovers.

Parent-of-origin effects on quantitative phenotypes in a large Hutterite pedigree

The impact of the parental origin of associated alleles in GWAS has been largely ignored. Yet sequence variants could affect traits differently depending on whether they are inherited from the mother or the father, as in imprinted regions, where identical inherited DNA sequences can have different effects based on the parental origin. To explore parent-of-origin effects (POEs), we studied 21 quantitative phenotypes in a large Hutterite pedigree to identify variants with single parent (maternal-only or paternal-only) effects, and then variants with opposite parental effects. Here we show that POEs, which can be opposite in direction, are relatively common in humans, have potentially important clinical effects, and will be missed in traditional GWAS. We identified POEs with 11 phenotypes, most of which are risk factors for cardiovascular disease. Many of the loci identified are characteristic of imprinted regions and are associated with the expression of nearby genes.

A homozygous loss-of-function mutation leading to CYBC1 deficiency causes chronic granulomatous disease

Mutations in genes encoding subunits of the phagocyte NADPH oxidase complex are recognized to cause chronic granulomatous disease (CGD), a severe primary immunodeficiency. Here we describe how deficiency of CYBC1, a previously uncharacterized protein in humans (C17orf62), leads to reduced expression of NADPH oxidase’s main subunit (gp91^phox) and results in CGD. Analyzing two brothers diagnosed with CGD we identify a homozygous loss-of-function mutation, p.Tyr2Ter, in CYBC1. Imputation of p.Tyr2Ter into 155K chip-genotyped Icelanders reveals six additional homozygotes, all with signs of CGD, manifesting as colitis, rare infections, or a severely impaired PMA-induced neutrophil oxidative burst. Homozygosity for p.Tyr2Ter consequently associates with inflammatory bowel disease (IBD) in Iceland (P = 8.3 × 10⁻⁸; OR = 67.6), as well as reduced height (P = 3.3 × 10⁻⁴; −8.5 cm). Overall, we find that CYBC1 deficiency results in CGD characterized by colitis and a distinct profile of infections indicative of macrophage dysfunction.

Mutations in genes encoding NAPDH oxidase subunits are known to be causative for the primary immunodeficiency chronic granulomatous disease (CGD). Here, the authors identify CYBC1 mutations in patients with CGD and show that CYBC1 is important for formation of the NADPH complex and respiratory burst.

Parent of origin gene expression in a founder population identifies two new candidate imprinted genes at known imprinted regions

Genomic imprinting is the phenomena that leads to silencing of one copy of a gene inherited from a specific parent. Mutations in imprinted regions have been involved in diseases showing parent of origin effects. Identifying genes with evidence of parent of origin expression patterns in family studies allows the detection of more subtle imprinting. Here, we use allele specific expression in lymphoblastoid cell lines from 306 Hutterites related in a single pedigree to provide formal evidence for parent of origin effects. We take advantage of phased genotype data to assign parent of origin to RNA-seq reads in individuals with gene expression data. Our approach identified known imprinted genes, two putative novel imprinted genes, PXDC1 and PWAR6, and 14 genes with asymmetrical parent of origin gene expression. We used gene expression in peripheral blood leukocytes (PBL) to validate our findings, and then confirmed imprinting control regions (ICRs) using DNA methylation levels in the PBLs.

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<P < 10-300). Using an independent dataset from the Brisbane Systems Genetic Study, we replicated 76 out of the 82 identified loci. POEs were remarkably consistent across time points and were so strong at some loci that methylation levels enabled good discrimination of parental transmissions at these and surrounding genomic regions. The implication is that parental allelic transmissions could be modelled at many imprinted (and linked) loci in GWAS of unrelated individuals given a combination of genetic and methylation data. Novel regions showing parent of origin effects on methylation will require replication using a different technology and further functional experiments to confirm that such effects arise through a genomic imprinting mechanism.

Biobank-driven genomic discovery yields new insight into atrial fibrillation biology

To identify genetic variation underlying atrial fibrillation, the most common cardiac arrhythmia, we performed a genome-wide association study of >1,000,000 people, including 60,620 atrial fibrillation cases and 970,216 controls. We identified 142 independent risk variants at 111 loci and prioritized 151 functional candidate genes likely to be involved in atrial fibrillation. Many of the identified risk variants fall near genes where more deleterious mutations have been reported to cause serious heart defects in humans (GATA4, MYH6, NKX2-5, PITX2, TBX5)¹, or near genes important for striated muscle function and integrity (for example, CFL2, MYH7, PKP2, RBM20, SGCG, SSPN). Pathway and functional enrichment analyses also suggested that many of the putative atrial fibrillation genes act via cardiac structural remodeling, potentially in the form of an 'atrial cardiomyopathy'², either during fetal heart development or as a response to stress in the adult heart.

Coding variants in RPL3L and MYZAP increase risk of atrial fibrillation

Most sequence variants identified hitherto in genome-wide association studies (GWAS) of atrial fibrillation are common, non-coding variants associated with risk through unknown mechanisms. We performed a meta-analysis of GWAS of atrial fibrillation among 29,502 cases and 767,760 controls from Iceland and the UK Biobank with follow-up in samples from Norway and the US, focusing on low-frequency coding and splice variants aiming to identify causal genes. We observe associations with one missense (OR = 1.20) and one splice-donor variant (OR = 1.50) in RPL3L, the first ribosomal gene implicated in atrial fibrillation to our knowledge. Analysis of 167 RNA samples from the right atrium reveals that the splice-donor variant in RPL3L results in exon skipping. We also observe an association with a missense variant in MYZAP (OR = 1.38), encoding a component of the intercalated discs of cardiomyocytes. Both discoveries emphasize the close relationship between the mechanical and electrical function of the heart.

The multiple myeloma risk allele at 5q15 lowers ELL2 expression and increases ribosomal gene expression

Recently, we identified ELL2 as a susceptibility gene for multiple myeloma (MM). To understand its mechanism of action, we performed expression quantitative trait locus analysis in CD138⁺ plasma cells from 1630 MM patients from four populations. We show that the MM risk allele lowers ELL2 expression in these cells (P_combined = 2.5 × 10⁻²⁷; β_combined = −0.24 SD), but not in peripheral blood or other tissues. Consistent with this, several variants representing the MM risk allele map to regulatory genomic regions, and three yield reduced transcriptional activity in plasmocytoma cell lines. One of these (rs3777189-C) co-locates with the best-supported lead variants for ELL2 expression and MM risk, and reduces binding of MAFF/G/K family transcription factors. Moreover, further analysis reveals that the MM risk allele associates with upregulation of gene sets related to ribosome biogenesis, and knockout/knockdown and rescue experiments in plasmocytoma cell lines support a cause–effect relationship. Our results provide mechanistic insight into MM predisposition.

ELL2 was recently discovered as a susceptibility gene for multiple myeloma (MM). Here, they show that the MM risk allele lowers ELL2 expression in plasma cells, that it also upregulates gene sets related to ribosome biogenesis, and that one of the linked variants reduces binding of MAFF/G/K family transcription factors.

Genetic short stature

Adult height and growth patterns are largely genetically programmed. Studies in twins have indicated that the heritability of height is high (>80%), suggesting that genetic variation is the main determinant of stature. Height exhibits a normal (Gaussian) distribution according to sex, age, and ancestry. Short stature is usually defined as a height which is 2 standard deviations (S.D.) less than the mean height of a specific population. This definition includes 2.3% of the population and usually includes healthy individuals. In this group of short stature non-syndromic conditions, the genetic influence occurs polygenically or oligogenically. As a rule, each common genetic variant accounts for a small effect (1mm) on individual height variation. Recently, several studies demonstrated that some rare variants can cause greater effect on height, without causing a syndromic condition. In more extreme cases, height SDS below 2.5 or 3 (which would comprise approximately 0.6 and 0.1% of the population, respectively) is frequently associated with syndromic conditions and are usually caused by a monogenic defect. More than 1,000 inherited/genetic diseases have growth disorder as an important phenotype. These conditions are usually responsible for syndromic short stature. In the coming years, we expect to discover several genetic causes of short stature, thereby explaining the phenotype of what we currently classify as short stature of unknown cause. These discoveries will have a profound impact on the follow-up and treatment of these children.

Association between Heights of Dialysis Patients and Outcomes: Results from a Retrospective Cohort Study of the International MONitoring Dialysis Outcomes (MONDO) Database Initiative

BACKGROUND

Tall people have improved metabolic profiles and better cardiovascular outcomes, a relationship inverted in dialysis patients. We investigated the relationship between height and outcomes in incident hemodialysis (HD) patients commencing treatment in an analysis of the international Monitoring Dialysis Outcomes (MONDO) database.

METHODS

In this retrospective cohort study, we included incident HD patients commencing treatment between -January 1, 2006 and December 31, 2010 and investigated the association between height and mortality using the MONDO database. A 6-months baseline period preceded 2.5 years of follow-up, during which we recorded patient mortality. Patients were stratified in region-specific deciles of the respective database's population (Asia Pacific, North and South America, and Europe) and we developed Cox-proportional hazard models (additionally adjusted for age, gender, post-dialysis weight, eKt/V, albumin, interdialytic weight gain, phosphorus, and predialysis systolic blood pressure) for each database.

RESULTS

We studied 23,353 patients (62 ± 15 years old, 42% female, body mass index 26 ± 6 kg/m2, height 165 ± 10 cm). We found a trend of increasing hazard ratio of death (HR) with increasing height for Asia Pacific, Europe, and South America. In the fully adjusted models, for South America, we found a trend of increasing HR without significance among deciles >5. In Europe, deciles 8-10 had significantly increased HR. No clear trend was found in North America.

CONCLUSION

We found an increasing risk of death with increasing height in all regions, except North America. While the reasons remain unclear, further research may be warranted.

Associations between height and blood pressure in the United States population

The mechanisms linking short stature with an increase in cardiovascular and cerebrovascular disease risk remain elusive. This study tested the hypothesis that significant associations are present between height and blood pressure in a representative sample of the US adult population.Participants were 12,988 men and women from a multiethnic sample (age ≥ 18 years) evaluated in the 1999 to 2006 National Health and Nutrition Examination Survey who were not taking antihypertensive medications and who had complete height, weight, % body fat, and systolic and diastolic arterial blood pressure (SBP and DBP) measurements; mean arterial blood pressure and pulse pressure (MBP and PP) were calculated. Multiple regression models for men and women were developed with each blood pressure as dependent variable and height, age, race/ethnicity, body mass index, % body fat, socioeconomic status, activity level, and smoking history as potential independent variables.Greater height was associated with significantly lower SBP and PP, and higher DBP (all P < .001) in combined race/ethnic-sex group models beginning in the 4th decade. Predicted blood pressure differences between people who are short and tall increased thereafter with greater age except for MBP. Socioeconomic status, activity level, and smoking history did not consistently contribute to blood pressure prediction models.Height-associated blood pressure effects were present in US adults who appeared in the 4th decade and increased in magnitude with greater age thereafter. These observations, in the largest and most diverse population sample evaluated to date, provide support for postulated mechanisms linking adult stature with cardiovascular and cerebrovascular disease risk.

rs10732516 polymorphism at the IGF2/H19 locus associates with a genotype-specific trend in placental DNA methylation and head circumference of prenatally alcohol-exposed newborns

STUDY QUESTION

Does prenatal alcohol exposure (PAE) affect regulation of the insulin-like growth factor 2 (IGF2)/H19 locus in placenta and the growth-restricted phenotype of newborns?

SUMMARY ANSWER

PAE results in genotype-specific trends in both placental DNA methylation at the IGF2/H19 locus and head circumference (HC) of newborns.

WHAT IS KNOWN ALREADY

PAE can disturb development of the nervous system and lead to restricted growth of the head, even microcephaly. To clarify the etiology of alcohol-induced growth restriction, we focused on the imprinted IGF2/H19 locus known to be important for normal placental and embryonic growth. The expression of IGF2 and a negative growth controller H19 are regulated by the H19 imprinting control region (H19 ICR) with seven-binding sites for the methylation-sensitive zinc-finger regulatory protein CTCF. A single nucleotide polymorphism rs10732516 G/A in the sixth-binding site has shown to associate with genotype-specific DNA methylation profiles at the H19 ICR.

STUDY DESIGN SIZE DURATION

By grouping 39 alcohol-exposed and 100 control samples according to rs10732516 polymorphism we explored alcohol-induced, genotype-specific changes in DNA methylation at the H19 ICR and the promoter region of H19 (H19 differentially methylated region). Also, IGF2 and H19 mRNA expression level in placenta as well as the phenotypes of newborns were examined.

PARTICIPANTS/MATERIALS SETTING METHODS

We explored alcohol-induced, genotype-specific changes in placental DNA methylation by MassARRAY EpiTYPER and allele-specific changes by bisulphite sequencing. IGF2 and H19 expression in placenta were analyzed by quantitative PCR and the HC, birthweight and birth length of newborns were examined using national growth charts.

MAIN RESULTS AND THE ROLE OF CHANCE

We observed a consistent trend in genotype-specific changes in DNA methylation at H19 ICR in alcohol-exposed placentas. DNA methylation level in the normally highly methylated paternal allele of rs10732516 paternal A/maternal G genotype was decreased in alcohol-exposed placentas. In addition to decreased IGF2 mRNA expression in alcohol-exposed placentas of this specific genotype (P = 0.03), we observed significantly increased expression of H19 in relation to IGF2 when comparing all alcohol-exposed placentas to unexposed controls (P = 0.006). Furthermore, phenotypic examination showed a significant genotype-specific association between the alcohol exposure and HC of newborns (P = 0.001).

LIMITATIONS REASONS FOR CAUTION

Owing to the exceptional character of the alcohol-exposed human samples collected in this study, the sample size is restricted. An increased sample size and functional studies are needed to confirm these data and clarify the biological significance or causality of the observed associations.

WIDER IMPLICATIONS OF THE FINDINGS

Our results suggest that the rs10732516 polymorphism associates with the alcohol-induced alterations in DNA methylation profiles and head growth in a parent-of-origin manner. We also introduce a novel genotype-specific approach for exploring environmental effects on the IGF2/H19 locus and ultimately on embryonic growth.

STUDY FUNDING/COMPETING INTERESTS

This work was supported by the Academy of Finland (258304), The Finnish Foundation for Alcohol Studies, Finnish Cultural Foundation, Juho Vainio Foundation, Yrjö Jahnsson Foundation and Arvo and Lea Ylppö Foundation. No competing interests are declared.