Hundreds of variants clustered in genomic loci and biological pathways affect human height

The GWAS candidate far upstream element binding protein 3 (FUBP3) is required for normal skeletal growth, and adult bone mass and strength in mice

Bone mineral density (BMD) and height are highly heritable traits for which hundreds of genetic loci have been linked through genome wide association studies (GWAS). FUBP3 is a DNA and RNA binding protein best characterised as a transcriptional regulator of c-Myc, but little is known about its role in vivo. Single nucleotide polymorphisms in FUBP3 at the 9q34.11 locus have been associated with BMD, fracture and height in multiple GWAS, but FUBP3 has no previously established role in the skeleton. We analysed Fubp3-deficient mice to determine the consequence of FUBP3 deficiency in vivo. Mice lacking Fubp3 had reduced survival to adulthood and impaired skeletal growth. Bone mass was decreased, most strikingly in the vertebrae, with altered trabecular micro-architecture. Fubp3 deficient bones were also weak. These data provide the first functional demonstration that Fubp3 is required for normal skeletal growth and development and maintenance of adult bone structure and strength, indicating that FUBP3 contributes to the GWAS association of 9q34.11 with variation in height, BMD and fracture.

Examining the association between fetal HLA-C, maternal KIR haplotypes and birth weight

Human birth weight is under stabilizing selection, seeking balance between extremes of high and low, thereby reducing fetal and maternal perinatal mortality risk. Certain combinations of maternal killer immunoglobulin-like receptor (KIR) and paternally derived fetal human leuokocyte antigen-C (HLA-C) alleles were previously associated with higher risk of high and low birth weight in a study with limited sample size (n=1,316). Using recently developed methods to impute HLA and KIR haplotypes using single nucleotide polymorphism (SNP) genotype data, we tested associations of fetal HLA and maternal KIR genotypes with offspring birth weight in a large sample. We imputed KIR haplotypes using the KIR*IMP imputation software in 10,602 mother-offspring pairs of European descent from singleton pregnancies from five studies. Using mixed linear regression models to account for mothers with multiple children, we tested associations between maternal KIR A vs B haplotypes (AA, AB/BA, BB genotypes) as well as copy number of activating receptor gene KIR2DS1 (0, 1, 2 copies of the gene) in the presence of fetal HLA C1/C2 alleles, and offspring birth weight. Associations were analyzed in each cohort before performing a meta-analysis to estimate the interaction effects between maternal KIR and fetal HLA-C2 on birth weight across the entire sample. The KIR haplotypes achieved imputation accuracy estimated at >95% in most of the cohorts. No interaction effects were observed between either the maternal A vs. B haplotype or the maternal KIR2DS1 locus and fetal HLA-C. When specifically trying to replicate the previously associated combination of maternal KIR2DS1 and paternally inherited fetal HLA-C2, there was a negligible change in offspring birth weight for each additional KIR2DS1 allele and HLA-C2 of paternal origin (7g lower birth weight per allele [95% CI: -54, 40], P = 0.78). We found little evidence of association between birth weight and maternal KIR haplotypes or fetal HLA-C2 and were unable to replicate previously reported findings. Our observations reinforce the importance of replication and the use of large sample sizes in the validation of genetic associations.

Evolution of long scalp hair in humans

The ability to grow long scalp hair is a distinct human characteristic. It probably originally evolved to aid in cooling the sun-exposed head, although the genetic determinants of long hair are largely unknown. Despite ancestral variations in hair growth, long scalp hair is common to all extant human populations, which suggests its emergence before or concurrently with the emergence of anatomically modern humans (AMHs), approximately 300 000 years ago. Long scalp hair in AMHs was also a trait that was selected because it conveyed essential signals related to an individual's age, sexual maturity, health and social status. Biologically, hair length is primarily determined by the amount of time that a hair follicle spends in the active growth phase (anagen). While anagen duration is typically tightly regulated in most mammals, the inherent ability of a hair follicle to continuously recruit new dividing progenitors to its base, where hair fibre is generated, theoretically removes limits on maximal anagen duration. We propose a model wherein hair cycle progression into and out of anagen is regulated by evolutionary malleable molecular checkpoints. Several animal species and domesticated animal breeds display long body hair, which suggests that extremely long scalp hair in humans emerged via attenuation of an existing out-of-anagen checkpoint mechanism rather than via a newly evolved molecular programme. Studying congenital and somatic mosaicism conditions featuring altered hair length could potentially unveil the currently unknown molecular basis underlying this human trait.

Genome-wide association study for resistance to Macrophomina phaseolina in maize (Zea mays L.)

Maize (Zea mays L.) is a frequently used food source in human and animal nutrition. Macrophomina phaseolina is a fungal pathogen causing charcoal rot disease in many plants, especially maize. This pathogen causes high yield losses in maize. The development of resistant maize genotypes is of great importance in controlling this disease. In this study, the population structure of 120 different maize genotypes with varying levels of disease resistance was determined and genome-wide association studies were performed. Each genotype was subjected to the pathogen under controlled conditions and their phenotypic responses to the disease were analyzed. Afterwards, single nucleotide polymorphisms were determined by DArT-seq sequencing. After filtering the SNP data, 37,470 clean SNPs were obtained. The population structure was analyzed with STRUCTURE software, and it was determined that the population was divided into two subgroups. The relationship between phenotypic and genotypic data was analyzed using the MLM (Q + K) model in TASSEL software. As a result, seven SNPs markers located on four different chromosomes were associated with disease resistance. The related markers can be used in the future for the development of maize varieties resistant to M. phaseolina by marker-assisted selection.

The structure and function of the ghrelin receptor coding for drug actions

Drugs targeting the ghrelin receptor hold therapeutic potential in anorexia, obesity and diabetes. However, developing effective drugs is challenging. To tackle this common issue across a broad drug target, this study aims to understand how anamorelin, the only approved drug targeting the ghrelin receptor, operates compared to other synthetic drugs. Our research elucidated the receptor's structure with anamorelin and miniGq, unveiling anamorelin's superagonistic activity. We demonstrated that ligands with distinct chemical structures uniquely bind to the receptor, resulting in diverse conformations and biasing signal transduction. Moreover, our study showcased the utility of structural information in effectively identifying natural genetic variations altering drug action and causing severe functional deficiencies, offering a basis for selecting the right medication on the basis of the individual's genomic sequence. Thus, by building on structural analysis, this study enhances the foundational framework for selecting therapeutic agents targeting the ghrelin receptor, by effectively leveraging signaling bias and genetic variations.

MARSweb: a fully automated web service for set-based association testing

BACKGROUND

Despite the successes in GWAS, there is still a large gap between the known heritability and the part explained by the SNPs identified by GWAS. Set-based analysis is one of the approaches that has tried to identify associations between multiple variants in a locus a trait, leveraging allelic heterogeneity to increase power in association testing. MARS is a set-based analysis method that integrates likelihood ratio test with a recently developed fine mapping technique to accurately account for causal status of variants in a risk locus. Unfortunately, due to its complex running process, time complexity, and the requirement of high-performance computing resources, it is not widely used.

RESULTS

