Discovering predisposing genes for hereditary breast cancer using deep learning.

Breast cancer (BC) is the most common malignancy affecting Western women today. It is estimated that as many as 10% of BC cases can be attributed to germline variants. However, the genetic basis of the majority of familial BC cases has yet to be identified. Discovering predisposing genes contributing to familial BC is challenging due to their presumed rarity, low penetrance, and complex biological mechanisms. Here, we focused on an analysis of rare missense variants in a cohort of 12 families of Middle Eastern origins characterized by a high incidence of BC cases. We devised a novel, high-throughput, variant analysis pipeline adapted for family studies, which aims to analyze variants at the protein level by employing state-of-the-art machine learning models and three-dimensional protein structural analysis. Using our pipeline, we analyzed 1218 rare missense variants that are shared between affected family members and classified 80 genes as candidate pathogenic. Among these genes, we found significant functional enrichment in peroxisomal and mitochondrial biological pathways which segregated across seven families in the study and covered diverse ethnic groups. We present multiple evidence that peroxisomal and mitochondrial pathways play an important, yet underappreciated, role in both germline BC predisposition and BC survival.

Differentiation shifts from a reversible to an irreversible heterochromatin state at the DM1 locus.

Epigenetic defects caused by hereditary or de novo mutations are implicated in various human diseases. It remains uncertain whether correcting the underlying mutation can reverse these defects in patient cells. Here we show by the analysis of myotonic dystrophy type 1 (DM1)-related locus that in mutant human embryonic stem cells (hESCs), DNA methylation and H3K9me3 enrichments are completely abolished by repeat excision (CTG2000 expansion), whereas in patient myoblasts (CTG2600 expansion), repeat deletion fails to do so. This distinction between undifferentiated and differentiated cells arises during cell differentiation, and can be reversed by reprogramming of gene-edited myoblasts. We demonstrate that abnormal methylation in DM1 is distinctively maintained in the undifferentiated state by the activity of the de novo DNMTs (DNMT3b in tandem with DNMT3a). Overall, the findings highlight a crucial difference in heterochromatin maintenance between undifferentiated (sequence-dependent) and differentiated (sequence-independent) cells, thus underscoring the role of differentiation as a locking mechanism for repressive epigenetic modifications at the DM1 locus.

Using multi-scale genomics to associate poorly annotated genes with rare diseases.

BACKGROUND: Next-generation sequencing (NGS) has significantly transformed the landscape of identifying disease-causing genes associated with genetic disorders. However, a substantial portion of sequenced patients remains undiagnosed. This may be attributed not only to the challenges posed by harder-to-detect variants, such as non-coding and structural variations but also to the existence of variants in genes not previously associated with the patient's clinical phenotype. This study introduces EvORanker, an algorithm that integrates unbiased data from 1,028 eukaryotic genomes to link mutated genes to clinical phenotypes. METHODS: EvORanker utilizes clinical data, multi-scale phylogenetic profiling, and other omics data to prioritize disease-associated genes. It was evaluated on solved exomes and simulated genomes, compared with existing methods, and applied to 6260 knockout genes with mouse phenotypes lacking human associations. Additionally, EvORanker was made accessible as a user-friendly web tool. RESULTS: In the analyzed exomic cohort, EvORanker accurately identified the "true" disease gene as the top candidate in 69% of cases and within the top 5 candidates in 95% of cases, consistent with results from the simulated dataset. Notably, EvORanker outperformed existing methods, particularly for poorly annotated genes. In the case of the 6260 knockout genes with mouse phenotypes, EvORanker linked 41% of these genes to observed human disease phenotypes. Furthermore, in two unsolved cases, EvORanker successfully identified DLGAP2 and LPCAT3 as disease candidates for previously uncharacterized genetic syndromes. CONCLUSIONS: We highlight clade-based phylogenetic profiling as a powerful systematic approach for prioritizing potential disease genes. Our study showcases the efficacy of EvORanker in associating poorly annotated genes to disease phenotypes observed in patients. The EvORanker server is freely available at https://ccanavati.shinyapps.io/EvORanker/ .

SHaploseek is a sequencing-only, high-resolution method for comprehensive preimplantation genetic testing.

