Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity

Causal associations between childhood obesity and delayed puberty or height: a bidirectional two-sample Mendelian randomization study

OBJECTIVES

Childhood obesity is thought to influence pubertal development, according to observational studies. However, the exact causal relationship remains unclear due to the complexity of factors affecting pubertal development.

METHODS

To explore the association between exposure (childhood obesity) and outcome (delayed puberty, height), we utilized various methods, including inverse-variance weighted (IVW), weighted median, weighted mode, and MR Egger regression. Additionally, sensitivity analyses were conducted using MR-Egger, MR-PRESSO, Cochran's Q, and leave-one-out techniques to ensure the robustness of the results. Additionally, reverse MR analysis was conducted to explore potential reverse causation.

RESULTS

The IVW analysis revealed no significant genetic causal link between childhood obesity and delayed puberty or height (all p>0.05). In the reverse analysis, height had a causal association with childhood obesity (OR=0.85, 95 % CI=0.76-0.96). The Cochran's Q test highlighted heterogeneity in the results concerning childhood obesity and height (p<0.05). But the MR-Egger intercept and MR-PRESSO test confirmed no impact the results pleiotropic bias, supported by leave-one-out sensitivity analysis.

CONCLUSIONS

Our study found no significant genetic causal association between childhood obesity and delayed puberty or height. However, height was causally associated with childhood obesity. Future research should utilize advanced analytical methods to better understand the determinants of pubertal development.

The hippo pathway: a molecular bridge between environmental cues and pace of life

The pace of life (POL) is shaped by a complex interplay between genetic and environmental factors, influencing growth, maturation, and lifespan across species. The Hippo signaling pathway, a key regulator of organ size and cellular homeostasis, has emerged as a central integrator of environmental cues that modulate POL traits. In this review, we explore how the Hippo pathway links environmental factors-such as temperature fluctuations and dietary energy availability-to molecular mechanisms governing metabolic balance, hormonal signaling, and reproductive timing. Specifically, we highlight the regulatory interactions between the Hippo pathway and metabolic sensors (AMPK, mTOR, SIRT1 and DLK1-Notch), as well as hormonal signals (IGF-1, kisspeptin, leptin, cortisol, thyroid and sex steroids), which together orchestrate key life-history traits, including growth rates, lifespan and sexual maturation, with a particular emphasis on their role in reproductive timing. Furthermore, we consider its role as a potential coordinator of POL-related molecular processes, such as telomere dynamics and epigenetic mechanisms, within a broader regulatory network. By integrating insights from molecular biology and eco-evolutionary perspectives, we propose future directions to dissect the Hippo pathway's role in POL regulation across taxa. Understanding these interactions will provide new perspectives on how organisms adaptively adjust life-history strategies in response to environmental variability.

Effects of childhood and adult height on later life cardiovascular disease risk estimated through Mendelian randomization

Taller individuals are at elevated and protected risk of various cardiovascular disease endpoints. Whether this is due to a direct consequence of their height during childhood, a long-term effect of remaining tall throughout the lifecourse, or confounding by other factors, is unknown. We sought to address this by harnessing human genetic data from the UK Biobank to separate the independent effects of childhood and adulthood height using an approach known as lifecourse Mendelian randomization (MR). Protective effects of taller childhood height on risk of later life coronary artery disease (OR = 0.78 per change in height category, 95% CI = 0.70 to 0.86, P = 4 × 10- 10) and stroke (OR = 0.93, 95% CI = 0.86 to 1.00, P = 0.03) using data from large-scale consortia were found using a univariable model, although evidence of these effects attenuated in a multivariable setting upon accounting for adulthood height. In contrast, direct effects of taller childhood height on increased risk of later life atrial fibrillation (OR = 1.61, 95% CI = 1.42 to 1.79, P = 5 × 10- 7) and thoracic aortic aneurysm (OR = 1.55, 95% CI = 1.16 to 1.95, P = 0.03) were found even after accounting for adulthood height. Evidence for both of these direct effects was replicated in the Million Veterans Program. The protective effect of childhood height on risk of coronary artery disease and stroke can be largely explained by taller children typically becoming taller individuals in later life. Conversely, the independent effect of childhood height on increased risk of atrial fibrillation and thoracic aortic aneurysm may point towards developmental mechanisms in early life which confer a lifelong risk on these disease outcomes.

Impact of Obesity on Pubertal Timing and Male Fertility

Childhood obesity has profound effects on puberty in boys and girls, altering its timing, progression, and associated hormonal changes. Also, later male fertility could be impaired by childhood and pubertal obesity in light of the impact of inflammatory markers on semen quality. The aim of this narrative review is to explore the intricate relationship between childhood obesity and its impact on pubertal development and fertility, with a specific focus on boys. Such a relationship between obesity and pubertal timing in males is highly influenced by metabolic, hormonal, genetic, epigenetic, and environmental factors. While many studies suggest that obesity accelerates pubertal onset in boys, some studies do not confirm these findings, especially in cases of severe obesity. In fact, delayed puberty has also been reported in certain instances. Obesity influences fertility through different central and peripheral processes, including an altered endocrine milieu, inflammatory environment, and epigenetic modifications that alter semen quality and vitality, leading to subfertility or infertility. The early identification and management of potential issues associated with obesity are crucial for ensuring optimal reproductive health in adulthood. Further research is essential to clarify these associations and to develop targeted interventions aimed at preventing the negative health outcomes associated with obesity-related disruptions in puberty and fertility.

Cross-trait multivariate GWAS confirms health implications of pubertal timing

Pubertal timing is highly variable and is associated with long-term health outcomes. Phenotypes associated with pubertal timing include age at menarche, age at voice break, age at first facial hair and growth spurt, and pubertal timing seems to have a shared genetic architecture between the sexes. However, puberty phenotypes have primarily been assessed separately, failing to account for shared genetics, which limits the reliability of the purported health implications. Here, we model the common genetic architecture for puberty timing using a multivariate GWAS, with an effective population of 514,750 European participants. We find 266 independent variants in 197 loci, including 18 novel variants. Transcriptomic, proteome imputation and fine-mapping analyses reveal genes causal for pubertal timing, including KDM4C, LEPR, CCNC, ACP1, and PCSK1. Linkage disequilibrium score regression and Mendelian randomisation analysis establish causal associations between earlier puberty and both accelerated ageing and the risk of developing cardiovascular disease and osteoporosis. We find that alanine aminotransferase, glycated haemoglobin, high-density lipoprotein cholesterol and Parabacteroides levels are mediators of these relationships, and establish that controlling oily fish and retinol intake may be beneficial for promoting healthy pubertal development.

Transcriptome and DNA methylation analysis of the goat pineal gland during puberty

Previous studies have confirmed that methylation regulates gene transcription in the hypothalamus-pituitary-gonadal axis during puberty initiation, but little is known about the regulation of DNA methylation on gene expression in the pineal gland. To screen pineal gland candidate genes related to the onset of goat puberty and regulated by genome methylation, we collected pineal glands from prepubertal and pubertal female goats, then, determined the DNA methylation profile by whole genome bisulfite sequencing and the transcriptome by RNA sequencing on Illumina HiSeqTM2500. We analyzed differentially expressed genes between the Pre group and Pub group using the DESeq2 software (version 1.20.0), and applied the Benjamini and Hochberg method for adjusting P-values. Genes with a P-value less than 0.05 and an absolute log2 fold change greater than 0 were considered differentially expressed genes. Results showed that there was no significant difference in the whole-genome methylation level of the pineal gland between prepubertal and pubertal goats, but the methylation pattern changed significantly, indicating that genomic DNA methylation of the pineal gland might play a role in regulating the initiation of goat puberty. Changes in DNA methylation patterns affected some pineal gland transcriptomes, while the transcriptional level of most genes remained unaffected by DNA methylation differences. Genes regulated by DNA methylation regulates genes primarily involved in metabolic processes, oxidative phosphorylation, and signaling pathways related to thermogenesis. Methylation significantly regulated the expression of genes such as ATP5F1D, CACNB2, and PTEN, while genes like LIN28B, GIP, OPN1SW, and DCC showed the most notable fold changes, which may indicate their involvement in the onset of puberty.

