PreMode predicts mode-of-action of missense variants by deep graph representation learning of protein sequence and structural context.

Accurate prediction of the functional impact of missense variants is important for disease gene discovery, clinical genetic diagnostics, therapeutic strategies, and protein engineering. Previous efforts have focused on predicting a binary pathogenicity classification, but the functional impact of missense variants is multi-dimensional. Pathogenic missense variants in the same gene may act through different modes of action (i.e., gain/loss-of-function) by affecting different aspects of protein function. They may result in distinct clinical conditions that require different treatments. We developed a new method, PreMode, to perform gene-specific mode-of-action predictions. PreMode models effects of coding sequence variants using SE(3)-equivariant graph neural networks on protein sequences and structures. Using the largest-to-date set of missense variants with known modes of action, we showed that PreMode reached state-of-the-art performance in multiple types of mode-of-action predictions by efficient transfer-learning. Additionally, PreMode's prediction of G/LoF variants in a kinase is consistent with inactive-active conformation transition energy changes. Finally, we show that PreMode enables efficient study design of deep mutational scans and optimization in protein engineering.

Genome Sequencing is Critical for Forecasting Outcomes following Congenital Cardiac Surgery.

While genome sequencing has transformed medicine by elucidating the genetic underpinnings of both rare and common complex disorders, its utility to predict clinical outcomes remains understudied. Here, we used artificial intelligence (AI) technologies to explore the predictive value of genome sequencing in forecasting clinical outcomes following surgery for congenital heart defects (CHD). We report results for a cohort of 2,253 CHD patients from the Pediatric Cardiac Genomics Consortium with a broad range of complex heart defects, pre- and post-operative clinical variables and exome sequencing. Damaging genotypes in chromatin-modifying and cilia-related genes were associated with an elevated risk of adverse post-operative outcomes, including mortality, cardiac arrest and prolonged mechanical ventilation. The impact of damaging genotypes was further amplified in the context of specific CHD phenotypes, surgical complexity and extra-cardiac anomalies. The absence of a damaging genotype in chromatin-modifying and cilia-related genes was also informative, reducing the risk for adverse postoperative outcomes. Thus, genome sequencing enriches the ability to forecast outcomes following congenital cardiac surgery.

Homozygous missense variants in YKT6 result in loss of function and are associated with developmental delay, with or without severe infantile liver disease and risk for hepatocellular carcinoma.

PURPOSE: YKT6 plays important roles in multiple intracellular vesicle trafficking events but has not been associated with Mendelian diseases. METHODS: We report three unrelated individuals with rare homozygous missense variants in YKT6 who exhibited neurological disease with or without a progressive infantile liver disease. We modeled the variants in Drosophila. We generated wild-type and variant genomic rescue constructs (GRs) of the fly ortholog dYkt6 and compared their ability in rescuing the loss-of-function phenotypes in mutant flies. We also generated a dYkt6KozakGAL4 allele to assess the expression pattern of dYkt6. RESULTS: Two individuals are homozygous for YKT6 [NM_006555.3:c.554A>G p.(Tyr185Cys)] and exhibited normal prenatal course followed by failure to thrive, developmental delay and progressive liver disease. Haplotype analysis identified a shared homozygous region flanking the variant, suggesting a common ancestry. The third individual is homozygous for YKT6 [NM_006555.3:c.191A>G p.(Tyr64Cys)] and exhibited neurodevelopmental disorders and optic atrophy. Fly dYkt6 is essential and is expressed in the fat body (analogous to liver) and central nervous system. Wild-type GR can rescue the lethality and autophagic flux defects whereas the variants are less efficient in rescuing the phenotypes. CONCLUSION: The YKT6 variants are partial loss-of-function alleles and the p.(Tyr185Cys) is more severe than p.(Tyr64Cys).

Functional dissection of human cardiac enhancers and noncoding de novo variants in congenital heart disease.

