Synthesizing genome-wide association studies and expression microarray reveals novel genes that act in the human growth plate to modulate height

Oxidative stress and mitochondrial impairment: Key drivers in neurodegenerative disorders

Mitochondrial dysfunction and oxidative stress are critical factors in the pathogenesis of neurodegenerative diseases. The complex interplay between these factors exacerbates neuronal damage and accelerates disease progression. In neurodegenerative diseases, mitochondrial dysfunction impairs ATP production and promotes the generation of reactive oxygen species (ROS). The accumulation of ROS further damages mitochondrial DNA, proteins, and lipids, creating a vicious cycle of oxidative stress and mitochondrial impairment. This review aims to elucidate the mechanisms by which mitochondrial dysfunction and oxidative stress lead to neurodegeneration, and to highlight potential therapeutic targets to mitigate their harmful effects.

Nutrition, Other Environmental Influences, and Genetics in the Determination of Human Stature

Linear growth during three distinct stages of life determines attained stature in adulthood: namely, in utero, early postnatal life, and puberty and the adolescent period. Individual host factors, genetics, and the environment, including nutrition, influence attained human stature. Each period of physical growth has its specific biological and environmental considerations. Recent epidemiologic investigations reveal a strong influence of prenatal factors on linear size at birth that in turn influence the postnatal growth trajectory. Although average population height changes have been documented in high-income regions, stature as a complex human trait is not well understood or easily modified. This review summarizes the biology of linear growth and its major drivers, including nutrition from a life-course perspective, the genetics of programmed growth patterns or height, and gene-environment interactions that determine human stature in toto over the life span. Implications for public health interventions and knowledge gaps are discussed.

Growth disorders caused by variants in epigenetic regulators: progress and prospects

Epigenetic modifications play an important role in regulation of transcription and gene expression. The molecular machinery governing epigenetic modifications, also known as epigenetic regulators, include non-coding RNA, chromatin remodelers, and enzymes or proteins responsible for binding, reading, writing and erasing DNA and histone modifications. Recent advancement in human genetics and high throughput sequencing technology have allowed the identification of causative variants, many of which are epigenetic regulators, for a wide variety of childhood growth disorders that include skeletal dysplasias, idiopathic short stature, and generalized overgrowth syndromes. In this review, we highlight the connection between epigenetic modifications, genetic variants in epigenetic regulators and childhood growth disorders being established over the past decade, discuss their insights into skeletal biology, and the potential of epidrugs as a new type of therapeutic intervention.

Identification of novel genes including NAV2 associated with isolated tall stature

Very tall people attract much attention and represent a clinically and genetically heterogenous group of individuals. Identifying the genetic etiology can provide important insights into the molecular mechanisms regulating linear growth. We studied a three-generation pedigree with five isolated (non-syndromic) tall members and one individual with normal stature by whole exome sequencing; the tallest man had a height of 211 cm. Six heterozygous gene variants predicted as damaging were shared among the four genetically related tall individuals and not present in a family member with normal height. To gain insight into the putative role of these candidate genes in bone growth, we assessed the transcriptome of murine growth plate by microarray and RNA Seq. Two (Ift140, Nav2) of the six genes were well-expressed in the growth plate. Nav2 (p-value 1.91E-62) as well as Ift140 (p-value of 2.98E-06) showed significant downregulation of gene expression between the proliferative and hypertrophic zone, suggesting that these genes may be involved in the regulation of chondrocyte proliferation and/or hypertrophic differentiation. IFT140, NAV2 and SCAF11 have also significantly associated with height in GWAS studies. Pathway and network analysis indicated functional connections between IFT140, NAV2 and SCAF11 and previously associated (tall) stature genes. Knockout of the all-trans retinoic acid responsive gene, neuron navigator 2 NAV2, in Xenopus supports its functional role as a growth promotor. Collectively, our data expand the spectrum of genes with a putative role in tall stature phenotypes and, among other genes, highlight NAV2 as an interesting gene to this phenotype.

Boosting the power of genome-wide association studies within and across ancestries by using polygenic scores

Genome-wide association studies (GWASs) have been mostly conducted in populations of European ancestry, which currently limits the transferability of their findings to other populations. Here, we show, through theory, simulations and applications to real data, that adjustment of GWAS analyses for polygenic scores (PGSs) increases the statistical power for discovery across all ancestries. We applied this method to analyze seven traits available in three large biobanks with participants of East Asian ancestry (n = 340,000 in total) and report 139 additional associations across traits. We also present a two-stage meta-analysis strategy whereby, in contributing cohorts, a PGS-adjusted GWAS is rerun using PGSs derived from a first round of a standard meta-analysis. On average, across traits, this approach yields a 1.26-fold increase in the number of detected associations (range 1.07- to 1.76-fold increase). Altogether, our study demonstrates the value of using PGSs to increase the power of GWASs in underrepresented populations and promotes such an analytical strategy for future GWAS meta-analyses.

Biobank-scale methods and projections for sparse polygenic prediction from machine learning

In this paper we characterize the performance of linear models trained via widely-used sparse machine learning algorithms. We build polygenic scores and examine performance as a function of training set size, genetic ancestral background, and training method. We show that predictor performance is most strongly dependent on size of training data, with smaller gains from algorithmic improvements. We find that LASSO generally performs as well as the best methods, judged by a variety of metrics. We also investigate performance characteristics of predictors trained on one genetic ancestry group when applied to another. Using LASSO, we develop a novel method for projecting AUC and correlation as a function of data size (i.e., for new biobanks) and characterize the asymptotic limit of performance. Additionally, for LASSO (compressed sensing) we show that performance metrics and predictor sparsity are in agreement with theoretical predictions from the Donoho-Tanner phase transition. Specifically, a future predictor trained in the Taiwan Precision Medicine Initiative for asthma can achieve an AUC of [Formula: see text] and for height a correlation of [Formula: see text] for a Taiwanese population. This is above the measured values of [Formula: see text] and [Formula: see text], respectively, for UK Biobank trained predictors applied to a European population.

Genetics of skeletal proportions in two different populations

Human height can be divided into sitting height and leg length, reflecting growth of different parts of the skeleton whose relative proportions are captured by the ratio of sitting to total height (as sitting height ratio, SHR). Height is a highly heritable trait, and its genetic basis has been well-studied. However, the genetic determinants of skeletal proportion are much less well-characterized. Expanding substantially on past work, we performed a genome-wide association study (GWAS) of SHR in ∼450,000 individuals with European ancestry and ∼100,000 individuals with East Asian ancestry from the UK and China Kadoorie Biobanks. We identified 565 loci independently associated with SHR, including all genomic regions implicated in prior GWAS in these ancestries. While SHR loci largely overlap height-associated loci (P < 0.001), the fine-mapped SHR signals were often distinct from height. We additionally used fine-mapped signals to identify 36 credible sets with heterogeneous effects across ancestries. Lastly, we used SHR, sitting height, and leg length to identify genetic variation acting on specific body regions rather than on overall human height.

