Search

Abstract Bronchopulmonary dysplasia (BPD) is a common complication of preterm birth. Despite this, genetic drivers of BPD are poorly understood. The objective of this study is to better understand the impact of single nucleotide polymorphisms (SNPs) previously associated with BPD by examining associations with other phenotypes. We drew pediatric subjects from the biorepository at the Center for Applied Genomics to identify associations between these SNPs and 2,146 imputed phenotypes. Methylation data, external cohorts, and in silico validation methods were used to corroborate significant associations. We identified 60 SNPs that were previously associated with BPD. We found a significant association between rs3771150 and rs3771171 and mean eosinophil percentage in a European cohort of 6,999 patients and replicated this in external cohorts. Both SNPs were also associated with asthma, COPD and FEV1/FVC ratio. These SNPs displayed associations with methylation probes and were functionally linked to ST2 (IL1RL1) levels in blood and lung tissue. Our findings support a genetic justification for the epidemiological link between BPD and asthma. Given the well-established link between ST2 and type 2 inflammation in asthma, these findings provide a rationale for future studies exploring the role of type 2 inflammation in the pathogenesis of BPD.

Abstract This study aims to assess the effect of familial structures on the still-missing heritability estimate and prediction accuracy of Type 2 Diabetes (T2D) using pedigree estimated risk values (ERV) and genomic ERV. We used 11,818 individuals (T2D cases: 2,210) with genotype (649,932 SNPs) and pedigree information from the ongoing periodic cohort study of the Iranian population project. We considered three different familial structure scenarios, including (i) all families, (ii) all families with ≥ 1 generation, and (iii) families with ≥ 1 generation in which both case and control individuals are presented. Comprehensive simulation strategies were implemented to quantify the difference between estimates of $$\:{\text{h}}^{2}$$ and $$\:{\text{h}}_{\text{S}\text{N}\text{P}}^{2}$$ . A proportion of still-missing heritability in T2D could be explained by overestimation of pedigree-based heritability due to the presence of families with individuals having only one of the two disease statuses. Our research findings underscore the significance of including families with only case/control individuals in cohort studies. The presence of such family structures (as observed in scenarios i and ii) contributes to a more accurate estimation of disease heritability, addressing the underestimation that was previously overlooked in prior research. However, when predicting disease risk, the absence of these families (as seen in scenario iii) can yield the highest prediction accuracy and the strongest correlation with Polygenic Risk Scores. Our findings represent the first evidence of the important contribution of familial structure for heritability estimations and genomic prediction studies in T2D.

Abstract Sickle cell disease (SCD) is an inherited blood disorder marked by homozygosity of hemoglobin S, which is a defective hemoglobin caused by a missense mutation in the β-globin gene. However, clinical phenotypes of SCD vary among patients. To investigate genetic variants associated with various clinical phenotypes of SCD, we genotyped DNA samples from 520 SCD subjects and used a genome-wide association study (GWAS) approach to identify genetic variants associated with phenotypic features of SCD. For HbF levels, the previously reported 2p16.1 locus (BCL11A) reached genome significance (rs1427407, P = 8.58 × 10–10) in our GWAS as expected. In addition, we found a new genome-wide significance locus at 15q14 (rs8182015, P = 2.07 × 10–8) near gene EMC7. GWAS of acute chest syndrome (ACS) detected a locus (rs79915189, P = 3.70 × 10–8) near gene IDH2 at 15q26.1. The SNP, rs79915189, is also an expression quantitative trait locus (eQTL) of IDH2 in multiple tissues. For vasoocclusive episode (VOE), GWAS detected multiple significant signals at 2p25.1 (rs62118798, P = 4.27 × 10–8), 15q26.1 (rs62020555, P = 2.04 × 10–9) and 15q26.3 (rs117797325, P = 4.63 × 10–8). Our findings provide novel insights into the genetic mechanisms of SCD suggesting that common genetic variants play an important role in the presentation of the clinical phenotypes of patients with SCD.

Abstract Jaundice affects almost all neonates in their first days of life and is caused by the accumulation of bilirubin. Although the core biochemistry of bilirubin metabolism is well understood, it is not clear why some neonates experience more severe jaundice and require treatment with phototherapy. Here, we present the first genome-wide association study of neonatal jaundice to date in nearly 30,000 parent-offspring trios from Norway (cases ≈ 2000). The alternate allele of a common missense variant affecting the sequence of UGT1A4 reduces the susceptibility to jaundice five-fold, which replicated in separate cohorts of neonates of African American and European ancestries. eQTL colocalization analyses indicate that the association may be driven by regulation of UGT1A1 in the intestines, but not in the liver. Our results reveal marked differences in the genetic variants involved in neonatal jaundice compared to those regulating bilirubin levels in adults, suggesting distinct genetic mechanisms for the same biological pathways.

Abstract Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is uncharacterized. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio = 0.55 per standard deviation increase in PGSWBC [95%CI, 0.30−0.94], p = 0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n = 1724, hazard ratio [HR] = 0.78 [0.69−0.88], p = 4.0 × 10−5) or immunosuppressant (n = 354, HR = 0.61 [0.38–0.99], p = 0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n = 1,466, HR = 0.62 [0.44−0.87], p = 0.006). Collectively, these findings suggest that there are genetically predisposed individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