To address these issues, we proposed a fully automated web-based analysis service, MARSweb. By providing a web service, we minimized the effort required for initial configuration. Additionally, users can perform analyses by simply uploading their data without needing to familiarize themselves with intricate analysis procedures. Furthermore, it facilitates easier interpretation of results by integrating advanced visualization tools. We confirmed the performance of MARSweb by detecting eGenes and performing pathway analysis of the genes using a Yeast Dataset.

CONCLUSIONS

MARSweb is a web-based analysis service that fully automates set-based analysis. It offers an intuitive user interface, making complex analyses more accessible while significantly reducing processing time for enhanced efficiency. MARSweb is available for use at http://cblab.dongguk.edu/MARSweb and its source code is available at https://github.com/DGU-CBLAB/MARSweb .

Genome wide association study reveals novel associations with face morphology

Genome-wide association studies (GWAS) on the Middle Eastern population, including the United Arab Emirates (UAE), have been relatively limited. The present study aims to investigate genotype-face morphology associations in the UAE population through Genome Wide Association Studies (GWAS). Phenotypic data (44 face measurements) from 172 Emiratis was obtained through three-dimensional (3D) scanning technology and an automatic face landmarking technique. GWAS analysis revealed associations of 19 genetic loci with six face features, 14 of which are novel. The GWAS analysis revealed 11 significant relationships between 44 face parameters and 242 SNPs, exceeding the GWAS significance threshold. These phenotypes were previously associated with body height, craniofacial defects, and facial characters. The most significant associations of these genetic variations were related to six main facial features which were facial convexity, left orbital protrusion, mandibular contour, nasolabial angle D, inferior facial angle B, and inferior facial angle A. To the best of our knowledge, this is the first GWAS study to investigate the association of SNP variations with face morphology in the Middle Eastern population.

An activating calcium-sensing receptor variant with biased signaling reveals a critical residue for Gα11 coupling

Autosomal dominant hypocalcemia (ADH) is due to enhanced calcium-dependent signaling caused by heterozygous gain-of-function (GOF) variants in the CASR gene (ADH1) or in the GNA11 gene, encoding Gα11 (ADH2). Both ADH1 and ADH2 are associated with hypocalcemia and normal or inappropriately low levels of circulating PTH. ADH1 patients typically manifest hypercalciuria, while ADH2 is associated with short stature in approximately 42% of cases. We evaluated a 10-yr-old boy with hypoparathyroidism and short stature. Biochemical analyses revealed hypocalcemia, hyperphosphatemia, and inconsistent hypercalciuria. Genetic analyses revealed a de novo heterozygous p.Leu723Arg variant in CASR. We characterized the expression of recombinant WT and Leu723Arg calcium-sensing receptor (CaSR) proteins in HEK293 cells and assessed G protein activation in vitro by CaSR using bioluminescence resonance energy transfer. Transient expression studies showed the Leu723Arg variant was normally expressed but resulted in a significantly lower EC50 for extracellular calcium activation of G11 but not other G proteins (ie, Gi, Gq, Gs). The Leu723Arg substitution has a novel GOF phenotype that leads to biased CaSR activation of G11 signaling, suggesting that residue 723 specifies activation of G11 but not other G proteins. Similar studies of a previously described CaSR variant associated with hypoparathyroidism and short stature, Leu616Val, showed no changes in any G protein pathways, indicating it is likely to be a benign variant. Given the preferential activation of G11 by the Leu723Arg CaSR variant, we propose that the patient's short stature shares a similar basis to that in patients with ADH2 due to GOF variants in GNA11.

Beneficial death: A substantial element of evolution?

If a shortened lifespan is evolutionarily advantageous, it becomes more likely that nature will strive to change it accordingly, affecting how we understand aging. Premature mortality because of aging would seem detrimental to the individual, but under what circumstances can it be of value? Based on a relative incremental increase in fitness, simulations were performed to reveal the benefit of death. This modification allows for continuous evolution in the model and establishes an optimal lifespan even under challenging conditions. As a result, shorter-lived individuals achieve faster adaptation through more frequent generational turnover, displacing longer-lived ones and likely providing a competitive advantage between species. Contrary to previous assumptions, this work proposes a mechanism by which early death, e.g., due to aging, may contribute to evolution.

Functional genomics of human skeletal development and the patterning of height heritability

Underlying variation in height are regulatory changes to chondrocytes, cartilage cells comprising long-bone growth plates. Currently, we lack knowledge on epigenetic regulation and gene expression of chondrocytes sampled across the human skeleton, and therefore we cannot understand basic regulatory mechanisms controlling height biology. We first rectify this issue by generating extensive epigenetic and transcriptomic maps from chondrocytes sampled from different growth plates across developing human skeletons, discovering novel regulatory networks shaping human bone/joint development. Next, using these maps in tandem with height genome-wide association study (GWAS) signals, we disentangle the regulatory impacts that skeletal element-specific versus global-acting variants have on skeletal growth, revealing the prime importance of regulatory pleiotropy in controlling height variation. Finally, as height is highly heritable, and thus often the test case for complex-trait genetics, we leverage these datasets within a testable omnigenic model framework to discover novel chondrocyte developmental modules and peripheral-acting factors shaping height biology and skeletal growth.

Functional genomics implicates natural killer cells in the pathogenesis of ankylosing spondylitis

Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased functional genomics approach. We integrated genome-wide association study (GWAS) data with epigenomic and transcriptomic datasets of human immune cells. To quantify enrichment of cell type-specific open chromatin or gene expression in AS risk loci, we used three published methods-LDSC-SEG, SNPsea, and scDRS-that have successfully identified relevant cell types in other diseases. Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA sequencing data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Using the human Space-Time Gut Cell Atlas, we also found significant upregulation of AS-associated genes predominantly in NK cells. We performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes. This revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two understudied loci, ENTR1 (SDCCAG3) and B3GNT2. Our findings suggest that NK cells may play a crucial role in AS development and highlight four putative target genes for functional follow-up in NK cells.

Variants in WASHC3, a component of the WASH complex, cause short stature, variable neurodevelopmental abnormalities, and distinctive facial dysmorphism

Purpose

Genetic defects that impair growth plate chondrogenesis cause a phenotype that varies from skeletal dysplasia to mild short stature with or without other syndromic features. In many individuals with impaired skeletal growth, the genetic causes remain unknown.

Method

Exome sequence was performed in 3 unrelated families with short stature, distinctive facies, and neurodevelopmental abnormalities. The impact of identified variants was studied in vitro.

Results

Exome sequencing identified variants in WASHC3, a component of the WASH complex. In the first family, a de-novo-dominant missense variant (p.L69F) impaired WASHC3 participation in the WASH complex, altered PTH1R endosomal trafficking, diminished PTH1R signaling, and affected growth plate chondrocyte hypertrophic differentiation, providing a likely explanation for the short stature. Knockdown of other WASH complex components also diminished PTH1R signaling. In the second and third families, a homozygous variant in the start codon (p.M1?) markedly reduced WASHC3 protein expression.

Conclusion

In combination with prior studies of WASH complex proteins, our findings provide evidence that the WASH complex is required for normal skeletal growth and that, consequently, genetic abnormalities impairing the function of the WASH complex (WASHopathy) cause short stature, as well as distinctive facies and variable neurodevelopmental abnormalities.