Recent advances in genomic technologies expand the scope and efficiency of preimplantation genetic testing (PGT). We previously developed Haploseek, a clinically-validated, variant-agnostic comprehensive PGT solution. Haploseek is based on microarray genotyping of the embryo's parents and relatives, combined with low-pass sequencing of the embryos. Here, to increase throughput and versatility, we aimed to develop a sequencing-only implementation of Haploseek. Accordingly, we developed SHaploseek, a universal PGT method to determine genome-wide haplotypes of each embryo based on low-pass (≤ 5x) sequencing of the parents and relative(s) along with ultra-low-pass (0.2-0.4x) sequencing of the embryos. We used SHaploseek to analyze five single lymphoblast cells and 31 embryos. We validated the genome-wide haplotype predictions against either bulk DNA, Haploseek, or, at focal genomic sites, PCR-based PGT results. SHaploseek achieved > 99% concordance with bulk DNA in two families from which single cells were derived from grown-up children. In embryos from 12 PGT families, all of SHaploseek's focal site haplotype predictions were concordant with clinical PCR-based PGT results. Genome-wide, there was > 99% median concordance between Haploseek and SHaploseek's haplotype predictions. Concordance remained high at all assayed sequencing depths ≥ 2x, as well as with only 1ng of parental DNA input. In subtelomeric regions, significantly more haplotype predictions were high-confidence in SHaploseek compared to Haploseek. In summary, SHaploseek constitutes a single-platform, accurate, and cost-effective comprehensive PGT solution.

Loss of function of FIGNL1, a DNA damage response gene, causes human ovarian dysgenesis.

Ovarian dysgenesis (OD), an XX disorder of sex development, presents with primary amenorrhea, hypergonadotrophic hypogonadism, and infertility. In an Ashkenazi Jewish patient with OD, whole exome sequencing identified compound heterozygous frameshifts in FIGNL1, a DNA damage response (DDR) gene: c.189del and c.1519_1523del. Chromosomal breakage was significantly increased in patient cells, both spontaneously, and following mitomycin C exposure. Transfection of DYK-tagged FIGNL1 constructs in HEK293 cells showed no detectable protein in FIGNL1c.189del and truncation with reduced expression in FIGNL1c.1519_1523del (64% of wild-type [WT], P = .003). FIGNL1 forms nuclear foci increased by phleomycin treatment (20.6 ± 1.6 vs 14.8 ± 2.4, P = .02). However, mutant constructs showed reduced DYK-FIGNL1 foci formation in non-treated cells (0.8 ± 0.9 and 5.6 ± 1.5 vs 14.8 ± 2.4 in DYK-FIGNL1WT, P < .001) and no increase with phleomycin treatment. In conclusion, FIGNL1 loss of function is a newly characterized OD gene, highlighting the DDR pathway's role in ovarian development and maintenance and suggesting chromosomal breakage as an assessment tool in XX-DSD patients.

Expanded clinical validation of Haploseek for comprehensive preimplantation genetic testing.

PURPOSE: We previously developed Haploseek, a method for comprehensive preimplantation genetic testing (PGT). However, some key features were missing, and the method has not yet been systematically validated. METHODS: We extended Haploseek to incorporate DNA from embryo grandparents and to allow testing of variants on chromosome X or in regions where parents share common haplotypes. We then validated Haploseek on 151 embryo biopsies from 27 clinical PGT cases. We sequenced all biopsies to low coverage (0.2×), and performed single-nucleotide polymorphism (SNP) microarray genotyping on the embryos' parents and siblings/grandparents. We used the extended Haploseek to predict chromosome copy-number variants (CNVs) and relevant variant-flanking haplotypes in each embryo. We validated haplotype predictions for each clinical sample against polymerase chain reaction (PCR)-based PGT case results, and CNV predictions against established commercial kits. RESULTS: For each of the 151 embryo biopsies, all Haploseek-derived haplotypes and CNVs were concordant with clinical PGT results. The cases included 17 autosomal dominant, 5 autosomal recessive, and 3 X-linked monogenic disorders. In addition, we evaluated 1 Robertsonian and 2 reciprocal translocations, and 17 cases of chromosome copy-number counting were performed. CONCLUSION: Our results demonstrate that Haploseek is clinically accurate and fit for all standard clinical PGT applications.