ObjectiveThe study delved into the epigenetic factors associated with periodontal disease in two lineages of mice, namely C57bl/6 and Balb/c. Its primary objective was to elucidate alterations in the methylome of mice with distinct genetic backgrounds following systemic microbial challenge, employing high-throughput DNA methylation analysis as the investigative tool.MethodsPorphyromonas gingivalis (Pg)was orally administered to induce periodontitis in both Balb/c and C57bl/6 lineage. After euthanasia, genomic DNA from both maxilla and blood were subjected to bisulfite conversion, PCR amplification and genome-wide DNA methylation analysis using the Ovation RRBS Methyl-Seq System coupled with the Illumina Infinium Mouse Methylation BeadChip.ResultsOf particular significance was the distinct methylation profile observed within the Pg-induced group of the Balb/c lineage, contrasting with both the control and Pg-induced groups of the C57bl/6 lineage. Utilizing rigorous filtering criteria, we successfully identified a substantial number of differentially methylated regions (DMRs) across various tissues and comparison groups, shedding light on the prevailing hypermethylation in non-induced cohorts and hypomethylation in induced groups. The comparison between blood and maxilla samples underscored the unique methylation patterns specific to the jaw tissue. Our comprehensive methylome analysis further unveiled statistically significant disparities, particularly within promoter regions, in several comparison groups.ConclusionThe differential DNA methylation patterns observed between C57bl/6 and Balb/c mouse lines suggest that epigenetic factors contribute to the variations in disease susceptibility. The identified differentially methylated regions associated with immune regulation and inflammatory response provide potential targets for further investigation. These findings emphasize the importance of considering epigenetic mechanisms in the development and progression of periodontitis.

Abstract Dyslipidemia, as a metabolic risk factor, with the strongest and most heritable independent cause of cardiovascular diseases worldwide. We investigated the familial transmission patterns of dyslipidemia through a longitudinal family-based cohort, the Tehran Cardiometabolic Genetic Study (TCGS) in Iran. We enrolled 18,729 individuals (45% were males) aged > 18 years (mean: 38.15 (15.82)) and observed them over five 3-year follow-up periods. We evaluated the serum concentrations of total cholesterol, triglyceride, high-density lipoprotein cholesterol, and low-density lipoprotein cholesterol with the first measurement among longitudinal measures and the average measurements (AM) of the five periods. Heritability analysis was conducted using a mixed-effect framework with likelihood-based and Bayesian approaches. The periodic prevalence and heritability of dyslipidemia were estimated to be 65.7 and 42%, respectively. The likelihood of an individual having at least one dyslipidemic parent reveals an OR = 6.94 (CI 5.28–9.30) compared to those who do not have dyslipidemic parents. The most considerable intraclass correlation of family members was for the same-sex siblings, with ICC ~ 25.5%. For serum concentrations, heritability ranged from 33.64 to 60.95%. Taken together, these findings demonstrate that familial transmission of dyslipidemia in the Tehran population is strong, especially within the same-gender siblings. According to previous reports, the heritability of dyslipidemia in this population is considerably higher than the global average.

Abstract 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.

Abstract Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10–14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.

Abstract Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Postural Orthostatic Tachycardia Syndrome (POTS) affects up to 1% of the US population, predominantly women, and is characterized by a complex, elusive etiology and heterogeneous phenotypes. This review delves into the intricate physiology and etiology of POTS, decoding the roles of the sinoatrial node, the autonomic nervous system, fluid dynamics, and the interplay between the immune and endocrine systems. It further examines key contributing factors such as dysautonomia, thoracic hypovolemia, autonomic neuropathies, sympathetic denervation, autoimmune responses, and associations with conditions such as small-fiber neuropathy and mast cell activation syndrome. Given the numerous mysteries surrounding POTS, we also cautiously bring attention to sinoatrial node and myocardial function, particularly in how the heart responds to stress despite exhibiting a normal cardiac phenotype at rest. The potential of genomic research in elucidating the underlying mechanisms of POTS is emphasized, suggesting this as a valuable approach that is likely to improve our understanding of the genetic underpinnings of POTS. The review introduces a tentative classification system for the etiological factors in POTS, which seeks to capture the condition’s diverse aspects by categorizing various etiological factors and acknowledging co-occurring conditions. This classification, while aiming to enhance understanding and optimize treatment targets, is presented as a preliminary model needing further study and refinement. This review underscores the ongoing need for research to unravel the complexities of POTS and to develop targeted therapies that can improve patient outcomes.