Genome-wide CRISPR screening of chondrocyte maturation newly implicates genes in skeletal growth and height-associated GWAS loci

Alterations in the growth and maturation of chondrocytes can lead to variation in human height, including monogenic disorders of skeletal growth. We aimed to identify genes and pathways relevant to human growth by pairing human height genome-wide association studies (GWASs) with genome-wide knockout (KO) screens of growth-plate chondrocyte proliferation and maturation in vitro. We identified 145 genes that alter chondrocyte proliferation and maturation at early and/or late time points in culture, with 90% of genes validating in secondary screening. These genes are enriched in monogenic growth disorder genes and in KEGG pathways critical for skeletal growth and endochondral ossification. Further, common variants near these genes capture height heritability independent of genes computationally prioritized from GWASs. Our study emphasizes the value of functional studies in biologically relevant tissues as orthogonal datasets to refine likely causal genes from GWASs and implicates new genetic regulators of chondrocyte proliferation and maturation.

A Genetic Approach in the Evaluation of Short Stature

Short stature is considered a condition in which the height is 2 standard deviations below the mean height of a given age, sex, and population group. Human height is a polygenic and heterogeneous characteristic, and its heritability is reported to be approximately 80%. More than 600 variants associated with human growth were detected in the genome-wide association studies. Rare and common variants concurrently affect human height. The rare variations that play a role in human height determination and have a strong impact on protein functions lead to monogenic short stature phenotypes, which are a highly heterogeneous group. With rapidly developing technologies in the last decade, molecular genetic tests have begun to be used widely in clinical genetics, and thus, the genetic etiology of several rare diseases has been elucidated. Identifying the genetic etiology underlying idiopathic short stature which represents phenotypically heterogeneous group of diseases ranging from isolated short stature to severe and syndromic short stature has promoted the understanding of the genetic regulation of growth plate and longitudinal bone growth. In cases of short stature, definite molecular diagnosis based on genetic evaluation enables the patient and family to receive genetic counseling on the natural course of the disease, prognosis, genetic basis, and recurrence risk. The determination of the genetic etiology in growth disorders is essential for the development of novel targeted therapies and crucial in the development of mutation-specific treatments in the future.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Genes with specificity for expression in the round cell layer of the growth plate are enriched in genomewide association study (GWAS) of human height

Human adult height reflects the outcome of childhood skeletal growth. Growth plate (epiphyseal) chondrocytes are key determinants of height. As epiphyseal chondrocytes mature and proliferate, they pass through three developmental stages, which are organized into three distinct layers in the growth plate: (i) resting (round), (ii) proliferative (flat), and (iii) hypertrophic. Recent genomewide association studies (GWASs) of human height identified numerous associated loci, which are enriched for genes expressed in growth plate chondrocytes. However, it remains unclear which specific genes expressed in which layers of the growth plate regulate skeletal growth and human height. To connect the genetics of height and growth plate biology, we analyzed GWAS data through the lens of gene expression in the three dissected layers of murine newborn tibial growth plate. For each gene, we derived a specificity score for each growth plate layer and regressed these scores against gene-level p values from recent height GWAS data. We found that specificity for expression in the round cell layer, which contains chondrocytes early in maturation, is significantly associated with height GWAS p values (p = 8.5 × 10-9 ); this association remains after conditioning on specificity for the other cell layers. The association also remains after conditioning on membership in an "Online Mendelian Inheritance in Man (OMIM) gene set" (genes known to cause monogenic skeletal growth disorders, p < 9.7 × 10-6 ). We replicated the association in RNA-sequencing (RNA-seq) data from maturing chondrocytes sampled at early and late time points during differentiation in vitro: we found that expression early in differentiation is significantly associated with p values from height GWASs (p = 6.1 × 10-10 ) and that this association remains after conditioning on expression at 10 days in culture and on the OMIM gene set (p < 0.006). These findings newly implicate genes highlighted by GWASs of height and specifically expressed in the round cell layer as being potentially important regulators of skeletal biology. © 2021 American Society for Bone and Mineral Research (ASBMR).

Genetic variant of TTLL11 gene and subsequent ciliary defects are associated with idiopathic scoliosis in a 5-generation UK family

Idiopathic scoliosis (IS) is a complex 3D deformation of the spine with a strong genetic component, most commonly found in adolescent girls. Adolescent idiopathic scoliosis (AIS) affects around 3% of the general population. In a 5-generation UK family, linkage analysis identified the locus 9q31.2-q34.2 as a candidate region for AIS; however, the causative gene remained unidentified. Here, using exome sequencing we identified a rare insertion c.1569_1570insTT in the tubulin tyrosine ligase like gene, member 11 (TTLL11) within that locus, as the IS causative gene in this British family. Two other TTLL11 mutations were also identified in two additional AIS cases in the same cohort. Analyses of primary cells of individuals carrying the c.1569_1570insTT (NM_194252) mutation reveal a defect at the primary cilia level, which is less present, smaller and less polyglutamylated compared to control. Further, in a zebrafish, the knock down of ttll11, and the mutated ttll11 confirmed its role in spine development and ciliary function in the fish retina. These findings provide evidence that mutations in TTLL11, a ciliary gene, contribute to the pathogenesis of IS.

Growth plate gene involment and isolated short stature

PURPOSE

Short stature is a common clinical presentation, thus it is widely accepted that it is a polygenic trait. However, genome wide association and next generation sequencing studies have recently challenged this view, suggesting that many of the children classified as idiopathic short stature could instead have monogenic defects. Linear growth is determined primarily by chondrogenesis at the growth plate. This process results from chondrocyte proliferation, hypertrophy, and extracellular matrix secretion, and it is perfectly coordinated by complex networks of local paracrine and endocrine factors. Alterations in genes which control growth plate development can explain a large number of cases of isolated short stature, allowing an etiological diagnosis.

METHODS/RESULTS

We reviewed recent data on the genetic alterations in fundamental cellular processes, paracrine signaling, and cartilage matrix formation associated with impaired growth plate chondrogenesis. In particular we focused on growth plate gene involvement in nonsyndromic short stature.

CONCLUSIONS