Helicase-inactivating BRIP1 mutation yields Fanconi anemia with microcephaly and other congenital abnormalities.

Fanconi anemia is a genetically and phenotypically heterogeneous disorder characterized by congenital anomalies, bone marrow failure, cancer, and sensitivity of chromosomes to DNA cross-linking agents. One of the 22 genes responsible for Fanconi anemia is BRIP1, in which biallelic truncating mutations lead to Fanconi anemia group J and monoallelic truncating mutations predispose to certain cancers. However, of the more than 1000 reported missense mutations in BRIP1, very few have been functionally characterized. We evaluated the functional consequence of BRIP1 p.R848H (c.2543G > A), which was homozygous in two cousins with low birth weight, microcephaly, upper limb abnormalities, and imperforate anus and for whom chromosome breakage analysis of patient cells revealed increased mitomycin C sensitivity. BRIP1 p.R848H alters a highly conserved residue in the catalytic DNA helicase domain. We show that BRIP1 p.R848H leads to a defect in helicase activity. Heterozygosity at this missense has been reported in multiple cancer patients but, in the absence of functional studies, classified as of unknown significance. Our results support that this mutation is pathogenic for Fanconi anemia in homozygotes and for increased cancer susceptibility in heterozygous carriers.

Genomic analysis of inherited hearing loss in the Palestinian population.

The genetic characterization of a common phenotype for an entire population reveals both the causes of that phenotype for that place and the power of family-based, population-wide genomic analysis for gene and mutation discovery. We characterized the genetics of hearing loss throughout the Palestinian population, enrolling 2,198 participants from 491 families from all parts of the West Bank and Gaza. In Palestinian families with no prior history of hearing loss, we estimate that 56% of hearing loss is genetic and 44% is not genetic. For the great majority (87%) of families with inherited hearing loss, panel-based genomic DNA sequencing, followed by segregation analysis of large kindreds and transcriptional analysis of participant RNA, enabled identification of the causal genes and mutations, including at distant noncoding sites. Genetic heterogeneity of hearing loss was striking with respect to both genes and alleles: The 337 solved families harbored 143 different mutations in 48 different genes. For one in four solved families, a transcription-altering mutation was the responsible allele. Many of these mutations were cryptic, either exonic alterations of splice enhancers or silencers or deeply intronic events. Experimentally calibrated in silico analysis of transcriptional effects yielded inferences of high confidence for effects on splicing even of mutations in genes not expressed in accessible tissue. Most (58%) of all hearing loss in the population was attributable to consanguinity. Given the ongoing decline in consanguineous marriage, inherited hearing loss will likely be much rarer in the next generation.

The novel founder homozygous V225M mutation in the HSD17B3 gene causes aberrant splicing and XY-DSD.

PURPOSE: Mutations in the gene HSD17B3 encoding the 17-beta hydroxysteroid dehydrogenase 3 enzyme cause testosterone insufficiency leading to XY disorders of sex development. In this study the clinical and molecular characteristics of three patients from consanguineous families are elucidated. METHODS: We identified three patients from two unrelated families with XY DSD and a novel homozygous HSD17B3:c. 673G>A mutation. The effect of the mutation on splicing was determined in RNA extracted from the testis of one patient. RESULTS: Three patients presented at ages 0.1, 8 and 0.7 years with ambiguous genitalia and an XY Karyotype. Endocrine workup showed normal cortisol and mineralocorticoid levels with a low testosterone/androstenedione ratio. Whole-exome sequencing, carried out in the first family, revealed a homozygous novel mutation in the HSD17B3 gene: c. 673G>A, p. V225M. The same mutation was found by Sanger sequencing in the third unrelated patient. Haplotype analysis of a 4 Mb region surrounding the HSD17B3 gene on chromosome 9 revealed that the mutation resides on the same allele in all three patients. The mutation, being the first nucleic acid on exon 10, affects splicing and causes exon 10 skipping in one of our patients' testes. CONCLUSION: The novel homozygous c. 673G>A, p. V225M mutation in the 17HSDB3 gene is likely a founder mutation and causes severe XY-DSD. It changes a conserved amino acid residue, and also alters 17HSDB3 gene transcription by causing skipping of exon 10, thereby contributing to an imbalance in the relevant protein isoforms and consequently, significant decreased 17HDSB3 enzymatic activity.

