Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

VarChat: the generative AI assistant for the interpretation of human genomic variations.

MOTIVATION: In the modern era of genomic research, the scientific community is witnessing an explosive growth in the volume of published findings. While this abundance of data offers invaluable insights, it also places a pressing responsibility on genetic professionals and researchers to stay informed about the latest findings and their clinical significance. Genomic variant interpretation is currently facing a challenge in identifying the most up-to-date and relevant scientific papers, while also extracting meaningful information to accelerate the process from clinical assessment to reporting. Computer-aided literature search and summarization can play a pivotal role in this context. By synthesizing complex genomic findings into concise, interpretable summaries, this approach facilitates the translation of extensive genomic datasets into clinically relevant insights. RESULTS: To bridge this gap, we present VarChat (varchat.engenome.com), an innovative tool based on generative AI, developed to find and summarize the fragmented scientific literature associated with genomic variants into brief yet informative texts. VarChat provides users with a concise description of specific genetic variants, detailing their impact on related proteins and possible effects on human health. In addition, VarChat offers direct links to related scientific trustable sources, and encourages deeper research. AVAILABILITY AND IMPLEMENTATION: varchat.engenome.com.

Molecular characterization of diffuse large B-cell lymphomas associated with hepatitis C virus infection.

Hepatitis C virus (HCV)-associated diffuse large B-cell lymphoma (DLBCL) displays peculiar clinicopathological characteristics, but its molecular landscape is not fully elucidated. In this study, we investigated the clinicopathological and molecular features of 54 patients with HCV-associated DLBCL. The median age was 71 years. An underlying marginal zone lymphoma component was detected in 14.8% of cases. FISH analysis showed rearrangements involving BCL6 in 50.9% of cases, MYC in 11.3% and BCL2 in 3.7%. Lymph2Cx-based assay was successful in 38 cases, recognizing 16 cases (42.1%) as ABC and 16 cases as GCB subtypes, while six resulted unclassified. ABC cases exhibited a higher lymphoma-related mortality (LRM). Next-generation sequencing analysis showed mutations in 158/184 evaluated genes. The most frequently mutated genes were KMT2D (42.6%), SETD1B (33.3%), RERE (29.4%), FAS and PIM1 (27.8%) and TBL1XR1 (25.9%). A mutation in the NOTCH pathway was detected in 25.9% of cases and was associated with worst LRM. Cluster analysis by LymphGen classified 29/54 cases within definite groups, including BN2 in 14 (48.2%), ST2 in seven (24.2%) and MCD and EZB in four each (13.8%). Overall, these results indicate a preferential marginal zone origin for a consistent subgroup of HCV-associated DLBCL cases and suggest potential implications for molecularly targeted therapies.

Cardiovascular Disease Burden, Mortality, and Sudden Death Risk in Epilepsy: A UK Biobank Study.

BACKGROUND: Sudden death is the leading cause of mortality in medically refractory epilepsy. Middle-aged persons with epilepsy (PWE) are under investigated regarding their mortality risk and burden of cardiovascular disease (CVD). METHODS: Using UK Biobank, we identified 7786 (1.6%) participants with diagnoses of epilepsy and 6,171,803 person-years of follow-up (mean 12.30 years, standard deviation 1.74); 566 patients with previous histories of stroke were excluded. The 7220 PWE comprised the study cohort with the remaining 494,676 without epilepsy as the comparator group. Prevalence of CVD was determined using validated diagnostic codes. Cox proportional hazards regression was used to assess all-cause mortality and sudden death risk. RESULTS: Hypertension, coronary artery disease, heart failure, valvular heart disease, and congenital heart disease were more prevalent in PWE. Arrhythmias including atrial fibrillation/flutter (12.2% vs 6.9%; P < 0.01), bradyarrhythmias (7.7% vs 3.5%; P < 0.01), conduction defects (6.1% vs 2.6%; P < 0.01), and ventricular arrhythmias (2.3% vs 1.0%; P < 0.01), as well as cardiac implantable electric devices (4.6% vs 2.0%; P < 0.01) were more prevalent in PWE. PWE had higher adjusted all-cause mortality (hazard ratio [HR], 3.9; 95% confidence interval [CI], 3.01-3.39), and sudden death-specific mortality (HR, 6.65; 95% CI, 4.53-9.77); and were almost 2 years younger at death (68.1 vs 69.8; P < 0.001). CONCLUSIONS: Middle-aged PWE have increased all-cause and sudden death-specific mortality and higher burden of CVD including arrhythmias and heart failure. Further work is required to elucidate mechanisms underlying all-cause mortality and sudden death risk in PWE of middle age, to identify prognostic biomarkers and develop preventative therapies in PWE.

