Functional interrogation of DNA damage response variants with base editing screens.

Mutations in DNA damage response (DDR) genes endanger genome integrity and predispose to cancer and genetic disorders. Here, using CRISPR-dependent cytosine base editing screens, we identify > 2,000 sgRNAs that generate nucleotide variants in 86 DDR genes, resulting in altered cellular fitness upon DNA damage. Among those variants, we discover loss- and gain-of-function mutants in the Tudor domain of the DDR regulator 53BP1 that define a non-canonical surface required for binding the deubiquitinase USP28. Moreover, we characterize variants of the TRAIP ubiquitin ligase that define a domain, whose loss renders cells resistant to topoisomerase I inhibition. Finally, we identify mutations in the ATM kinase with opposing genome stability phenotypes and loss-of-function mutations in the CHK2 kinase previously categorized as variants of uncertain significance for breast cancer. We anticipate that this resource will enable the discovery of additional DDR gene functions and expedite studies of DDR variants in human disease.

Functional analysis of germline VANGL2 variants using rescue assays of vangl2 knockout zebrafish.

Developmental studies have shown that the evolutionarily conserved Wnt Planar Cell Polarity (PCP) pathway is essential for the development of a diverse range of tissues and organs including the brain, spinal cord, heart and sensory organs, as well as establishment of the left-right body axis. Germline mutations in the highly conserved PCP gene VANGL2 in humans have only been associated with central nervous system malformations, and functional testing to understand variant impact has not been performed. Here we report three new families with missense variants in VANGL2 associated with heterotaxy and congenital heart disease p.(Arg169His), non-syndromic hearing loss p.(Glu465Ala) and congenital heart disease with brain defects p.(Arg135Trp). To test the in vivo impact of these and previously described variants, we have established clinically-relevant assays using mRNA rescue of the vangl2 mutant zebrafish. We show that all variants disrupt Vangl2 function, although to different extents and depending on the developmental process. We also begin to identify that different VANGL2 missense variants may be haploinsufficient and discuss evidence in support of pathogenicity. Together, this study demonstrates that zebrafish present a suitable pipeline to investigate variants of unknown significance and suggests new avenues for investigation of the different developmental contexts of VANGL2 function that are clinically meaningful.

Genetic insights into the age-specific biological mechanisms governing human ovarian aging.

There is currently little evidence that the genetic basis of human phenotype varies significantly across the lifespan. However, time-to-event phenotypes are understudied and can be thought of as reflecting an underlying hazard, which is unlikely to be constant through life when values take a broad range. Here, we find that 74% of 245 genome-wide significant genetic associations with age at natural menopause (ANM) in the UK Biobank show a form of age-specific effect. Nineteen of these replicated discoveries are identified only by our modeling framework, which determines the time dependency of DNA-variant age-at-onset associations without a significant multiple-testing burden. Across the range of early to late menopause, we find evidence for significantly different underlying biological pathways, changes in the signs of genetic correlations of ANM to health indicators and outcomes, and differences in inferred causal relationships. We find that DNA damage response processes only act to shape ovarian reserve and depletion for women of early ANM. Genetically mediated delays in ANM were associated with increased relative risk of breast cancer and leiomyoma at all ages and with high cholesterol and heart failure for late-ANM women. These findings suggest that a better understanding of the age dependency of genetic risk factor relationships among health indicators and outcomes is achievable through appropriate statistical modeling of large-scale biobank data.

Estimating clinical risk in gene regions from population sequencing cohort data.

While pathogenic variants can significantly increase disease risk, it is still challenging to estimate the clinical impact of rare missense variants more generally. Even in genes such as BRCA2 or PALB2, large cohort studies find no significant association between breast cancer and rare missense variants collectively. Here, we introduce REGatta, a method to estimate clinical risk from variants in smaller segments of individual genes. We first define these regions by using the density of pathogenic diagnostic reports and then calculate the relative risk in each region by using over 200,000 exome sequences in the UK Biobank. We apply this method in 13 genes with established roles across several monogenic disorders. In genes with no significant difference at the gene level, this approach significantly separates disease risk for individuals with rare missense variants at higher or lower risk (BRCA2 regional model OR = 1.46 [1.12, 1.79], p = 0.0036 vs. BRCA2 gene model OR = 0.96 [0.85, 1.07] p = 0.4171). We find high concordance between these regional risk estimates and high-throughput functional assays of variant impact. We compare our method with existing methods and the use of protein domains (Pfam) as regions and find REGatta better identifies individuals at elevated or reduced risk. These regions provide useful priors and are potentially useful for improving risk assessment for genes associated with monogenic diseases.