Abstract Children with birth defects (BD) express distinct clinical features that often have various medical consequences, one of which is predisposition to the development of cancers. Identification of the underlying genetic mechanisms related to the development of cancer in BD patients would allow for preventive measures. We performed a whole genome sequencing (WGS) study on blood-derived DNA samples from 1566 individuals without chromosomal anomalies, including 454 BD probands with at least one type of malignant tumors, 767 cancer-free BD probands, and 345 healthy individuals. Exclusive recurrent variants were identified in BD-cancer and BD-only patients and mapped to their corresponding genomic regions. We observed statistically significant overlaps for protein-coding/ncRNA with exclusive variants in exons, introns, ncRNAs, and 3’UTR regions. Exclusive exonic variants, especially synonymous variants, tend to occur in prior exons locus in BD-cancer children. Intronic variants close to splicing site (

Abstract Genetic variants in chromatin regulators are frequently found in neurodevelopmental disorders, but their effect in disease etiology is rarely determined. Here, we uncover and functionally define pathogenic variants in the chromatin modifier EZH1 as the cause of dominant and recessive neurodevelopmental disorders in 19 individuals. EZH1 encodes one of the two alternative histone H3 lysine 27 methyltransferases of the PRC2 complex. Unlike the other PRC2 subunits, which are involved in cancers and developmental syndromes, the implication of EZH1 in human development and disease is largely unknown. Using cellular and biochemical studies, we demonstrate that recessive variants impair EZH1 expression causing loss of function effects, while dominant variants are missense mutations that affect evolutionarily conserved aminoacids, likely impacting EZH1 structure or function. Accordingly, we found increased methyltransferase activity leading to gain of function of two EZH1 missense variants. Furthermore, we show that EZH1 is necessary and sufficient for differentiation of neural progenitor cells in the developing chick embryo neural tube. Finally, using human pluripotent stem cell-derived neural cultures and forebrain organoids, we demonstrate that EZH1 variants perturb cortical neuron differentiation. Overall, our work reveals a critical role of EZH1 in neurogenesis regulation and provides molecular diagnosis for previously undefined neurodevelopmental disorders.

Application of polyethylene glycol (PEG) to a peripheral nerve injury at the time of primary neurorrhaphy is thought to prevent Wallerian degeneration via direct axolemma fusion. The molecular mechanisms of nerve fusion and recovery are unclear. Our study tested the hypothesis that PEG alters gene expression in neural and muscular environments as part of its restorative properties. Lewis rats underwent unilateral sciatic nerve transection with immediate primary repair. Subjects were randomly assigned to receive either PEG treatment or standard repair at the time of neurorrhaphy. Samples of sciatic nerve distal to the injury and tibialis muscle at the site of innervation were harvested at 24 hours and 4 weeks postoperatively. Total RNA sequencing and subsequent bioinformatics analyses were used to identify significant differences in differentially expressed genes (DEGs) and their related biological pathways (p

Abstract Background Neurodevelopmental disorders (NDDs), such as attention deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD), are examples of complex and partially overlapping phenotypes that often lack definitive corroborating genetic information. ADHD and ASD have complex genetic associations implicated by rare recurrent copy number variations (CNVs). Both of these NDDs have been shown to share similar biological etiologies as well as genetic pleiotropy. Methods Platforms aimed at investigating genetic-based associations, such as high-density microarray technologies, have been groundbreaking techniques in the field of complex diseases, aimed at elucidating the underlying disease biology. Previous studies have uncovered CNVs associated with genes within shared candidate genomic networks, including glutamate receptor genes, across multiple different NDDs. To examine shared biological pathways across two of the most common NDDs, we investigated CNVs across 15,689 individuals with ADHD (n = 7920), ASD (n = 4318), or both (n = 3,416), as well as 19,993 controls. Cases and controls were matched by genotype array (i.e., Illumina array versions). Three case–control association studies each calculated and compared the observed vs. expected frequency of CNVs across individual genes, loci, pathways, and gene networks. Quality control measures of confidence in CNV-calling, prior to association analyses, included visual inspection of genotype and hybridization intensity. Results Here, we report results from CNV analysis in search for individual genes, loci, pathways, and gene networks. To extend our previous observations implicating a key role of the metabotropic glutamate receptor (mGluR) network in both ADHD and autism, we exhaustively queried patients with ASD and/or ADHD for CNVs associated with the 273 genomic regions of interest within the mGluR gene network (genes with one or two degrees protein–protein interaction with mGluR 1–8 genes). Among CNVs in mGluR network genes, we uncovered CNTN4 deletions enriched in NDD cases (P = 3.22E − 26, OR = 2.49). Additionally, we uncovered PRLHR deletions in 40 ADHD cases and 12 controls (P = 5.26E − 13, OR = 8.45) as well as clinically diagnostic relevant 22q11.2 duplications and 16p11.2 duplications in 23 ADHD + ASD cases and 9 controls (P = 4.08E − 13, OR = 15.05) and 22q11.2 duplications in 34 ADHD + ASD cases and 51 controls (P = 9.21E − 9, OR = 3.93); those control samples were not with previous 22qDS diagnosis in their EHR records. Conclusion Together, these results suggest that disruption in neuronal cell-adhesion pathways confers significant risk to NDDs and showcase that rare recurrent CNVs in CNTN4, 22q11.2, and 16p11.2 are overrepresented in NDDs that constitute patients predominantly suffering from ADHD and ASD. Trial registration ClinicalTrials.gov Identifier: NCT02286817 First Posted: 10 November 14, ClinicalTrials.gov Identifier: NCT02777931 first posted: 19 May 2016, ClinicalTrials.gov Identifier: NCT03006367 first posted: 30 December 2016, ClinicalTrials.gov Identifier: NCT02895906 first posted: 12 September 2016.

Autoimmune diseases arise from atypical immune responses that attack self-tissue epitopes, and their development is intricately connected to the disruption of the JAK-STAT signaling pathway, where SOCS proteins play crucial roles. Conditions such as autoimmune uveitis, psoriasis, lupus, and autoimmune encephalitis exhibit immune system dysfunctions associated with JAK-STAT signaling dysregulation. Emerging therapeutic strategies utilize JAK-STAT inhibitors and SOCS mimetics to modulate immune responses and alleviate autoimmune manifestations. Although more research and clinical studies are required to assess their effectiveness, safety profiles, and potential for personalized therapeutic approaches in autoimmune conditions, JAK-STAT inhibitors and SOCS mimetics show promise as potential treatment options. This review explores the action, effectiveness, safety profiles, and future prospects of JAK inhibitors and SOCS mimetics as therapeutic agents for psoriasis, autoimmune uveitis, systemic lupus erythematosus, and autoimmune encephalitis. The findings underscore the importance of investigating these targeted therapies to advance treatment options for individuals suffering from autoimmune diseases.

Abstract Thoracic insufficiency syndromes are a genetically and phenotypically heterogeneous group of disorders characterized by congenital abnormalities or progressive deformation of the chest wall and/or vertebrae that result in restrictive lung disease and compromised respiratory capacity. We performed whole exome sequencing on a cohort of 42 children with thoracic insufficiency to elucidate the underlying molecular etiologies of syndromic and non-syndromic thoracic insufficiency and predict extra-skeletal manifestations and disease progression. Molecular diagnosis was established in 24/42 probands (57%), with 18/24 (75%) probands having definitive diagnoses as defined by laboratory and clinical criteria and 6/24 (25%) probands having strong candidate genes. Gene identified in cohort patients most commonly encoded components of the primary cilium, connective tissue, and extracellular matrix. A novel association between KIF7 and USP9X variants and thoracic insufficiency was identified. We report and expand the genetic and phenotypic spectrum of a cohort of children with thoracic insufficiency, reinforce the prevalence of extra-skeletal manifestations in thoracic insufficiency syndromes, and expand the phenotype of KIF7 and USP9X-related disease to include thoracic insufficiency.

Abstract Background Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. Results To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3–5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. Conclusions Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Single-cell multimodal sequencing technologies are developed to simultaneously profile different modalities of data in the same cell. Here the authors develops a multimodal deep clustering method for the analysis of single-cell multi-omics data that supports clustering different types of multi-omics data and multi-batch data, as well as downstream differential expression analysis.

Heart failure is a major cause of cardiovascular morbidity and mortality. Here, the authors report results of a genome-wide association study meta-analysis, characterizing the role of common genetic variants in heart failure, finding overlap with common cardiovascular risk factors and imaging measures of cardiac structure/function.

Abstract Background Children with birth defects (BD) are more likely to develop cancer and the increased risk of cancer persists into adulthood. Prior population-based assessments have demonstrated that even non-chromosomal BDs are associated with at least two-fold increase of cancer risk. Identification of variants that are associated with malignant tumor in BD patients without chromosomal anomalies may improve our understanding of the underlying molecular mechanisms and provide clues for early cancer detection in children with BD. Methods In this study, whole genome sequencing (WGS) data of blood-derived DNA for 1653 individuals without chromosomal anomalies were acquired from the Kids First Data Resource Center (DRC), including 541 BD probands with at least one type of malignant tumors, 767 BD probands without malignant tumor, and 345 healthy family members who are the parents or siblings of the probands. Recurrent variants exclusively seen in cancer patients were selected and mapped to their corresponding genomic regions. The targeted genes/non-coding RNAs were further reduced using random forest and forward feature selection (ffs) models. Results The filtered genes/non-coding RNAs, including variants in non-coding areas, showed enrichment in cancer-related pathways. To further support the validity of these variants, blood WGS data of additional 40 independent BD probands, including 25 patients with at least one type of cancers from unrelated projects, were acquired. The counts of variants of interest identified in the Kid First data showed clear deviation in the validation dataset between BD patients with cancer and without cancer. Furthermore, a deep learning model was built to assess the predictive abilities in the 40 patients using variants of interest identified in the Kids First cohort as feature vectors. The accuracies are ~ 75%, with the noteworthy observation that variants mapped to non-coding regions provided the highest accuracy (31 out of 40 patients were labeled correctly). Conclusion We present for the first time a panorama of genetic variants that are associated with cancers in non-chromosomal BD patients, implying that our approach may potentially serve for the early detection of malignant tumors in patients with BD.

Abstract Background Obesity is a complex trait caused by a combination of genetic, environmental and lifestyle factors that contributes to the risks of numerous serious diseases. Predictive measures of body mass index (BMI) hold significant promise, with implications for the prevention and early intervention of obesity, promoting overall improvement in health. Main body An effective BMI polygenic risk score (PRS) model can assist with the prediction and early detection of obesity at the individual level, aligning with the objectives of precision medicine in the management of obesity. However, potential health disparities may emerge among under‐represented populations with a high prevalence of obesity, primarily due to the lack of genomic data available for these populations. The development of a trans‐ethnic (TE) PRS for BMI necessitates collective action from the research community, and requires genomic data from diverse populations. Conclusion The current BMI‐PRS model exhibits moderate performance, which could be improved in several key areas: integrating genomic information through TE GWAS studies, including admixed populations, considering gene–environment interactions, implementing advanced machine learning techniques, and incorporating genomic‐informed risk assessment (GIRA) with the inclusion of family history, environmental and lifestyle factors for the risk prediction.

Abstract Background While polygenic risk scores hold significant promise in estimating an individual's risk of developing a complex trait such as obesity, their application in the clinic has, to date, been limited by a lack of data from non‐European populations. As a collaboration model of the International Hundred K+ Cohorts Consortium (IHCC), we endeavored to develop a globally applicable trans‐ethnic PRS for body mass index (BMI) through this relatively new international effort. Methods The polygenic risk score (PRS) model was developed, trained and tested at the Center for Applied Genomics (CAG) of The Children's Hospital of Philadelphia (CHOP) based on a BMI meta‐analysis from the GIANT consortium. The validated PRS models were subsequently disseminated to the participating sites. Scores were generated by each site locally on their cohorts and summary statistics returned to CAG for final analysis. Results We show that in the absence of a well powered trans‐ethnic GWAS from which to derive marker SNPs and effect estimates for PRS, trans‐ethnic scores can be generated from European ancestry GWAS using Bayesian approaches such as LDpred, by adjusting the summary statistics using trans‐ethnic linkage disequilibrium reference panels. The ported trans‐ethnic scores outperform population specific‐PRS across all non‐European ancestry populations investigated including East Asians and three‐way admixed Brazilian cohort. Conclusions Here we show that for a truly polygenic trait such as BMI adjusting the summary statistics of a well powered European ancestry study using trans‐ethnic LD reference results in a score that is predictive across a range of ancestries including East Asians and three‐way admixed Brazilians.

Abstract Complex chromosomal rearrangements involve the restructuring of genetic material within a single chromosome or across multiple chromosomes. These events can cause serious human disease by disrupting coding DNA and gene regulatory elements via deletions, duplications, and structural rearrangements. Here we describe a 5-year-old female with severe developmental delay, dysmorphic features, multi-suture craniosynostosis, and growth failure found to have a complex series of balanced intra- and inter-chromosomal rearrangements involving chromosomes 4, 11, 13, and X. Initial clinical studies were performed by karyotype, chromosomal microarray, and FISH with research-based short-read genome sequencing coupled with sanger sequencing to precisely map her breakpoints to the base pair resolution to understand the molecular basis of her phenotype. Genome analysis revealed two pathogenic deletions at 4p16.1-p15.32 and 4q31.1, accounting for her developmental delay and dysmorphism. We identified over 60 breakpoints, many with blunt ends and limited homology, supporting a role for non-homologous end joining in restructuring and resolution of the seminal chromoplexy event. We propose that the complexity of our patient’s genomic rearrangements with a high number of breakpoints causes dysregulation of gene expression by three-dimensional chromatin interactions or topologically associating domains leading to growth failure and craniosynostosis. Our work supports an important role for genome sequencing in understanding the molecular basis of complex chromosomal rearrangements in human disease.

Abstract Variability in response to short-acting β2-agonists (e.g., albuterol) among patients with asthma from diverse racial/ethnic groups may contribute to asthma disparities. We sought to identify genetic variants associated with bronchodilator response (BDR) to identify potential mechanisms of drug response and risk factors for worse asthma outcomes. Genome-wide association studies of bronchodilator response (BDR) were performed using TOPMed Whole Genome Sequencing data of the Asthma Translational Genomic Collaboration (ATGC), which corresponded to 1136 Puerto Rican, 656 Mexican and 4337 African American patients with asthma. With the population-specific GWAS results, a trans-ethnic meta-analysis was performed to identify BDR-associated variants shared across the three populations. Replication analysis was carried out in three pediatric asthma cohorts, including CAMP (Childhood Asthma Management Program; n = 560), GACRS (Genetics of Asthma in Costa Rica Study; n = 967) and HPR (Hartford-Puerto Rico; n = 417). A genome-wide significant locus (rs35661809; P = 3.61 × 10–8) in LINC02220, a non-coding RNA gene, was identified in Puerto Ricans. While this region was devoid of protein-coding genes, capture Hi-C data showed a distal interaction with the promoter of the DNAH5 gene in lung tissue. In replication analysis, the GACRS cohort yielded a nominal association (1-tailed P

Abstract Background In over half of pediatric cases, ADHD presents with comorbidities, and often, it is unclear whether the symptoms causing impairment are due to the comorbidity or the underlying ADHD. Comorbid conditions increase the likelihood for a more severe and persistent course and complicate treatment decisions. Therefore, it is highly important to establish an algorithm that identifies ADHD and comorbidities in order to improve research on ADHD using biorepository and other electronic record data. Methods It is feasible to accurately distinguish between ADHD in isolation from ADHD with comorbidities using an electronic algorithm designed to include other psychiatric disorders. We sought to develop an EHR phenotype algorithm to discriminate cases with ADHD in isolation from cases with ADHD with comorbidities more effectively for efficient future searches in large biorepositories. We developed a multi-source algorithm allowing for a more complete view of the patient’s EHR, leveraging the biobank of the Center for Applied Genomics (CAG) at Children’s Hospital of Philadelphia (CHOP). We mined EHRs from 2009 to 2016 using International Statistical Classification of Diseases and Related Health Problems (ICD) codes, medication history and keywords specific to ADHD, and comorbid psychiatric disorders to facilitate genotype-phenotype correlation efforts. Chart abstractions and behavioral surveys added evidence in support of the psychiatric diagnoses. Most notably, the algorithm did not exclude other psychiatric disorders, as is the case in many previous algorithms. Controls lacked psychiatric and other neurological disorders. Participants enrolled in various CAG studies at CHOP and completed a broad informed consent, including consent for prospective analyses of EHRs. We created and validated an EHR-based algorithm to classify ADHD and comorbid psychiatric status in a pediatric healthcare network to be used in future genetic analyses and discovery-based studies. Results In this retrospective case-control study that included data from 51,293 subjects, 5840 ADHD cases were discovered of which 46.1% had ADHD alone and 53.9% had ADHD with psychiatric comorbidities. Our primary study outcome was to examine whether the algorithm could identify and distinguish ADHD exclusive cases from ADHD comorbid cases. The results indicate ICD codes coupled with medication searches revealed the most cases. We discovered ADHD-related keywords did not increase yield. However, we found including ADHD-specific medications increased our number of cases by 21%. Positive predictive values (PPVs) were 95% for ADHD cases and 93% for controls. Conclusion We established a new algorithm and demonstrated the feasibility of the electronic algorithm approach to accurately diagnose ADHD and comorbid conditions, verifying the efficiency of our large biorepository for further genetic discovery-based analyses. Trial registration ClinicalTrials.gov, NCT02286817 . First posted on 10 November 2014. ClinicalTrials.gov, NCT02777931 . First posted on 19 May 2016. ClinicalTrials.gov, NCT03006367 . First posted on 30 December 2016. ClinicalTrials.gov, NCT02895906 . First posted on 12 September 2016.

Central conducting lymphatic anomaly (CCLA) due to congenital maldevelopment of the lymphatics can result in debilitating and life-threatening disease with limited treatment options. We identified 4 individuals with CCLA, lymphedema, and microcystic lymphatic malformation due to pathogenic, mosaic variants in KRAS. To determine the functional impact of these variants and identify a targeted therapy for these individuals, we used primary human dermal lymphatic endothelial cells (HDLECs) and zebrafish larvae to model the lymphatic dysplasia. Expression of the p.Gly12Asp and p.Gly13Asp variants in HDLECs in a 2‑dimensional (2D) model and 3D organoid model led to increased ERK phosphorylation, demonstrating these variants activate the RAS/MAPK pathway. Expression of activating KRAS variants in the venous and lymphatic endothelium in zebrafish resulted in lymphatic dysplasia and edema similar to the individuals in the study. Treatment with MEK inhibition significantly reduced the phenotypes in both the organoid and the zebrafish model systems. In conclusion, we present the molecular characterization of the observed lymphatic anomalies due to pathogenic, somatic, activating KRAS variants in humans. Our preclinical studies suggest that MEK inhibition should be studied in future clinical trials for CCLA due to activating KRAS pathogenic variants.

Hallux valgus is a common form of foot deformity, and genetic factors contribute substantially to the pathogenesis of hallux valgus deformity. We conducted a genetic study on the structural variants underlying familial hallux valgus using whole exome sequencing approach. Twenty individuals from five hallux valgus families and two sporadic cases were included in this study. A total of 372 copy number variations were found and passed quality control filtering. Among them, 43 were only present in cases but not in controls or healthy individuals in the database of genomic variants. The genes covered by these copy number variations were enriched in gene sets related to immune signaling pathway, and cytochrome P450 metabolism. The hereditary CNVs demonstrate a dominant inheritance pattern. Two candidate pathogenic CNVs were further validated by quantitative-PCR. This study suggests that hallux valgus is a degenerative joint disease involving the dysregulation of immune and metabolism signaling pathways.

Background Cardiometabolic diseases are highly comorbid, but their relationship with female‐specific or overwhelmingly female‐predominant health conditions (breast cancer, endometriosis, pregnancy complications) is understudied. This study aimed to estimate the cross‐trait genetic overlap and influence of genetic burden of cardiometabolic traits on health conditions unique to women. Methods and Results Using electronic health record data from 71 008 ancestrally diverse women, we examined relationships between 23 obstetrical/gynecological conditions and 4 cardiometabolic phenotypes (body mass index, coronary artery disease, type 2 diabetes, and hypertension) by performing 4 analyses: (1) cross‐trait genetic correlation analyses to compare genetic architecture, (2) polygenic risk score–based association tests to characterize shared genetic effects on disease risk, (3) Mendelian randomization for significant associations to assess cross‐trait causal relationships, and (4) chronology analyses to visualize the timeline of events unique to groups of women with high and low genetic burden for cardiometabolic traits and highlight the disease prevalence in risk groups by age. We observed 27 significant associations between cardiometabolic polygenic scores and obstetrical/gynecological conditions (body mass index and endometrial cancer, body mass index and polycystic ovarian syndrome, type 2 diabetes and gestational diabetes, type 2 diabetes and polycystic ovarian syndrome). Mendelian randomization analysis provided additional evidence of independent causal effects. We also identified an inverse association between coronary artery disease and breast cancer. High cardiometabolic polygenic scores were associated with early development of polycystic ovarian syndrome and gestational hypertension. Conclusions We conclude that polygenic susceptibility to cardiometabolic traits is associated with elevated risk of certain female‐specific health conditions.

Abstract Background Asthma is a complex condition largely attributed to the interactions among genes and environments as a heterogeneous phenotype. Obesity is significantly associated with asthma development, and genetic studies on obese vs. non-obese asthma are warranted. Methods To investigate asthma in the minority African American (AA) population with or without obesity, we performed a whole genome sequencing (WGS) study on blood-derived DNA of 4289 AA individuals, included 2226 asthma patients (1364 with obesity and 862 without obesity) and 2006 controls without asthma. The burden analysis of functional rare coding variants was performed by comparing asthma vs. controls and by stratified analysis of obese vs. non-obese asthma, respectively. Results Among the top 66 genes with P

Abstract Variants of uncertain significance (VUS) are commonly found following genomic sequencing, particularly in ethnically diverse populations that are underrepresented in large population databases. Functional characterization of VUS may assist in variant reclassification, however these studies are not readily available and often rely on research funding and good will. We present four individuals from three families at different stages of their diagnostic trajectory with recurrent acute liver failure (RALF) and biallelic NBAS variants, confirmed by either trio analysis or cDNA studies. Functional characterization was undertaken, measuring NBAS and p31 levels by Western blotting, demonstrating reduced NBAS levels in two of three families, and reduced p31 levels in all three families. These results provided functional characterization of the molecular impact of a missense VUS, allowing reclassification of the variant and molecular confirmation of NBAS‐associated RALF. Importantly, p31 was decreased in all individuals, including an individual with two missense variants where NBAS protein levels were preserved. These results highlight the importance of access to timely functional studies after identification of putative variants, and the importance of considering a range of assays to validate variants whose pathogenicity is uncertain. We suggest that funding models for genomic sequencing should consider incorporating capabilities for adjunct RNA, protein, biochemical, and other specialized tests to increase the diagnostic yield which will lead to improved medical care, increased equity, and access to molecular diagnoses for all patients.

Abstract Introduction Currently, one of the commonly used methods for disseminating electronic health record (EHR)-based phenotype algorithms is providing a narrative description of the algorithm logic, often accompanied by flowcharts. A challenge with this mode of dissemination is the potential for under-specification in the algorithm definition, which leads to ambiguity and vagueness. Methods This study examines incidents of under-specification that occurred during the implementation of 34 narrative phenotyping algorithms in the electronic Medical Record and Genomics (eMERGE) network. We reviewed the online communication history between algorithm developers and implementers within the Phenotype Knowledge Base (PheKB) platform, where questions could be raised and answered regarding the intended implementation of a phenotype algorithm. Results We developed a taxonomy of under-specification categories via an iterative review process between two groups of annotators. Under-specifications that lead to ambiguity and vagueness were consistently found across narrative phenotype algorithms developed by all involved eMERGE sites. Discussion and conclusion Our findings highlight that under-specification is an impediment to the accuracy and efficiency of the implementation of current narrative phenotyping algorithms, and we propose approaches for mitigating these issues and improved methods for disseminating EHR phenotyping algorithms.

IntroductionType 1 diabetes, a disorder caused by autoimmune destruction of pancreatic insulin-producing cells, is more difficult to manage when it presents at a younger age. We sought to identify genetic correlates of the age of onset by conducting the first genome-wide association study (GWAS) treating the age of first diagnosis as a quantitative trait.MethodsWe performed GWAS with a discovery cohort of 4,014 cases and a replication cohort of 493 independent cases. Genome-wide significant SNPs were mapped to a causal variant by Bayesian conditional analysis and gel shift assay. The causal protein-coding gene was identified and characterized by RNA interference treatment of primary human pan-CD4+ T cells with RNA-seq of the transcriptome. The candidate gene was evaluated functionally in primary cells by CD69 staining and proliferation assays.ResultsOur GWAS replicated the known association of the age of diagnosis with the human leukocyte antigen complex (HLA-DQB1). The second signal identified was in an intron of the NELL1 gene on chromosome 11 and fine-mapped to variant rs10833518 (P < 1.54 × 10−9). Homozygosity for the risk allele leads to average age of onset one year earlier. Knock-down of HIV TAT-interacting protein 2 (HTATIP2), but not other genes in the locus, resulted in alterations to gene expression in signal transduction pathways including MAP kinases and PI3-kinase. Higher levels of HTATIP2 expression are associated with increased viability, proliferation, and activation of T cells in the presence of signals from antigen and cytokine receptors.DiscussionThis study implicates HTATIP2 as a new type 1 diabetes gene acting via T cell regulation. Larger population sample sizes are expected to reveal additional loci.

Abstract Background SYNJ1 encodes Synaptojanin‐1, a dual‐function poly‐phosphoinositide phosphatase that is expressed in the brain to regulate neuronal synaptic vesicle dynamics. Biallelic SYNJ1 variants cause a spectrum of clinical manifestations, from early onset parkinsonism to developmental and epileptic encephalopathy. Methods Proband‐only exome sequencing was used to identify a homozygous SYNJ1 pathogenic variant in an individual with epileptic encephalopathy. Sanger sequencing was used to confirm the variant. Results We present an Afro‐Caribbean female who developed uncontrollable seizures shortly after birth, accompanied by developmental delay and severe generalized dystonia. She had homozygosity for a novel c.242‐2A > G variant in SYNJ1 with both parents being heterozygous carriers. An older sister was reported to have had a similar presentation but was not examined. Both siblings died at an approximate age of 16 years. Conclusions We report a novel pathogenic variant in SYNJ1 present in homozygosity leading to developmental and epileptic encephalopathy. Currently, there are only 4 reports describing 10 individuals with SYNJ1‐related developmental and epileptic encephalopathy. This case expands the clinical knowledge and the allelic heterogeneity associated with SYNJ1 variants.

Summary: WDR5 is a broadly studied, highly conserved key protein involved in a wide array of biological functions. Among these functions, WDR5 is a part of several protein complexes that affect gene regulation via post-translational modification of histones. We collected data from 11 unrelated individuals with six different rare de novo germline missense variants in WDR5; one identical variant was found in five individuals and another variant in two individuals. All individuals had neurodevelopmental disorders including speech/language delays (n = 11), intellectual disability (n = 9), epilepsy (n = 7), and autism spectrum disorder (n = 4). Additional phenotypic features included abnormal growth parameters (n = 7), heart anomalies (n = 2), and hearing loss (n = 2). Three-dimensional protein structures indicate that all the residues affected by these variants are located at the surface of one side of the WDR5 protein. It is predicted that five out of the six amino acid substitutions disrupt interactions of WDR5 with RbBP5 and/or KMT2A/C, as part of the COMPASS (complex proteins associated with Set1) family complexes. Our experimental approaches in Drosophila melanogaster and human cell lines show normal protein expression, localization, and protein-protein interactions for all tested variants. These results, together with the clustering of variants in a specific region of WDR5 and the absence of truncating variants so far, suggest that dominant-negative or gain-of-function mechanisms might be at play. All in all, we define a neurodevelopmental disorder associated with missense variants in WDR5 and a broad range of features. This finding highlights the important role of genes encoding COMPASS family proteins in neurodevelopmental disorders.

Nonischemic dilated cardiomyopathy is a common form of heart muscle disease in which genetic factors play a critical etiological role. In this regard, both rare disease‐causing mutations and common disease‐susceptible variants, in the Bcl‐2–associated athanogene 3 (BAG3) gene have been reported, highlighting the critical role of BAG3 in cardiomyocytes and in the development of dilated cardiomyopathy. The phenotypic effects of the BAG3 mutations help investigators understand the structure and function of the BAG3 gene. Indeed, we report herein that all of the known pathogenic/likely pathogenic variants affect at least 1 of 3 protein functional domains, ie, the WW domain, the second IPV (Ile‐Pro‐Val) domain, or the BAG domain, whereas none of the missense nontruncating pathogenic/likely pathogenic variants affect the proline‐rich repeat (PXXP) domain. A common variant, p.Cys151Arg, associated with reduced susceptibility to dilated cardiomyopathy demonstrated a significant difference in allele frequencies among diverse human populations, suggesting evolutionary selective pressure. As BAG3‐related therapies for heart failure move from the laboratory to the clinic, the ability to provide precision medicine will depend in large part on having a thorough understanding of the potential effects of both common and uncommon genetic variants on these target proteins. The current review article provides a roadmap that investigators can utilize to determine the potential interactions between a patient's genotype, their phenotype, and their response to therapeutic interventions with both gene delivery and small molecules.

Genetic studies of eczema to date have mostly explored common genetic variation. Here, the authors perform a large meta-analysis for common and rare variants and discover 8 loci associated with eczema. Over 20% of the heritability of the condition is attributable to rare variants.

Globally, asthma prevention and treatment remain a challenge. Ambient air pollution (AAP) is an environmental risk factor of special interest in asthma research. AAP is poorly defined and has been subdivided either by the origin of the air pollution or by the specific bioactive compounds. The link between AAP exposure and asthma exacerbations is well established and has been extensively reviewed. In this narrative review, we discuss the specific genetic variants that have been associated with increased AAP sensitivity and impact in paediatric asthma. We highlight the relative importance of variants associated with genes with a role in oxidant defences and the nuclear factor-κB pathway supporting a potential central role for these pathways in AAP sensitivity.

We report a 20-year-old, female, adopted Indian patient with over 662 Mb regions of homozy-gosity who presented with intellectual disability, ataxia, schizophrenia, retinal dystrophy, moder-ate-to-severe progressive sensorineural hearing loss (SNHL), congenital hypothyroidism, cleft mi-tral valve with mild mitral valve regurgitation, and dysmorphic features. Exome analysis first on a clinical basis and subsequently on research reanalysis uncovered pathogenic variants in three nu-clear genes following two modes of inheritance that were causal to her complex phenotype. These included (1) compound heterozygous variants in BBS6 potentially causative for Bardet–Biedl syn-drome 6; (2) a homozygous, known pathogenic variant in the stereocilin (STRC) gene associated with nonsyndromic deafness; and (3) a homozygous variant in dual oxidase 2 (DUOX2) gene asso-ciated with congenital hypothyroidism. A variant of uncertain significance was identified in a fourth gene, troponin T2 (TNNT2), associated with cardiomyopathy but not the cleft mitral valve, with mild mitral regurgitation seen in this case. This patient was the product of an apparent first-degree relationship, explaining the multiple independent inherited findings. This case high-lights the need to carefully evaluate multiple independent genetic etiologies for complex pheno-types, particularly in the case of consanguinity, rather than presuming unexplained features are expansions of known gene disorders.

Abstract With polygenic risk score (PRS) for autoimmune type 1 diabetes (T1D), this study identified T1D cases with low T1D PRS and searched for susceptibility loci in these cases. Our hypothesis is that genetic effects (likely mediated by relatively rare genetic variants) of non-mainstream (or non-autoimmune) T1D might have been diluted in the previous studies on T1D cases in general. Two cohorts for the PRS modeling and testing respectively were included. The first cohort consisted of 3302 T1D cases and 6181 controls, and the independent second cohort consisted of 3297 T1D cases and 6169 controls. Cases with low T1D PRS were identified using PRSice-2 and compared to controls with low T1D PRS by genome-wide association (GWA) test. Thirteen novel genetic loci with high imputation quality (Quality Score r2 > 0.91) were identified of SNPs/SNVs associated with low PRS T1D at genome-wide significance (P ≤ 5.0 × E−08), in addition to 4 established T1D loci, 3 reported loci by our previous study, as well as 9 potential novel loci represented by rare SNVs, but with relatively low imputation quality (Quality Score r2

Hui-Qi Qu et al. present a case-control association study to examine the genetic architecture of European patients with type 1 diabetes (T1D) of low genetic risk score. They identify 41 loci associated with T1D of low genetic risk score, including regions previously linked to viral infection and autism spectrum disorder, altogether providing further insight into the pathogenesis of T1D.