Rare coding mutations cause ∼45% of congenital heart disease (CHD). Noncoding mutations that perturb cis-regulatory elements (CREs) likely contribute to the remaining cases, but their identification has been problematic. Using a lentiviral massively parallel reporter assay (lentiMPRA) in human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), we functionally evaluated 6,590 noncoding de novo variants (ncDNVs) prioritized from the whole-genome sequencing of 750 CHD trios. A total of 403 ncDNVs substantially affected cardiac CRE activity. A majority increased enhancer activity, often at regions with undetectable reference sequence activity. Of ten DNVs tested by introduction into their native genomic context, four altered the expression of neighboring genes and iPSC-CM transcriptional state. To prioritize future DNVs for functional testing, we used the MPRA data to develop a regression model, EpiCard. Analysis of an independent CHD cohort by EpiCard found enrichment of DNVs. Together, we developed a scalable system to measure the effect of ncDNVs on CRE activity and deployed it to systematically assess the contribution of ncDNVs to CHD.

Dynamic establishment of recipient resident memory T cell repertoire after human intestinal transplantation.

BACKGROUND: Understanding formation of the human tissue resident memory T cell (TRM) repertoire requires longitudinal access to human non-lymphoid tissues. METHODS: By applying flow cytometry and next generation sequencing to serial blood, lymphoid tissue, and gut samples from 16 intestinal transplantation (ITx) patients, we assessed the origin, distribution, and specificity of human TRMs at phenotypic and clonal levels. FINDINGS: Donor age ≥1 year and blood T cell macrochimerism (peak level ≥4%) were associated with delayed establishment of stable recipient TRM repertoires in the transplanted ileum. T cell receptor (TCR) overlap between paired gut and blood repertoires from ITx patients was significantly greater than that in healthy controls, demonstrating increased gut-blood crosstalk after ITx. Crosstalk with the circulating pool remained high for years of follow-up. TCR sequences identifiable in pre-Tx recipient gut but not those in lymphoid tissues alone were more likely to populate post-Tx ileal allografts. Clones detected in both pre-Tx gut and lymphoid tissue had distinct transcriptional profiles from those identifiable in only one tissue. Recipient T cells were distributed widely throughout the gut, including allograft and native colon, which had substantial repertoire overlap. Both alloreactive and microbe-reactive recipient T cells persisted in transplanted ileum, contributing to the TRM repertoire. INTERPRETATION: Our studies reveal human intestinal TRM repertoire establishment from the circulation, preferentially involving lymphoid tissue counterparts of recipient intestinal T cell clones, including TRMs. We have described the temporal and spatial dynamics of this active crosstalk between the circulating pool and the intestinal TRM pool. FUNDING: This study was funded by the National Institute of Allergy and Infectious Diseases (NIAID) P01 grant AI106697.

The osteoclastic activity in apical distal region of molar mesial roots affects orthodontic tooth movement and root resorption in rats.

The utilization of optimal orthodontic force is crucial to prevent undesirable side effects and ensure efficient tooth movement during orthodontic treatment. However, the sensitivity of existing detection techniques is not sufficient, and the criteria for evaluating optimal force have not been yet established. Here, by employing 3D finite element analysis methodology, we found that the apical distal region (A-D region) of mesial roots is particularly sensitive to orthodontic force in rats. Tartrate-resistant acidic phosphatase (TRAP)-positive osteoclasts began accumulating in the A-D region under the force of 40 grams (g), leading to alveolar bone resorption and tooth movement. When the force reached 80 g, TRAP-positive osteoclasts started appearing on the root surface in the A-D region. Additionally, micro-computed tomography revealed a significant root resorption at 80 g. Notably, the A-D region was identified as a major contributor to whole root resorption. It was determined that 40 g is the minimum effective force for tooth movement with minimal side effects according to the analysis of tooth movement, inclination, and hyalinization. These findings suggest that the A-D region with its changes on the root surface is an important consideration and sensitive indicator when evaluating orthodontic forces for a rat model. Collectively, our investigations into this region would aid in offering valuable implications for preventing and minimizing root resorption during patients' orthodontic treatment.

Rescuing lung development through embryonic inhibition of histone acetylation.