Temporal ablation of the ciliary protein IFT88 alters normal brainwave patterns

The primary cilium is a hair-like organelle that hosts molecular machinery for various developmental and homeostatic signaling pathways. Its alteration can cause rare ciliopathies such as the Bardet-Biedl and Joubert syndromes, but is also linked to Alzheimer's disease, clinical depression, and autism spectrum disorder. These afflictions are caused by disturbances in a wide variety of genes but a common phenotype amongst them is cognitive impairment. While cilia-mediated neural function has been widely examined in early neurodevelopment, their function in the adult brain is not well understood. To help elucidate the role of cilia in neural activity, we temporally induced the ablation of IFT88, a gene encoding the intraflagellar transport 88 protein which is neccessary for ciliogenesis, in adult mice before performing memory-related behavioral assays and electroencephalogram/electromyogram (EEG/EMG) recordings. Inducible IFT88 KO mice exhibited severe learning deficits in trace fear conditioning and Morris water maze tests. They had strongly affected brainwave activity both under isoflurane induced anesthesia and during normal activity. And additionally, inducible IFT88 KO mice had altered sleep architecture and attenuated phase-amplitude coupling, a process that underlies learning and memory formation. These results highlight the growing significance of primary cilia for healthy neural function in the adult brain.

Height estimation in children and adolescents using body composition big data: Machine-learning and explainable artificial intelligence approach

Objective

To develop an accurate and interpretable height estimation model for children and adolescents using body composition variables and explainable artificial intelligence approaches.

Methods

A light gradient boosting method was employed on a dataset of 278,301 measurements from 54,374 children and adolescents aged 6-18 years. The model incorporated anthropometric and body composition measures. Model interpretability was enhanced through feature importance analysis, Shapley additive explanations, partial dependence plots, and accumulated local effects.

Results

The models achieved high accuracy with mean absolute percentage errors of 1.64% and 1.63% for boys and girls, respectively. Soft lean mass (SLM), body fat mass percentage (BFMP), skeletal muscle mass, and skeletal muscle mass percentage were consistently identified as key factors influencing height estimation. Analysis revealed a positive correlation between SLM and estimated height, while BFMP exhibited an inverse relationship with height projections.

Conclusion

These findings provide valuable insights into the relationship between body composition and height, underlining the potential of body composition variables as accurate height predictors in children and adolescents. The model's interpretability and accuracy make it a promising tool for pediatric growth assessment and monitoring.

Temporal Ablation of the Ciliary Protein IFT88 Alters Normal Brainwave Patterns

The primary cilium is a hair-like organelle that hosts molecular machinery for various developmental and homeostatic signaling pathways. Its alteration can cause rare ciliopathies such as the Bardet-Biedl and Joubert syndromes, but is also linked to Alzheimer's disease, clinical depression, and autism spectrum disorder. These afflictions are caused by disturbances in a wide variety of genes but a common phenotype amongst them is cognitive impairment. While cilia-mediated neural function has been widely examined in early neurodevelopment, their function in the adult brain is not well understood. To help elucidate the role of cilia in neural activity, we temporally induced the ablation of IFT88, a gene encoding the intraflagellar transport 88 protein which is neccessary for ciliogenesis, in adult mice before performing memory-related behavioral assays and electroencephalogram/electromyogram (EEG/EMG) recordings. Inducible IFT88 KO mice exhibited severe learning deficits in trace fear conditioning and Morris water maze tests. They had strongly affected brainwave activity both under isoflurane induced anesthesia and during normal activity. And additionally, inducible IFT88 KO mice had altered sleep architecture and attenuated phase-amplitude coupling, a process that underlies learning and memory formation. These results highlight the growing significance of primary cilia for healthy neural function in the adult brain.

A genome-wide association study of neonatal metabolites

Genetic factors significantly influence the concentration of metabolites in adults. Nevertheless, the genetic influence on neonatal metabolites remains uncertain. To bridge this gap, we employed genotype imputation techniques on large-scale low-pass genome data obtained from non-invasive prenatal testing. Subsequently, we conducted association studies on a total of 75 metabolic components in neonates. The study identified 19 previously reported associations and 11 novel associations between single-nucleotide polymorphisms and metabolic components. These associations were initially found in the discovery cohort (8,744 participants) and subsequently confirmed in a replication cohort (19,041 participants). The average heritability of metabolic components was estimated to be 76.2%, with a range of 69%-78.8%. These findings offer valuable insights into the genetic architecture of neonatal metabolism.

Genetically inferred birthweight, height, and puberty timing and risk of osteosarcoma

INTRODUCTION

Several studies have linked increased risk of osteosarcoma with tall stature, high birthweight, and early puberty, although evidence is inconsistent. We used genetic risk scores (GRS) based on established genetic loci for these traits and evaluated associations between genetically inferred birthweight, height, and puberty timing with osteosarcoma.

METHODS

Using genotype data from two genome-wide association studies, totaling 1039 cases and 2923 controls of European ancestry, association analyses were conducted using logistic regression for each study and meta-analyzed to estimate pooled odds ratios (ORs) and 95% confidence intervals (CIs). Subgroup analyses were conducted by case diagnosis age, metastasis status, tumor location, tumor histology, and presence of a known pathogenic variant in a cancer susceptibility gene.

RESULTS

Genetically inferred higher birthweight was associated with an increased risk of osteosarcoma (OR =1.59, 95% CI 1.07-2.38, P = 0.02). This association was strongest in cases without metastatic disease (OR =2.46, 95% CI 1.44-4.19, P = 9.5 ×10-04). Although there was no overall association between osteosarcoma and genetically inferred taller stature (OR=1.06, 95% CI 0.96-1.17, P = 0.28), the GRS for taller stature was associated with an increased risk of osteosarcoma in 154 cases with a known pathogenic cancer susceptibility gene variant (OR=1.29, 95% CI 1.03-1.63, P = 0.03). There were no significant associations between the GRS for puberty timing and osteosarcoma.

CONCLUSION

A genetic propensity to higher birthweight was associated with increased osteosarcoma risk, suggesting that shared genetic factors or biological pathways that affect birthweight may contribute to osteosarcoma pathogenesis.

Mendelian randomization analysis on the impacts of age at menarche on adult height: A Taiwanese population study

BACKGROUNDS

Ample evidence supports potential influence of age at menarche (AM) on adult height (AH), but multiple confounders may affect causal estimates. To address this issue, the Mendelian randomization (MR) analysis was used to explore the causal impacts of AM on AH.

METHODS

Using data (n = 57,349) from the publicly accessible Taiwan Biobank and randomly splitting them into 2 equal-size subsets, we identified single nucleotide polymorphisms (SNPs) significantly associated with AM in the exploration subset and used these SNPs as instrumental variables to estimate the effects of instruments on AH in the validation subset based on two stage least squares (2SLS) regression. In addition, three more summary statistics-based approaches, namely inverse variance weighted (IVW), MR-Egger, and weighted median (WM) analyses, were used to verify the findings. We also performed heterogeneity and sensitivity analyses to evaluate the robustness of the results.

RESULTS

We identified 4 leading SNPs associated with AM at the genome-wide significant level, whereas rs9409082 may exert some pleiotropic effects on AH. After eliminating rs9409082, the 2SLS analysis indicated that one year delay in genetically determined AM predicted 1.5 cm height gain in adulthood (β = 1.508, 95% confidence interval [0.852, 2.163]). The causal relationship was also supported by WM (β = 1.183, [0.329, 2.038]) and IVW (β = 1.493, [0.523, 2.463]) methods.