Syndrome-informed phenotyping identifies a polygenic background for achondroplasia-like facial variation in the general population

Human craniofacial shape is highly variable yet highly heritable with numerous genetic variants interacting through multiple layers of development. Here, we hypothesize that Mendelian phenotypes represent the extremes of a phenotypic spectrum and, using achondroplasia as an example, we introduce a syndrome-informed phenotyping approach to identify genomic loci associated with achondroplasia-like facial variation in the general population. We compare three-dimensional facial scans from 43 individuals with achondroplasia and 8246 controls to calculate achondroplasia-like facial scores. Multivariate GWAS of the control scores reveals a polygenic basis for facial variation along an achondroplasia-specific shape axis, identifying genes primarily involved in skeletal development. Jointly modeling these genes in two independent control samples, both human and mouse, shows craniofacial effects approximating the characteristic achondroplasia phenotype. These findings suggest that both complex and Mendelian genetic variation act on the same developmentally determined axes of facial variation, providing insights into the genetic intersection of complex traits and Mendelian disorders.

M2 Macrophages are Major Mediators of Germline Risk of Endometriosis and Explain Pleiotropy with Comorbid Traits

Endometriosis is a common gynecologic condition that causes chronic life-altering symptoms including pain, infertility, and elevated cancer risk. There is an urgent need for new non-hormonal targeted therapeutics to treat endometriosis, but until very recently, the cellular and molecular signatures of endometriotic lesions were undefined, severely hindering the development of clinical advances. Integrating inherited risk data from analyses of >450,000 individuals with ∼350,000 single cell transcriptomes from 21 patients, we uncover M2-macrophages as candidate drivers of disease susceptibility, and nominate IL1 signaling as a central hub impacted by germline genetic variation associated with endometriosis. Extensive functional follow-up confirmed these associations and revealed a pleiotropic role for this pathway in endometriosis. Population-scale expression quantitative trail locus analysis demonstrated that genetic variation controlling IL1A expression is also associated with endometriosis risk variants. Manipulation of IL1 signaling in state-of-the-art in vitro decidualized assembloids impacted epithelial differentiation, and in an in vivo endometriosis model, treatment with anakinra (an interleukin-1 receptor antagonist) resulted in a significant, dose-dependent reduction in both spontaneous pain and evoked pain. Together these studies highlight non-diagnostic cell types as central to endometriosis susceptibility and support IL1 signaling as an important actionable pathway for this disease.

Genomic prediction and genome-wide association studies of morphological traits and distraction index in Korean Sapsaree dogs

The Korean Sapsaree dog is a native breed known for its distinctive appearance and historical significance in Korean culture. The accurate estimation of breeding values is essential for the genetic improvement and conservation of such indigenous breeds. This study aimed to evaluate the accuracy of breeding values for body height, body length, chest width, hair length, and distraction index (DI) traits in Korean Sapsaree dogs. Additionally, a genome-wide association study (GWAS) was conducted to identify the genomic regions and nearby candidate genes influencing these traits. Phenotypic data were collected from 378 Korean Sapsaree dogs, and of these, 234 individuals were genotyped using the 170k Illumina CanineHD BeadChip. The accuracy of genomic predictions was evaluated using the traditional BLUP method with phenotypes only on genotyped animals (PBLUP-G), another traditional BLUP method using a pedigree-based relationship matrix (PBLUP) for all individuals, a GBLUP method based on a genomic relationship matrix, and a single-step GBLUP (ssGBLUP) method. Heritability estimates for body height, body length, chest width, hair length, and DI were 0.45, 0.39, 0.32, 0.55, and 0.50, respectively. Accuracy values varied across methods, with ranges of 0.22 to 0.31 for PBLUP-G, 0.30 to 0.57 for PBLUP, 0.31 to 0.54 for GBLUP, and 0.39 to 0.67 for ssGBLUP. Through GWAS, 194 genome-wide significant SNPs associated with studied Sapsaree traits were identified. The selection of the most promising candidate genes was based on gene ontology (GO) terms and functions previously identified to influence traits. Notable genes included CCKAR and DCAF16 for body height, PDZRN3 and CNTN1 for body length, TRIM63, KDELR2, and SUPT3H for chest width, RSPO2, EIF3E, PKHD1L1, TRPS1, and EXT1 for hair length, and DDHD1, BMP4, SEMA3C, and FOXP1 for the DI. These findings suggest that significant QTL, combined with functional candidate genes, can be leveraged to improve the genetic quality of the Sapsaree population. This study provides a foundation for more effective breeding strategies aimed at preserving and enhancing the unique traits of this Korean dog breed.

Exome-wide comparative analyses revealed differentiating genomic regions for performance traits in Indian native buffaloes

In India, 20 breeds of buffalo have been identified and registered, yet limited studies have been conducted to explore the performance potential of these breeds, especially in the Indian native breeds. This study is a maiden attempt to delineate the important variants and unique genes through exome sequencing for milk yield, milk composition, fertility, and adaptation traits in Indian local breeds of buffalo. In the present study, whole exome sequencing was performed on Chhattisgarhi (n = 3), Chilika (n = 4), Gojri (n = 3), and Murrah (n = 4) buffalo breeds and after stringent quality control, 4333, 6829, 4130, and 4854 InDels were revealed, respectively. Exome-wide FST along 100-kb sliding windows detected 27, 98, 38, and 35 outlier windows in Chhattisgarhi, Chilika, Gojri, and Murrah, respectively. The comparative exome analysis of InDels and subsequent gene ontology revealed unique breed specific genes for milk yield (CAMSAP3), milk composition (CLCN1, NUDT3), fertility (PTGER3) and adaptation (KCNA3, TH) traits. Study provides insight into mechanism of how these breeds have evolved under natural selection, the impact of these events on their respective genomes, and their importance in maintaining purity of these breeds for the traits under study. Additionally, this result will underwrite to the genetic acquaintance of these breeds for breeding application, and in understanding of evolution of these Indian local breeds.

Methodologies underpinning polygenic risk scores estimation: a comprehensive overview

Polygenic risk scores (PRS) have emerged as a promising tool for predicting disease risk and treatment outcomes using genomic data. Thousands of genome-wide association studies (GWAS), primarily involving populations of European ancestry, have supported the development of PRS models. However, these models have not been adequately evaluated in non-European populations, raising concerns about their clinical validity and predictive power across diverse groups. Addressing this issue requires developing novel risk prediction frameworks that leverage genetic characteristics across diverse populations, considering host-microbiome interactions and a broad range of health measures. One of the key aspects in evaluating PRS is understanding the strengths and limitations of various methods for constructing them. In this review, we analyze strengths and limitations of different methods for constructing PRS, including traditional weighted approaches and new methods such as Bayesian and Frequentist penalized regression approaches. Finally, we summarize recent advances in PRS calculation methods development, and highlight key areas for future research, including development of models robust across diverse populations by underlining the complex interplay between genetic variants across diverse ancestral backgrounds in disease risk as well as treatment response prediction. PRS hold great promise for improving disease risk prediction and personalized medicine; therefore, their implementation must be guided by careful consideration of their limitations, biases, and ethical implications to ensure that they are used in a fair, equitable, and responsible manner.