Critical assessment of variant prioritization methods for rare disease diagnosis within the Rare Genomes Project.

Background: A major obstacle faced by rare disease families is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years, and causal variants are identified in under 50%. The Rare Genomes Project (RGP) is a direct-to-participant research study on the utility of genome sequencing (GS) for diagnosis and gene discovery. Families are consented for sharing of sequence and phenotype data with researchers, allowing development of a Critical Assessment of Genome Interpretation (CAGI) community challenge, placing variant prioritization models head-to-head in a real-life clinical diagnostic setting. Methods: Predictors were provided a dataset of phenotype terms and variant calls from GS of 175 RGP individuals (65 families), including 35 solved training set families, with causal variants specified, and 30 test set families (14 solved, 16 unsolved). The challenge tasked teams with identifying the causal variants in as many test set families as possible. Ranked variant predictions were submitted with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on rank position of true positive causal variants and maximum F-measure, based on precision and recall of causal variants across EPCR thresholds. Results: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performing teams recalled the causal variants in up to 13 of 14 solved families by prioritizing high quality variant calls that were rare, predicted deleterious, segregating correctly, and consistent with reported phenotype. In unsolved families, newly discovered diagnostic variants were returned to two families following confirmatory RNA sequencing, and two prioritized novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant, in an unsolved proband with phenotype overlap with asparagine synthetase deficiency. Conclusions: By objective assessment of variant predictions, we provide insights into current state-of-the-art algorithms and platforms for genome sequencing analysis for rare disease diagnosis and explore areas for future optimization. Identification of diagnostic variants in unsolved families promotes synergy between researchers with clinical and computational expertise as a means of advancing the field of clinical genome interpretation.

Selected memory T cells infused post-haploidentical hematopoietic stem cell transplantation persist and hyperexpand.

