Genetic influence on within-person longitudinal change in anthropometric traits in the UK Biobank.

The causes of temporal fluctuations in adult traits are poorly understood. Here, we investigate the genetic determinants of within-person trait variability of 8 repeatedly measured anthropometric traits in 50,117 individuals from the UK Biobank. We found that within-person (non-directional) variability had a SNP-based heritability of 2-5% for height, sitting height, body mass index (BMI) and weight (P ≤ 2.4 × 10-3). We also analysed longitudinal trait change and show a loss of both average height and weight beyond about 70 years of age. A variant tracking the Alzheimer's risk APOE- E 4 allele (rs429358) was significantly associated with weight loss ( ß = -0.047 kg per yr, s.e. 0.007, P = 2.2 × 10-11), and using 2-sample Mendelian Randomisation we detected a relationship consistent with causality between decreased lumbar spine bone mineral density and height loss (bxy = 0.011, s.e. 0.003, P = 3.5 × 10-4). Finally, population-level variance quantitative trait loci (vQTL) were consistent with within-person variability for several traits, indicating an overlap between trait variability assessed at the population or individual level. Our findings help elucidate the genetic influence on trait-change within an individual and highlight disease risks associated with these changes.

Y chromosome haplogroups are associated with birth size in Japanese men.

The Y chromosome is classified into haplogroups (A-T) based on a combination of several DNA polymorphisms. Japanese men are mainly classified into haplogroups C, D, and O, which have been further subdivided. The distribution of Y-chromosome haplogroups varies by ethnicity. The phylogenetic age, origin, and migration also differ. I hypothesized that Y chromosome haplogroups may be associated with height and/or weight at birth. An association analysis of height and weight at birth with Y chromosome haplogroups was performed in 288 Japanese men. Men belonging to haplogroup O1b2 were significantly associated with short stature at birth (beta = －1.88, standard error (SE) = 0.55, P = 0.00076), and those belonging to D1a2a-12f2b were significantly associated with increased birth weight (beta = 174, SE = 64, P = 0.0069). Y chromosome haplogroups are associated with physical birth characteristics in modern Japanese men. J. Med. Invest. 71 : 129-133, February, 2024.

Differences in Childhood Growth Parameters Between Patients With Somatic and Heritable Retinoblastoma.

Purpose: Little is known regarding differences in childhood growth between somatic and heritable retinoblastoma (Rb) populations. We aimed to compare childhood growth parameters between somatic and heritable Rb cohorts at birth and at time of diagnosis with Rb. Methods: A multinational, longitudinal cohort study was conducted with patients from 11 centers in 10 countries who presented with treatment naïve Rb from January to December 2019. Variables of interest included age, sex, and size characteristics at birth and at time of presentation, as well as germline mutation status. After Bonferroni correction, results were statistically significant if the P value was less than 0.005. Results: We enrolled 696 patients, with 253 analyzed after exclusion criteria applied. Between somatic (n = 39) and heritable (n = 214) Rb cohorts, with males and females analyzed separately, there was no significant difference in birth weight percentile, weight percentile at time of diagnosis, length percentile at time of diagnosis, weight-for-length percentile at time of diagnosis, or change of weight percentile from birth to time of diagnosis. Patients with heritable Rb had a smaller mean weight percentile at birth and smaller mean weight and length percentiles at time of diagnosis with Rb, although this difference was not statistically significant. All cohorts experienced a slight negative change of weight percentile from birth to time of diagnosis. No cohort mean percentiles met criteria for failure to thrive, defined as less than the 5th percentile. Conclusions: Children with Rb seem to have normal birth and childhood growth patterns. There is no definitive evidence that somatic or heritable Rb has a biological or environmental impact on childhood growth parameters.

Correlation-based tests for the formal comparison of polygenic scores in multiple populations.

Polygenic scores (PGS) are measures of genetic risk, derived from the results of genome wide association studies (GWAS). Previous work has proposed the coefficient of determination (R2) as an appropriate measure by which to compare PGS performance in a validation dataset. Here we propose correlation-based methods for evaluating PGS performance by adapting previous work which produced a statistical framework and robust test statistics for the comparison of multiple correlation measures in multiple populations. This flexible framework can be extended to a wider variety of hypothesis tests than currently available methods. We assess our proposed method in simulation and demonstrate its utility with two examples, assessing previously developed PGS for low-density lipoprotein cholesterol and height in multiple populations in the All of Us cohort. Finally, we provide an R package 'coranova' with both parametric and nonparametric implementations of the described methods.