Analyses of GWAS signal using GRIN identify additional genes contributing to suicidal behavior

Genome-wide association studies (GWAS) identify genetic variants underlying complex traits but are limited by stringent genome-wide significance thresholds. We present GRIN (Gene set Refinement through Interacting Networks), which increases confidence in the expanded gene set by retaining genes strongly connected by biological networks when GWAS thresholds are relaxed. GRIN was validated on both simulated interrelated gene sets as well as multiple GWAS traits. From multiple GWAS summary statistics of suicide attempt, a complex phenotype, GRIN identified additional genes that replicated across independent cohorts and retained biologically interrelated genes despite a relaxed significance threshold. We present a conceptual model of how these retained genes interact through neurobiological pathways that may influence suicidal behavior, and identify existing drugs associated with these pathways that would not have been identified under traditional GWAS thresholds. We demonstrate GRIN's utility in boosting GWAS results by increasing the number of true positive genes identified from GWAS results.

Divergent gene expression patterns in alcohol and opioid use disorders lead to consistent alterations in functional networks within the dorsolateral prefrontal cortex

Substance Use Disorders (SUDs) manifest as persistent drug-seeking behavior despite adverse consequences, with Alcohol Use Disorder (AUD) and Opioid Use Disorder (OUD) representing prevalent forms associated with significant mortality rates and economic burdens. The co-occurrence of AUD and OUD is common, necessitating a deeper comprehension of their intricate interactions. While the causal link between these disorders remains elusive, shared genetic factors are hypothesized. Leveraging public datasets, we employed genomic and transcriptomic analyses to explore conserved and distinct molecular pathways within the dorsolateral prefrontal cortex associated with AUD and OUD. Our findings unveil modest transcriptomic overlap at the gene level between the two disorders but substantial convergence on shared biological pathways. Notably, these pathways predominantly involve inflammatory processes, synaptic plasticity, and key intracellular signaling regulators. Integration of transcriptomic data with the latest genome-wide association studies (GWAS) for problematic alcohol use (PAU) and OUD not only corroborated our transcriptomic findings but also confirmed the limited shared heritability between the disorders. Overall, our study indicates that while alcohol and opioids induce diverse transcriptional alterations at the gene level, they converge on select biological pathways, offering promising avenues for novel therapeutic targets aimed at addressing both disorders simultaneously.

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.

Genome-Wide Association Study on Body Conformation Traits in Xinjiang Brown Cattle

Body conformation traits are linked to the health, longevity, reproductivity, and production performance of cattle. These traits are also crucial for herd selection and developing new breeds. This study utilized pedigree information and phenotypic (1185 records) and genomic (The resequencing of 496 Xinjiang Brown cattle generated approximately 74.9 billion reads.) data of Xinjiang Brown cattle to estimate the genetic parameters, perform factor analysis, and conduct a genome-wide association study (GWAS) for these traits. Our results indicated that most traits exhibit moderate to high heritability. The principal factors, which explained 59.12% of the total variance, effectively represented body frame, muscularity, rump, feet and legs, and mammary system traits. Their heritability estimates range from 0.17 to 0.73, with genetic correlations ranging from -0.53 to 0.33. The GWAS identified 102 significant SNPs associated with 12 body conformation traits. A few of the SNPs were located near previously reported genes and quantitative trait loci (QTLs), while others were novel. The key candidate genes such as LCORL, NCAPG, and FAM184B were annotated within 500 Kb upstream and downstream of the significant SNPs. Therefore, factor analysis can be used to simplify multidimensional conformation traits into new variables, thus reducing the computational burden. The identified candidate genes from GWAS can be incorporated into the genomic selection of Xinjiang Brown cattle, enhancing the reliability of breeding programs.

QSOX2 Deficiency-induced short stature, gastrointestinal dysmotility and immune dysfunction

Postnatal growth failure is often attributed to dysregulated somatotropin action, however marked genetic and phenotypic heterogeneity exist. We report five patients from three families who present with short stature, immune dysfunction, atopic eczema and gastrointestinal pathology associated with recessive variants in QSOX2. QSOX2 encodes a nuclear membrane protein linked to disulphide isomerase and oxidoreductase activity. Loss of QSOX2 disrupts Growth hormone-mediated STAT5B nuclear translocation despite enhanced Growth hormone-induced STAT5B phosphorylation. Moreover, patient-derived dermal fibroblasts demonstrate Growth hormone-induced mitochondriopathy and reduced mitochondrial membrane potential. Located at the nuclear membrane, QSOX2 acts as a gatekeeper for regulating stabilisation and import of phosphorylated-STAT5B. Altogether, QSOX2 deficiency modulates human growth by impairing Growth hormone-STAT5B downstream activities and mitochondrial dynamics, which contribute to multi-system dysfunction. Furthermore, our work suggests that therapeutic recombinant insulin-like growth factor-1 may circumvent the Growth hormone-STAT5B dysregulation induced by pathological QSOX2 variants and potentially alleviate organ specific disease.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic risk to childhood-onset asthma

Asthma is a complex disease caused by genetic and environmental factors. Studies show that wheezing during rhinovirus infection correlates with childhood asthma development. Over 150 non-coding risk variants for asthma have been identified, many affecting gene regulation in T cells, but the effects of most risk variants remain unknown. We hypothesized that airway epithelial cells could also mediate genetic susceptibility to asthma given they are the first line of defense against respiratory viruses and allergens. We integrated genetic data with transcriptomics of airway epithelial cells subject to different stimuli. We demonstrate that rhinovirus infection significantly upregulates childhood-onset asthma-associated genes, particularly in non-ciliated cells. This enrichment is also observed with influenza infection but not with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or cytokine activation. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Functional Characterization of Accessible Chromatin in Common Wheat

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

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

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

Natural variation in starvation sensitivity maps to a point mutation in phospholipase IPLA2-VIA in Drosophila melanogaster

Resistance to starvation is a classic complex trait, where genetic and environmental variables can greatly modify an animal's ability to survive without nutrients. In this study, we investigate the genetic basis of starvation resistance using complementary quantitative and classical genetic mapping in Drosophila melanogaster. Using the Drosophila Genetics Reference Panel (DGRP) as a starting point, we queried the genetic basis of starvation sensitivity in one of the most sensitive DGRP lines. We localize a major effect locus modifying starvation resistance to the phospholipase iPLA2-VIA. This finding is consistent with the work of others which demonstrate the importance of lipid regulation in starvation stress. Furthermore, we demonstrate that iPLA2-VIA plays a role in the maintenance of sugar reserves post-starvation, which highlights a key dynamic between lipid remodeling, sugar metabolism and resistance to starvation stress.

You are What Your Parents Expect: Height and Local Reference Points

Recent estimates are that about 150 million children under five years of age are stunted, with substantial negative consequences for their schooling, cognitive skills, health, and economic productivity. Therefore, understanding what determines such growth retardation is significant for designing public policies that aim to address this issue. We build a model for nutritional choices and health with reference-dependent preferences. Parents care about the health of their children relative to some reference population. In our empirical model, we use height as the health outcome that parents target. Reference height is an equilibrium object determined by earlier cohorts' parents' nutritional choices in the same village. We explore the exogenous variation in reference height produced by a protein-supplementation experiment in Guatemala to estimate our model's parameters. We use our model to decompose the impact of the protein intervention on height into price and reference-point effects. We find that the changes in reference points account for 65% of the height difference between two-year-old children in experimental and control villages in the sixth annual cohort born after the initiation of the intervention.