A major barrier to the impact of genomic diagnosis in patients with congenital malformations is the lack of understanding regarding how sequence variants contribute to disease pathogenesis and whether this information could be used to generate patient-specific therapies. Congenital diaphragmatic hernia (CDH) is among the most common and severe of all structural malformations; however, its underlying mechanisms are unclear. We identified loss-of-function sequence variants in the epigenomic regulator gene SIN3A in two patients with complex CDH. Tissue-specific deletion of Sin3a in mice resulted in defects in diaphragm development, lung hypoplasia, and pulmonary hypertension, the cardinal features of CDH and major causes of CDH-associated mortality. Loss of SIN3A in the lung mesenchyme resulted in reduced cellular differentiation, impaired cell proliferation, and increased DNA damage. Treatment of embryonic Sin3a mutant mice with anacardic acid, an inhibitor of histone acetyltransferase, reduced DNA damage, increased cell proliferation and differentiation, improved lung and pulmonary vascular development, and reduced pulmonary hypertension. These findings demonstrate that restoring the balance of histone acetylation can improve lung development in the Sin3a mouse model of CDH.

Plasticity of intragraft alloreactive T cell clones in human gut correlates with transplant outcomes.

The site of transition between tissue-resident memory (TRM) and circulating phenotypes of T cells is unknown. We integrated clonotype, alloreactivity, and gene expression profiles of graft-repopulating recipient T cells in the intestinal mucosa at the single-cell level after human intestinal transplantation. Host-versus-graft (HvG)-reactive T cells were mainly distributed to TRM, effector T (Teff)/TRM, and T follicular helper compartments. RNA velocity analysis demonstrated a trajectory from TRM to Teff/TRM clusters in association with rejection. By integrating pre- and post-transplantation (Tx) mixed lymphocyte reaction-determined alloreactive repertoires, we observed that pre-existing HvG-reactive T cells that demonstrated tolerance in the circulation were dominated by TRM profiles in quiescent allografts. Putative de novo HvG-reactive clones showed a transcriptional profile skewed to cytotoxic effectors in rejecting grafts. Inferred protein regulon network analysis revealed upstream regulators that accounted for the effector and tolerant T cell states. We demonstrate Teff/TRM interchangeability for individual T cell clones with known (allo)recognition in the human gut, providing novel insight into TRM biology.

Thalidomide for Recurrent Bleeding Due to Small-Intestinal Angiodysplasia.

BACKGROUND: Recurrent bleeding from the small intestine accounts for 5 to 10% of cases of gastrointestinal bleeding and remains a therapeutic challenge. Thalidomide has been evaluated for the treatment of recurrent bleeding due to small-intestinal angiodysplasia (SIA), but confirmatory trials are lacking. METHODS: We conducted a multicenter, double-blind, randomized, placebo-controlled trial to investigate the efficacy and safety of thalidomide for the treatment of recurrent bleeding due to SIA. Eligible patients with recurrent bleeding (at least four episodes of bleeding during the previous year) due to SIA were randomly assigned to receive thalidomide at an oral daily dose of 100 mg or 50 mg or placebo for 4 months. Patients were followed for at least 1 year after the end of the 4-month treatment period. The primary end point was effective response, which was defined as a reduction of at least 50% in the number of bleeding episodes that occurred during the year after the end of thalidomide treatment as compared with the number that occurred during the year before treatment. Key secondary end points were cessation of bleeding without rebleeding, blood transfusion, hospitalization because of bleeding, duration of bleeding, and hemoglobin levels. RESULTS: Overall, 150 patients underwent randomization: 51 to the 100-mg thalidomide group, 49 to the 50-mg thalidomide group, and 50 to the placebo group. The percentages of patients with an effective response in the 100-mg thalidomide group, 50-mg thalidomide group, and placebo group were 68.6%, 51.0%, and 16.0%, respectively (P<0.001 for simultaneous comparison across the three groups). The results of the analyses of the secondary end points supported those of the primary end point. Adverse events were more common in the thalidomide groups than in the placebo group overall; specific events included constipation, somnolence, limb numbness, peripheral edema, dizziness, and elevated liver-enzyme levels. CONCLUSIONS: In this placebo-controlled trial, treatment with thalidomide resulted in a reduction in bleeding in patients with recurrent bleeding due to SIA. (Funded by the National Natural Science Foundation of China and the Shanghai Municipal Education Commission, Gaofeng Clinical Medicine; ClinicalTrials.gov number, NCT02707484.).