Significance tests for R2 of out-of-sample prediction using polygenic scores.

The coefficient of determination (R2) is a well-established measure to indicate the predictive ability of polygenic scores (PGSs). However, the sampling variance of R2 is rarely considered so that 95% confidence intervals (CI) are not usually reported. Moreover, when comparisons are made between PGSs based on different discovery samples, the sampling covariance of R2 is required to test the difference between them. Here, we show how to estimate the variance and covariance of R2 values to assess the 95% CI and p value of the R2 difference. We apply this approach to real data calculating PGSs in 28,880 European participants derived from UK Biobank (UKBB) and Biobank Japan (BBJ) GWAS summary statistics for cholesterol and BMI. We quantify the significantly higher predictive ability of UKBB PGSs compared to BBJ PGSs (p value 7.6e-31 for cholesterol and 1.4e-50 for BMI). A joint model of UKBB and BBJ PGSs significantly improves the predictive ability, compared to a model of UKBB PGS only (p value 3.5e-05 for cholesterol and 1.3e-28 for BMI). We also show that the predictive ability of regulatory SNPs is significantly enriched over non-regulatory SNPs for cholesterol (p value 8.9e-26 for UKBB and 3.8e-17 for BBJ). We suggest that the proposed approach (available in R package r2redux) should be used to test the statistical significance of difference between pairs of PGSs, which may help to draw a correct conclusion about the comparative predictive ability of PGSs.

Identification of hidden associations among eukaryotic genes through statistical analysis of coevolutionary transitions.

Coevolution at the gene level, as reflected by correlated events of gene loss or gain, can be revealed by phylogenetic profile analysis. The optimal method and metric for comparing phylogenetic profiles, especially in eukaryotic genomes, are not yet established. Here, we describe a procedure suitable for large-scale analysis, which can reveal coevolution based on the assessment of the statistical significance of correlated presence/absence transitions between gene pairs. This metric can identify coevolution in profiles with low overall similarities and is not affected by similarities lacking coevolutionary information. We applied the procedure to a large collection of 60,912 orthologous gene groups (orthogroups) in 1,264 eukaryotic genomes extracted from OrthoDB. We found significant cotransition scores for 7,825 orthogroups associated in 2,401 coevolving modules linking known and unknown genes in protein complexes and biological pathways. To demonstrate the ability of the method to predict hidden gene associations, we validated through experiments the involvement of vertebrate malate synthase-like genes in the conversion of (S)-ureidoglycolate into glyoxylate and urea, the last step of purine catabolism. This identification explains the presence of glyoxylate cycle genes in metazoa and suggests an anaplerotic role of purine degradation in early eukaryotes.

Bayes factor functions for reporting outcomes of hypothesis tests.

Bayes factors represent a useful alternative to P-values for reporting outcomes of hypothesis tests by providing direct measures of the relative support that data provide to competing hypotheses. Unfortunately, the competing hypotheses have to be specified, and the calculation of Bayes factors in high-dimensional settings can be difficult. To address these problems, we define Bayes factor functions (BFFs) directly from common test statistics. BFFs depend on a single noncentrality parameter that can be expressed as a function of standardized effects, and plots of BFFs versus effect size provide informative summaries of hypothesis tests that can be easily aggregated across studies. Such summaries eliminate the need for arbitrary P-value thresholds to define "statistical significance." Because BFFs are defined using nonlocal alternative prior densities, they provide more rapid accumulation of evidence in favor of true null hypotheses without sacrificing efficiency in supporting true alternative hypotheses. BFFs can be expressed in closed form and can be computed easily from z, t, χ2, and F statistics.