DeepRisk: A deep learning approach for genome-wide assessment of common disease risk

The potential for being able to identify individuals at high disease risk solely based on genotype data has garnered significant interest. Although widely applied, traditional polygenic risk scoring methods fall short, as they are built on additive models that fail to capture the intricate associations among single nucleotide polymorphisms (SNPs). This presents a limitation, as genetic diseases often arise from complex interactions between multiple SNPs. To address this challenge, we developed DeepRisk, a biological knowledge-driven deep learning method for modeling these complex, nonlinear associations among SNPs, to provide a more effective method for scoring the risk of common diseases with genome-wide genotype data. Evaluations demonstrated that DeepRisk outperforms existing PRS-based methods in identifying individuals at high risk for four common diseases: Alzheimer's disease, inflammatory bowel disease, type 2 diabetes, and breast cancer.

Testing for differences in polygenic scores in the presence of confounding

Polygenic scores have become an important tool in human genetics, enabling the prediction of individuals' phenotypes from their genotypes. Understanding how the pattern of differences in polygenic score predictions across individuals intersects with variation in ancestry can provide insights into the evolutionary forces acting on the trait in question, and is important for understanding health disparities. However, because most polygenic scores are computed using effect estimates from population samples, they are susceptible to confounding by both genetic and environmental effects that are correlated with ancestry. The extent to which this confounding drives patterns in the distribution of polygenic scores depends on patterns of population structure in both the original estimation panel and in the prediction/test panel. Here, we use theory from population and statistical genetics, together with simulations, to study the procedure of testing for an association between polygenic scores and axes of ancestry variation in the presence of confounding. We use a general model of genetic relatedness to describe how confounding in the estimation panel biases the distribution of polygenic scores in a way that depends on the degree of overlap in population structure between panels. We then show how this confounding can bias tests for associations between polygenic scores and important axes of ancestry variation in the test panel. Specifically, for any given test, there exists a single axis of population structure in the GWAS panel that needs to be controlled for in order to protect the test. Based on this result, we propose a new approach for directly estimating this axis of population structure in the GWAS panel. We then use simulations to compare the performance of this approach to the standard approach in which the principal components of the GWAS panel genotypes are used to control for stratification.

Functional genomics implicates natural killer cells in the pathogenesis of ankylosing spondylitis

Objective

Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased functional genomics approach.

Methods

We integrated genome-wide association study (GWAS) data with epigenomic and transcriptomic datasets of human immune cells. To quantify enrichment of cell type-specific open chromatin or gene expression in AS risk loci, we used three published methods that have successfully identified relevant cell types in other diseases. We performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes.

Results

Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA-seq data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Using the human Space-Time Gut Cell Atlas, we also found significant upregulation of AS-associated genes predominantly in NK cells. Co-localization analysis revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two under-studied loci, ENTR1 (aka SDCCAG3) and B3GNT2.

Conclusion

Our findings suggest that NK cells may play a crucial role in AS development and highlight four putative target genes for functional follow-up in NK cells.

Divergent gene expression patterns in alcohol and opioid use disorders lead to consistent alterations in functional networks within the Dorsolateral Prefrontal Cortex

Substance Use Disorders (SUDs) manifest as persistent drug-seeking behavior despite adverse consequences, with Alcohol Use Disorder (AUD) and Opioid Use Disorder (OUD) representing prevalent forms associated with significant mortality rates and economic burdens. The co-occurrence of AUD and OUD is common, necessitating a deeper comprehension of their intricate interactions. While the causal link between these disorders remains elusive, shared genetic factors are hypothesized. Leveraging public datasets, we employed genomic and transcriptomic analyses to explore conserved and distinct molecular pathways within the dorsolateral prefrontal cortex associated with AUD and OUD. Our findings unveil modest transcriptomic overlap at the gene level between the two disorders but substantial convergence on shared biological pathways. Notably, these pathways predominantly involve inflammatory processes, synaptic plasticity, and key intracellular signaling regulators. Integration of transcriptomic data with the latest genome-wide association studies (GWAS) for problematic alcohol use (PAU) and OUD not only corroborated our transcriptomic findings but also confirmed the limited shared heritability between the disorders. Overall, our study indicates that while alcohol and opioids induce diverse transcriptional alterations at the gene level, they converge on select biological pathways, offering promising avenues for novel therapeutic targets aimed at addressing both disorders simultaneously.

Polymorphisms of CYP7A1 and HADHB Genes and Their Effects on Milk Production Traits in Chinese Holstein Cows

Our preliminary research proposed the cytochrome P450 family 7 subfamily A member 1 (CYP7A1) and hydroxyacyl-coenzyme A dehydrogenase trifunctional multienzyme complex beta subunit (HADHB) genes as candidates for association with milk-production traits in dairy cattle because of their differential expression across different lactation stages in the liver tissues of Chinese Holstein cows and their potential roles in lipid metabolism. Hence, we identified single-nucleotide polymorphisms (SNPs) of the CYP7A1 and HADHB genes and validated their genetic effects on milk-production traits in a Chinese Holstein population with the goal of providing valuable genetic markers for genomic selection (GS) in dairy cattle, This study identified five SNPs, 14:g.24676921A>G, 14:g.24676224G>A, 14:g.24675708G>T, 14:g.24665961C>T, and 14:g.24664026A>G, in the CYP7A1 gene and three SNPs, 11:g.73256269T>C, 11:g.73256227A>C, and 11:g.73242290C>T, in HADHB. The single-SNP association analysis revealed significant associations (p value ≤ 0.0461) between the eight SNPs of CYP7A1 and HADHB genes and 305-day milk, fat and protein yields. Additionally, using Haploview 4.2, we found that the five SNPs of CYP7A1 formed two haplotype blocks and that the two SNPs of HADHB formed one haplotype block; notably, all three haplotype blocks were also significantly associated with milk, fat and protein yields (p value ≤ 0.0315). Further prediction of transcription factor binding sites (TFBSs) based on Jaspar software (version 2023) showed that the 14:g.24676921A>G, 14:g.24675708G>T, 11:g.73256269T>C, and 11:g.73256227A>C SNPs could alter the 5' terminal TFBS of the CYP7A1 and HADHB genes. The 14:g.24665961C>T SNP caused changes in the structural stability of the mRNA for the CYP7A1 gene. These alterations have the potential to influence gene expression and, consequently, the phenotype associated with milk-production traits. In summary, we have confirmed the genetic effects of CYP7A1 and HADHB genes on milk-production traits in dairy cattle and identified potential functional mutations that we suggest could be used for GS of dairy cattle and in-depth mechanistic studies of animals.

Estimating the overall fraction of phenotypic variance attributed to high-dimensional predictors measured with error