The association between DNA methylation and human height and a prospective model of DNA methylation-based height prediction.

As a vital anthropometric characteristic, human height information not only helps to understand overall developmental status and genetic risk factors, but is also important for forensic DNA phenotyping. We utilized linear regression analysis to test the association between each CpG probe and the height phenotype. Next, we designed a methylation sequencing panel targeting 959 CpGs and subsequent height inference models were constructed for the Chinese population. A total of 11,730 height-associated sites were identified. By employing KPCA and deep neural networks, a prediction model was developed, of which the cross-validation RMSE, MAE and R2 were 5.62 cm, 4.45 cm and 0.64, respectively. Genetic factors could explain 39.4% of the methylation level variance of sites used in the height inference models. Collectively, we demonstrated an association between height and DNA methylation status through an EWAS analysis. Targeted methylation sequencing of only 959 CpGs combined with deep learning techniques could provide a model to estimate human height with higher accuracy than SNP-based prediction models.

Novel signals and polygenic score for height are associated with pubertal growth traits in Southwestern American Indians.

Most genetic variants associated with adult height have been identified through large genome-wide association studies (GWASs) in European-ancestry cohorts. However, it is unclear how these variants influence linear growth during adolescence. This study uses anthropometric and genotypic data from a longitudinal study conducted in an American Indian community in Arizona between 1965-2007. Growth parameters (i.e. height, velocity, and timing of growth spurt) were derived from the Preece-Baines growth model, a parametric growth curve fitted to longitudinal height data, in 787 participants with height measurements spanning the whole period of growth. Heritability estimates suggested that genetic factors could explain 25% to 71% of the variance of pubertal growth traits. We performed a GWAS of growth parameters, testing their associations with 5 077 595 imputed or directly genotyped variants. Six variants associated with height at peak velocity (P < 5 × 10-8, adjusted for sex, birth year and principal components). Implicated genes include NUDT3, previously associated with adult height, and PACSIN1. Two novel variants associated with duration of growth spurt (P < 5 × 10-8) in LOC105375344, an uncharacterized gene with unknown function. We finally examined the association of growth parameters with a polygenic score for height derived from 9557 single nucleotide polymorphisms (SNPs) identified in the GIANT meta-analysis for which genotypic data were available for the American Indian study population. Height polygenic score was correlated with the magnitude and velocity of height growth that occurred before and at the peak of the adolescent growth spurt, indicating overlapping genetic architecture, with no influence on the timing of adolescent growth.

Trans-ancestral genome-wide association study of longitudinal pubertal height growth and shared heritability with adult health outcomes.

BACKGROUND: Pubertal growth patterns correlate with future health outcomes. However, the genetic mechanisms mediating growth trajectories remain largely unknown. Here, we modeled longitudinal height growth with Super-Imposition by Translation And Rotation (SITAR) growth curve analysis on ~ 56,000 trans-ancestry samples with repeated height measurements from age 5 years to adulthood. We performed genetic analysis on six phenotypes representing the magnitude, timing, and intensity of the pubertal growth spurt. To investigate the lifelong impact of genetic variants associated with pubertal growth trajectories, we performed genetic correlation analyses and phenome-wide association studies in the Penn Medicine BioBank and the UK Biobank. RESULTS: Large-scale growth modeling enables an unprecedented view of adolescent growth across contemporary and 20th-century pediatric cohorts. We identify 26 genome-wide significant loci and leverage trans-ancestry data to perform fine-mapping. Our data reveals genetic relationships between pediatric height growth and health across the life course, with different growth trajectories correlated with different outcomes. For instance, a faster tempo of pubertal growth correlates with higher bone mineral density, HOMA-IR, fasting insulin, type 2 diabetes, and lung cancer, whereas being taller at early puberty, taller across puberty, and having quicker pubertal growth were associated with higher risk for atrial fibrillation. CONCLUSION: We report novel genetic associations with the tempo of pubertal growth and find that genetic determinants of growth are correlated with reproductive, glycemic, respiratory, and cardiac traits in adulthood. These results aid in identifying specific growth trajectories impacting lifelong health and show that there may not be a single "optimal" pubertal growth pattern.

Molecular diagnosis is an important indicator for response to growth hormone therapy in children with short stature.