The identification of genetic basis of growth failure will have a significant impact on the care of children affected with short stature.

Evidence That Non-Syndromic Familial Tall Stature Has an Oligogenic Origin Including Ciliary Genes

Human growth is a complex trait. A considerable number of gene defects have been shown to cause short stature, but there are only few examples of genetic causes of non-syndromic tall stature. Besides rare variants with large effects and common risk alleles with small effect size, oligogenic effects may contribute to this phenotype. Exome sequencing was carried out in a tall male (height 3.5 SDS) and his parents. Filtered damaging variants with high CADD scores were validated by Sanger sequencing in the trio and three other affected and one unaffected family members. Network analysis was carried out to assess links between the candidate genes, and the transcriptome of murine growth plate was analyzed by microarray as well as RNA Seq. Heterozygous gene variants in CEP104, CROCC, NEK1, TOM1L2, and TSTD2 predicted as damaging were found to be shared between the four tall family members. Three of the five genes (CEP104, CROCC, and NEK1) belong to the ciliary gene family. All genes are expressed in mouse growth plate. Pathway and network analyses indicated close functional connections. Together, these data expand the spectrum of genes with a role in linear growth and tall stature phenotypes.

Physiological regulation of bone length and skeletal proportion in mammals

NEW FINDINGS

What is the topic of this review? Mechanisms regulating bone length and skeletal proportions What advances does it highlight? The study of differential bone length between leg and finger bones, metatarsals of the Egyptian jerboa and genomic analysis of giraffes.

ABSTRACT

Among mammalian species, skeletal structures vary greatly in size and shape, leading to a dramatic variety of body sizes and proportions. How different bones grow to different lengths, whether among different species, different individuals of the same species, or even in different anatomical parts of our the body, has always been a fascinating subject of research in biology and physiology. In the current review, we focus on some of the recent advances in the field and discuss how these provided important new insights into the mechanisms regulating bone length and skeletal proportions.

Short Stature is Progressive in Patients with Heterozygous NPR2 Mutations

BACKGROUND

NPR2 encodes atrial natriuretic peptide receptor B (ANPRB), a regulator of skeletal growth. Biallelic loss-of-function mutations in NPR2 result in acromesomelic dysplasia Maroteaux type (AMDM; OMIM 602875), while heterozygous mutations may account for 2% to 6% of idiopathic short stature (ISS).

OBJECTIVE

Describe the physical proportions and growth characteristics of an extended family with novel NPR2 mutations including members with AMDM, ISS, or normal stature.

DESIGN AND PARTICIPANTS

We performed whole exome sequencing in 2 healthy parents and 2 children with AMDM. Detailed genotyping and phenotyping were performed on members of a multigenerational family in an academic medical center. We expressed mutant proteins in mammalian cells and characterized expression and function.

RESULTS

The sisters with AMDM were compound heterozygotes for missense mutations in the NPR2 gene, a novel p.P93S (maternal) and the previously reported p.R989L (paternal). Both mutant ANPRB proteins were normally expressed in HEK293T cells and exhibited dominant negative effects on wild-type ANPRB catalytic activity. Heterozygous relatives had proportionate short stature (height z-scores -2.06 ± 0.97, median ± SD) compared with their wild-type siblings (-1.37 ± 0.59). Height z-scores progressively and significantly decreased as NPR2-heterozygous children matured, while remaining constant in their wild-type siblings.

CONCLUSIONS

Biallelic NPR2 mutations cause severe skeletal dysplasia (AMDM), whereas heterozygous mutations lead to a subtler phenotype characterized by progressive short stature with by increasing loss of height potential with age.

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.

Genetic regulation of linear growth

Linear growth occurs at the growth plate. Therefore, genetic defects that interfere with the normal function of the growth plate can cause linear growth disorders. Many genetic causes of growth disorders have already been identified in humans. However, recent genome-wide approaches have broadened our knowledge of the mechanisms of linear growth, not only providing novel monogenic causes of growth disorders but also revealing single nucleotide polymorphisms in genes that affect height in the general population. The genes identified as causative of linear growth disorders are heterogeneous, playing a role in various growth-regulating mechanisms including those involving the extracellular matrix, intracellular signaling, paracrine signaling, endocrine signaling, and epigenetic regulation. Understanding the underlying genetic defects in linear growth is important for clinicians and researchers in order to provide proper diagnoses, management, and genetic counseling, as well as to develop better treatment approaches for children with growth disorders.

Polygenic adaptation and convergent evolution on growth and cardiac genetic pathways in African and Asian rainforest hunter-gatherers

Different human populations facing similar environmental challenges have sometimes evolved convergent biological adaptations, for example, hypoxia resistance at high altitudes and depigmented skin in northern latitudes on separate continents. The "pygmy" phenotype (small adult body size), characteristic of hunter-gatherer populations inhabiting both African and Asian tropical rainforests, is often highlighted as another case of convergent adaptation in humans. However, the degree to which phenotypic convergence in this polygenic trait is due to convergent versus population-specific genetic changes is unknown. To address this question, we analyzed high-coverage sequence data from the protein-coding portion of the genomes of two pairs of populations: Batwa rainforest hunter-gatherers and neighboring Bakiga agriculturalists from Uganda and Andamanese rainforest hunter-gatherers and Brahmin agriculturalists from India. We observed signatures of convergent positive selection between the rainforest hunter-gatherers across the set of genes with "growth factor binding" functions ([Formula: see text]). Unexpectedly, for the rainforest groups, we also observed convergent and population-specific signatures of positive selection in pathways related to cardiac development (e.g., "cardiac muscle tissue development"; [Formula: see text]). We hypothesize that the growth hormone subresponsiveness likely underlying the adult small body-size phenotype may have led to compensatory changes in cardiac pathways, in which this hormone also plays an essential role. Importantly, in the agriculturalist populations, we did not observe similar patterns of positive selection on sets of genes associated with growth or cardiac development, indicating our results most likely reflect a history of convergent adaptation to the similar ecology of rainforests rather than a more general evolutionary pattern.

Insights and Implications of Genome-Wide Association Studies of Height

CONTEXT

In the last decade, genome-wide association studies (GWASs) have catalyzed our understanding of the genetics of height and have identified hundreds of regions of the genome associated with adult height and other height-related body measurements.

EVIDENCE ACQUISITION

GWASs related to height were identified via PubMed search and a review of the GWAS catalog.

EVIDENCE SYNTHESIS