CONCLUSIONS

Evidence from the present MR study supports a causal relationship between later AM and taller AH.

The complex web of obesity: from genetics to precision medicine

INTRODUCTION

Obesity is a growing public health concern affecting both children and adults. Since it involves both genetic and environmental components, the management of obesity requires both, an understanding of the underlying genetics and changes in lifestyle. The knowledge of obesity genetics will enable the possibility of precision medicine in anti-obesity medications.

AREAS COVERED

Here, we explore health complications and the prevalence of obesity. We discuss disruptions in energy balance as a symptom of obesity, examining evolutionary theories, its multi-factorial origins, and heritability. Additionally, we discuss monogenic and polygenic obesity, the converging biological pathways, potential pharmacogenomics applications, and existing anti-obesity medications - specifically focussing on the leptin-melanocortin and incretin pathways. Comparisons between childhood and adult obesity genetics are made, along with insights into structural variants, epigenetic changes, and environmental influences on epigenetic signatures.

EXPERT OPINION

With recent advancements in anti-obesity drugs, genetic studies pinpoint new targets and allow for repurposing existing drugs. This creates opportunities for genotype-informed treatment options. Also, lifestyle interventions can help in the prevention and treatment of obesity by altering the epigenetic signatures. The comparison of genetic architecture in adults and children revealed a significant overlap. However, more robust studies with diverse ethnic representation is required in childhood obesity.

Benchmarking Mendelian randomization methods for causal inference using genome-wide association study summary statistics

Mendelian randomization (MR), which utilizes genetic variants as instrumental variables (IVs), has gained popularity as a method for causal inference between phenotypes using genetic data. While efforts have been made to relax IV assumptions and develop new methods for causal inference in the presence of invalid IVs due to confounding, the reliability of MR methods in real-world applications remains uncertain. Instead of using simulated datasets, we conducted a benchmark study evaluating 16 two-sample summary-level MR methods using real-world genetic datasets to provide guidelines for the best practices. Our study focused on the following crucial aspects: type I error control in the presence of various confounding scenarios (e.g., population stratification, pleiotropy, and family-level confounders like assortative mating), the accuracy of causal effect estimates, replicability, and power. By comprehensively evaluating the performance of compared methods over one thousand exposure-outcome trait pairs, our study not only provides valuable insights into the performance and limitations of the compared methods but also offers practical guidance for researchers to choose appropriate MR methods for causal inference.

Ontogenetic variation in the marine foraging of Atlantic salmon functionally links genomic diversity with a major life history polymorphism

The ecological role of heritable phenotypic variation in free-living populations remains largely unknown. Knowledge of the genetic basis of functional ecological processes can link genomic and phenotypic diversity, providing insight into polymorphism evolution and how populations respond to environmental changes. By quantifying the marine diet of Atlantic salmon, we assessed how foraging behaviour changes along the ontogeny, and in relation to genetic variation in two loci with major effects on age at maturity (six6 and vgll3). We used a two-component, zero-inflated negative binomial model to simultaneously quantify foraging frequency and foraging outcome, separately for fish and crustaceans diets. We found that older salmon forage for both prey types more actively (as evidenced by increased foraging frequency), but with a decreased efficiency (as evidenced by fewer prey in the diet), suggesting an age-dependent shift in foraging dynamics. The vgll3 locus was linked to age-dependent changes in foraging behaviour: Younger salmon with vgll3LL (the genotype associated with late maturation) tended to forage crustaceans more often than those with vgll3EE (the genotype associated with early maturation), whereas the pattern was reversed in older salmon. Vgll3 LL genotype was also linked to a marginal increase in fish acquisition, especially in younger salmon, while six6 was not a factor explaining the diet variation. Our results suggest a functional role for marine feeding behaviour linking genomic diversity at vgll3 with age at maturity among salmon, with potential age-dependent trade-offs maintaining the genetic variation. A shared genetic basis between dietary ecology and age at maturity likely subjects Atlantic salmon populations to evolution induced by bottom-up changes in marine productivity.

The association between dietary, physical activity and the DNA methylation of PPARGC1A, HLA-DQA1 and ADCY3 in pregnant women with gestational diabetes mellitus: a nest case-control study

BACKGROUND

Gestational diabetes mellitus (GDM) is associated with DNA methylation and lifestyle. The effects of DNA methylation on GDM, and the interaction between DNA methylation and lifestyle factors are not well elucidated. The objective of this study was to explore the association between GDM, DNA methylation and lifestyle factors.

METHODS

A nest case-control design was performed. Sociodemographic data, dietary intake and daily physical activity information of pregnant women were collected. Bisulfate pyrosequencing was used to detect the DNA methylation level of PPARGC1A, HLA-DQA1, and ADCY3 genes. The differences of DNA methylation levels between the GDM group and the control group were compared. The correlation between clinical characteristics, dietary, physical activity and DNA methylation level was analyzed.

RESULTS

A total of 253 pregnant women were enrolled, of which, 60 participants (GDM: 30; control: 30) were included in the final analysis. There were no significant differences in DNA methylation levels of six methylated sites between the two groups in this study (P > 0.05). Daily intake of potato and poultry were associated with DNA methylation level of the CpG 1 site of the ADCY3 gene in all participants and the control group (P < 0.05). Duration of folic acid intake before pregnancy was correlated with the methylation level of the CpG 1 site of the ADCY3 gene in all participants (r = 0.341, P = 0.04) and the control group (r = 0.431, P = 0.025). Daily oil intake was correlated with the methylation level of CpG 2 (r = 0.627, P = 0.016) and CpG 3 (r = 0.563, P = 0.036) of PPARGC1A in the GDM group.

CONCLUSION

The association between the DNA methylation levels and GDM wasn't validated. There were associations between dietary and DNA methylation in pregnant women. A large-sample-sized and longitudinal study is warranted to further investigate the impacts of lifestyle on DNA methylation.

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.

Longitudinal method comparison: modeling polygenic risk for post-traumatic stress disorder over time in individuals of African and European ancestry

Cross-sectional data allow the investigation of how genetics influence health at a single time point, but to understand how the genome impacts phenotype development, one must use repeated measures data. Ignoring the dependency inherent in repeated measures can exacerbate false positives and requires the utilization of methods other than general or generalized linear models. Many methods can accommodate longitudinal data, including the commonly used linear mixed model and generalized estimating equation, as well as the less popular fixed-effects model, cluster-robust standard error adjustment, and aggregate regression. We simulated longitudinal data and applied these five methods alongside naïve linear regression, which ignored the dependency and served as a baseline, to compare their power, false positive rate, estimation accuracy, and precision. The results showed that the naïve linear regression and fixed-effects models incurred high false positive rates when analyzing a predictor that is fixed over time, making them unviable for studying time-invariant genetic effects. The linear mixed models maintained low false positive rates and unbiased estimation. The generalized estimating equation was similar to the former in terms of power and estimation, but it had increased false positives when the sample size was low, as did cluster-robust standard error adjustment. Aggregate regression produced biased estimates when predictor effects varied over time. To show how the method choice affects downstream results, we performed longitudinal analyses in an adolescent cohort of African and European ancestry. We examined how developing post-traumatic stress symptoms were predicted by polygenic risk, traumatic events, exposure to sexual abuse, and income using four approaches-linear mixed models, generalized estimating equations, cluster-robust standard error adjustment, and aggregate regression. While the directions of effect were generally consistent, coefficient magnitudes and statistical significance differed across methods. Our in-depth comparison of longitudinal methods showed that linear mixed models and generalized estimating equations were applicable in most scenarios requiring longitudinal modeling, but no approach produced identical results even if fit to the same data. Since result discrepancies can result from methodological choices, it is crucial that researchers determine their model a priori, refrain from testing multiple approaches to obtain favorable results, and utilize as similar as possible methods when seeking to replicate results.