Evolutionary history of MEK1 illuminates the nature of deleterious mutations.

Mutations in signal transduction pathways lead to various diseases including cancers. MEK1 kinase, encoded by the human MAP2K1 gene, is one of the central components of the MAPK pathway and more than a hundred somatic mutations in the MAP2K1 gene were identified in various tumors. Germline mutations deregulating MEK1 also lead to congenital abnormalities, such as the cardiofaciocutaneous syndrome and arteriovenous malformation. Evaluating variants associated with a disease is a challenge, and computational genomic approaches aid in this process. Establishing evolutionary history of a gene improves computational prediction of disease-causing mutations; however, the evolutionary history of MEK1 is not well understood. Here, by revealing a precise evolutionary history of MEK1, we construct a well-defined dataset of MEK1 metazoan orthologs, which provides sufficient depth to distinguish between conserved and variable amino acid positions. We matched known and predicted disease-causing and benign mutations to evolutionary changes observed in corresponding amino acid positions and found that all known and many suspected disease-causing mutations are evolutionarily intolerable. We selected several variants that cannot be unambiguously assessed by automated prediction tools but that are confidently identified as "damaging" by our approach, for experimental validation in Drosophila. In all cases, evolutionary intolerant variants caused increased mortality and severe defects in fruit fly embryos confirming their damaging nature. We anticipate that our analysis will serve as a blueprint to help evaluate known and novel missense variants in MEK1 and that our approach will contribute to improving automated tools for disease-associated variant interpretation.

In vivo versus in silico assessment of potentially pathogenic missense variants in human reproductive genes.

Infertility is a heterogeneous condition, with genetic causes thought to underlie a substantial fraction of cases. Genome sequencing is becoming increasingly important for genetic diagnosis of diseases including idiopathic infertility; however, most rare or minor alleles identified in patients are variants of uncertain significance (VUS). Interpreting the functional impacts of VUS is challenging but profoundly important for clinical management and genetic counseling. To determine the consequences of these variants in key fertility genes, we functionally evaluated 11 missense variants in the genes ANKRD31, BRDT, DMC1, EXO1, FKBP6, MCM9, M1AP, MEI1, MSH4 and SEPT12 by generating genome-edited mouse models. Nine variants were classified as deleterious by most functional prediction algorithms, and two disrupted a protein-protein interaction (PPI) in the yeast two hybrid (Y2H) assay. Though these genes are essential for normal meiosis or spermiogenesis in mice, only one variant, observed in the MCM9 gene of a male infertility patient, compromised fertility or gametogenesis in the mouse models. To explore the disconnect between predictions and outcomes, we compared pathogenicity calls of missense variants made by ten widely used algorithms to 1) those annotated in ClinVar and 2) those evaluated in mice. All the algorithms performed poorly in terms of predicting the effects of human missense variants modeled in mice. These studies emphasize caution in the genetic diagnoses of infertile patients based primarily on pathogenicity prediction algorithms and emphasize the need for alternative and efficient in vitro or in vivo functional validation models for more effective and accurate VUS description to either pathogenic or benign categories.

The functional impact of BRCA1 BRCT domain variants using multiplexed DNA double-strand break repair assays.

Pathogenic variants in BRCA1 are associated with a greatly increased risk of hereditary breast and ovarian cancer (HBOC). With the increased availability and affordability of genetic testing, many individuals have been identified with BRCA1 variants of uncertain significance (VUSs), which are individually detected in the population too infrequently to ascertain a clinical risk. Functional assays can be used to experimentally assess the effects of these variants. In this study, we used multiplexed DNA repair assays of variants in the BRCA1 carboxyl terminus to functionally characterize 2,271 variants for homology-directed repair function (HDR) and 1,427 variants for cisplatin resistance (CR). We found a high level of consistent results (Pearson's r = 0.74) in the two multiplexed functional assays with non-functional variants located within regions of the BRCA1 protein necessary for its tumor suppression activity. In addition, functional categorizations of variants tested in the multiplex HDR and CR assays correlated with known clinical significance and with other functional assays for BRCA1 (Pearson's r = 0.53 to 0.71). The results of the multiplex HDR and CR assays are useful resources for characterizing large numbers of BRCA1 VUSs.