Nucleotide composition affects codon usage toward the 3'-end.

The 3'-end of the coding sequence in several species is known to show specific codon usage bias. Several factors have been suggested to underlie this phenomenon, including selection against translation efficiency, selection for translation accuracy, and selection against RNA folding. All are supported by some evidence, but there is no general agreement as to which factors are the main determinants. Nor is it known how universal this phenomenon is, and whether the same factors explain it in different species. To answer these questions, we developed a measure that quantifies the codon usage bias at the gene end, and used it to compute this bias for 91 species that span the three domains of life. In addition, we characterized the codons in each species by features that allow discrimination between the different factors. Combining all these data, we were able to show that there is a universal trend to favor AT-rich codons toward the gene end. Moreover, we suggest that this trend is explained by avoidance from forming RNA secondary structures around the stop codon, which may interfere with normal translation termination.

Inclusion of hemimegalencephaly into the phenotypic spectrum of NPRL3 pathogenic variants in familial focal epilepsy with variable foci.

Despite tremendous progress through next generation sequencing technologies, familial focal epilepsies are insufficiently understood. We sought to identify the genetic basis in multiplex Palestinian families with familial focal epilepsy with variable foci (FFEVF). Family I with 10 affected individuals and Family II with five affected individuals underwent detailed phenotyping over three generations. The phenotypic spectrum of the two families varied from nonlesional focal epilepsy including nocturnal frontal lobe epilepsy to severe structural epilepsy due to hemimegalencephaly. Whole-exome sequencing and single nucleotide polymorphism array analysis revealed pathogenic variants in NPRL3 in each family, a partial ~38-kb deletion encompassing eight exons (exons 8-15) and the 3'-untranslated region of the NPRL3 gene in Family I, and a de novo nonsense variant c.1063C>T, p.Gln355* in Family II. Furthermore, we identified a truncating variant in the PDCD10 gene in addition to the NPRL3 variant in a patient with focal epilepsy from Family I. The individual also had developmental delay and multiple cerebral cavernomas, possibly demonstrating a digenic contribution to the individual's phenotype. Our results implicate the association of NPRL3 with hemimegalencephaly, expanding the phenotypic spectrum of NPRL3 in FFEVF and underlining that partial deletions are part of the genotypic spectrum of NPRL3 variants.

Off the street phasing (OTSP): no hassle haplotype phasing for molecular PGD applications.

PURPOSE: Pre-implantation genetic diagnosis (PGD) for molecular disorders requires the construction of parental haplotypes. Classically, haplotype resolution ("phasing") is obtained by genotyping multiple polymorphic markers in both parents and at least one additional relative. However, this process is time-consuming, and immediate family members are not always available. The recent availability of massive genomic data for many populations promises to eliminate the needs for developing family-specific assays and for recruiting additional family members. In this study, we aimed to validate population-assisted haplotype phasing for PGD. METHODS: Targeted sequencing of CFTR gene variants and ~ 1700 flanking polymorphic SNPs (± 2 Mb) was performed on 54 individuals from 12 PGD families of (a) Full Ashkenazi (FA; n = 16), (b) mixed Ashkenazi (MA; n = 23 individuals with at least one Ashkenazi and one non-Ashkenazi grandparents), or (c) non-Ashkenazi (NA; n = 15) descent. Heterozygous genotype calls in each individual were phased using various whole genome reference panels and appropriate computational models. All computationally derived haplotype predictions were benchmarked against trio-based phasing. RESULTS: Using the Ashkenazi reference panel, phasing of FA was highly accurate (99.4% ± 0.2% accuracy); phasing of MA was less accurate (95.4% ± 4.5% accuracy); and phasing of NA was predictably low (83.4% ± 6.6% accuracy). Strikingly, for founder mutation carriers, our haplotyping approach facilitated near perfect phasing accuracy (99.9% ± 0.1% and 98.2% ± 2.8% accuracy for W1282X and delF508 carriers, respectively). CONCLUSIONS: Our results demonstrate the feasibility of replacing classical haplotype phasing with population-based phasing with uncompromised accuracy.

Haploseek: a 24-hour all-in-one method for preimplantation genetic diagnosis (PGD) of monogenic disease and aneuploidy.