Variable parallelism in the genomic basis of age at maturity across spatial scales in Atlantic Salmon

Complex traits often exhibit complex underlying genetic architectures resulting from a combination of evolution from standing variation, hard and soft sweeps, and alleles of varying effect size. Increasingly, studies implicate both large-effect loci and polygenic patterns underpinning adaptation, but the extent that common genetic architectures are utilized during repeated adaptation is not well understood. Sea age or age at maturation represents a significant life history trait in Atlantic Salmon (Salmo salar), the genetic basis of which has been studied extensively in European Atlantic populations, with repeated identification of large-effect loci. However, the genetic basis of sea age within North American Atlantic Salmon populations remains unclear, as does the potential for a parallel trans-Atlantic genomic basis to sea age. Here, we used a large single-nucleotide polymorphism (SNP) array and low-coverage whole-genome resequencing to explore the genomic basis of sea age variation in North American Atlantic Salmon. We found significant associations at the gene and SNP level with a large-effect locus (vgll3) previously identified in European populations, indicating genetic parallelism, but found that this pattern varied based on both sex and geographic region. We also identified nonrepeated sets of highly predictive loci associated with sea age among populations and sexes within North America, indicating polygenicity and low rates of genomic parallelism. Despite low genome-wide parallelism, we uncovered a set of conserved molecular pathways associated with sea age that were consistently enriched among comparisons, including calcium signaling, MapK signaling, focal adhesion, and phosphatidylinositol signaling. Together, our results indicate parallelism of the molecular basis of sea age in North American Atlantic Salmon across large-effect genes and molecular pathways despite population-specific patterns of polygenicity. These findings reveal roles for both contingency and repeated adaptation at the molecular level in the evolution of life history variation.

Non-linear development of EEG coherence in adolescents and young adults shown by the analysis of neurophysiological trajectories and their covariance

To contribute to the understanding of changes in the factors governing the development of neural connectivity, the developmental structure of EEG coherence in adolescents and young adults was analyzed using the means, variances, and covariances of high alpha frequency band coherence measures from a set of 27 coherence pairs obtained from a sample of 1426 participants from the COGA study with 5006 observations over ages 12 through 31. Means and covariances were calculated at 96 age centers by a LOESS method. In the current study, trajectories of covariance matrices considered as individual units were determined by tensorial analysis: calculation of Riemannian geodesic (non-Euclidean) distances between matrices and application of both linear and non-linear dimension reduction techniques to these distances. Results were evaluated by bootstrap methods. Mean coherence trajectories for males and females were very similar, showing a steady upward trend from ages 12 to 20 which diminishes gradually from 20 to 25 and reaches stability from 25 to 31. In contrast, the individual covariance trajectories of males and female differed, with the male covariance levels becoming greater than that of females during the developmental process. Tensorial determination of the distances from the initial covariance matrix of subsequent covariance matrices to age 20 had the same trajectory as the mean coherence values. Tensorial determination of the trajectories of the covariance matrices of males and females based on their all pairs geodesic distances revealed a non-linear pattern in the multi-dimensional space of each of the trajectories: A steady increase in one dimension is accompanied by deviations from it peaking at age 20 which have both transient and lasting effects. There is a precise temporal parallelism of this pattern of covariance in males and females, while there is a consistent distance between the male and female covariance structures throughout the developmental process. Between region differences (anterior-posterior) within each sex are greater than between sex differences within regions. Examining development using multiple methods provides unique insight into the developmental process.

Interactions Between Genetic, Prenatal Substance Use, Puberty, and Parenting are Less Important for Understanding Adolescents' Internalizing, Externalizing, and Substance Use than Developmental Cascades in Multifactorial Models

This study tested interactions among puberty-related genetic risk, prenatal substance use, harsh discipline, and pubertal timing for the severity and directionality (i.e., differentiation) of externalizing and internalizing problems and adolescent substance use. This is a companion paper to Marceau et al. (2021) which examined the same influences in developmental cascade models. Data were from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (n = 4504 White boys, n = 4287 White girls assessed from the prenatal period through 18.5 years). We hypothesized generally that later predictors would strengthen the influence of puberty-related genetic risk, prenatal substance use exposure, and pubertal risk on psychopathology and substance use (two-way interactions), and that later predictors would strengthen the interactions of earlier influences on psychopathology and substance use (three-way interactions). Interactions were sparse. Although all fourteen interactions showed that later influences can exacerbate or trigger the effects of earlier ones, they often were not in the expected direction. The most robust moderator was parental discipline, and differing and synergistic effects of biological and socially-relevant aspects of puberty were found. In all, the influences examined here operate more robustly in developmental cascades than in interaction with each other for the development of psychopathology and transitions to substance use.

Maternal Age at Menarche Genes Determines Fetal Growth Restriction Risk

We aimed to explore the potential link of maternal age at menarche (mAAM) gene polymorphisms with risk of the fetal growth restriction (FGR). This case (FGR)-control (FGR free) study included 904 women (273 FGR and 631 control) in the third trimester of gestation examined/treated in the Departments of Obstetrics. For single nucleotide polymorphism (SNP) multiplex genotyping, 50 candidate loci of mAAM were chosen. The relationship of mAAM SNPs and FGR was appreciated by regression procedures (logistic/model-based multifactor dimensionality reduction [MB-MDR]) with subsequent in silico assessment of the assumed functionality pithy of FGR-related loci. Three mAAM-appertain loci were FGR-linked to genes such as KISS1 (rs7538038) (effect allele G-odds ratio (OR)allelic = 0.63/pperm = 0.0003; ORadditive = 0.61/pperm = 0.001; ORdominant = 0.56/pperm = 0.001), NKX2-1 (rs999460) (effect allele A-ORallelic = 1.37/pperm = 0.003; ORadditive = 1.45/pperm = 0.002; ORrecessive = 2.41/pperm = 0.0002), GPRC5B (rs12444979) (effect allele T-ORallelic = 1.67/pperm = 0.0003; ORdominant = 1.59/pperm = 0.011; ORadditive = 1.56/pperm = 0.009). The haplotype ACA FSHB gene (rs555621*rs11031010*rs1782507) was FRG-correlated (OR = 0.71/pperm = 0.05). Ten FGR-implicated interworking models were founded for 13 SNPs (pperm ≤ 0.001). The rs999460 NKX2-1 and rs12444979 GPRC5B interplays significantly influenced the FGR risk (these SNPs were present in 50% of models). FGR-related mAAM-appertain 15 polymorphic variants and 350 linked SNPs were functionally momentous in relation to 39 genes participating in the regulation of hormone levels, the ovulation cycle process, male gonad development and vitamin D metabolism. Thus, this study showed, for the first time, that the mAAM-appertain genes determine FGR risk.

Whole-Genome Scanning for Selection Signatures Reveals Candidate Genes Associated with Growth and Tail Length in Sheep

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

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.

Interrogating Causal Effects of Body Composition and Puberty-Related Risk Factors on Adolescent Idiopathic Scoliosis: A Two-Sample Mendelian Randomization Study

Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty-related exposures and AIS risk in Europeans and Asians. For our two-sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist-hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25-hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome-wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance-weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy-sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Investigating the role of genetic variation in vgll3 and six6 in the domestication of gilthead seabream (Sparus aurata Linnaeus) and European seabass (Dicentrarchus labrax Linnaeus)

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age-at-maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth-related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene-level genotyping'). The second dataset examined 300-600 k SNPs located in the chromosomes of the two genes, derived from a meta-analysis of a Pool-Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome-level genotyping'). The gene-level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome-level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Maternal Age at Menarche Gene Polymorphisms Are Associated with Offspring Birth Weight