Haploidentical hematopoietic stem cell transplantation (haplo-HSCT) with post-transplant cyclophosphamide is a curative treatment for many hematological malignancies, yet a majority of patients still suffers from recurrent infections. Post-transplant infusion of memory T-cells could potentially enhance immunological protection without increasing the risk of eliciting acute graft-versus-host disease, which is mainly induced by naïve T-cells. Here, we performed longitudinal analysis of the lymphocyte compartment in 19 patients who underwent haplo-HSCT previously enrolled in a phase II prospective clinical trial (www.clinicaltrials.gov as #NCT04687982), in which they received post-transplant CD45RA-depleted donor lymphocyte infusions (DLI). T-cell receptor sequencing analysis showed that, surprisingly, CD45RA-depleted DLI do not increase T-cell clonal diversity, but lead to prominent expansion of a selected number of infused memory T-cell clones, suggesting recruitment of these cells in the immune response. Pathogen-specific memory T-cells, including cytomegalovirus (CMV)-specific cells, were engrafted and were able to persist for at least 1 month. Deep immunophenotyping revealed strong polyfunctional effector CMV-specific T-cell responses in the majority of patients, with their expansion correlating with the frequency of CMV-specific cells in the donor. These findings provide a rationale behind the suggested improved protection against viral infections in patients receiving CD45RA-depleted DLI.

Phenotypic Variation in Two Siblings Affected with Shwachman-Diamond Syndrome: The Use of Expert Variant Interpreter (eVai) Suggests Clinical Relevance of a Variant in the KMT2A Gene.

Introduction. Shwachman-Diamond Syndrome (SDS) is an autosomal-recessive disorder characterized by neutropenia, pancreatic exocrine insufficiency, skeletal dysplasia, and an increased risk for leukemic transformation. Biallelic mutations in the SBDS gene have been found in about 90% of patients. The clinical spectrum of SDS in patients is wide, and variability has been noticed between different patients, siblings, and even within the same patient over time. Herein, we present two SDS siblings (UPN42 and UPN43) carrying the same SBDS mutations and showing relevant differences in their phenotypic presentation. Study aim. We attempted to understand whether other germline variants, in addition to SBDS, could explain some of the clinical variability noticed between the siblings. Methods. Whole-exome sequencing (WES) was performed. Human Phenotype Ontology (HPO) terms were defined for each patient, and the WES data were analyzed using the eVai and DIVAs platforms. Results. In UPN43, we found and confirmed, using Sanger sequencing, a novel de novo variant (c.10663G > A, p.Gly3555Ser) in the KMT2A gene that is associated with autosomal-dominant Wiedemann−Steiner Syndrome. The variant is classified as pathogenic according to different in silico prediction tools. Interestingly, it was found to be related to some of the HPO terms that describe UPN43. Conclusions. We postulate that the KMT2A variant found in UPN43 has a concomitant and co-occurring clinical effect, in addition to SBDS mutation. This dual molecular effect, supported by in silico prediction, could help to understand some of the clinical variations found among the siblings. In the future, these new data are likely to be useful for personalized medicine and therapy for selected cases.

Single-cell profiling identifies impaired adaptive NK cells expanded after HCMV reactivation in haploidentical HSCT.

Haploidentical hematopoietic stem cell transplantation (h-HSCT) represents an efficient curative approach for patients affected by hematologic malignancies in which the reduced intensity conditioning induces a state of immunologic tolerance between donor and recipient. However, opportunistic viral infections greatly affect h-HSCT clinical outcomes. NK cells are the first lymphocytes that recover after transplant and provide a prompt defense against human cytomegalovirus (HCMV) infection/reactivation. By undertaking a longitudinal single-cell computational profiling of multiparametric flow cytometry, we show that HCMV accelerates NK cell immune reconstitution together with the expansion of CD158b1b2jpos/NKG2Aneg/NKG2Cpos/NKp30lo NK cells. The frequency of this subset correlates with HCMV viremia, further increases in recipients experiencing multiple episodes of viral reactivations, and persists for months after the infection. The transcriptional profile of FACS-sorted CD158b1b2jpos NK cells confirmed the ability of HCMV to deregulate NKG2C, NKG2A, and NKp30 gene expression, thus inducing the expansion of NK cells with adaptive traits. These NK cells are characterized by the downmodulation of several gene pathways associated with cell migration, the cell cycle, and effector-functions, as well as by a state of metabolic/cellular exhaustion. This profile reflects the functional impairments of adaptive NK cells to produce IFN-γ, a phenomenon also due to the viral-induced expression of lymphocyte-activation gene 3 (LAG-3) and programmed cell death protein 1 (PD-1) checkpoint inhibitors.

Feasibility and Efficacy of CD45RA+ Depleted Donor Lymphocytes Infusion After Haploidentical Transplantation With Post-Transplantation Cyclophosphamide in Patients With Hematological Malignancies.

Allogeneic stem cell transplantation from haploidentical donor using post-transplantation cyclophosphamide has been used to cure hematological diseases. Because of slow immunological reconstitution, there is an increased incidence of viral infection. The aim of our study was to prospectively evaluate the efficacy and the feasibility of a CD45RA+ depleted donor lymphocytes infusion (DLI) in terms of reduction of viral infection early after haploidentical transplantation. This a prospective single-center study. We enrolled 23 patients, of whom 19 were evaluable. Graft-versus-host disease (GVHD) prophylaxis was the same for all patients. The primary endpoint was 100-day cumulative incidence of viral infections. The primary endpoint was met, because the 100-day cumulative incidence of viral infection was 32%. The median time from transplantation to first CD45RA+ depleted DLI was 55 days (range, 46-63). 28% of patients had cytomegalovirus reactivation, no patients reactivated human herpesvirus-6; 1 patient developed BK virus related hemorrhagic cystitis. Most of the patients received the planned 3 infusions. Only 1 patient had development of grade 2 acute GVHD, and 2 patients had moderate chronic GVHD. All evaluable patients were off immunosuppressive therapy at last follow-up. The median follow-up was 12 months (range, 3-23), the 1-year overall survival and progression-free survival were 79% and 75%, respectively; the 100-day and 1-year non-relapse mortality were 5% and 12%, respectively. CD45RA+ depleted DLI are feasible in patients treated with haploidentical transplantation. The toxic profile is good with a low risk for development of both acute and chronic GVHD.

Two subsets of stem-like CD8+ memory T cell progenitors with distinct fate commitments in humans.

T cell memory relies on the generation of antigen-specific progenitors with stem-like properties. However, the identity of these progenitors has remained unclear, precluding a full understanding of the differentiation trajectories that underpin the heterogeneity of antigen-experienced T cells. We used a systematic approach guided by single-cell RNA-sequencing data to map the organizational structure of the human CD8+ memory T cell pool under physiological conditions. We identified two previously unrecognized subsets of clonally, epigenetically, functionally, phenotypically and transcriptionally distinct stem-like CD8+ memory T cells. Progenitors lacking the inhibitory receptors programmed death-1 (PD-1) and T cell immunoreceptor with Ig and ITIM domains (TIGIT) were committed to a functional lineage, whereas progenitors expressing PD-1 and TIGIT were committed to a dysfunctional, exhausted-like lineage. Collectively, these data reveal the existence of parallel differentiation programs in the human CD8+ memory T cell pool, with potentially broad implications for the development of immunotherapies and vaccines.

CXCR3 Identifies Human Naive CD8+ T Cells with Enhanced Effector Differentiation Potential.

In mice, the ability of naive T (TN) cells to mount an effector response correlates with TCR sensitivity for self-derived Ags, which can be quantified indirectly by measuring surface expression levels of CD5. Equivalent findings have not been reported previously in humans. We identified two discrete subsets of human CD8+ TN cells, defined by the absence or presence of the chemokine receptor CXCR3. The more abundant CXCR3+ TN cell subset displayed an effector-like transcriptional profile and expressed TCRs with physicochemical characteristics indicative of enhanced interactions with peptide-HLA class I Ags. Moreover, CXCR3+ TN cells frequently produced IL-2 and TNF in response to nonspecific activation directly ex vivo and differentiated readily into Ag-specific effector cells in vitro. Comparative analyses further revealed that human CXCR3+ TN cells were transcriptionally equivalent to murine CXCR3+ TN cells, which expressed high levels of CD5. These findings provide support for the notion that effector differentiation is shaped by heterogeneity in the preimmune repertoire of human CD8+ T cells.

Antioxidant metabolism regulates CD8+ T memory stem cell formation and antitumor immunity.

Adoptive T cell transfer (ACT) immunotherapy benefits from early differentiated stem cell memory T (Tscm) cells capable of persisting in the long term and generating potent antitumor effectors. Due to their paucity ex vivo, Tscm cells can be derived from naive precursors, but the molecular signals at the basis of Tscm cell generation are ill-defined. We found that less differentiated human circulating CD8+ T cells display substantial antioxidant capacity ex vivo compared with more differentiated central and effector memory T cells. Limiting ROS metabolism with antioxidants during naive T cell activation hindered terminal differentiation, while allowing expansion and generation of Tscm cells. N-acetylcysteine (NAC), the most effective molecule in this regard, induced transcriptional and metabolic programs characteristic of self-renewing memory T cells. Upon ACT, NAC-generated Tscm cells established long-term memory in vivo and exerted more potent antitumor immunity in a xenogeneic model when redirected with CD19-specific CAR, highlighting the translational relevance of NAC as a simple and inexpensive method to improve ACT.

Background fluorescence and spreading error are major contributors of variability in high-dimensional flow cytometry data visualization by t-distributed stochastic neighboring embedding.

Multidimensional single-cell analysis requires approaches to visualize complex data in intuitive 2D graphs. In this regard, t-distributed stochastic neighboring embedding (tSNE) is the most popular algorithm for single-cell RNA sequencing and cytometry by time-of-flight (CyTOF), but its application to polychromatic flow cytometry, including the recently developed 30-parameter platform, is still under investigation. We identified differential distribution of background values between samples, generated by either background calculation or spreading error (SE), as a major source of variability in polychromatic flow cytometry data representation by tSNE, ultimately resulting in the identification of erroneous heterogeneity among cell populations. Biexponential transformation of raw data and limiting SE during panel development dramatically improved data visualization. These aspects must be taken into consideration when using computational approaches as discovery tools in large sets of samples from independent experiments or immunomonitoring in clinical trials.

Curtailed T-cell activation curbs effector differentiation and generates CD8+ T cells with a naturally-occurring memory stem cell phenotype.

Human T memory stem (TSCM ) cells with superior persistence capacity and effector functions are emerging as important players in the maintenance of long-lived T-cell memory and are thus considered an attractive population to be used in adoptive transfer-based immunotherapy of cancer. However, the molecular signals regulating their generation remain poorly defined. Here we show that curtailed T-cell receptor stimulation curbs human effector CD8+ T-cell differentiation and allows the generation of CD45RO- CD45RA+ CCR7+ CD27+ CD95+ -phenotype cells from highly purified naïve T-cell precursors, resembling naturally-occurring human TSCM . These cells proliferate extensively in vitro and in vivo, express low amounts of effector-associated genes and transcription factors and undergo considerable self-renewal in response to IL-15 while retaining effector differentiation potential. Such a phenotype is associated with a lower number of mitochondria compared to highly-activated effector T cells committed to terminal differentiation. These results shed light on the molecular signals that are required to generate long-lived memory T cells with potential application in adoptive cell transfer immunotherapy.

Transducin-like enhancer of split-1 is expressed and functional in human macrophages.

Macrophages display heterogeneous phenotypes, including the classical M1 proinflammatory and the alternative M2 anti-inflammatory polarization states. The transducin-like enhancer of split-1 (TLE1) is a transcriptional corepressor whose functions in macrophages have not been studied yet. We report that TLE1 is highly expressed in human alternative macrophages in vitro and in atherosclerotic plaques as well as in adipose tissue M1/M2 mixed macrophages. TLE1 silencing in alternative macrophages decreases the expression of the M2 markers IL-1Ra and IL-10, while it exacerbates TNFα and CCL3 induction by lipopolysaccharide. Hence, TLE1 is expressed in human macrophages where it has potential anti-inflammatory and alternative phenotype promoting properties.

Macrophage phenotypes and their modulation in atherosclerosis.

Atherosclerosis is the result of a chronic inflammatory response in the arterial wall related to uptake of low-density lipoprotein by macrophages and their subsequent transformation in foam cells. Monocyte-derived macrophages are the principal mediators of tissue homeostasis and repair, response to pathogens and inflammation. However, macrophages are a homogeneous cell population presenting a continuum phenotypic spectrum with, at the extremes, the classically Th-1 polarized M1 and alternatively Th-2 polarized M2 macrophage phenotypes, which have been well described. Moreover, M2 macrophages also present several subtypes often termed M2a, b, c and d, each of them expressing specific markers and exhibiting specialized properties. Macrophage plasticity is mirrored also in the atherosclerotic lesions, where different stimuli can influence the phenotype giving rise to a complex system of subpopulations, such as Mox, Mhem, M(Hb) and M4 macrophages. An abundant literature has described the potential modulators of the reciprocal skewing between pro-inflammatory M1 and anti-inflammatory M2 macrophages including lesion stage and localization, miRNA, transcription factors such as PPARγ, KLF4 and NR4A family members, high-density lipoproteins and plaque lipid content, pathways such as the rapamycin-mTOR1 pathway, molecules such as thioredoxin-1, infection by helminths and irradiation. We hope to provide an overview of the macrophage phenotype complexity in cardiovascular diseases, particularly atherosclerosis.

miR-206 controls LXRα expression and promotes LXR-mediated cholesterol efflux in macrophages.

Liver X receptors (LXRα and LXRß) are key transcription factors in cholesterol metabolism that regulate cholesterol biosynthesis/efflux and bile acid metabolism/excretion in the liver and numerous organs. In macrophages, LXR signaling modulates cholesterol handling and the inflammatory response, pathways involved in atherosclerosis. Since regulatory pathways of LXR transcription control are well understood, in the present study we aimed at identifying post-transcriptional regulators of LXR activity. MicroRNAs (miRs) are such post-transcriptional regulators of genes that in the canonical pathway mediate mRNA inactivation. In silico analysis identified miR-206 as a putative regulator of LXRα but not LXRß. Indeed, as recently shown, we found that miR-206 represses LXRα activity and expression of LXRα and its target genes in hepatic cells. Interestingly, miR-206 regulates LXRα differently in macrophages. Stably overexpressing miR-206 in THP-1 human macrophages revealed an up-regulation and miR-206 knockdown led to a down-regulation of LXRα and its target genes. In support of these results, bone marrow-derived macrophages (BMDMs) from miR-206 KO mice also exhibited lower expression of LXRα target genes. The physiological relevance of these findings was proven by gain- and loss-of-function of miR-206; overexpression of miR-206 enhanced cholesterol efflux in human macrophages and knocking out miR-206 decreased cholesterol efflux from MPMs. Moreover, we show that miR-206 expression in macrophages is repressed by LXRα activation, while oxidized LDL and inflammatory stimuli profoundly induced miR-206 expression. We therefore propose a feed-back loop between miR-206 and LXRα that might be part of an LXR auto-regulatory mechanism to fine tune LXR activity.