Disrupted endosomal trafficking of the Vangl-Celsr polarity complex underlies congenital anomalies in trachea-esophageal morphogenesis.

Disruptions in foregut morphogenesis can result in life-threatening conditions where the trachea and esophagus fail to separate properly, such as esophageal atresia (EA) and tracheoesophageal fistulas (TEF). The developmental basis of these congenital anomalies is poorly understood, but recent genome sequencing reveals that de novo variants in intracellular trafficking genes are enriched in EA/TEF patients. Here we show that mutation of orthologous genes in Xenopus disrupts trachea-esophageal separation similar to EA/TEF patients. We show that the Rab11a recycling endosome pathway is required to localize Vangl-Celsr polarity complexes at the cell surface where opposite sides of the common foregut tube fuse. Partial loss of endosome trafficking or the Vangl/Celsr complex disrupts epithelial polarity and cell division orientation. Mutant cells accumulate at the fusion point, fail to downregulate cadherin, and do not separate into distinct trachea and esophagus. These data provide new insights into the mechanisms of congenital anomalies and general paradigms of tissue fusion during organogenesis.

Rare predicted loss of function alleles in Bassoon (BSN) are associated with obesity.

Bassoon (BSN) is a component of a hetero-dimeric presynaptic cytomatrix protein that orchestrates neurotransmitter release with Piccolo (PCLO) from glutamatergic neurons throughout the brain. Heterozygous missense variants in BSN have previously been associated with neurodegenerative disorders in humans. We performed an exome-wide association analysis of ultra-rare variants in about 140,000 unrelated individuals from the UK Biobank to search for new genes associated with obesity. We found that rare heterozygous predicted loss of function (pLoF) variants in BSN are associated with higher BMI with p-value of 3.6e-12 in the UK biobank cohort. Additionally, we identified two individuals (one of whom has a de novo variant) with a heterozygous pLoF variant in a cohort of early onset or extreme obesity and report the clinical histories of these individuals with non-syndromic obesity with no history of neurobehavioral or cognitive disability. The BMI association was replicated in the All of Us whole genome sequencing data. Heterozygous pLoF BSN variants constitute a new etiology for obesity.

PLS3 missense variants affecting the actin-binding domains cause X-linked congenital diaphragmatic hernia and body-wall defects.

Congenital diaphragmatic hernia (CDH) is a relatively common and genetically heterogeneous structural birth defect associated with high mortality and morbidity. We describe eight unrelated families with an X-linked condition characterized by diaphragm defects, variable anterior body-wall anomalies, and/or facial dysmorphism. Using linkage analysis and exome or genome sequencing, we found that missense variants in plastin 3 (PLS3), a gene encoding an actin bundling protein, co-segregate with disease in all families. Loss-of-function variants in PLS3 have been previously associated with X-linked osteoporosis (MIM: 300910), so we used in silico protein modeling and a mouse model to address these seemingly disparate clinical phenotypes. The missense variants in individuals with CDH are located within the actin-binding domains of the protein but are not predicted to affect protein structure, whereas the variants in individuals with osteoporosis are predicted to result in loss of function. A mouse knockin model of a variant identified in one of the CDH-affected families, c.1497G>C (p.Trp499Cys), shows partial perinatal lethality and recapitulates the key findings of the human phenotype, including diaphragm and abdominal-wall defects. Both the mouse model and one adult human male with a CDH-associated PLS3 variant were observed to have increased rather than decreased bone mineral density. Together, these clinical and functional data in humans and mice reveal that specific missense variants affecting the actin-binding domains of PLS3 might have a gain-of-function effect and cause a Mendelian congenital disorder.