BACKGROUND: Significant differences have been observed in the efficacy of recombinant human growth hormone (rhGH) treatment for short children. The present study aimed to identify the genetic etiology of short stature and to assess the role of molecular diagnosis in predicting responses to rhGH treatment. METHODS: A total of 407 short children were included in the present study, 226 of whom received rhGH treatment. Whole-exome sequencing (WES) was conducted on short children to identify the underlying genetic etiology. Correlations between molecular diagnosis and the efficacy of rhGH treatment were examined. RESULTS: Pathogenic or likely pathogenic mutations were identified in 86 of the 407 patients (21.1%), including 36 (41.9%) novel variants. Among the multiple pathways affecting short stature, genes involved in fundamental cellular processes (38.7%) play a larger role, especially the RAS-MAPK pathway. In general, patients without pathogenic mutations responded better to rhGH than those with mutations. Furthermore, patients with hormone signaling pathway mutations had a better response to rhGH, while those with paracrine factor mutations had a worse response to rhGH. CONCLUSIONS: This study highlights the utility of WES in identifying genetic etiology in children with short stature. Identifying likely causal mutations is an important factor in predicting rhGH response.

The Genetic Landscape of Children Born Small for Gestational Age with Persistent Short Stature.

INTRODUCTION: Among children born small for gestational age, 10-15% fail to catch up and remain short (SGA-SS). The underlying mechanisms are mostly unknown. We aimed to decipher genetic aetiologies of SGA-SS within a large single-centre cohort. METHODS: Out of 820 patients treated with growth hormone (GH), 256 were classified as SGA-SS (birth length and/or birth weight <-2 SD for gestational age and life-minimum height <-2.5 SD). Those with the DNA triplet available (child and both parents) were included in the study (176/256). Targeted testing (karyotype/FISH/MLPA/specific Sanger sequencing) was performed if a specific genetic disorder was clinically suggestive. All remaining patients underwent MS-MLPA to identify Silver-Russell syndrome, and those with unknown genetic aetiology were subsequently examined using whole-exome sequencing or targeted panel of 398 growth-related genes. Genetic variants were classified using ACMG guidelines. RESULTS: The genetic aetiology was elucidated in 74/176 (42%) children. Of these, 12/74 (16%) had pathogenic or likely pathogenic (P/LP) gene variants affecting pituitary development (LHX4, OTX2, PROKR2, PTCH1, POU1F1), the GH-IGF-1 or IGF-2 axis (GHSR, IGFALS, IGF1R, STAT3, HMGA2), 2/74 (3%) the thyroid axis (TRHR, THRA), 17/74 (23%) the cartilaginous matrix (ACAN, various collagens, FLNB, MATN3), and 7/74 (9%) the paracrine chondrocyte regulation (FGFR3, FGFR2, NPR2). In 12/74 (16%), we revealed P/LP affecting fundamental intracellular/intranuclear processes (CDC42, KMT2D, LMNA, NSD1, PTPN11, SRCAP, SON, SOS1, SOX9, TLK2). SHOX deficiency was found in 7/74 (9%), Silver-Russell syndrome in 12/74 (16%) (11p15, UPD7), and miscellaneous chromosomal aberrations in 5/74 (7%) children. CONCLUSIONS: The high diagnostic yield sheds a new light on the genetic landscape of SGA-SS, with a central role for the growth plate with substantial contributions from the GH-IGF-1 and thyroid axes and intracellular regulation and signalling.

Growth charts in DYRK1A syndrome.