A calibrated functional patch-clamp assay to enhance clinical variant interpretation in KCNH2-related long QT syndrome.

Modern sequencing technologies have revolutionized our detection of gene variants. However, in most genes, including KCNH2, the majority of missense variants are currently classified as variants of uncertain significance (VUSs). The aim of this study was to investigate the utility of an automated patch-clamp assay for aiding clinical variant classification in KCNH2. The assay was designed according to recommendations proposed by the Clinical Genome Sequence Variant Interpretation Working Group. Thirty-one variants (17 pathogenic/likely pathogenic, 14 benign/likely benign) were classified internally as variant controls. They were heterozygously expressed in Flp-In HEK293 cells for assessing the effects of variants on current density and channel gating in order to determine the sensitivity and specificity of the assay. All 17 pathogenic variant controls had reduced current density, and 13 of 14 benign variant controls had normal current density, which enabled determination of normal and abnormal ranges for applying evidence of moderate or supporting strength for VUS reclassification. Inclusion of functional assay evidence enabled us to reclassify 6 out of 44 KCNH2 VUSs as likely pathogenic. The high-throughput patch-clamp assay can provide moderate-strength evidence for clinical interpretation of clinical KCNH2 variants and demonstrates the value of developing automated patch-clamp assays for functional characterization of ion channel gene variants.

Identifying local associations in biological time series: algorithms, statistical significance, and applications.

Local associations refer to spatial-temporal correlations that emerge from the biological realm, such as time-dependent gene co-expression or seasonal interactions between microbes. One can reveal the intricate dynamics and inherent interactions of biological systems by examining the biological time series data for these associations. To accomplish this goal, local similarity analysis algorithms and statistical methods that facilitate the local alignment of time series and assess the significance of the resulting alignments have been developed. Although these algorithms were initially devised for gene expression analysis from microarrays, they have been adapted and accelerated for multi-omics next generation sequencing datasets, achieving high scientific impact. In this review, we present an overview of the historical developments and recent advances for local similarity analysis algorithms, their statistical properties, and real applications in analyzing biological time series data. The benchmark data and analysis scripts used in this review are freely available at http://github.com/labxscut/lsareview.

Classification of PTEN missense VUS through exascale simulations.

Phosphatase and tensin homolog (PTEN), a tumor suppressor with dual phosphatase properties, is a key factor in PI3K/AKT signaling pathway. Pathogenic germline variation in PTEN can abrogate its ability to dephosphorylate, causing high cancer risk. Lack of functional evidence lets numerous PTEN variants be classified as variants of uncertain significance (VUS). Utilizing Molecular Dynamics (MD) simulations, we performed a thorough evaluation for 147 PTEN missense VUS, sorting them into 66 deleterious and 81 tolerated variants. Utilizing replica exchange molecular dynamic (REMD) simulations, we further assessed the variants situated in the catalytic core of PTEN's phosphatase domain and uncovered conformational alterations influencing the structural stability of the phosphatase domain. There was a high degree of agreement between our results and the variants classified by Variant Abundance by Massively Parallel Sequencing, saturation mutagenesis, multiplexed functional data and experimental assays. Our extensive analysis of PTEN missense VUS should benefit their clinical applications in PTEN-related cancer. SIGNIFICANCE STATEMENT: Classification of PTEN variants affecting its lipid phosphatase activity is important for understanding the roles of PTEN variation in the pathogenesis of hereditary and sporadic malignancies. Of the 3000 variants identified in PTEN, 1296 (43%) were assigned as VUS. Here, we applied MD and REMD simulations to investigate the effects of PTEN missense VUS on the structural integrity of the PTEN phosphatase domain consisting the WPD, P and TI active sites. We classified a total of 147 missense VUS into 66 deleterious and 81 tolerated variants by referring to the control group comprising 54 pathogenic and 12 benign variants. The classification was largely in concordance with these classified by experimental approaches.