The GWAS results demonstrate that height is highly polygenic: that is, many thousands of genetic variants distributed across the genome each contribute to an individual's height. These height-associated regions of the genome are enriched for genes in known biological pathways involved in growth, such as fibroblast growth factor signaling, as well as for genes expressed in relevant tissues, such as the growth plate. GWASs can also uncover previously unappreciated biological pathways, such as the STC2/PAPPA/IGFBP4 pathway. The genes implicated by GWASs are often the same genes that are the genetic causes of Mendelian growth disorders or skeletal dysplasias, and GWAS results can provide complementary information about these disorders.

CONCLUSIONS

Here, we review the rationale behind GWASs and what we have learned from GWASs for height, including how it has enhanced our understanding of the underlying biology of human growth. We also highlight the implications of GWASs in terms of prediction of adult height and our understanding of Mendelian growth disorders.

Growth plate expression profiling: Large and small breed dogs provide new insights in endochondral bone formation

The difference in the adult height of mammals, and hence in endochondral bone formation, is not yet fully understood and may serve to identify targets for bone and cartilage regeneration. In line with this hypothesis, the intra-species disparity between the adult height of Great Danes and Miniature Poodles was investigated at a transcriptional level. Microarray analysis of the growth plate of five Great Danes and five Miniature Poodles revealed 2,981 unique genes that were differentially expressed, including many genes with an unknown role in skeletal development. A signaling pathway impact analysis indicated activation of the cell cycle, extracellular matrix receptor interaction and the tight junction pathway, and inhibition of pathways associated with inflammation and the complement cascade. In additional validation steps, the gene expression profile of the separate growth plate zones for both dog breeds were determined. Given that the BMP signaling is known for its crucial role in skeletal development and fracture healing, and BMP-2 is used in orthopaedic and spine procedures for bone augmentation, further investigations concentrated on the BMP pathway.The canonical BMP-2 and BMP-6 signaling pathway was activated in the Great Danes compared to Miniature Poodles. In conclusion, investigating the differential expression of genes involved in endochondral bone formation in small and large breed dogs, could be a game changing strategy to provide new insights in growth plate development and identify new targets for bone and cartilage regeneration. © 2017 The Authors. Journal of Orthopaedic Research® published by Wiley Periodicals, Inc. on behalf of the Orthopaedic Research Society. J Orthop Res 36:138-148, 2018.

Epigenetic profiling of growth plate chondrocytes sheds insight into regulatory genetic variation influencing height

GWAS have identified hundreds of height-associated loci. However, determining causal mechanisms is challenging, especially since height-relevant tissues (e.g. growth plates) are difficult to study. To uncover mechanisms by which height GWAS variants function, we performed epigenetic profiling of murine femoral growth plates. The profiled open chromatin regions recapitulate known chondrocyte and skeletal biology, are enriched at height GWAS loci, particularly near differentially expressed growth plate genes, and enriched for binding motifs of transcription factors with roles in chondrocyte biology. At specific loci, our analyses identified compelling mechanisms for GWAS variants. For example, at CHSY1, we identified a candidate causal variant (rs9920291) overlapping an open chromatin region. Reporter assays demonstrated that rs9920291 shows allelic regulatory activity, and CRISPR/Cas9 targeting of human chondrocytes demonstrates that the region regulates CHSY1 expression. Thus, integrating biologically relevant epigenetic information (here, from growth plates) with genetic association results can identify biological mechanisms important for human growth.

Genetics of Short Stature

Short stature is a common and heterogeneous condition that is often genetic in etiology. For most children with genetic short stature, the specific molecular causes remain unknown; but with advances in exome/genome sequencing and bioinformatics approaches, new genetic causes of growth disorders have been identified, contributing to the understanding of the underlying molecular mechanisms of longitudinal bone growth and growth failure. Identifying new genetic causes of growth disorders has the potential to improve diagnosis, prognostic accuracy, and individualized management, and help avoid unnecessary testing for endocrine and other disorders.

Functional interactors of three genome-wide association study genes are differentially expressed in severe chronic obstructive pulmonary disease lung tissue

In comparison to genome-wide association studies (GWAS), there has been poor replication of gene expression studies in chronic obstructive pulmonary disease (COPD). We performed microarray gene expression profiling on a large sample of resected lung tissues from subjects with severe COPD. Comparing 111 COPD cases and 40 control smokers, 204 genes were differentially expressed; none were at significant GWAS loci. The top differentially expressed gene was HMGB1, which interacts with AGER, a known COPD GWAS gene. Differentially expressed genes showed enrichment for putative interactors of the first three identified COPD GWAS genes IREB2, HHIP, and FAM13A, based on gene sets derived from protein and RNA binding studies, RNA-interference, a murine smoking model, and expression quantitative trait locus analyses. The gene module most highly associated for COPD in Weighted Gene Co-Expression Network Analysis (WGCNA) was enriched for B cell pathways, and shared seventeen genes with a mouse smoking model and twenty genes with previous emphysema studies. As in other common diseases, genes at COPD GWAS loci were not differentially expressed; however, using a combination of network methods, experimental studies and careful phenotype definition, we found differential expression of putative interactors of these genes, and we replicated previous human and mouse microarray results.

Growth in patients with type 1 diabetes

PURPOSE OF REVIEW

As the incidence of type 1 diabetes (T1DM) continues to rise, complications including impairment of childhood growth remain a major concern. This review provides an overview of alterations in growth patterns before and after the onset of T1DM.

RECENT FINDINGS

Recent advances in this field include several prospective investigations of height and weight trajectories in children leading up to the development of islet autoimmunity and T1DM as well as evaluations of larger cohorts of T1DM patients to better assess predictors of altered growth. In addition, genetic and metabolic investigations have improved our understanding of the more rare severe growth impairment of Mauriac Syndrome.

SUMMARY

Despite advances in medical care of children with T1DM, growth remains suboptimal in this population and likely reflects ongoing metabolic derangement linked with classic microvascular diabetic complications.

A systematic SNP selection approach to identify mechanisms underlying disease aetiology: linking height to post-menopausal breast and colorectal cancer risk