In prospective genomic studies (e.g., DNA methylation, metagenomics, and transcriptomics), it is crucial to estimate the overall fraction of phenotypic variance (OFPV) attributed to the high-dimensional genomic variables, a concept similar to heritability analyses in genome-wide association studies (GWAS). Unlike genetic variants in GWAS, these genomic variables are typically measured with error due to technical limitation and temporal instability. While the existing methods developed for GWAS can be used, ignoring measurement error may severely underestimate OFPV and mislead the design of future studies. Assuming that measurement error variances are distributed similarly between causal and noncausal variables, we show that the asymptotic attenuation factor equals to the average intraclass correlation coefficients of all genomic variables, which can be estimated based on a pilot study with repeated measurements. We illustrate the method by estimating the contribution of microbiome taxa to body mass index and multiple allergy traits in the American Gut Project. Finally, we show that measurement error does not cause meaningful bias when estimating the correlation of effect sizes for two traits.

Trans-eQTL mapping in gene sets identifies network effects of genetic variants

Nearly all trait-associated variants identified in genome-wide association studies (GWASs) are noncoding. The cis regulatory effects of these variants have been extensively characterized, but how they affect gene regulation in trans has been the subject of fewer studies because of the difficulty in detecting trans-expression quantitative loci (eQTLs). We developed trans-PCO for detecting trans effects of genetic variants on gene networks. Our simulations demonstrate that trans-PCO substantially outperforms existing trans-eQTL mapping methods. We applied trans-PCO to two gene expression datasets from whole blood, DGN (N = 913) and eQTLGen (N = 31,684), and identified 14,985 high-quality trans-eSNP-module pairs associated with 197 co-expression gene modules and biological processes. We performed colocalization analyses between GWAS loci of 46 complex traits and the trans-eQTLs. We demonstrated that the identified trans effects can help us understand how trait-associated variants affect gene regulatory networks and biological pathways.

Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits

BACKGROUND

Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected.

RESULTS

All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count.

CONCLUSIONS

Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Epigenetic regulators controlling osteogenic lineage commitment and bone formation

Bone formation and homeostasis are controlled by environmental factors and endocrine regulatory cues that initiate intracellular signaling pathways capable of modulating gene expression in the nucleus. Bone-related gene expression is controlled by nucleosome-based chromatin architecture that limits the accessibility of lineage-specific gene regulatory DNA sequences and sequence-specific transcription factors. From a developmental perspective, bone-specific gene expression must be suppressed during the early stages of embryogenesis to prevent the premature mineralization of skeletal elements during fetal growth in utero. Hence, bone formation is initially inhibited by gene suppressive epigenetic regulators, while other epigenetic regulators actively support osteoblast differentiation. Prominent epigenetic regulators that stimulate or attenuate osteogenesis include lysine methyl transferases (e.g., EZH2, SMYD2, SUV420H2), lysine deacetylases (e.g., HDAC1, HDAC3, HDAC4, HDAC7, SIRT1, SIRT3), arginine methyl transferases (e.g., PRMT1, PRMT4/CARM1, PRMT5), dioxygenases (e.g., TET2), bromodomain proteins (e.g., BRD2, BRD4) and chromodomain proteins (e.g., CBX1, CBX2, CBX5). This narrative review provides a broad overview of the covalent modifications of DNA and histone proteins that involve hundreds of enzymes that add, read, or delete these epigenetic modifications that are relevant for self-renewal and differentiation of mesenchymal stem cells, skeletal stem cells and osteoblasts during osteogenesis.

Advancements in elucidating the pathogenesis of actinic keratosis: present state and future prospects

Solar keratosis, also known as actinic keratosis (AK), is becoming increasingly prevalent. It is a benign tumor that develops in the epidermis. Individuals with AK typically exhibit irregular, red, scaly bumps or patches as a result of prolonged exposure to UV rays. These growths primarily appear on sun-exposed areas of the skin such as the face, scalp, and hands. Presently, dermatologists are actively studying AK due to its rising incidence rate in the United States. However, the underlying causes of AK remain poorly understood. Previous research has indicated that the onset of AK involves various mechanisms including UV ray-induced inflammation, oxidative stress, complex mutagenesis, resulting immunosuppression, inhibited apoptosis, dysregulated cell cycle, altered cell proliferation, tissue remodeling, and human papillomavirus (HPV) infection. AK can develop in three ways: spontaneous regression, persistence, or progression into invasive cutaneous squamous cell carcinoma (cSCC). Multiple risk factors and diverse signaling pathways collectively contribute to its complex pathogenesis. To mitigate the risk of cancerous changes associated with long-term UV radiation exposure, prompt identification, management, and prevention of AK are crucial. The objective of this review is to elucidate the primary mechanisms underlying AK malignancy and identify potential treatment targets for dermatologists in clinical settings.

Exploration of the Polymorphism Distribution of Bovine HMGA2 Gene in Worldwide Breeds and Its Associations with Ovarian Traits

The high-mobility group AT-hook 2(HMGA2) gene has been widely studied in the context of cancer and animal growth. However, recently, several studies have uncovered its critical role in cell proliferation. A genome-wide association study (GWAS) further suggests that the HMGA2 gene is a candidate gene in fertility, indicating its connection not only to growth traits but also to reproduction, specifically ovarian traits. Thus, this study aimed to analyze the distribution of the HMGA2 gene in 54 bovine breeds worldwide, identify important short fragment variants (indels), and investigate the relationship between HMGA2 and ovarian development. The dataset included genotypic information from a bovine population of 634 individuals (n = 634). After genotyping and analyzing four selected loci, we found that one out of four loci, rs133750033 (P4-D22-bp), was polymorphic. Our results also reveal that this indel of HMGA2 is significantly associated with certain ovarian traits (p < 0.05). Specifically, it has connection with ovarian length (p = 0.004) and ovarian height (p = 0.026) during diestrus. Additionally, we discovered a higher expression of the HMGA2 gene in Asian cattle breeds. In summary, this study suggests that HMGA2 has the potential to serve as an animal fertility testing marker gene. Moreover, these findings contribute to a more promising outlook for the bovine industry.

Tall stature and gigantism in adult patients with acromegaly

OBJECTIVES

Increased height in patients with acromegaly could be a manifestation of growth hormone (GH) excess before epiphysis closure. The aim of this study was to evaluate the relationship between the height of adult patients with GH excess related to mid-parental height (MPH) and population mean and to find whether taller patients with acromegaly come from tall families.

METHODS

This is a single-centre, observational study involving 135 consecutive patients with acromegaly diagnosed as adults and no family history of GH excess. We established three categories for height for patients with acromegaly: normal stature, tall stature (TS, height above the 97th percentile (1.88 standard deviations (SD)) to <3 SD for gender- and country-specific data or as a height which was greater than 1.5 SD but less than 2 SD above the MPH) and gigantism (height which was greater than 3 SD) above the gender- and country-specific mean or greater than 2 SD above MPH).

RESULTS

Thirteen percent (17/135) of patients (53% females) met the criteria for gigantism, 10% (14/135) fulfilled the criteria for TS (57% females). Parents and adult siblings were not taller than the population mean.

CONCLUSION

In a group of 135 consecutive adult patients with acromegaly, 23% had increased height based on country-specific and MPH data: 13% presented with gigantism while 10% had TS. The frequency of gigantism and TS in patients diagnosed with GH excess as adults is not higher in males than in females. Patients with acromegaly come from normal-stature families.