In this study, the association between maternal age at menarche (AAM)-related polymorphisms and offspring birth weight (BW) was studied. The work was performed on a sample of 716 pregnant women and their newborns. All pregnant women underwent genotyping of 50 SNPs of AAM candidate genes. Regression methods (linear and Model-Based Multifactor Dimensionality Reduction (MB-MDR)) with permutation procedures (the indicator pperm was calculated) were used to identify the correlation between SNPs and newborn weight (transformed BW values were analyzed) and in silico bioinformatic examination was applied to assess the intended functionality of BW-associated loci. Four AAM-related genetic variants were BW-associated including genes such as POMC (rs7589318) (βadditive = 0.202/pperm = 0.015), KDM3B (rs757647) (βrecessive = 0.323/pperm = 0.005), INHBA (rs1079866) (βadditive = 0.110/pperm = 0.014) and NKX2-1 (rs999460) (βrecessive = -0.176/pperm = 0.015). Ten BW-significant models of interSNPs interactions (pperm ≤ 0.001) were identified for 20 polymorphisms. SNPs rs7538038 KISS1, rs713586 RBJ, rs12324955 FTO and rs713586 RBJ-rs12324955 FTO two-locus interaction were included in the largest number of BW-associated models (30% models each). BW-associated AAM-linked 22 SNPs and 350 proxy loci were functionally related to 49 genes relevant to pathways such as the hormone biosynthesis/process and female/male gonad development. In conclusion, maternal AMM-related genes polymorphism is associated with the offspring BW.

Genetic influences vary by age and sex: Trajectories of the association of cholinergic system variants and theta band event related oscillations

To characterize systemic changes in genetic effects on brain development, the age variation of the associations of cholinergic genetic variants and theta band event-related oscillations (EROs) was studied in a sample of 2140 adolescents and young adults, ages 12 to 25 from the COGA prospective study. The theta band EROs were elicited in visual and auditory oddball (target detection) tasks and measured by EEG recording. Associations were found to vary with age, sex, task modality (auditory or visual), and scalp locality. Seven of the twenty-one muscarinic and nicotinic cholinergic SNPs studied in the analysis, from CHRM2, CHRNA3, CHRNA5, and CHRNB4, had significant effects on theta band EROs with considerable age spans for some sex-modality combination. No SNP-age-modality combination had significant effects in the same direction for males and females. Results suggest that nicotinic receptor associations are stronger before age 18, while muscarinic receptor associations are stronger after age 18.

Support Interval for Two-Sample Summary Data-Based Mendelian Randomization

The summary-data-based Mendelian randomization (SMR) method is gaining popularity in estimating the causal effect of an exposure on an outcome. In practice, the instrument SNP is often selected from the genome-wide association study (GWAS) on the exposure but no correction is made for such selection in downstream analysis, leading to a biased estimate of the effect size and invalid inference. We address this issue by using the likelihood derived from the sampling distribution of the estimated SNP effects in the exposure GWAS and the outcome GWAS. This likelihood takes into account how the instrument SNPs are selected. Since the effective sample size is 1, the asymptotic theory does not apply. We use a support for a profile likelihood as an interval estimate of the causal effect. Simulation studies indicate that this support has robust coverage while the confidence interval implied by the SMR method has lower-than-nominal coverage. Furthermore, the variance of the two-stage least squares estimate of the causal effect is shown to be the same as the variance used for SMR for one-sample data when there is no selection.

gBLUP-GWAS identifies candidate genes, signaling pathways, and putative functional polymorphisms for age at puberty in gilts

Successful development of replacement gilts determines their reproductive longevity and lifetime productivity. Selection for reproductive longevity is challenging due to low heritability and expression late in life. In pigs, age at puberty is the earliest known indicator for reproductive longevity and gilts that reach puberty earlier have a greater probability of producing more lifetime litters. Failure of gilts to reach puberty and display a pubertal estrus is a major reason for early removal of replacement gilts. To identify genomic sources of variation in age at puberty for improving genetic selection for early age at puberty and related traits, gilts (n = 4,986) from a multigeneration population representing commercially available maternal genetic lines were used for a genomic best linear unbiased prediction-based genome-wide association. Twenty-one genome-wide significant single nucleotide polymorphisms (SNP) located on Sus scrofa chromosomes (SSC) 1, 2, 9, and 14 were identified with additive effects ranging from -1.61 to 1.92 d (P < 0.0001 to 0.0671). Novel candidate genes and signaling pathways were identified for age at puberty. The locus on SSC9 (83.7 to 86.7 Mb) was characterized by long range linkage disequilibrium and harbors the AHR transcription factor gene. A second candidate gene on SSC2 (82.7 Mb), ANKRA2, is a corepressor for AHR, suggesting a possible involvement of AHR signaling in regulating pubertal onset in pigs. Putative functional SNP associated with age at puberty in the AHR and ANKRA2 genes were identified. Combined analysis of these SNP showed that an increase in the number of favorable alleles reduced pubertal age by 5.84 ± 1.65 d (P < 0.001). Candidate genes for age at puberty showed pleiotropic effects with other fertility functions such as gonadotropin secretion (FOXD1), follicular development (BMP4), pregnancy (LIF), and litter size (MEF2C). Several candidate genes and signaling pathways identified in this study play a physiological role in the hypothalamic-pituitary-gonadal axis and mechanisms permitting puberty onset. Variants located in or near these genes require further characterization to identify their impact on pubertal onset in gilts. Because age at puberty is an indicator of future reproductive success, these SNP are expected to improve genomic predictions for component traits of sow fertility and lifetime productivity expressed later in life.

Nutrition in school-age children: a rationale for revisiting priorities

Middle childhood and early adolescence have received disproportionately low levels of scientific attention relative to other life stages, especially as related to nutrition and health. This is partly due to the justified emphasis on the first 1000 days of life, and the idea that early deficits and consequences may not be fully reversible. In addition, these stages of life may superficially appear less "eventful" than infancy or late adolescence. Finally, there has been historical ambiguity and inconsistency in terminology, depending on whether viewing "childhood" through physiologic, social, legal, or other lenses. Nevertheless, this age bracket, which encompasses most of the primary education and basic schooling years for most individuals, is marked by significant changes, inflection points, and sexually driven divergence in somatic and brain growth and development trajectories. These constitute transformative changes, and thus middle childhood and early adolescence represents a major and last opportunity to influence long-term health and productivity. This review highlights the specificities of growth and development in school age, with a focus on middle childhood and early adolescence (5 years-15 years of age, for the purposes of this review), the role of nutrition, the short- and long-term consequences of inadequate nutrition, and the current global status of nutrition in this age group. Adequate attention and emphasis on nutrition in the school-age years is critical: (a) for maintaining an adequate course of somatic and cognitive development, (b) for taking advantage of this last major opportunity to correct deficits of undernutrition and "catch-up" to normal life course development, and (c) for addressing the nutritional inadequacies and mitigating the longer-term consequences of overnutrition. This review summarizes and provides a rationale for prioritizing nutrition in school-age children, and for the need to revisit priorities and focus on this part of the life cycle to maximize individuals' potential and their contribution to society.

Life-history genotype explains variation in migration activity in Atlantic salmon (Salmo salar)

One of the most well-known life-history continuums is the fast-slow axis, where 'fast' individuals mature earlier than 'slow' individuals. 'Fast' individuals are predicted to be more active than 'slow' individuals because high activity is required to maintain a fast life-history strategy. Recent meta-analyses revealed mixed evidence for such integration. Here, we test whether known life-history genotypes differ in activity expression by using Atlantic salmon (Salmo salar) as a model. In salmon, variation in Vgll3, a transcription cofactor, explains approximately 40% of variation in maturation timing. We predicted that the allele related to early maturation (vgll3*E) would be associated with higher activity. We used an automated surveillance system to follow approximately 1900 juveniles including both migrants and non-migrants (i.e. smolt and parr fish, respectively) in semi-natural conditions over 31 days (approx. 580 000 activity measurements). In migrants, but not in non-migrants, vgll3 explained variation in activity according to our prediction in a sex-dependent manner. Specifically, in females the vgll3*E allele was related to increasing activity, whereas in males the vgll3*L allele (later maturation allele) was related to increasing activity. These sex-dependent effects might be a mechanism maintaining within-population genetic life-history variation.