Data from GWAS suggest that SNPs associated with complex diseases or traits tend to co-segregate in regions of low recombination, harbouring functionally linked gene clusters. This phenomenon allows for selecting a limited number of SNPs from GWAS repositories for large-scale studies investigating shared mechanisms between diseases. For example, we were interested in shared mechanisms between adult-attained height and post-menopausal breast cancer (BC) and colorectal cancer (CRC) risk, because height is a risk factor for these cancers, though likely not a causal factor. Using SNPs from public GWAS repositories at p-values < 1 × 10⁻⁵ and a genomic sliding window of 1 mega base pair, we identified SNP clusters including at least one SNP associated with height and one SNP associated with either post-menopausal BC or CRC risk (or both). SNPs were annotated to genes using HapMap and GRAIL and analysed for significantly overrepresented pathways using ConsensuspathDB. Twelve clusters including 56 SNPs annotated to 26 genes were prioritised because these included at least one height- and one BC risk- or CRC risk-associated SNP annotated to the same gene. Annotated genes were involved in Indian hedgehog signalling (p-value = 7.78 × 10⁻⁷) and several cancer site-specific pathways. This systematic approach identified a limited number of clustered SNPs, which pinpoint potential shared mechanisms linking together the complex phenotypes height, post-menopausal BC and CRC.

A Weighted SNP Correlation Network Method for Estimating Polygenic Risk Scores

Polygenic scores are useful for examining the joint associations of genetic markers. However, because traditional methods involve summing weighted allele counts, they may fail to capture the complex nature of biology. Here we describe a network-based method, which we call weighted SNP correlation network analysis (WSCNA), and demonstrate how it could be used to generate meaningful polygenic scores. Using data on human height in a US population of non-Hispanic whites, we illustrate how this method can be used to identify SNP networks from GWAS data, create network-specific polygenic scores, examine network topology to identify hub SNPs, and gain biological insights into complex traits. In our example, we show that this method explains a larger proportion of the variance in human height than traditional polygenic score methods. We also identify hub genes and pathways that have previously been identified as influencing human height. In moving forward, this method may be useful for generating genetic susceptibility measures for other health related traits, examining genetic pleiotropy, identifying at-risk individuals, examining gene score by environmental effects, and gaining a deeper understanding of the underlying biology of complex traits.

EZH1 and EZH2 promote skeletal growth by repressing inhibitors of chondrocyte proliferation and hypertrophy

Histone methyltransferases EZH1 and EZH2 catalyse the trimethylation of histone H3 at lysine 27 (H3K27), which serves as an epigenetic signal for chromatin condensation and transcriptional repression. Genome-wide associated studies have implicated EZH2 in the control of height and mutations in EZH2 cause Weaver syndrome, which includes skeletal overgrowth. Here we show that the combined loss of Ezh1 and Ezh2 in chondrocytes severely impairs skeletal growth in mice. Both of the principal processes underlying growth plate chondrogenesis, chondrocyte proliferation and hypertrophy, are compromised. The decrease in chondrocyte proliferation is due in part to derepression of cyclin-dependent kinase inhibitors Ink4a/b, while ineffective chondrocyte hypertrophy is due to the suppression of IGF signalling by the increased expression of IGF-binding proteins. Collectively, our findings reveal a critical role for H3K27 methylation in the regulation of chondrocyte proliferation and hypertrophy in the growth plate, which are the central determinants of skeletal growth.

Chondrocyte FGFR3 Regulates Bone Mass by Inhibiting Osteogenesis

Chondrogenesis can regulate bone formation. Fibroblast growth factor receptor 3, highly expressed in chondrocytes, is a negative regulator of bone growth. To investigate whether chondrocyte FGFR3 regulates osteogenesis, thereby contributing to postnatal bone formation and bone remodeling, mice with conditional knock-out of Fgfr3 in chondrocytes (mutant (MUT)) were generated. MUT mice displayed overgrowth of bone with lengthened growth plates. Bone mass of MUT mice was significantly increased at both 1 month and 4 months of age. Histological analysis showed that osteoblast number and bone formation were remarkably enhanced after deletion of Fgfr3 in chondrocytes. Chondrocyte-osteoblast co-culture assay further revealed that Fgfr3 deficiency in chondrocytes promoted differentiation and mineralization of osteoblasts by up-regulating the expressions of Ihh, Bmp2, Bmp4, Bmp7, Wnt4, and Tgf-β1, as well as down-regulating Nog expression. In addition, osteoclastogenesis was also impaired in MUT mice with decreased number of osteoclasts lining trabecular bone, which may be related to the reduced ratio of Rankl to Opg in Fgfr3-deficient chondrocytes. This study reveals that chondrocyte FGFR3 is involved in the regulation of bone formation and bone remodeling by a paracrine mechanism.

A Track Record on SHOX: From Basic Research to Complex Models and Therapy

SHOX deficiency is the most frequent genetic growth disorder associated with isolated and syndromic forms of short stature. Caused by mutations in the homeobox gene SHOX, its varied clinical manifestations include isolated short stature, Léri-Weill dyschondrosteosis, and Langer mesomelic dysplasia. In addition, SHOX deficiency contributes to the skeletal features in Turner syndrome. Causative SHOX mutations have allowed downstream pathology to be linked to defined molecular lesions. Expression levels of SHOX are tightly regulated, and almost half of the pathogenic mutations have affected enhancers. Clinical severity of SHOX deficiency varies between genders and ranges from normal stature to profound mesomelic skeletal dysplasia. Treatment options for children with SHOX deficiency are available. Two decades of research support the concept of SHOX as a transcription factor that integrates diverse aspects of bone development, growth plate biology, and apoptosis. Due to its absence in mouse, the animal models of choice have become chicken and zebrafish. These models, therefore, together with micromass cultures and primary cell lines, have been used to address SHOX function. Pathway and network analyses have identified interactors, target genes, and regulators. Here, we summarize recent data and give insight into the critical molecular and cellular functions of SHOX in the etiopathogenesis of short stature and limb development.

MECHANISMS IN ENDOCRINOLOGY: Novel genetic causes of short stature

The fast technological development, particularly single nucleotide polymorphism array, array-comparative genomic hybridization, and whole exome sequencing, has led to the discovery of many novel genetic causes of growth failure. In this review we discuss a selection of these, according to a diagnostic classification centred on the epiphyseal growth plate. We successively discuss disorders in hormone signalling, paracrine factors, matrix molecules, intracellular pathways, and fundamental cellular processes, followed by chromosomal aberrations including copy number variants (CNVs) and imprinting disorders associated with short stature. Many novel causes of GH deficiency (GHD) as part of combined pituitary hormone deficiency have been uncovered. The most frequent genetic causes of isolated GHD are GH1 and GHRHR defects, but several novel causes have recently been found, such as GHSR, RNPC3, and IFT172 mutations. Besides well-defined causes of GH insensitivity (GHR, STAT5B, IGFALS, IGF1 defects), disorders of NFκB signalling, STAT3 and IGF2 have recently been discovered. Heterozygous IGF1R defects are a relatively frequent cause of prenatal and postnatal growth retardation. TRHA mutations cause a syndromic form of short stature with elevated T3/T4 ratio. Disorders of signalling of various paracrine factors (FGFs, BMPs, WNTs, PTHrP/IHH, and CNP/NPR2) or genetic defects affecting cartilage extracellular matrix usually cause disproportionate short stature. Heterozygous NPR2 or SHOX defects may be found in ∼3% of short children, and also rasopathies (e.g., Noonan syndrome) can be found in children without clear syndromic appearance. Numerous other syndromes associated with short stature are caused by genetic defects in fundamental cellular processes, chromosomal abnormalities, CNVs, and imprinting disorders.