Machine learning nominates the inositol pathway and novel genes in Parkinson's disease

There are 78 loci associated with Parkinson's disease in the most recent genome-wide association study (GWAS), yet the specific genes driving these associations are mostly unknown. Herein, we aimed to nominate the top candidate gene from each Parkinson's disease locus and identify variants and pathways potentially involved in Parkinson's disease. We trained a machine learning model to predict Parkinson's disease-associated genes from GWAS loci using genomic, transcriptomic and epigenomic data from brain tissues and dopaminergic neurons. We nominated candidate genes in each locus and identified novel pathways potentially involved in Parkinson's disease, such as the inositol phosphate biosynthetic pathway (INPP5F, IP6K2, ITPKB and PPIP5K2). Specific common coding variants in SPNS1 and MLX may be involved in Parkinson's disease, and burden tests of rare variants further support that CNIP3, LSM7, NUCKS1 and the polyol/inositol phosphate biosynthetic pathway are associated with the disease. Functional studies are needed to further analyse the involvements of these genes and pathways in Parkinson's disease.

Characterization of HMGA2 variants expands the spectrum of Silver-Russell syndrome

Silver-Russell syndrome (SRS) is a heterogeneous disorder characterized by intrauterine and postnatal growth retardation. HMGA2 variants are a rare cause of SRS and its functional role in human linear growth is unclear. Patients with suspected SRS negative for 11p15LOM/mUPD7 underwent whole-exome and/or targeted-genome sequencing. Mutant HMGA2 protein expression and nuclear localization were assessed. Two Hmga2-knockin mouse models were generated. Five clinical SRS patients harbored HMGA2 variants with differing functional impacts: 2 stop-gain nonsense variants (c.49G>T, c.52C>T), c.166A>G missense variant, and 2 frameshift variants (c.144delC, c.145delA) leading to an identical, extended-length protein. Phenotypic features were highly variable. Nuclear localization was reduced/absent for all variants except c.166A>G. Homozygous knockin mice recapitulating the c.166A>G variant (Hmga2K56E) exhibited a growth-restricted phenotype. An Hmga2Ter76-knockin mouse model lacked detectable full-length Hmga2 protein, similarly to patient 3 and 5 variants. These mice were infertile, with a pygmy phenotype. We report a heterogeneous group of individuals with SRS harboring variants in HMGA2 and describe the first Hmga2 missense knockin mouse model (Hmga2K56E) to our knowledge causing a growth-restricted phenotype. In patients with clinical features of SRS but negative genetic screening, HMGA2 should be included in next-generation sequencing testing approaches.

Emerging roles for the ADAMTS-like family of matricellular proteins in cardiovascular disease through regulation of the extracellular microenvironment

Dysregulation of the extracellular matrix (ECM) occurs widely across cardiovascular pathologies. Recent work has revealed important roles for the «a disintegrin-like and metalloprotease domain with thrombospondin-type 1 motifs like" (ADAMTSL) family of secreted glycoproteins in cardiovascular tissues during development and disease. Key insights in this regard have come from naturally occurring gene mutations in humans and animals that result in severe diseases with cardiovascular manifestations or aortopathies. Expression of ADAMTSL genes is greatly increased in the myocardium during heart failure. Genetically modified mice recapitulate phenotypes of patients with ADAMTSL mutations and demonstrate important functions in the ECM. The novel functions thus disclosed are intriguing because, while these proteins are neither structural, nor proteases like the related ADAMTS proteases, they appear to act as regulatory, i.e., matricellular proteins. Evidence from genetic variants, genetically engineered mouse mutants, and in vitro investigations have revealed regulatory functions of ADAMTSLs related to fibrillin microfibrils and growth factor signaling. Interestingly, the ability to regulate transforming growth factor (TGF)β signaling may be a shared characteristic of some ADAMTSLs. TGFβ signaling is important in cardiovascular development, health and disease and a central driver of ECM remodeling and cardiac fibrosis. New strategies to target dysregulated TGFβ signaling are warranted in aortopathies and cardiac fibrosis. With their emerging roles in cardiovascular tissues, the ADAMTSL proteins may provide causative genes, diagnostic biomarkers and novel treatment targets in cardiovascular disease. Here, we discuss the relevance of ADAMTSLs to cardiovascular medicine.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic susceptibility to childhood-onset asthma

Asthma is a complex disease caused by genetic and environmental factors. Epidemiological studies have shown that in children, wheezing during rhinovirus infection (a cause of the common cold) is associated with asthma development during childhood. This has led scientists to hypothesize there could be a causal relationship between rhinovirus infection and asthma or that RV-induced wheezing identifies individuals at increased risk for asthma development. However, not all children who wheeze when they have a cold develop asthma. Genome-wide association studies (GWAS) have identified hundreds of genetic variants contributing to asthma susceptibility, with the vast majority of likely causal variants being non-coding. Integrative analyses with transcriptomic and epigenomic datasets have indicated that T cells drive asthma risk, which has been supported by mouse studies. However, the datasets ascertained in these integrative analyses lack airway epithelial cells. Furthermore, large-scale transcriptomic T cell studies have not identified the regulatory effects of most non-coding risk variants in asthma GWAS, indicating there could be additional cell types harboring these "missing regulatory effects". Given that airway epithelial cells are the first line of defense against rhinovirus, we hypothesized they could be mediators of genetic susceptibility to asthma. Here we integrate GWAS data with transcriptomic datasets of airway epithelial cells subject to stimuli that could induce activation states relevant to asthma. We demonstrate that epithelial cultures infected with rhinovirus significantly upregulate childhood-onset asthma-associated genes. We show that this upregulation occurs specifically in non-ciliated epithelial cells. This enrichment for genes in asthma risk loci, or 'asthma heritability enrichment' is also significant for epithelial genes upregulated with influenza infection, but not with SARS-CoV-2 infection or cytokine activation. Additionally, cells from patients with asthma showed a stronger heritability enrichment compared to cells from healthy individuals. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Genomic Determinants of Knee Joint Biomechanics: An Exploration into the Molecular Basis of Locomotor Function, a Narrative Review

In recent years, the nexus between genetics and biomechanics has garnered significant attention, elucidating the role of genomic determinants in shaping the biomechanical attributes of human joints, specifically the knee. This review seeks to provide a comprehensive exploration of the molecular basis underlying knee joint locomotor function. Leveraging advancements in genomic sequencing, we identified specific genetic markers and polymorphisms tied to key biomechanical features of the knee, such as ligament elasticity, meniscal resilience, and cartilage health. Particular attention was devoted to collagen genes like COL1A1 and COL5A1 and their influence on ligamentous strength and injury susceptibility. We further investigated the genetic underpinnings of knee osteoarthritis onset and progression, as well as the potential for personalized rehabilitation strategies tailored to an individual's genetic profile. We reviewed the impact of genetic factors on knee biomechanics and highlighted the importance of personalized orthopedic interventions. The results hold significant implications for injury prevention, treatment optimization, and the future of regenerative medicine, targeting not only knee joint health but joint health in general.

Pregnancy associated plasma protein-A2 (PAPP-A2) and stanniocalcin-2 (STC2) but not PAPP-A are associated with circulating total IGF-1 in a human adult population