DYRK1A Syndrome (OMIM #614104) is caused by pathogenic variations in the DYRK1A gene located on 21q22. Haploinsufficiency of DYRK1A causes a syndrome with global psychomotor delay and intellectual disability. Low birth weight, growth restriction with feeding difficulties, stature insufficiency, and microcephaly are frequently reported. This study aims to create specific growth charts for individuals with DYRK1A Syndrome and identify parameters for size prognosis. Growth parameters were obtained for 92 individuals with DYRK1A Syndrome (49 males vs. 43 females). The data were obtained from pediatric records, parent reporting, and scientific literature. Growth charts for height, weight, body mass index (BMI), and occipitofrontal circumference (OFC) were generated using generalized additive models through R package gamlss. The growth curves include height, weight, and OFC measurements for patients aged 0-5 years. In accordance with the literature, the charts show that individuals are more likely to present intrauterine growth restriction with low birth weight and microcephaly. The growth is then characterized by severe microcephaly, low weight, and short stature. This study proposes growth charts for widespread use in the management of patients with DYRK1A syndrome.

Growth reference charts for children with hypochondroplasia.

Hypochondroplasia (HCH) is a rare skeletal dysplasia causing mild short stature. There is a paucity of growth reference charts for this population. Anthropometric data were collected to generate height, weight, and head circumference (HC) growth reference charts for children with a diagnosis of HCH. Mixed longitudinal anthropometric data and genetic analysis results were collected from 14 European specialized skeletal dysplasia centers. Growth charts were generated using Generalized Additive Models for Location, Scale, and Shape. Measurements for height (983), weight (896), and HC (389) were collected from 188 (79 female) children with a diagnosis of HCH aged 0-18 years. Of the 84 children who underwent genetic testing, a pathogenic variant in FGFR3 was identified in 92% (77). The data were used to generate growth references for height, weight, and HC, plotted as charts with seven centiles from 2nd to 98th, for ages 0-4 and 0-16 years. HCH-specific growth charts are important in the clinical care of these children. They help to identify if other comorbidities are present that affect growth and development and serve as an important benchmark for any prospective interventional research studies and trials.

Height of non-Hispanic white adults with homeostatic iron regulator HFE genotypes p.C282Y/p.C282Y and wt/wt.

BACKGROUND: We sought to evaluate height in white adults with hemochromatosis. METHODS: We analyzed the height of (1) post-screening examination participants with HFE p.C282Y/p.C282Y (rs1800562) and wt/wt (absence of p.C282Y and p.H63D (rs1799945)) and (2) referred hemochromatosis probands with p.C282Y/p.C282Y. RESULTS: There were 762 participants (270 p.C282Y/p.C282Y, 492 wt/wt; 343 men, 419 women) and 180 probands (104 men, 76 women). Median height of male participants with p.C282Y/p.C282Y or wt/wt was 177.8 cm. Median height of female participants was greater in those with p.C282Y/p.C282Y than wt/wt (165.1 cm vs 162.6 cm, respectively; p = 0.0298). Median height of p.C282Y/p.C282Y participants and probands was the same (men 177.8 cm; women 165.1 cm). Regressions on height of male and female participants revealed no associations with HFE genotype and inverse and positive associations with age and weight, respectively. Height of female participants was positively and inversely associated with transferrin saturation and serum ferritin, respectively. Regressions on height of male and female probands revealed positive associations with weight. CONCLUSIONS: The height of men with HFE p.C282Y/p.C282Y and wt/wt does not differ significantly. The height of female participants was greater in those with p.C282Y/p.C282Y than wt/wt. We found no independent association of HFE genotype with the height of men or women.

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.

Novel heterozygous mutation in the SHOX gene leading to familial idiopathic short stature: A case report and literature review.

BACKGROUND: The pathogenic mutation of short stature homeobox (SHOX) gene is one of the main genetic causes of short stature in children, with an incidence rate of 1/1000~1/2000 and the main clinical manifestations are short stature and (or) limb skeletal abnormalities. SHOX gene mutations are mostly large deletions of regulatory sequence genes, while exon mutations are relatively rare. The pathogenic rate of mutations occurring in exon 5 is only 1/50 000~1/100 000. This study reviewed the clinical data of a child with SHOX gene mutation in exon 5, and analyzed the clinical phenotype, pathogenesis, diagnosis, treatment and prognosis of SHOX gene mutation in combination with relevant literature at home and abroad. CASE PRESENTATION: The patient was an 8-year-old girl with a height of 105.2 cm (-4.31 standard deviations). Her sitting height/height ratio was 56.8% (>55.5%), and she exhibited high-arched palate, irregular dentition, micrognathia, short fingers, and a normal growth hormone stimulation test. Whole-exome sequencing was performed, and Sanger sequencing was used for site validation. The sequencing results revealed a heterozygous mutation of c.577G > A in exon 5 of the SHOX gene, inherited from the father. The clinical symptoms of the proband were consistent with the phenotype of short stature idiopathic familial associated with SHOX gene mutations. The father, grandfather, uncle, and sister of the proband all had the c.577G > A heterozygous mutation. Therefore, the clinical diagnosis was childhood short stature caused by SHOX gene defects. The SHOX: c.577G > A mutation is likely to be the genetic etiology of familial idiopathic short stature in this family, and this novel mutation enriches the mutation spectrum of the SHOX gene. CONCLUSION: This is the first case report of familial idiopathic dwarfism caused by mutation at the c.577G > A locus of exon 5 of SHOX gene in the world. This novel mutation enriches the mutation spectrum of the SHOX gene. It is important to emphasize genetic testing, including the SHOX gene, in patients with familial idiopathic short stature and to provide timely growth hormone therapy to individuals with short stature caused by SHOX gene mutations in order to improve their adult height.

Mexican Biobank advances population and medical genomics of diverse ancestries.

Latin America continues to be severely underrepresented in genomics research, and fine-scale genetic histories and complex trait architectures remain hidden owing to insufficient data1. To fill this gap, the Mexican Biobank project genotyped 6,057 individuals from 898 rural and urban localities across all 32 states in Mexico at a resolution of 1.8 million genome-wide markers with linked complex trait and disease information creating a valuable nationwide genotype-phenotype database. Here, using ancestry deconvolution and inference of identity-by-descent segments, we inferred ancestral population sizes across Mesoamerican regions over time, unravelling Indigenous, colonial and postcolonial demographic dynamics2-6. We observed variation in runs of homozygosity among genomic regions with different ancestries reflecting distinct demographic histories and, in turn, different distributions of rare deleterious variants. We conducted genome-wide association studies (GWAS) for 22 complex traits and found that several traits are better predicted using the Mexican Biobank GWAS compared to the UK Biobank GWAS7,8. We identified genetic and environmental factors associating with trait variation, such as the length of the genome in runs of homozygosity as a predictor for body mass index, triglycerides, glucose and height. This study provides insights into the genetic histories of individuals in Mexico and dissects their complex trait architectures, both crucial for making precision and preventive medicine initiatives accessible worldwide.

Longitudinal assessment of bone health index as a measure of bone health in short-statured children before and during treatment with recombinant growth hormone.

OBJECTIVES: The aim of our study was the longitudinal assessment of bone health index (BHI) in short-statured children during growth hormone (GH) treatment to estimate changes in their bone health. METHODS: 256 short-statured children (isolated GH deficiency (IGHD) n=121, multiple pituitary hormone deficiency (MPHD) n=49, intrauterine growth retardation (small for gestational age (SGA)) n=52, SHOX (short stature homeobox gene) deficiency n=9, Ullrich Turner syndrome (UTS) n=25) who started with GH between 2010 and 2018 were included. Annual bone ages (Greulich and Pyle, GP) and BHI were, retrospectively, analysed in consecutive radiographs of the left hand (BoneXpert software) from GH therapy start (T0) up to 10 years (T10) thereafter, with T max indicating the individual time point of the last available radiograph. The results are presented as the median (25 %/75 % interquartile ranges, IQR) and statistical analyses were performed using non-parametric tests as appropriate. RESULTS: The BHI standard deviation scores (SDS) were reduced (-0.97, -1.8/-0.3) as bone ages were retarded (-1.6 years, -2.31/-0.97) in all patients before start of GH and were significantly lower in patients with growth hormone deficiency (GHD) (-1.04, -1.85/-0.56; n=170) compared to non-GHD patients (-0.79, -1.56/-0.01; n=86; p=0.022). BHI SDS increased to -0.17 (-1/0.58) after 1 year of GH (T1, 0.5-1.49, p<0.001) and to -0.20 (-1/-0.50, p<0.001) after 5.3 years (T max, 3.45/7.25). CONCLUSIONS: BHI SDS are reduced in treatment-naive short-statured children regardless of their GH status, increase initially with GH treatment while plateauing thereafter, suggesting sustained improved bone health.

Nonlinear Catch-Up Growth in Height, Weight, and Head Circumference from Birth to Adolescence: A Longitudinal Twin Study.

Owing to high rates of prenatal complications, twins are, on average, substantially smaller than population norms on physical measurements including height, weight, and head circumference at birth. By early childhood, twins are physically average. This study is the first to explore the process of catch-up growth by fitting asymptotic growth models to age-standardized height, weight, and head circumference measurements in a community sample of twins (n = 1281, 52.3% female) followed at up to 17 time points from birth to 15 years. Catch-up growth was rapid over the first year and plateaued around the population mean by early childhood. Shared environmental factors accounted for the majority of individual differences in initial physical size (57.7-65.5%), whereas additive genetic factors accounted for the majority of individual differences in the upper asymptotes of height, weight, and head circumference (73.4-92.6%). Both additive genetic and shared environmental factors were associated with variance in how quickly twins caught up. Gestational age and family SES emerged as important environmental correlates of physical catch-up growth.