PURPOSE: To develop an economical, user-friendly, and accurate all-in-one next-generation sequencing (NGS)-based workflow for single-cell gene variant detection combined with comprehensive chromosome screening in a 24-hour workflow protocol. METHODS: We subjected single lymphoblast cells or blastomere/blastocyst biopsies from four different families to low coverage (0.3×-1.4×) genome sequencing. We combined copy-number variant (CNV) detection and whole-genome haplotype phase prediction via Haploseek, a novel, user-friendly analysis pipeline. We validated haplotype predictions for each sample by comparing with clinical preimplantation genetic diagnosis (PGD) case results or by single-nucleotide polymorphism (SNP) microarray analysis of bulk DNA from each respective lymphoblast culture donor. CNV predictions were validated by established commercial kits for single-cell CNV prediction. RESULTS: Haplotype phasing of the single lymphoblast/embryo biopsy sequencing data was highly concordant with relevant ground truth haplotypes in all samples/biopsies from all four families. In addition, whole-genome copy-number assessments were concordant with the results of a commercial kit. CONCLUSION: Our results demonstrate the establishment of a reliable method for all-in-one molecular and chromosomal diagnosis of single cells. Important features of the Haploseek pipeline include rapid sample processing, rapid sequencing, streamlined analysis, and user-friendly reporting, so as to expedite clinical PGD implementation.

Noninvasive paternal exclusion testing for cystic fibrosis in the first five to eight weeks of gestation.

Prenatal genetic testing is not generally applicable to the very early stages of pregnancy (prior to week 8 gestation), a time period that is crucial to pregnant couples with high risk for transmission of genetic disease to their fetus. Therefore, we developed a new ultra-sensitive targeted next generation sequencing method for noninvasive haplotype-based paternal allele exclusion testing of the cystic fibrosis-associated gene, CFTR. This new method was compared to a conventional library prep and sequencing analysis method and all test results were validated by amniotic fluid testing at later stages of pregnancy. Out of 7 enrolled couples, who provided at least two blood samples (at least one week apart) for noninvasive CFTR testing, a result was obtained for 6 fetuses. Using the new hypersensitive method, all six couples (100%) received a correct diagnosis for the paternal allele as opposed to 3/6 (50%) when tested with the conventional strategy. Among 4 couples who provided just one early pregnancy blood draw for analysis, diagnosis was possible in one fetus, but only using the ultra-sensitive method. Thus, we describe a novel noninvasive CFTR screening method which demonstrates unprecedented fetal allele typing accuracy in the earliest stages of pregnancy.

Mitochondrial PITRM1 peptidase loss-of-function in childhood cerebellar atrophy.

OBJECTIVE: To identify the genetic basis of a childhood-onset syndrome of variable severity characterised by progressive spinocerebellar ataxia, mental retardation, psychotic episodes and cerebellar atrophy. METHODS: Identification of the underlying mutations by whole exome and whole genome sequencing. Consequences were examined in patients' cells and in yeast. RESULTS: Two brothers from a consanguineous Palestinian family presented with progressive spinocerebellar ataxia, mental retardation and psychotic episodes. Serial brain imaging showed severe progressive cerebellar atrophy. Whole exome sequencing revealed a novel mutation: pitrilysin metallopeptidase 1 (PITRM1) c.2795C>T, p.T931M, homozygous in the affected children and resulting in 95% reduction in PITRM1 protein. Whole genome sequencing revealed a chromosome X structural rearrangement that also segregated with the disease. Independently, two siblings from a second Palestinian family presented with similar, somewhat milder symptoms and the same PITRM1 mutation on a shared haplotype. PITRM1T931M carrier frequency was 0.027 (3/110) in the village of the first family evaluated, and 0/300 among Palestinians from other locales. PITRM1 is a mitochondrial matrix enzyme that degrades 10-65 amino acid oligopeptides, including the mitochondrial fraction of amyloid-beta peptide. Analysis of peptide cleavage activity by the PITRM1T931M protein revealed a significant decrease in the degradation capacity specifically of peptides ≥40 amino acids. CONCLUSION: PITRM1T931M results in childhood-onset recessive cerebellar pathology. Severity of PITRM1-related disease may be affected by the degree of impairment in cleavage of mitochondrial long peptides. Disruption and deletion of X linked regulatory segments may also contribute to severity.