The pappalysins pregnancy associated plasma protein-A (PAPP-A) and -A2 (PAPP-A2) act as proteinases of insulin-like growth factor-1 (IGF-1) binding proteins, while stanniocalcin-2 (STC2) was identified as a pappalysin inhibitor. While there is some evidence from studies in children and adolescents, it is unclear whether these molecules are related to concentrations of IGF-1 and its binding proteins in adults. We investigated cross-sectionally the association of circulating PAPP-A, PAPP-A2 and STC2 with IGF-1 and its binding proteins (IGFBPs) in 394 adult pretest participants (20-69 years) of the German National Cohort Berlin North study center. Plasma PAPP-A, PAPP-A2, total and free IGF-1, IGFBP-1, IGFBP-2, IGFBP-3, IGFBP-5 and STC2 were measured by ELISAs. The associations of PAPP-A, PAPP-A2 and STC2 with IGF-1 or IGFBPs were investigated using multivariable linear regression analyses adjusting for age, sex, body mass index and pretest phase. We observed significant inverse associations of PAPP-A2 (difference in concentrations per 0.5 ng/mL higher PAPP-A2 levels) with total IGF-1 (- 4.3 ng/mL; 95% CI - 7.0; - 1.6), the IGF-1:IGFBP-3 molar ratio (- 0.34%; 95%-CI - 0.59; - 0.09), but not free IGF-1 and a positive association with IGFBP-2 (11.9 ng/mL; 95% CI 5.0; 18.8). PAPP-A was not related to total or free IGF-1, but positively associated with IGFBP-5. STC2 was inversely related to total IGF-1, IGFBP-2 and IGFBP-3 and positively to IGFBP-1. This first investigation of these associations in a general adult population supports the hypothesis that PAPP-A2 as well as STC2 play a role for IGF-1 and its binding proteins, especially for total IGF-1. The role of PAPP-A2 and STC2 for health and disease in adults warrants further investigation.

Genome-wide association studies for economically important traits in mink using copy number variation

Copy number variations (CNVs) are structural variants consisting of duplications and deletions of DNA segments, which are known to play important roles in the genetics of complex traits in livestock species. However, CNV-based genome-wide association studies (GWAS) have remained unexplored in American mink. Therefore, the purpose of the current study was to investigate the association between CNVs and complex traits in American mink. A CNV-based GWAS was performed with the ParseCNV2 software program using deregressed estimated breeding values of 27 traits as pseudophenotypes, categorized into traits of growth and feed efficiency, reproduction, pelt quality, and Aleutian disease tests. The study identified a total of 10,137 CNVs (6968 duplications and 3169 deletions) using the Affymetrix Mink 70K single nucleotide polymorphism (SNP) array in 2986 American mink. The association analyses identified 250 CNV regions (CNVRs) associated with at least one of the studied traits. These CNVRs overlapped with a total of 320 potential candidate genes, and among them, several genes have been known to be related to the traits such as ARID1B, APPL1, TOX, and GPC5 (growth and feed efficiency traits); GRM1, RNASE10, WNT3, WNT3A, and WNT9B (reproduction traits); MYO10, and LIMS1 (pelt quality traits); and IFNGR2, APEX1, UBE3A, and STX11 (Aleutian disease tests). Overall, the results of the study provide potential candidate genes that may regulate economically important traits and therefore may be used as genetic markers in mink genomic breeding programs.

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.

Concurrent impact of de novo mutations on cranial and cortical development in nonsyndromic craniosynostosis

OBJECTIVE

Nonsyndromic craniosynostosis (nsCS), characterized by premature cranial suture fusion, is considered a primary skull disorder in which impact on neurodevelopment, if present, results from the mechanical hindrance of brain growth. Despite surgical repair of the cranial defect, neurocognitive deficits persist in nearly half of affected children. Therefore, the authors performed a functional genomics analysis of nsCS to determine when, where, and in what cell types nsCS-associated genes converge during development.

METHODS

The authors integrated whole-exome sequencing data from 291 nsCS proband-parent trios with 29,803 single-cell transcriptomes of the prenatal and postnatal neurocranial complex to inform when, where, and in what cell types nsCS-mutated genes might exert their pathophysiological effects.

RESULTS

The authors found that nsCS-mutated genes converged in cranial osteoprogenitors and pial fibroblasts and their transcriptional networks that regulate both skull ossification and cerebral neurogenesis. Nonsyndromic CS-mutated genes also converged in inhibitory neurons and gene coexpression modules that overlapped with autism and other developmental disorders. Ligand-receptor cell-cell communication analysis uncovered crosstalk between suture osteoblasts and neurons via the nsCS-associated BMP, FGF, and noncanonical WNT signaling pathways.

CONCLUSIONS

These data implicate a concurrent impact of nsCS-associated de novo mutations on cranial morphogenesis and cortical development via cell- and non-cell-autonomous mechanisms in a developmental nexus of fetal osteoblasts, pial fibroblasts, and neurons. These results suggest that neurodevelopmental outcomes in nsCS patients may be driven more by mutational status than surgical technique.

Syndrome-informed phenotyping identifies a polygenic background for achondroplasia-like facial variation in the general population

Human craniofacial shape is highly variable yet highly heritable with genetic variants interacting through multiple layers of development. Here, we hypothesize that Mendelian phenotypes represent the extremes of a phenotypic spectrum and, using achondroplasia as an example, we introduce a syndrome-informed phenotyping approach to identify genomic loci associated with achondroplasia-like facial variation in the normal population. We compared three-dimensional facial scans from 43 individuals with achondroplasia and 8246 controls to calculate achondroplasia-like facial scores. Multivariate GWAS of the control scores revealed a polygenic basis for normal facial variation along an achondroplasia-specific shape axis, identifying genes primarily involved in skeletal development. Jointly modeling these genes in two independent control samples showed craniofacial effects approximating the characteristic achondroplasia phenotype. These findings suggest that both complex and Mendelian genetic variation act on the same developmentally determined axes of facial variation, providing new insights into the genetic intersection of complex traits and Mendelian disorders.

High-throughput screening of glucocorticoid-induced enhancer activity reveals mechanisms of stress-related psychiatric disorders

Exposure to stressful life events increases the risk for psychiatric disorders. Mechanistic insight into the genetic factors moderating the impact of stress can increase our understanding of disease processes. Here, we test 3,662 single nucleotide polymorphisms (SNPs) from preselected expression quantitative trait loci in massively parallel reporter assays to identify genetic variants that modulate the activity of regulatory elements sensitive to glucocorticoids, important mediators of the stress response. Of the tested SNP sequences, 547 were located in glucocorticoid-responsive regulatory elements of which 233 showed allele-dependent activity. Transcripts regulated by these functional variants were enriched for those differentially expressed in psychiatric disorders in the postmortem brain. Phenome-wide Mendelian randomization analysis in 4,439 phenotypes revealed potentially causal associations specifically in neurobehavioral traits, including major depression and other psychiatric disorders. Finally, a functional gene score derived from these variants was significantly associated with differences in the physiological stress response, suggesting that these variants may alter disease risk by moderating the individual set point of the stress response.