Homologous recombination-deficient mutation cluster in tumor suppressor RAD51C identified by comprehensive analysis of cancer variants.

Mutations in homologous recombination (HR) genes, including BRCA1, BRCA2, and the RAD51 paralog RAD51C, predispose to tumorigenesis and sensitize cancers to DNA-damaging agents and poly(ADP ribose) polymerase inhibitors. However, ∼800 missense variants of unknown significance have been identified for RAD51C alone, impairing cancer risk assessment and therapeutic strategies. Here, we interrogated >50 RAD51C missense variants, finding that mutations in residues conserved with RAD51 strongly predicted HR deficiency and disrupted interactions with other RAD51 paralogs. A cluster of mutations was identified in and around the Walker A box that led to impairments in HR, interactions with three other RAD51 paralogs, binding to single-stranded DNA, and ATP hydrolysis. We generated structural models of the two RAD51 paralog complexes containing RAD51C, RAD51B-RAD51C-RAD51D-XRCC2 and RAD51C-XRCC3. Together with our functional and biochemical analyses, the structural models predict ATP binding at the interface of RAD51C interactions with other RAD51 paralogs, similar to interactions between monomers in RAD51 filaments, and explain the failure of RAD51C variants in binding multiple paralogs. Ovarian cancer patients with variants in this cluster showed exceptionally long survival, which may be relevant to the reversion potential of the variants. This comprehensive analysis provides a framework for RAD51C variant classification. Importantly, it also provides insight into the functioning of the RAD51 paralog complexes.

PLCG2-associated immune dysregulation (PLAID) comprises broad and distinct clinical presentations related to functional classes of genetic variants.

BACKGROUND: Pathogenic variants of phospholipase C gamma 2 (PLCG2) cause 2 related forms of autosomal-dominant immune dysregulation (ID), PLCγ2-associated antibody deficiency and immune dysregulation (PLAID) and autoinflammatory PLAID (APLAID). Since describing these conditions, many PLCG2 variants of uncertain significance have been identified by clinical sequencing of patients with diverse features of ID. OBJECTIVE: We sought to functionally classify PLCG2 variants and explore known and novel genotype-function-phenotype relationships. METHODS: Clinical data from patients with PLCG2 variants were obtained via standardized questionnaire. PLCG2 variants were generated by mutagenesis of enhanced green fluorescent protein (EGFP)-PLCG2 plasmid, which was overexpressed in Plcg2-deficient DT-40 B cells. B-cell receptor-induced calcium flux and extracellular signal-regulated kinase phosphorylation were assayed by flow cytometry. In some cases, stimulation-induced calcium flux was also measured in primary patient cells. RESULTS: Three-fourths of PLCG2 variants produced functional alteration of B-cell activation, in vitro. Thirteen variants led to gain of function (GOF); however, most functional variants defined a new class of PLCG2 mutation, monoallelic loss of function (LOF). Susceptibility to infection and autoinflammation were common with both GOF and LOF variants, whereas a new phenotypic cluster consisting of humoral immune deficiency, autoinflammation, susceptibility to herpesvirus infection, and natural killer cell dysfunction was observed in association with multiple heterozygous LOF variants detected in both familial and sporadic cases. In some cases, PLCG2 variants produced greater effects in natural killer cells than in B cells. CONCLUSIONS: This work expands the genotypic and phenotypic associations with functional variation in PLCG2, including a novel form of ID in carriers of heterozygous loss of PLCG2 function. It also demonstrates the need for more diverse assays for assessing the impact of PLCG2 variants on human disease.

A pan-cancer analysis of prognostic significance and immunological role of lysosomal-associated membrane protein 3.