Copy number variations in 119 Chinese children with idiopathic short stature identified by the custom genome-wide microarray

Background

Idiopathic short stature (ISS) refers to short stature with no evident etiologies. The custom genome-wide microarray specifically designed to cover height-related genes may be helpful to detect copy number variations (CNVs) in ISS patients, which may be missed by the general microarray. The aim of the study was to validate the applicability of the custom microarray and to analyze CNVs in Chinese ISS children.

Results

Sixty non-polymorphic CNVs were identified in 119 ISS patients. There were 13 small CNVs with a size below 50 kb, accounting for 21.7 % of all the CNVs (13/60). Five pathogenic or possibly pathogenic CNVs were detected in five patients, including deletions at 22q11.21, duplications at 4q11-q13.1, 4q12 and Yp11.32-p11.2. Taking only the pathogenic variants into account, the diagnostic yield was 2.5 % (3/119). The TMEM165, POLR2B and PDGFRA genes were analyzed as candidate genes. A 15 kb deletion in the RASA2 gene was of interest for further investigation.

Conclusions

This study showed that the custom microarray is applicable to detect CNVs in patients with short stature. Candidate genes and CNVs detected in ISS patients may be helpful for CNV analysis of short stature, especially in East Asian population.

Electronic supplementary material

The online version of this article (doi:10.1186/s13039-016-0225-0) contains supplementary material, which is available to authorized users.

Integrating Diverse Types of Genomic Data to Identify Genes that Underlie Adverse Pregnancy Phenotypes

Progress in understanding complex genetic diseases has been bolstered by synthetic approaches that overlay diverse data types and analyses to identify functionally important genes. Pre-term birth (PTB), a major complication of pregnancy, is a leading cause of infant mortality worldwide. A major obstacle in addressing PTB is that the mechanisms controlling parturition and birth timing remain poorly understood. Integrative approaches that overlay datasets derived from comparative genomics with function-derived ones have potential to advance our understanding of the genetics of birth timing, and thus provide insights into the genes that may contribute to PTB. We intersected data from fast evolving coding and non-coding gene regions in the human and primate lineage with data from genes expressed in the placenta, from genes that show enriched expression only in the placenta, as well as from genes that are differentially expressed in four distinct PTB clinical subtypes. A large fraction of genes that are expressed in placenta, and differentially expressed in PTB clinical subtypes (23–34%) are fast evolving, and are associated with functions that include adhesion neurodevelopmental and immune processes. Functional categories of genes that express fast evolution in coding regions differ from those linked to fast evolution in non-coding regions. Finally, there is a surprising lack of overlap between fast evolving genes that are differentially expressed in four PTB clinical subtypes. Integrative approaches, especially those that incorporate evolutionary perspectives, can be successful in identifying potential genetic contributions to complex genetic diseases, such as PTB.

Short and tall stature: a new paradigm emerges

In the past, the growth hormone (GH)-insulin-like growth factor 1 (IGF-1) axis was often considered to be the main system that regulated childhood growth and, therefore, determined short stature and tall stature. However, findings have now revealed that the GH-IGF-1 axis is just one of many regulatory systems that control chondrogenesis in the growth plate, which is the biological process that drives height gain. Consequently, normal growth in children depends not only on GH and IGF-1 but also on multiple hormones, paracrine factors, extracellular matrix molecules and intracellular proteins that regulate the activity of growth plate chondrocytes. Mutations in the genes that encode many of these local proteins cause short stature or tall stature. Similarly, genome-wide association studies have revealed that the normal variation in height seems to be largely due to genes outside the GH-IGF-1 axis that affect growth at the growth plate through a wide variety of mechanisms. These findings point to a new conceptual framework for understanding short and tall stature that is centred not on two particular hormones but rather on the growth plate, which is the structure responsible for height gain.

Regulation of body growth

PURPOSE OF REVIEW

Recent basic studies have yielded important new insights into the molecular mechanisms that regulate growth locally. Simultaneously, clinical studies have identified new molecular defects that cause growth failure and overgrowth, and genome-wide association studies have elucidated the genetic basis for normal human height variation.

RECENT FINDINGS

The Hippo pathway has emerged as one of the major mechanisms controlling organ size. In addition, an extensive genetic program has been described that allows rapid body growth in the fetus and infant but then causes growth to slow with age in multiple tissues. In human genome-wide association studies, hundreds of loci associated with adult stature have been identified; many appear to involve genes that function locally in the growth plate. Clinical genetic studies have identified a new genetic abnormality, microduplication of Xq26.3, that is responsible for growth hormone excess, and a gene, DNMT3A, in which mutations cause an overgrowth syndrome through epigenetic mechanisms.

SUMMARY

These recent advances in our understanding of somatic growth not only provide insight into childhood growth disorders but also have broader medical applications because disruption of these regulatory systems contributes to oncogenesis.

Chondrocyte hypertrophy in skeletal development, growth, and disease

Most of our bones form through the process of endochondral ossification, which is tightly regulated by the activity of the cartilage growth plate. Chondrocyte maturation through the various stages of growth plate physiology ultimately results in hypertrophy. Chondrocyte hypertrophy is an essential contributor to longitudinal bone growth, but recent data suggest that these cells also play fundamental roles in signaling to other skeletal cells, thus coordinating endochondral ossification. On the other hand, ectopic hypertrophy of articular chondrocytes has been implicated in the pathogenesis of osteoarthritis. Thus, a better understanding of the processes that control chondrocyte hypertrophy in the growth plate as well as in articular cartilage is required for improved management of both skeletal growth disorders and osteoarthritis. This review summarizes recent findings on the regulation of hypertrophic chondrocyte differentiation, the cellular mechanisms involved in hypertrophy, and the role of chondrocyte hypertrophy in skeletal physiology and pathophysiology.

Human monoclonal antibody fragments targeting matrilin-3 in growth plate cartilage

PURPOSE