Significantly Enhancing Mechanical and Thermal Properties of Cellulose-Based Composites by Adding Small Amounts of Lysozyme-Modified Graphene Nanoplatelets via Forming Strong Double-Cross-Linked Interface Interactions.

Thermally conductive cellulose-based composites have great application potential in the thermal management of portable and wearable electronic devices. In this work, cellulose-based composites with excellent mechanical and thermal properties were developed by using lysozyme-modified graphene nanoplatelets (LmGNP), epichlorohydrin (ECH), and hydrolyzed cellulose via forming strong double-cross-linked interface interactions, including the hydrogen bond network generated between LmGNP and cellulose and the chemical cross-link of ECH. As for the composites containing 8 wt % LmGNP, the in-plane thermal conductivity was 3.341 W·m-1K-1, while the tensile stress was 114.60 MPa, which increased by 297.3 and 146.2%, respectively, compared to pure cellulose. Along with the good stability, insulation, and lightweight properties, the fabricated composites have the potential to become a promising heat dissipation material for wearable electronic devices.

VBASS enables integration of single cell gene expression data in Bayesian association analysis of rare variants.

Rare or de novo variants have substantial contribution to human diseases, but the statistical power to identify risk genes by rare variants is generally low due to rarity of genotype data. Previous studies have shown that risk genes usually have high expression in relevant cell types, although for many conditions the identity of these cell types are largely unknown. Recent efforts in single cell atlas in human and model organisms produced large amount of gene expression data. Here we present VBASS, a Bayesian method that integrates single-cell expression and de novo variant (DNV) data to improve power of disease risk gene discovery. VBASS models disease risk prior as a function of expression profiles, approximated by deep neural networks. It learns the weights of neural networks and parameters of Gamma-Poisson likelihood models of DNV counts jointly from expression and genetics data. On simulated data, VBASS shows proper error rate control and better power than state-of-the-art methods. We applied VBASS to published datasets and identified more candidate risk genes with supports from literature or data from independent cohorts. VBASS can be generalized to integrate other types of functional genomics data in statistical genetics analysis.

Site-specific development and progressive maturation of human tissue-resident memory T cells over infancy and childhood.

Infancy and childhood are critical life stages for generating immune memory to protect against pathogens; however, the timing, location, and pathways for memory development in humans remain elusive. Here, we investigated T cells in mucosal sites, lymphoid tissues, and blood from 96 pediatric donors aged 0-10 years using phenotypic, functional, and transcriptomic profiling. Our results revealed that memory T cells preferentially localized in the intestines and lungs during infancy and accumulated more rapidly in mucosal sites compared with blood and lymphoid organs, consistent with site-specific antigen exposure. Early life mucosal memory T cells exhibit distinct functional capacities and stem-like transcriptional profiles. In later childhood, they progressively adopt proinflammatory functions and tissue-resident signatures, coincident with increased T cell receptor (TCR) clonal expansion in mucosal and lymphoid sites. Together, our findings identify staged development of memory T cells targeted to tissues during the formative years, informing how we might promote and monitor immunity in children.

Contribution of Previously Unrecognized RNA Splice-Altering Variants to Congenital Heart Disease.

BACKGROUND: Known genetic causes of congenital heart disease (CHD) explain <40% of CHD cases, and interpreting the clinical significance of variants with uncertain functional impact remains challenging. We aim to improve diagnostic classification of variants in patients with CHD by assessing the impact of noncanonical splice region variants on RNA splicing. METHODS: We tested de novo variants from trio studies of 2649 CHD probands and their parents, as well as rare (allele frequency, <2×10-6) variants from 4472 CHD probands in the Pediatric Cardiac Genetics Consortium through a combined computational and in vitro approach. RESULTS: We identified 53 de novo and 74 rare variants in CHD cases that alter splicing and thus are loss of function. Of these, 77 variants are in known dominant, recessive, and candidate CHD genes, including KMT2D and RBFOX2. In 1 case, we confirmed the variant's predicted impact on RNA splicing in RNA transcripts from the proband's cardiac tissue. Two probands were found to have 2 loss-of-function variants for recessive CHD genes HECTD1 and DYNC2H1. In addition, SpliceAI-a predictive algorithm for altered RNA splicing-has a positive predictive value of ≈93% in our cohort. CONCLUSIONS: Through assessment of RNA splicing, we identified a new loss-of-function variant within a CHD gene in 78 probands, of whom 69 (1.5%; n=4472) did not have a previously established genetic explanation for CHD. Identification of splice-altering variants improves diagnostic classification and genetic diagnoses for CHD. REGISTRATION: URL: https://clinicaltrials.gov; Unique identifier: NCT01196182.