Lysosomal dysfunction can drive carcinogenesis. Lysosomal-associated membrane protein 3 (LAMP3), is a member of the Lysosome Associated Membrane Proteins and is involved in the malignant phenotype such as tumour metastasis and drug resistance, while the mechanisms that regulate the malignant progression of tumour remain vague. Our study aims to provide a more systematic and comprehensive understanding of the role of LAMP3 in the progression of various cancers by various databases.We explored the role of LAMP3 in pan-cancer using The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) database. Multiple online web platforms and software were used for data analysis, including HPA, TIMER, TISIDB, GEPIA, UALCAN, Kaplan-Meier plotter, DAVID and TIGER. The immunohistochemistry was used to quantify the LAMP3 and PD-L1 expression levels in cancer.High LAMP3 expression was found in most cancers and differentially expressed across molecular and immune subtypes. The expression of LAMP3 was involved in the immune-associated processes of Antigen processing and presentation, Th17 cell differentiation, Th1 and Th2 cell differentiation, and the immune-associated pathways of T cell receptor and B cell receptor signalling pathways in most cancers. It also correlated with genetic markers of immunomodulators in various cancers. LAMP3 and PD-L1 expression in BRCA and HNSC tissues was higher than that in corresponding adjacent normal tissues by immunohistochemistry. There is a significant correlation between the expression of LAMP3 and PD-L1.Our study elucidates that LAMP3 has different expression patterns and genetic alteration patterns in different tumours. It is a potential biomarker for immune-related cancer diagnosis, prognosis and efficacy prediction.

Structural implications of amyloidogenic rare variants Ser282Leu and Gln356Arg identified in h-BRCA1.

Preliminary studies have shown BRCA1 (170-1600) residues to be intrinsically disordered with unknown structural details. However, thousands of clinically reported variants have been identified in this central region of BRCA1. Therefore, we aimed to characterize h-BRCA1(260-553) to assess the structural basis for pathogenicity of two rare missense variants Ser282Leu, Gln356Arg identified from the Indian and Russian populations respectively. Small-angle X-ray scattering analysis revealed WT scores Rg -32 Å, Dmax -93 Å, and Rflex-51% which are partially disordered, whereas Ser282Leu variant displayed a higher degree of disorderedness and Gln356Arg was observed to be aggregated. WT protein also possesses an inherent propensity to undergo a disorder-to-order transition in the presence of cruciform DNA and 2,2,2-Trifluoroethanol (TFE). An increased alpha-helical pattern was observed with increasing concentration of TFE for the Gln356Arg mutant whereas Ser282Leu mutant showed significant differences only at the highest TFE concentration. Furthermore, higher thermal shift was observed for WT-DNA complex compared to the Gln356Arg and Ser282Leu protein-DNA complex. Moreover, mature amyloid-like fibrils were observed with 30 µM thioflavin T (ThT) at 37°C for Ser282Leu and Gln356Arg proteins while the WT protein exists in a protofibril state as observed by TEM. Gln356Arg formed higher-order aggregates with amyloidogenesis over time as monitored by ThT fluorescence. In addition, computational analyses confirmed larger conformational fluctuations for Ser282Leu and Gln356Arg mutants than for the WT. The global structural alterations caused by these variants provide a mechanistic approach for further classification of the variants of uncertain clinical significance in BRCA1 into amyloidogenic variants which may have a significant role in disease pathogenesis.

Decoding HiPSC-CM's Response to SARS-CoV-2: mapping the molecular landscape of cardiac injury.