Many genetic disorders, including chondrodysplasias, and acquired disorders impair growth plate function, resulting in short and sometimes malformed bones. There are multiple endocrine and paracrine factors that promote chondrogenesis at the growth plate, which could potentially be used to treat these disorders. Targeting these growth factors specifically to the growth plate might augment the therapeutic skeletal effect while diminishing undesirable effects on non-target tissues.

METHODS

Using yeast display technology, we selected single-chain variable antibody fragments that bound to human and mouse matrilin-3, an extracellular matrix protein specifically expressed in cartilage tissue. The ability of the selected antibody fragments to bind matrilin-3 and to bind cartilage tissue in vitro and in vivo was assessed by ELISA and immunohistochemistry.

RESULTS

We identified antibody fragments that bound matrilin-3 with high affinity and also bound with high tissue specificity to cartilage homogenates and to cartilage structures in mouse embryo sections. When injected intravenously in mice, the antibody fragments specifically homed to cartilage.

CONCLUSIONS

Yeast display successfully selected antibody fragments that are able to target cartilage tissue in vivo. Coupling these antibodies to chondrogenic endocrine and paracrine signaling molecules has the potential to open up new pharmacological approaches to treat childhood skeletal growth disorders.

Biological interpretation of genome-wide association studies using predicted gene functions

The main challenge for gaining biological insights from genetic associations is identifying which genes and pathways explain the associations. Here we present DEPICT, an integrative tool that employs predicted gene functions to systematically prioritize the most likely causal genes at associated loci, highlight enriched pathways and identify tissues/cell types where genes from associated loci are highly expressed. DEPICT is not limited to genes with established functions and prioritizes relevant gene sets for many phenotypes.

Environmental temperature impact on bone and cartilage growth

Environmental temperature can have a surprising impact on extremity growth in homeotherms, but the underlying mechanisms have remained elusive for over a century. Limbs of animals raised at warm ambient temperature are significantly and permanently longer than those of littermates housed at cooler temperature. These remarkably consistent lab results closely resemble the ecogeographical tenet described by Allen's "extremity size rule," that appendage length correlates with temperature and latitude. This phenotypic growth plasticity could have adaptive significance for thermal physiology. Shortened extremities help retain body heat in cold environments by decreasing surface area for potential heat loss. Homeotherms have evolved complex mechanisms to maintain tightly regulated internal temperatures in challenging environments, including "facultative extremity heterothermy" in which limb temperatures can parallel ambient. Environmental modulation of tissue temperature can have direct and immediate consequences on cell proliferation, metabolism, matrix production, and mineralization in cartilage. Temperature can also indirectly influence cartilage growth by modulating circulating levels and delivery routes of essential hormones and paracrine regulators. Using an integrated approach, this article synthesizes classic studies with new data that shed light on the basis and significance of this enigmatic growth phenomenon and its relevance for treating human bone elongation disorders. Discussion centers on the vasculature as a gateway to understanding the complex interconnection between direct (local) and indirect (systemic) mechanisms of temperature-enhanced bone lengthening. Recent advances in imaging modalities that enable the dynamic study of cartilage growth plates in vivo will be key to elucidating fundamental physiological mechanisms of long bone growth regulation.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Pharmacogenomics related to growth disorders

Growth disorders resulting in short stature are caused by a wide range of underlying pathophysiological processes. To improve height many of these conditions are treated with recombinant human growth hormone (rhGH). However, substantial inter-individual variability in growth response both in the short and long-term is recognised. Over the last decade, disease-specific growth prediction models have been developed that the clinician can use to define a child's potential response to rhGH and to optimise starting and maintenance doses of rhGH. These models, however, are not able to predict all the variations in treatment response. There has, therefore, been recent interest in using genetic information to contribute to the evaluation of responses to rhGH, including high-throughput technologies for assessing DNA markers (genome) and mRNA transcripts (transcriptome) as pharmacogenomic tools. This review will focus on how these pharmacogenomic approaches are being applied to growth disorders.

Recent research on the growth plate: Recent insights into the regulation of the growth plate

For most bones, elongation is driven primarily by chondrogenesis at the growth plates. This process results from chondrocyte proliferation, hypertrophy, and extracellular matrix secretion, and it is carefully orchestrated by complex networks of local paracrine factors and modulated by endocrine factors. We review here recent advances in the understanding of growth plate physiology. These advances include new approaches to study expression patterns of large numbers of genes in the growth plate, using microdissection followed by microarray. This approach has been combined with genome-wide association studies to provide insights into the regulation of the human growth plate. We also review recent studies elucidating the roles of bone morphogenetic proteins, fibroblast growth factors, C-type natriuretic peptide, and suppressor of cytokine signaling in the local regulation of growth plate chondrogenesis and longitudinal bone growth.

Short stature, accelerated bone maturation, and early growth cessation due to heterozygous aggrecan mutations

CONTEXT

Many children with idiopathic short stature have a delayed bone age. Idiopathic short stature with advanced bone age is far less common.

OBJECTIVE

The aim was to identify underlying genetic causes of short stature with advanced bone age.

SETTING AND DESIGN

We used whole-exome sequencing to study three families with autosomal-dominant short stature, advanced bone age, and premature growth cessation.

RESULTS

Affected individuals presented with short stature [adult heights -2.3 to -4.2 standard deviation scores (SDS)] with histories of early growth cessation or childhood short stature (height SDS -1.9 to -3.5 SDS), advancement of bone age, and normal endocrine evaluations. Whole-exome sequencing identified novel heterozygous variants in ACAN, which encodes aggrecan, a proteoglycan in the extracellular matrix of growth plate and other cartilaginous tissues. The variants were present in all affected, but in no unaffected, family members. In Family 1, a novel frameshift mutation in exon 3 (c.272delA) was identified, which is predicted to cause early truncation of the aggrecan protein. In Family 2, a base-pair substitution was found in a highly conserved location within a splice donor site (c.2026+1G>A), which is also likely to alter the amino acid sequence of a large portion of the protein. In Family 3, a missense variant (c.7064T>C) in exon 14 affects a highly conserved residue (L2355P) and is strongly predicted to perturb protein function.

CONCLUSIONS

Our study demonstrates that heterozygous mutations in ACAN can cause a mild skeletal dysplasia, which presents clinically as short stature with advanced bone age. The accelerating effect on skeletal maturation has not previously been noted in the few prior reports of human ACAN mutations. Our findings thus expand the spectrum of ACAN defects and provide a new molecular genetic etiology for the unusual child who presents with short stature and accelerated skeletal maturation.