Targeting RORα in macrophages to boost diabetic bone regeneration.

Diabetes mellitus (DM) has become a serious threat to human health. Bone regeneration deficiency and nonunion caused by DM is perceived as a worldwide epidemic, with a very high socioeconomic impact on public health. Here, we find that targeted activation of retinoic acid-related orphan receptor α (RORα) by SR1078 in the early stage of bone defect repair can significantly promote in situ bone regeneration of DM rats. Bone regeneration relies on the activation of macrophage RORα in the early bone repair, but RORα of DM rats fails to upregulation as hyperglycemic inflammatory microenvironment induced IGF1-AMPK signalling deficiency. Mechanistic investigations suggest that RORα is vital for macrophage-induced migration and proliferation of bone mesenchymal stem cells (BMSCs) via a CCL3/IL-6 depending manner. In summary, our study identifies RORα expressed in macrophages during the early stage of bone defect repair is crucial for in situ bone regeneration, and offers a novel strategy for bone regeneration therapy and fracture repair in DM patients.

Developmental growth plate cartilage formation suppressed by artificial light at night via inhibiting BMAL1-driven collagen hydroxylation.

Exposure to artificial light at night (LAN) can induce obesity, depressive disorder and osteoporosis, but the pernicious effects of excessive LAN exposure on tissue structure are poorly understood. Here, we demonstrated that artificial LAN can impair developmental growth plate cartilage extracellular matrix (ECM) formation and cause endoplasmic reticulum (ER) dilation, which in turn compromises bone formation. Excessive LAN exposure induces downregulation of the core circadian clock protein BMAL1, which leads to collagen accumulation in the ER. Further investigations suggest that BMAL1 is the direct transcriptional activator of prolyl 4-hydroxylase subunit alpha 1 (P4ha1) in chondrocytes, which orchestrates collagen prolyl hydroxylation and secretion. BMAL1 downregulation induced by LAN markedly inhibits proline hydroxylation and transport of collagen from ER to golgi, thereby inducing ER stress in chondrocytes. Restoration of BMAL1/P4HA1 signaling can effectively rescue the dysregulation of cartilage formation within the developmental growth plate induced by artificial LAN exposure. In summary, our investigations suggested that LAN is a significant risk factor in bone growth and development, and a proposed novel strategy targeting enhancement of BMAL1-mediated collagen hydroxylation could be a potential therapeutic approach to facilitate bone growth.

Predicted loss of function alleles in Bassoon (BSN) are associated with obesity.

Bassoon ( BSN ) is a component of a hetero-dimeric presynaptic cytomatrix protein that orchestrates neurotransmitter release with Piccolo ( PCLO ) from glutamatergic neurons throughout the brain. Heterozygous missense variants in BSN have previously been associated with neurodegenerative disorders in humans. We performed an exome-wide association analysis of ultra-rare variants in about 140,000 unrelated individuals from the UK Biobank to search for new genes associated with obesity. We found that rare heterozygous predicted loss of function (pLoF) variants in BSN are associated with higher BMI with log10-p value of 11.78 in the UK biobank cohort. The association was replicated in the All of Us whole genome sequencing data. Additionally, we have identified two individuals (one of whom has a de novo variant) with a heterozygous pLoF variant in a cohort of early onset or extreme obesity at Columbia University. Like the individuals identified in the UKBB and All of us Cohorts, these individuals have no history of neurobehavioral or cognitive disability. Heterozygosity for pLoF BSN variants constitutes a new etiology for obesity.