NUCKS1 is a highly modified, chromatin-associated protein involved in a diverse set of biological and pathophysiological processes

The Nuclear Casein and Cyclin-dependent Kinase Substrate 1 (NUCKS1) protein is highly conserved in vertebrates, predominantly localized to the nucleus and one of the most heavily modified proteins in the human proteome. NUCKS1 expression is high in stem cells and the brain, developmentally regulated in mice and associated with several diverse malignancies in humans, including cancer, metabolic syndrome and Parkinson's disease. NUCKS1 function has been linked to modulating chromatin architecture and transcription, DNA repair and cell cycle regulation. In this review, we summarize and discuss the published information on NUCKS1 and highlight the questions that remain to be addressed to better understand the complex biology of this multifaceted protein.

Proteomic Analysis of Hypothalamus and Pituitary Gland in Pre and Postpubertal Brahman Heifers

The hypothalamus and the pituitary gland are directly involved in the complex systemic changes that drive the onset of puberty in cattle. Here, we applied integrated bioinformatics to elucidate the critical proteins underlying puberty and uncover potential molecular mechanisms from the hypothalamus and pituitary gland of prepubertal (n = 6) and postpubertal (n = 6) cattle. Proteomic analysis in the hypothalamus and pituitary gland revealed 275 and 186 differentially abundant (DA) proteins, respectively (adjusted p-value < 0.01). The proteome profiles found herein were integrated with previously acquired transcriptome profiles. These transcriptomic studies used the same tissues harvested from the same heifers at pre- and post-puberty. This comparison detected a small number of matched transcripts and protein changes at puberty in each tissue, suggesting the need for multiple omics analyses for interpreting complex biological systems. In the hypothalamus, upregulated DA proteins at post-puberty were enriched in pathways related to puberty, including GnRH, calcium and oxytocin signalling pathways, whereas downregulated proteins were observed in the estrogen signalling pathway, axon guidance and GABAergic synapse. Additionally, this study revealed that ribosomal pathway proteins in the pituitary were involved in the pubertal development of mammals. The reported molecules and derived protein-protein networks are a starting point for future experimental approaches that might dissect with more detail the role of each molecule to provide new insights into the mechanisms of puberty onset in cattle.

The Association Between Puberty Timing and Body Mass Index in a Longitudinal Setting: The Contribution of Genetic Factors

We analyzed the contribution of genetic factors on the association between puberty timing and body mass index (BMI) using longitudinal data and two approaches: (i) genetic twin design and (ii) polygenic scores (PGS) of obesity indices. Our data were derived from Finnish cohorts: 9080 twins had information on puberty timing and BMI and 2468 twins also had genetic data. Early puberty timing was moderately associated with higher BMI in childhood in both boys and girls; in adulthood these correlations were weaker and largely disappeared after adjusting for childhood BMI. The largest proportion of these correlations was attributable to genetic factors. The higher PGSs of BMI and waist circumference were associated with earlier timing of puberty in girls, whereas weaker associations were found in boys. Early puberty is not an independent risk factor for adult obesity but rather reflects the association between puberty timing and childhood BMI contributed by genetic predisposition.

Maternal pre-pregnancy body mass index, gestational weight gain, and pubertal timing in daughters: A systematic review and meta-analysis of cohort studies

The timing of daughter's puberty onset is constantly earlier. It is still unclear about the maternal pre-pregnancy body mass index (BMI) and gestational weight gain (GWG) as important prenatal factors that may affect offspring's onset of puberty. Thus, we evaluated the association among maternal pre-pregnancy BMI, GWG, and daughters' early pubertal development based on the existing literature. Literature review was conducted in different databases, including Web of Science, Pubmed, Wiley, ScienceDirect, Web of Science, and Chinese National Knowledge Infrastructure databases up to June 2021. We selected random effects model or fixed effects model for meta-analysis according to the I² statistics value to obtain the summary measurement. A total of 12 cohort studies were included. Compared to maternal pre-pregnancy normal weight, maternal pre-pregnancy overall overweight/obesity (RR = 1.24; 95% CI 1.17 to 1.32), obesity (RR = 1.35; 95% CI 1.23 to 1.48), and overweight (RR = 1.17; 95% CI 1.09 to 1.26) were significantly associated with the increased risk of earlier timing of pubertal onset in daughters. Daughters born of mothers with pre-pregnancy overall overweight/obesity, obesity, and overweight had earlier pubertal onset compared to those born of mothers with normal weight ([mean difference = -3.03, 95% CI: -3.97 to -2.10], [mean difference = -3.50, 95% CI: -5.38 to -1.62], and [mean difference = -2.89, 95% CI: -4.07 to -1.71], respectively). The effects were also significant in the assessed three milestones (menarche, breast development, and pubic hair development). Maternal excessive GWG increased the risk of early pubertal timing in daughters (RR = 1.19; 95% CI 1.09 to 1.30).

Estimation of Causal Effect of Age at Menarche on Pubertal Height Growth Using Mendelian Randomization

We use Mendelian randomization to estimate the causal effect of age at menarche on late pubertal height growth and total pubertal height growth. The instrument SNPs selected from the exposure genome-wide association study (GWAS) are validated in additional population-matched exposure GWASs. Based on the inverse variance weighting method, there is a positive causal relationship of age at menarche on late pubertal growth (β^=0.56, 95% CI: (0.34, 0.78), p=3.16×10−7) and on total pubertal growth (β^=0.36, 95% CI: (0.14, 0.58), p=1.30×10−3). If the instrument SNPs are not validated in additional exposure GWASs, the estimated effect on late pubertal height growth increases by 3.6% to β^=0.58 (95% CI: (0.42, 0.73), p=4.38×10−13) while the estimates on total pubertal height growth increases by 41.7% to β^=0.51 (95% CI: (0.35, 0.67), p=2.96×10−11).

Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite their notable anthropometric variability, ancestry proportions, and high burden of growth stunting and overweight/obesity. To address this knowledge gap, we analyzed densely imputed genetic data in a sample of Hispanic/Latino adults to identify and fine-map genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (stage 1, n = 59,771) and generalized our findings in 9 additional studies (stage 2, n = 10,538). We conducted a trans-ancestral GWAS with summary statistics from HISLA stage 1 and existing consortia of European and African ancestries. In our HISLA stage 1 + 2 analyses, we discovered one BMI locus, as well as two BMI signals and another height signal each within established anthropometric loci. In our trans-ancestral meta-analysis, we discovered three BMI loci, one height locus, and one WHRadjBMI locus. We also identified 3 secondary signals for BMI, 28 for height, and 2 for WHRadjBMI in established loci. We show that 336 known BMI, 1,177 known height, and 143 known WHRadjBMI (combined) SNPs demonstrated suggestive transferability (nominal significance and effect estimate directional consistency) in Hispanic/Latino adults. Of these, 36 BMI, 124 height, and 11 WHRadjBMI SNPs were significant after trait-specific Bonferroni correction. Trans-ancestral meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our findings demonstrate that future studies may also benefit from leveraging diverse ancestries and differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

Ultrasonic bone age fractionates cognitive abilities in adolescence

Adolescent development is not only shaped by the mere passing of time and accumulating experience, but it also depends on pubertal timing and the cascade of maturational processes orchestrated by gonadal hormones. Although individual variability in puberty onset confounds adolescent studies, it has not been efficiently controlled for. Here we introduce ultrasonic bone age assessment to estimate biological maturity and disentangle the independent effects of chronological and biological age on adolescent cognitive abilities. Comparing cognitive performance of female participants with different skeletal maturity we uncover the impact of biological age on both IQ and specific abilities. We find that biological age has a selective effect on abilities: more mature individuals within the same age group have higher working memory capacity and processing speed, while those with higher chronological age have better verbal abilities, independently of their maturity. Based on our findings, bone age is a promising biomarker of adolescent maturity.