Recent research on the growth plate: Advances in fibroblast growth factor signaling in growth plate development and disorders

Skeletons are formed through two distinct developmental actions, intramembranous ossification and endochondral ossification. During embryonic development, most bone is formed by endochondral ossification. The growth plate is the developmental center for endochondral ossification. Multiple signaling pathways participate in the regulation of endochondral ossification. Fibroblast growth factor (FGF)/FGF receptor (FGFR) signaling has been found to play a vital role in the development and maintenance of growth plates. Missense mutations in FGFs and FGFRs can cause multiple genetic skeletal diseases with disordered endochondral ossification. Clarifying the molecular mechanisms of FGFs/FGFRs signaling in skeletal development and genetic skeletal diseases will have implications for the development of therapies for FGF-signaling-related skeletal dysplasias and growth plate injuries. In this review, we summarize the recent advances in elucidating the role of FGFs/FGFRs signaling in growth plate development, genetic skeletal disorders, and the promising therapies for those genetic skeletal diseases resulting from FGFs/FGFRs dysfunction. Finally, we also examine the potential important research in this field in the future.

Genome-wide association study of sexual maturation in males and females highlights a role for body mass and menarche loci in male puberty

Little is known about genes regulating male puberty. Further, while many identified pubertal timing variants associate with age at menarche, a late manifestation of puberty, and body mass, little is known about these variants' relationship to pubertal initiation or tempo. To address these questions, we performed genome-wide association meta-analysis in over 11 000 European samples with data on early pubertal traits, male genital and female breast development, measured by the Tanner scale. We report the first genome-wide significant locus for male sexual development upstream of myocardin-like 2 (MKL2) (P = 8.9 × 10(-9)), a menarche locus tagging a developmental pathway linking earlier puberty with reduced pubertal growth (P = 4.6 × 10(-5)) and short adult stature (p = 7.5 × 10(-6)) in both males and females. Furthermore, our results indicate that a proportion of menarche loci are important for pubertal initiation in both sexes. Consistent with epidemiological correlations between increased prepubertal body mass and earlier pubertal timing in girls, body mass index (BMI)-increasing alleles correlated with earlier breast development. In boys, some BMI-increasing alleles associated with earlier, and others with delayed, sexual development; these genetic results mimic the controversy in epidemiological studies, some of which show opposing correlations between prepubertal BMI and male puberty. Our results contribute to our understanding of the pubertal initiation program in both sexes and indicate that although mechanisms regulating pubertal onset in males and females may largely be shared, the relationship between body mass and pubertal timing in boys may be complex and requires further genetic studies.

Joint analysis of functional genomic data and genome-wide association studies of 18 human traits

Annotations of gene structures and regulatory elements can inform genome-wide association studies (GWASs). However, choosing the relevant annotations for interpreting an association study of a given trait remains challenging. I describe a statistical model that uses association statistics computed across the genome to identify classes of genomic elements that are enriched with or depleted of loci influencing a trait. The model naturally incorporates multiple types of annotations. I applied the model to GWASs of 18 human traits, including red blood cell traits, platelet traits, glucose levels, lipid levels, height, body mass index, and Crohn disease. For each trait, I used the model to evaluate the relevance of 450 different genomic annotations, including protein-coding genes, enhancers, and DNase-I hypersensitive sites in over 100 tissues and cell lines. The fraction of phenotype-associated SNPs influencing protein sequence ranged from around 2% (for platelet volume) up to around 20% (for low-density lipoprotein cholesterol), repressed chromatin was significantly depleted for SNPs associated with several traits, and cell-type-specific DNase-I hypersensitive sites were enriched with SNPs associated with several traits (for example, the spleen in platelet volume). Finally, reweighting each GWAS by using information from functional genomics increased the number of loci with high-confidence associations by around 5%.

Copy number variants in short children born small for gestational age

BACKGROUND/AIMS

In addition to genome-wide association studies (GWAS), height-associated genes may be uncovered by studying individuals with extreme short or tall stature.

METHODS

Genome-wide analysis for copy number variants (CNVs), using single nucleotide polymorphism (SNP) arrays, was performed in 49 index cases born small for gestational age with persistent short stature. Segregation analysis was performed, and genes in CNVs were compared with information from GWAS, gene expression in rodents' growth plates, and published information.

RESULTS

CNVs were detected in 13 cases. In 5 children a known cause of short stature was found: UPD7, UPD14, a duplication of the SHOX enhancer region, an IGF1R deletion, and a 22q11.21 deletion. In the remaining 8 cases, potential pathogenic CNVs were detected, either de novo (n = 1), segregating (n = 2), or not segregating with short stature (n = 5). Bioinformatic analysis of the de novo and segregating CNVs suggested that HOXD4, AGPS, PDE11A, OSBPL6, PRKRA and PLEKHA3, and possibly DGKB and TNFRSF11B are potential candidate genes. A SERPINA7 or NRK defect may be associated with an X-linked form of short stature.

CONCLUSION

SNP arrays detected 5 known causes of short stature with prenatal onset and suggested several potential candidate genes.

Towards identification of molecular mechanisms of short stature

Growth evaluations are among the most common referrals to pediatric endocrinologists. Although a number of pathologies, both primary endocrine and non-endocrine, can present with short stature, an estimated 80% of evaluations fail to identify a clear etiology, leaving a default designation of idiopathic short stature (ISS). As a group, several features among children with ISS are suggestive of pathophysiology of the GH–IGF-1 axis, including low serum levels of IGF-1 despite normal GH secretion. Candidate gene analysis of rare cases has demonstrated that severe mutations of genes of the GH–IGF-1 axis can present with a profound height phenotype, leading to speculation that a collection of mild mutations or polymorphisms of these genes can explain poor growth in a larger proportion of patients. Recent genome-wide association studies have identified ~180 genomic loci associated with height that together account for approximately 10% of height variation. With only modest representation of the GH–IGF-1 axis, there is little support for the long-held hypothesis that common genetic variants of the hormone pathway provide the molecular mechanism for poor growth in a substantial proportion of individuals. The height-associated common variants are not observed in the anticipated frequency in the shortest individuals, suggesting rare genetic factors with large effect are more plausible in this group. As we advance towards establishing a molecular mechanism for poor growth in a greater percentage of those currently labeled ISS, we highlight two strategies that will likely be offered with increasing frequency: (1) unbiased genetic technologies including array analysis for copy number variation and whole exome/genome sequencing and (2) epigenetic alterations of key genomic loci. Ultimately data from subsets with similar molecular etiologies may emerge that will allow tailored interventions to achieve the best clinical outcome.