The role of nucleotide composition in premature termination codon recognition.

BACKGROUND: It is not fully understood how a termination codon is recognized as premature (PTC) by the nonsense-mediated decay (NMD) machinery. This is particularly true for transcripts lacking an exon junction complex (EJC) along their 3' untranslated region (3'UTR), and thus degrade through the EJC-independent NMD pathway. RESULTS: Here, we analyzed data of transcript stability change following NMD repression and identified over 200 EJC-independent NMD-targets. We examined many features characterizing these transcripts, and compared them to NMD-insensitive transcripts, as well as to a group of transcripts that are destabilized following NMD repression (destabilized transcripts). CONCLUSIONS: We found that none of the known NMD-triggering features, such as the presence of upstream open reading frames, significantly characterizes EJC-independent NMD-targets. Instead, we saw that NMD-targets are strongly enriched with G nucleotides upstream of the termination codon, and even more so along their 3'UTR. We suggest that high G content around the termination codon impedes translation termination as a result of mRNA folding, thus triggering NMD. We also suggest that high G content in the 3'UTR helps to activate NMD by allowing for the accumulation of UPF1, or other NMD-promoting proteins, along the 3'UTR.

Proof-of-principle rapid noninvasive prenatal diagnosis of autosomal recessive founder mutations.

BACKGROUND: Noninvasive prenatal testing can be used to accurately detect chromosomal aneuploidies in circulating fetal DNA; however, the necessity of parental haplotype construction is a primary drawback to noninvasive prenatal diagnosis (NIPD) of monogenic disease. Family-specific haplotype assembly is essential for accurate diagnosis of minuscule amounts of circulating cell-free fetal DNA; however, current haplotyping techniques are too time-consuming and laborious to be carried out within the limited time constraints of prenatal testing, hampering practical application of NIPD in the clinic. Here, we have addressed this pitfall and devised a universal strategy for rapid NIPD of a prevalent mutation in the Ashkenazi Jewish (AJ) population. METHODS: Pregnant AJ couples, carrying mutation(s) in GBA, which encodes acid ß-glucosidase, were recruited at the SZMC Gaucher Clinic. Targeted next-generation sequencing of GBA-flanking SNPs was performed on peripheral blood samples from each couple, relevant mutation carrier family members, and unrelated individuals who are homozygotes for an AJ founder mutation. Allele-specific haplotypes were constructed based on linkage, and a consensus Gaucher disease-associated founder mutation-flanking haplotype was fine mapped. Together, these haplotypes were used for NIPD. All test results were validated by conventional prenatal or postnatal diagnostic methods. RESULTS: Ten parental alleles in eight unrelated fetuses were diagnosed successfully based on the noninvasive method developed in this study. The consensus mutation-flanking haplotype aided diagnosis for 6 of 9 founder mutation alleles. CONCLUSIONS: The founder NIPD method developed and described here is rapid, economical, and readily adaptable for prenatal testing of prevalent autosomal recessive disease-causing mutations in an assortment of worldwide populations. FUNDING: SZMC, Protalix Biotherapeutics Inc., and Centogene AG.

The unfolded protein response affects readthrough of premature termination codons.

One-third of monogenic inherited diseases result from premature termination codons (PTCs). Readthrough of in-frame PTCs enables synthesis of full-length functional proteins. However, extended variability in the response to readthrough treatment is found among patients, which correlates with the level of nonsense transcripts. Here, we aimed to reveal cellular pathways affecting this inter-patient variability. We show that activation of the unfolded protein response (UPR) governs the response to readthrough treatment by regulating the levels of transcripts carrying PTCs. Quantitative proteomic analyses showed substantial differences in UPR activation between patients carrying PTCs, correlating with their response. We further found a significant inverse correlation between the UPR and nonsense-mediated mRNA decay (NMD), suggesting a feedback loop between these homeostatic pathways. We uncovered and characterized the mechanism underlying this NMD-UPR feedback loop, which augments both UPR activation and NMD attenuation. Importantly, this feedback loop enhances the response to readthrough treatment, highlighting its clinical importance. Altogether, our study demonstrates the importance of the UPR and its regulatory network for genetic diseases caused by PTCs and for cell homeostasis under normal conditions.