BACKGROUND: Acute cardiac injury caused by coronavirus disease 2019 (COVID-19) increases mortality. Acute cardiac injury caused by COVID-19 requires understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) directly infects cardiomyocytes. This study provides a solid foundation for related studies by using a model of SARS-CoV-2 infection in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) at the transcriptome level, highlighting the relevance of this study to related studies. SARS-CoV-2 infection in hiPSC-CMs has previously been studied by bioinformatics without presenting the full molecular biological process. We present a unique bioinformatics view of the complete molecular biological process of SARS-CoV-2 infection in hiPSC-CMs. METHODS: To validate the RNA-seq datasets, we used GSE184715 and GSE150392 for the analytical studies, GSE193722 for validation at the cellular level, and GSE169241 for validation in heart tissue samples. GeneCards and MsigDB databases were used to find genes associated with the phenotype. In addition to differential expression analysis and principal component analysis (PCA), we also performed protein-protein interaction (PPI) analysis, functional enrichment analysis, hub gene analysis, upstream transcription factor prediction, and drug prediction. RESULTS: Differentially expressed genes (DEGs) were classified into four categories: cardiomyocyte cytoskeletal protein inhibition, proto-oncogene activation and inflammation, mitochondrial dysfunction, and intracellular cytoplasmic physiological function. Each of the hub genes showed good diagnostic prediction, which was well validated in other datasets. Inhibited biological functions included cardiomyocyte cytoskeletal proteins, adenosine triphosphate (ATP) synthesis and electron transport chain (ETC), glucose metabolism, amino acid metabolism, fatty acid metabolism, pyruvate metabolism, citric acid cycle, nucleic acid metabolism, replication, transcription, translation, ubiquitination, autophagy, and cellular transport. Proto-oncogenes, inflammation, nuclear factor-kappaB (NF-κB) pathways, and interferon signaling were activated, as well as inflammatory factors. Viral infection activates multiple pathways, including the interferon pathway, proto-oncogenes and mitochondrial oxidative stress, while inhibiting cardiomyocyte backbone proteins and energy metabolism. Infection limits intracellular synthesis and metabolism, as well as the raw materials for mitochondrial energy synthesis. Mitochondrial dysfunction and energy abnormalities are ultimately caused by proto-oncogene activation and SARS-CoV-2 infection. Activation of the interferon pathway, proto-oncogene up-regulation, and mitochondrial oxidative stress cause the inflammatory response and lead to diminished cardiomyocyte contraction. Replication, transcription, translation, ubiquitination, autophagy, and cellular transport are among the functions that decline physiologically. CONCLUSION: SARS-CoV-2 infection in hiPSC-CMs is fundamentally mediated via mitochondrial dysfunction. Therapeutic interventions targeting mitochondrial dysfunction may alleviate the cardiovascular complications associated with SARS-CoV-2 infection.

Deciphering pathogenicity of variants of uncertain significance with CRISPR-edited iPSCs.

Genetic variants play an important role in conferring risk for cardiovascular diseases (CVDs). With the rapid development of next-generation sequencing (NGS), thousands of genetic variants associated with CVDs have been identified by genome-wide association studies (GWAS), but the function of more than 40% of genetic variants is still unknown. This gap of knowledge is a barrier to the clinical application of the genetic information. However, determining the pathogenicity of a variant of uncertain significance (VUS) is challenging due to the lack of suitable model systems and accessible technologies. By combining clustered regularly interspaced short palindromic repeats (CRISPR) and human induced pluripotent stem cells (iPSCs), unprecedented advances are now possible in determining the pathogenicity of VUS in CVDs. Here, we summarize recent progress and new strategies in deciphering pathogenic variants for CVDs using CRISPR-edited human iPSCs.

Classification of 101 BRCA1 and BRCA2 variants of uncertain significance by cosegregation study: A powerful approach.

Up to 80% of BRCA1 and BRCA2 genetic variants remain of uncertain clinical significance (VUSs). Only variants classified as pathogenic or likely pathogenic can guide breast and ovarian cancer prevention measures and treatment by PARP inhibitors. We report the first results of the ongoing French national COVAR (cosegregation variant) study, the aim of which is to classify BRCA1/2 VUSs. The classification method was a multifactorial model combining different associations between VUSs and cancer, including cosegregation data. At this time, among the 653 variants selected, 101 (15%) distinct variants shared by 1,624 families were classified as pathogenic/likely pathogenic or benign/likely benign by the COVAR study. Sixty-six of the 101 (65%) variants classified by COVAR would have remained VUSs without cosegregation data. Of note, among the 34 variants classified as pathogenic by COVAR, 16 remained VUSs or likely pathogenic when following the ACMG/AMP variant classification guidelines. Although the initiation and organization of cosegregation analyses require a considerable effort, the growing number of available genetic tests results in an increasing number of families sharing a particular variant, and thereby increases the power of such analyses. Here we demonstrate that variant cosegregation analyses are a powerful tool for the classification of variants in the BRCA1/2 breast-ovarian cancer predisposition genes.