Gene Region Association Analysis of Longitudinal Quantitative Traits Based on a Function-On-Function Regression Model

In the process of growth and development in life, gene expressions that control quantitative traits will turn on or off with time. Studies of longitudinal traits are of great significance in revealing the genetic mechanism of biological development. With the development of ultra-high-density sequencing technology, the associated analysis has tremendous challenges to statistical methods. In this paper, a longitudinal functional data association test (LFDAT) method is proposed based on the function-on-function regression model. LFDAT can simultaneously treat phenotypic traits and marker information as continuum variables and analyze the association of longitudinal quantitative traits and gene regions. Simulation studies showed that: 1) LFDAT performs well for both linkage equilibrium simulation and linkage disequilibrium simulation, 2) LFDAT has better performance for gene regions (include common variants, low-frequency variants, rare variants and mixture), and 3) LFDAT can accurately identify gene switching in the growth and development stage. The longitudinal data of the Oryza sativa projected shoot area is analyzed by LFDAT. It showed that there is the advantage of quick calculations. Further, an association analysis was conducted between longitudinal traits and gene regions by integrating the micro effects of multiple related variants and using the information of the entire gene region. LFDAT provides a feasible method for studying the formation and expression of longitudinal traits.

Identifying Environmental Endocrine Disruptors Associated With the Age at Menarche by Integrating a Transcriptome-Wide Association Study With Chemical-Gene-Interaction Analysis

Menarche is the first occurrence of menstrual bleeding and one of the most important events of female puberty. Alarmingly, over the last several decades, the mean age at menarche (AAM) has decreased. Environmental endocrine disruptors (EEDs) are chemicals that may interfere with the endocrine system, resulting in adverse developmental, immunological, neurological, and reproductive effects in humans. Thus, the effects of EEDs on fertility and reproduction are growing concerns in modern societies. In this study, we aimed to determine the influence of genetic and environmental factors on AAM. We used data from an AAM genome-wide association study of 329,345 women to conduct a transcriptome-wide association study (TWAS) with FUSION software. As references, we determined the gene-expression levels in the hypothalamus, pituitary gland, ovaries, uterus, and whole blood. We performed Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses using the significantly dysregulated genes identified by the TWAS. Using the STRING database, we also generated a protein-protein-interaction network to analyze common AAM-specific genes identified by the TWAS with different tissues. We performed chemical-related gene set enrichment analysis (CGSEA) and identified significant TWAS genes to uncover relationships between different chemicals and AAM. The TWAS identified 9,848 genes; among these, 1580 genes were significant (P < 0.05), and 11 genes were significant among the hypothalamus, pituitary, ovary, uterus, and whole blood. CGSEA identified 1,634 chemicals, including 120 chemicals significantly correlated with AAM. In summary, we performed a TWAS (for genetic factors) and CGSEA (for environmental factors) focusing on AAM and identified several AAM-associated genes and EEDs. The results of this study expand our understanding of genetic and environmental factors related to the onset of female puberty.

Timing of Pubertal Development and Midlife Blood Pressure in Men and Women: A Mendelian Randomization Study

INTRODUCTION

Observational studies suggest earlier puberty is associated with higher adulthood blood pressure (BP), but these findings have not been replicated using Mendelian randomization (MR). We examined this question sex-specifically using larger genome-wide association studies (GWAS) with more extensive measures of pubertal timing.

METHODS

We obtained genetic instruments proxying pubertal maturation (age at menarche [AAM] or voice breaking [AVB]) from the largest published GWAS. We applied them to summary sex-specific genetic associations with systolic and diastolic BP z-scores, and self-reported hypertension in women (n = 194 174) and men (n = 167 020) from the UK Biobank, using inverse-variance weighted meta-analysis. We conducted sensitivity analyses using other MR methods, including multivariable MR adjusted for childhood obesity proxied by body mass index (BMI). We used late pubertal growth as a validation outcome.

RESULTS

AAM (beta per 1-year later = -0.030 [95% confidence interval, -0.055 to -0.005] and AVB (beta -0.058 [95% CI, -0.100 to -0.015]) were inversely associated with systolic BP independent of childhood BMI, as were diastolic BP (-0.035 [95% CI, -0.060 to -0.009] for AAM and -0.046 [95% CI, -0.089 to -0.004] for AVB) and self-reported hypertension (odds ratio 0.89 [95% CI, 0.84-0.95] for AAM and 0.87 [95% CI, 0.79-0.96] for AVB). AAM and AVB were positively associated with late pubertal growth, as expected. The results were robust to sensitivity analysis using other MR methods.

CONCLUSION

Timing of pubertal maturation was associated with adulthood BP independent of childhood BMI, highlighting the role of pubertal maturation timing in midlife BP.

Cis-regulatory architecture of human ESC-derived hypothalamic neuron differentiation aids in variant-to-gene mapping of relevant complex traits

The hypothalamus regulates metabolic homeostasis by influencing behavior and endocrine systems. Given its role governing key traits, such as body weight and reproductive timing, understanding the genetic regulation of hypothalamic development and function could yield insights into disease pathogenesis. However, given its inaccessibility, studying human hypothalamic gene regulation has proven challenging. To address this gap, we generate a high-resolution chromatin architecture atlas of an established embryonic stem cell derived hypothalamic-like neuron model across three stages of in vitro differentiation. We profile accessible chromatin and identify physical contacts between gene promoters and putative cis-regulatory elements to characterize global regulatory landscape changes during hypothalamic differentiation. Next, we integrate these data with GWAS loci for various complex traits, identifying multiple candidate effector genes. Our results reveal common target genes for these traits, potentially affecting core developmental pathways. Our atlas will enable future efforts to determine hypothalamic mechanisms influencing disease susceptibility.

CYP11B1 variants influence skeletal maturation via alternative splicing

We performed genome-wide association study meta-analysis to identify genetic determinants of skeletal age (SA) deviating in multiple growth disorders. The joint meta-analysis (N = 4557) in two multiethnic cohorts of school-aged children identified one locus, CYP11B1 (expression confined to the adrenal gland), robustly associated with SA (rs6471570-A; β = 0.14; P = 6.2 × 10-12). rs6410 (a synonymous variant in the first exon of CYP11B1 in high LD with rs6471570), was prioritized for functional follow-up being second most significant and the one closest to the first intron-exon boundary. In 208 adrenal RNA-seq samples from GTEx, C-allele of rs6410 was associated with intron 3 retention (P = 8.11 × 10-40), exon 4 inclusion (P = 4.29 × 10-34), and decreased exon 3 and 5 splicing (P = 7.85 × 10-43), replicated using RT-PCR in 15 adrenal samples. As CYP11B1 encodes 11-β-hydroxylase, involved in adrenal glucocorticoid and mineralocorticoid biosynthesis, our findings highlight the role of adrenal steroidogenesis in SA in healthy children, suggesting alternative splicing as a likely underlying mechanism.

Refining the genomic location of single nucleotide polymorphism variation affecting Atlantic salmon maturation timing at a key large-effect locus

Efforts to understand the genetic underpinnings of phenotypic variation are becoming more and more frequent in molecular ecology. Such efforts often lead to the identification of candidate regions showing signals of association and/or selection. These regions may contain multiple genes and therefore validation of which genes are actually responsible for the signal is required. In Atlantic salmon (Salmo salar), a large-effect locus for maturation timing, an ecologically important trait, occurs in a genomic region including two genes, vgll3 and akap11, but data for clearly determining which of the genes (or both) contribute to the association have been lacking. Here, we take advantage of natural recombination events detected between the two candidate genes in a salmon broodstock to reduce linkage disequilibrium at the locus, thus enabling delineation of the influence of variation at these two genes on early maturation. By rearing 5,895 males to maturation age, of which 81% had recombinant vgll3/akap11 allelic combinations, we found that vgll3 single nucleotide polymorphism (SNP) variation was strongly associated with early maturation, whereas there was little or no association between akap11 SNP variation and early maturation. These findings provide strong evidence supporting vgll3 as the primary candidate gene in the chromosome 25 locus for influencing early maturation. This will help guide future research for understanding the genetic processes controlling early maturation. This also exemplifies the utility of natural recombinants to more precisely map causal variation underlying ecologically important phenotypic diversity.

Predicting physiological aging rates from a range of quantitative traits using machine learning

It is widely thought that individuals age at different rates. A method that measures “physiological age” or physiological aging rate independent of chronological age could therefore help elucidate mechanisms of aging and inform an individual’s risk of morbidity and mortality. Here we present machine learning frameworks for inferring individual physiological age from a broad range of biochemical and physiological traits including blood phenotypes (e.g., high-density lipoprotein), cardiovascular functions (e.g., pulse wave velocity) and psychological traits (e.g., neuroticism) as main groups in two population cohorts SardiNIA (~6,100 participants) and InCHIANTI (~1,400 participants). The inferred physiological age was highly correlated with chronological age (R2 > 0.8). We further defined an individual’s physiological aging rate (PAR) as the ratio of the predicted physiological age to the chronological age. Notably, PAR was a significant predictor of survival, indicating an effect of aging rate on mortality. Our trait-based PAR was correlated with DNA methylation-based epigenetic aging score (r = 0.6), suggesting that both scores capture a common aging process. PAR was also substantially heritable (h2~0.3), and a subsequent genome-wide association study of PAR identified significant associations with two genetic loci, one of which is implicated in telomerase activity. Our findings support PAR as a proxy for an underlying whole-body aging mechanism. PAR may thus be useful to evaluate the efficacy of treatments that target aging-related deficits and controllable epidemiological factors.

Polygenic and major-locus contributions to sexual maturation timing in Atlantic salmon

Sexual maturation timing is a life-history trait central to the balance between mortality and reproduction. Maturation may be triggered when an underlying compound trait, called liability, exceeds a threshold. In many different species and especially fishes, this liability is approximated by growth and body condition. However, environmental vs. genetic contributions either directly or via growth and body condition to maturation timing remain unclear. Uncertainty exists also because the maturation process can reverse this causality and itself affect growth and body condition. In addition, disentangling the contributions of polygenic and major loci can be important. In many fishes, males mature before females, enabling the study of associations between male maturation and maturation-unbiased female liability traits. Using 40 Atlantic salmon families, longitudinal common-garden experimentation, and quantitative genetic analyses, we disentangled environmental from polygenic and major locus (vgll3) effects on male maturation, and sex-specific growth and condition. We detected polygenic heritabilities for maturation, growth, and body condition, and vgll3 effects on maturation and body condition but not on growth. Longitudinal patterns for sex-specific phenotypic liability, and for genetic variances and correlations between sexes suggested that early growth and condition indeed positively affected maturation initiation. However, towards spawning time, causality appeared reversed for males whereby maturation affected growth negatively and condition positively via both the environmental and genetic effects. Altogether, the results indicate that growth and condition are useful traits to study liability for maturation initiation, but only until maturation alters their expression, and that vgll3 contributes to maturation initiation via condition.

[Analysis on the law of height growth spurt in adolescence of children and adolescents in Zhongshan City]

OBJECTIVE

To analyze the characteristics of the age at peak height velocity and peak height velocity of primary and middle school students in Zhongshan City, Guangdong Province, and to explore the law of the sudden increase in adolescent height in this area, and to understand the law of height growth spurt in adolescence by longitudinal tracking of the height of children and adolescents in Zhong-shan City.

METHODS

Based on the physical examination database of primary and middle school students in Zhongshan City, Guangdong Province from 2005 to 2016, individuals who had been continuously tracked for more than 6 times were selected as research samples. SITAR model was used to fit the height data of the sample population, and the age at peak height velocity and peak height velocity were calcula-ted.

RESULTS

A total of 49 579 subjects were included in this study, including 26 524 boys and 26 008 urban students. The median follow-up ages of boys and girls were 7.74 and 7.72 years, respectively. The boy's height spurt peak age was (12.72±0.89) years, and later than the girls at the age of (10.98±0.95) years (t=207.639, P < 0.001), the boy's height spurt peak velocity of (10.12±1.49) cm/year, higher than the girls of (8.35±1.12) cm/year (t=150.826, P < 0.001). The gender differences of height spurt peak age and height spurt peak speed in urban and rural students were consistent with the whole sample. The height surge peak age of urban male students was earlier than that of rural male students, and the height surge peak speed of urban female students was lower than that of rural female students.

CONCLUSION

The peak age of the surge of girls was earlier than that of boys, but the peak rate of the surge of girls was lower than that of boys, the peak age of urban students was earlier than that of rural students, but the peak rate of urban boys was lower than that of rural boys in Guangdong Province.

[Analysis on the law of height growth spurt in adolescence of children and adolescents in Zhongshan City]

OBJECTIVE

To analyze the characteristics of the age at peak height velocity and peak height velocity of primary and middle school students in Zhongshan City, Guangdong Province, and to explore the law of the sudden increase in adolescent height in this area, and to understand the law of height growth spurt in adolescence by longitudinal tracking of the height of children and adolescents in Zhong-shan City.

METHODS

Based on the physical examination database of primary and middle school students in Zhongshan City, Guangdong Province from 2005 to 2016, individuals who had been continuously tracked for more than 6 times were selected as research samples. SITAR model was used to fit the height data of the sample population, and the age at peak height velocity and peak height velocity were calcula-ted.

RESULTS

A total of 49 579 subjects were included in this study, including 26 524 boys and 26 008 urban students. The median follow-up ages of boys and girls were 7.74 and 7.72 years, respectively. The boy's height spurt peak age was (12.72±0.89) years, and later than the girls at the age of (10.98±0.95) years (t=207.639, P < 0.001), the boy's height spurt peak velocity of (10.12±1.49) cm/year, higher than the girls of (8.35±1.12) cm/year (t=150.826, P < 0.001). The gender differences of height spurt peak age and height spurt peak speed in urban and rural students were consistent with the whole sample. The height surge peak age of urban male students was earlier than that of rural male students, and the height surge peak speed of urban female students was lower than that of rural female students.

CONCLUSION

The peak age of the surge of girls was earlier than that of boys, but the peak rate of the surge of girls was lower than that of boys, the peak age of urban students was earlier than that of rural students, but the peak rate of urban boys was lower than that of rural boys in Guangdong Province.

Interpreting type 1 diabetes risk with genetics and single-cell epigenomics

Genetic risk variants that have been identified in genome-wide association studies of complex diseases are primarily non-coding1. Translating these risk variants into mechanistic insights requires detailed maps of gene regulation in disease-relevant cell types2. Here we combined two approaches: a genome-wide association study of type 1 diabetes (T1D) using 520,580 samples, and the identification of candidate cis-regulatory elements (cCREs) in pancreas and peripheral blood mononuclear cells using single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) of 131,554 nuclei. Risk variants for T1D were enriched in cCREs that were active in T cells and other cell types, including acinar and ductal cells of the exocrine pancreas. Risk variants at multiple T1D signals overlapped with exocrine-specific cCREs that were linked to genes with exocrine-specific expression. At the CFTR locus, the T1D risk variant rs7795896 mapped to a ductal-specific cCRE that regulated CFTR; the risk allele reduced transcription factor binding, enhancer activity and CFTR expression in ductal cells. These findings support a role for the exocrine pancreas in the pathogenesis of T1D and highlight the power of large-scale genome-wide association studies and single-cell epigenomics for understanding the cellular origins of complex disease.