A unified analytic framework for prioritization of non-coding variants of uncertain significance in heritable breast and ovarian cancer

ATM Variant as a Cause of Hereditary Cutaneous Melanoma in a Spanish Family: Case Report

Introduction

Ataxia-Telangiectasia Mutated (ATM) is a cancer predisposition gene; carriers of germline pathogenic variants have an increased risk of developing malignancies, including breast, prostate, pancreatic, and ovarian cancer. Most ATM variants are of uncertain significance. Findings from genome-wide association studies (GWAS) suggest that ATM may be a low-risk melanoma susceptibility locus.

Case Report

We report the case of a Hispanic family whose members who have presented cutaneous melanoma have been found to be carriers for the ATM pathogenic variant c.3747-1G>C (rs730881364), one of whom was diagnosed at 24 years old.

Discussion

We describe for the first time the possible clinical association between ATM (c.3747-1G>C) and familial melanoma. In silico splice site analysis predicts that this alteration will weaken the native splice acceptor site and will result in the creation or strengthening of a novel splice acceptor site, assuming a variant that entails loss of functionality that is probably pathogenic and related to oncogenesis. However, we cannot exclude that cutaneous melanoma in both members and at an early age is the result of chance, environmental interaction, other uncontrolled external factors, or the interaction of other genetic alterations other than the ATM variant described in this study.

Deep learning-assisted genome-wide characterization of massively parallel reporter assays

Massively parallel reporter assay (MPRA) is a high-throughput method that enables the study of the regulatory activities of tens of thousands of DNA oligonucleotides in a single experiment. While MPRA experiments have grown in popularity, their small sample sizes compared to the scale of the human genome limits our understanding of the regulatory effects they detect. To address this, we develop a deep learning model, MpraNet, to distinguish potential MPRA targets from the background genome. This model achieves high discriminative performance (AUROC = 0.85) at differentiating MPRA positives from a set of control variants that mimic the background genome when applied to the lymphoblastoid cell line. We observe that existing functional scores represent very distinct functional effects, and most of them fail to characterize the regulatory effect that MPRA detects. Using MpraNet, we predict potential MPRA functional variants across the genome and identify the distributions of MPRA effect relative to other characteristics of genetic variation, including allele frequency, alternative functional annotations specified by FAVOR, and phenome-wide associations. We also observed that the predicted MPRA positives are not uniformly distributed across the genome; instead, they are clumped together in active regions comprising 9.95% of the genome and inactive regions comprising 89.07% of the genome. Furthermore, we propose our model as a screen to filter MPRA experiment candidates at genome-wide scale, enabling future experiments to be more cost-efficient by increasing precision relative to that observed from previous MPRAs.

PALB2 germline mutations in a multi-gene panel testing cohort of 1905 breast-ovarian cancer patients in Argentina

PURPOSE

PALB2 variants have been scarcely described in Argentinian and Latin-American reports. In this study, we describe molecular and clinical characteristics of PALB2 mutations found in multi-gene panels (MP) from breast-ovarian cancer (BOC) families in different institutions from Argentina.

METHODS

We retrospectively identified PALB2 pathogenic (PV) and likely pathogenic (LPV) variants from a cohort of 1905 MP results, provided by one local lab (Heritas) and SITHER (Hereditary Tumor Information System) public database. All patients met hereditary BOC clinical criteria for testing, according to current guidelines.

RESULTS

The frequency of PALB2 mutations is 2.78% (53/1905). Forty-eight (90.5%) are PV and five (9.5%) are LPV. Most of the 18 different mutations (89%) are nonsense and frameshift types and 2 variants are novel. One high-rate recurrent PV (Y551*) is present in 43% (23/53) of the unrelated index cases. From the 53 affected carriers, 94% have BC diagnosis with 14% of bilateral cases. BC phenotype is mainly invasive ductal (78%) with 62% of hormone-receptor positive and 22% of triple negative tumors. Self-reported ethnic background of the cohort is West European (66%) and native Latin-American (20%) which is representative of Buenos Aires and other big urban areas of the country.

CONCLUSION

This is the first report describing molecular and clinical characteristics of PALB2 carriers in Argentina. Frequency of PALB2 PV in Argentinian HBOC families is higher than in other reported populations. Y551* is a recurrent mutation that seems to be responsible for almost 50% of PALB2 cases.

Next step in molecular genetics of hereditary breast/ovarian cancer: Multigene panel testing in clinical actionably genes and prioritization algorithms in the study of variants of uncertain significance

INTRODUCTION

BRCA1 and BRCA2 are the two main genes causing hereditary breast and ovarian cancer (HBOC). However, thanks to the development of Next Generation Sequencing (NGS), other genes linked to this syndrome (CHEK2, BRIP1, ATM and PALB2 among others) can be analysed.

MATERIAL AND METHODS

an analysis by multigene panel testing was performed in 138 index cases (ICs) from HBOC Spanish families with a previous non-informative result for BRCA1/2. The BRCA Hereditary Cancer Master™ Plus kit, including 26 actionable and candidate genes related to HBOC was employed. Once classified, an algorithm was employed to prioritized those variants of unknown significance with a higher risk of having a deleterious effect. Moreover, a mRNA splicing assay was performed for the prioritized VUS c.3402+3A > C in ATM, located at intron 23.

RESULTS

A total of 82 variants were found: 70 VUS and 12 pathogenic or probably pathogenic variants. The diagnostic yield in actionable genes non-BRCA was 7.97% of the total tested ICs. Overall, 19 VUS were prioritized, which meant 27% of the 70 total VUS. RNA analysis of the variant 3402+3A > C confirmed a deleterious impact on splicing.

DISCUSSION

The implementation of a multigene panel in HBOC studied families improved the diagnostic yield, concordant with results obtained in previous publications. Due to the important number of VUS obtained in NGS, the application of a prioritization algorithm is needed in order to select those variants in which it is necessary to conduct further studies.

CDH1 Germline Variants in a Tunisian Cohort with Hereditary Diffuse Gastric Carcinoma

Mutational screening of the CDH1 gene is a standard treatment for patients who fulfill Hereditary Diffuse Gastric Cancer (HDGC) testing criteria. In this framework, the classification of variants found in this gene is a crucial step for the clinical management of patients at high risk for HDGC. The aim of our study was to identify CDH1 as well as CTNNA1 mutational profiles predisposing to HDGC in Tunisia. Thirty-four cases were included for this purpose. We performed Sanger sequencing for the entire coding region of both genes and MLPA (Multiplex Ligation Probe Amplification) assays to investigate large rearrangements of the CDH1 gene. As a result, three cases, all with the HDGC inclusion criteria (8.82% of the entire cohort), carried pathogenic and likely pathogenic variants of the CDH1 gene. These variants involve a novel splicing alteration, a missense c.2281G > A detected by Sanger sequencing, and a large rearrangement detected by MLPA. No pathogenic CTNNA1 variants were found. The large rearrangement is clearly pathogenic, implicating a large deletion of two exons. The novel splicing variant creates a cryptic site. The missense variant is a VUS (Variant with Uncertain Significance). With ACMG (American College of Medical Genetics and Genomics) classification and the evidence available, we thus suggest a revision of its status to likely pathogenic. Further functional studies or cosegregation analysis should be performed to confirm its pathogenicity. In addition, molecular exploration will be needed to understand the etiology of the other CDH1- and CTNNA1-negative cases fulfilling the HDGC inclusion criteria.

Unique somatic variants in DNA from urine exosomes of individuals with bladder cancer

Bladder cancer (BC), a heterogeneous disease characterized by high recurrence rates, is diagnosed and monitored by cystoscopy. Accurate clinical staging based on biopsy remains a challenge, and additional, objective diagnostic tools are needed urgently. We used exosomal DNA (exoDNA) as an analyte to examine cancer-associated mutations and compared the diagnostic utility of exoDNA from urine and serum of individuals with BC. In contrast to urine exosomes from healthy individuals, urine exosomes from individuals with BC contained significant amounts of DNA. Whole-exome sequencing of DNA from matched urine and serum exosomes, bladder tumors, and normal tissue (peripheral blood mononuclear cells) identified exonic and 3' UTR variants in frequently mutated genes in BC, detectable in urine exoDNA and matched tumor samples. Further analyses identified somatic variants in driver genes, unique to urine exoDNA, possibly because of the inherent intra-tumoral heterogeneity of BC, which is not fully represented in random small biopsies. Multiple variants were also found in untranslated portions of the genome, such as microRNA (miRNA)-binding regions of the KRAS gene. Gene network analyses revealed that exoDNA is associated with cancer, inflammation, and immunity in BC exosomes. Our findings show utility of exoDNA as an objective, non-invasive strategy to identify novel biomarkers and targets for BC.

Whole Exome Sequencing Identifies APCDD1 and HDAC5 Genes as Potentially Cancer Predisposing in Familial Colorectal Cancer

Germline mutations in predisposition genes account for only 20% of all familial colorectal cancers (CRC) and the remaining genetic burden may be due to rare high- to moderate-penetrance germline variants that are not explored. With the aim of identifying such potential cancer-predisposing variants, we performed whole exome sequencing on three CRC cases and three unaffected members of a Polish family and identified two novel heterozygous variants: a coding variant in APC downregulated 1 gene (APCDD1, p.R299H) and a non-coding variant in the 5′ untranslated region (UTR) of histone deacetylase 5 gene (HDAC5). Sanger sequencing confirmed the variants segregating with the disease and Taqman assays revealed 8 additional APCDD1 variants in a cohort of 1705 familial CRC patients and no further HDAC5 variants. Proliferation assays indicated an insignificant proliferative impact for the APCDD1 variant. Luciferase reporter assays using the HDAC5 variant resulted in an enhanced promoter activity. Targeting of transcription factor binding sites of SNAI-2 and TCF4 interrupted by the HDAC5 variant showed a significant impact of TCF4 on promoter activity of mutated HDAC5. Our findings contribute not only to the identification of unrecognized genetic causes of familial CRC but also underline the importance of 5’UTR variants affecting transcriptional regulation and the pathogenesis of complex disorders.

A proposed molecular mechanism for pathogenesis of severe RNA-viral pulmonary infections

Background: Certain riboviruses can cause severe pulmonary complications leading to death in some infected patients. We propose that DNA damage induced-apoptosis accelerates viral release, triggered by depletion of host RNA binding proteins (RBPs) from nuclear RNA bound to replicating viral sequences. Methods: Information theory-based analysis of interactions between RBPs and individual sequences in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2), Influenza A (H3N1), HIV-1, and Dengue genomes identifies strong RBP binding sites in these viral genomes. Replication and expression of viral sequences is expected to increasingly sequester RBPs - SRSF1 and RNPS1. Ordinarily, RBPs bound to nascent host transcripts prevents their annealing to complementary DNA. Their depletion induces destabilizing R-loops. Chromosomal breakage occurs when an excess of unresolved R-loops collide with incoming replication forks, overwhelming the DNA repair machinery. We estimated stoichiometry of inhibition of RBPs in host nuclear RNA by counting competing binding sites in replicating viral genomes and host RNA. Results: Host RBP binding sites are frequent and conserved among different strains of RNA viral genomes. Similar binding motifs of SRSF1 and RNPS1 explain why DNA damage resulting from SRSF1 depletion is complemented by expression of RNPS1. Clustering of strong RBP binding sites coincides with the distribution of RNA-DNA hybridization sites across the genome. SARS-CoV-2 replication is estimated to require 32.5-41.8 hours to effectively compete for binding of an equal proportion of SRSF1 binding sites in host encoded nuclear RNAs. Significant changes in expression of transcripts encoding DNA repair and apoptotic proteins were found in an analysis of influenza A and Dengue-infected cells in some individuals. Conclusions: R-loop-induced apoptosis indirectly resulting from viral replication could release significant quantities of membrane-associated virions into neighboring alveoli. These could infect adjacent pneumocytes and other tissues, rapidly compromising lung function, causing multiorgan system failure and other described symptoms.

A proposed molecular mechanism for pathogenesis of severe RNA-viral pulmonary infections

Background: Certain riboviruses can cause severe pulmonary complications leading to death in some infected patients. We propose that DNA damage induced-apoptosis accelerates viral release, triggered by depletion of host RNA binding proteins (RBPs) from nuclear RNA bound to replicating viral sequences. Methods: Information theory-based analysis of interactions between RBPs and individual sequences in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2), Influenza A (H3N2), HIV-1, and Dengue genomes identifies strong RBP binding sites in these viral genomes. Replication and expression of viral sequences is expected to increasingly sequester RBPs - SRSF1 and RNPS1. Ordinarily, RBPs bound to nascent host transcripts prevents their annealing to complementary DNA. Their depletion induces destabilizing R-loops. Chromosomal breakage occurs when an excess of unresolved R-loops collide with incoming replication forks, overwhelming the DNA repair machinery. We estimated stoichiometry of inhibition of RBPs in host nuclear RNA by counting competing binding sites in replicating viral genomes and host RNA. Results: Host RBP binding sites are frequent and conserved among different strains of RNA viral genomes. Similar binding motifs of SRSF1 and RNPS1 explain why DNA damage resulting from SRSF1 depletion is complemented by expression of RNPS1. Clustering of strong RBP binding sites coincides with the distribution of RNA-DNA hybridization sites across the genome. SARS-CoV-2 replication is estimated to require 32.5-41.8 hours to effectively compete for binding of an equal proportion of SRSF1 binding sites in host encoded nuclear RNAs. Significant changes in expression of transcripts encoding DNA repair and apoptotic proteins were found in an analysis of influenza A and Dengue-infected cells in some individuals. Conclusions: R-loop-induced apoptosis indirectly resulting from viral replication could release significant quantities of membrane-associated virions into neighboring alveoli. These could infect adjacent pneumocytes and other tissues, rapidly compromising lung function, causing multiorgan system failure and other described symptoms.

Expression Changes Confirm Genomic Variants Predicted to Result in Allele-Specific, Alternative mRNA Splicing

Splice isoform structure and abundance can be affected by either noncoding or masquerading coding variants that alter the structure or abundance of transcripts. When these variants are common in the population, these nonconstitutive transcripts are sufficiently frequent so as to resemble naturally occurring, alternative mRNA splicing. Prediction of the effects of such variants has been shown to be accurate using information theory-based methods. Single nucleotide polymorphisms (SNPs) predicted to significantly alter natural and/or cryptic splice site strength were shown to affect gene expression. Splicing changes for known SNP genotypes were confirmed in HapMap lymphoblastoid cell lines with gene expression microarrays and custom designed q-RT-PCR or TaqMan assays. The majority of these SNPs (15 of 22) as well as an independent set of 24 variants were then subjected to RNAseq analysis using the ValidSpliceMut web beacon (http://validsplicemut.cytognomix.com), which is based on data from the Cancer Genome Atlas and International Cancer Genome Consortium. SNPs from different genes analyzed with gene expression microarray and q-RT-PCR exhibited significant changes in affected splice site use. Thirteen SNPs directly affected exon inclusion and 10 altered cryptic site use. Homozygous SNP genotypes resulting in stronger splice sites exhibited higher levels of processed mRNA than alleles associated with weaker sites. Four SNPs exhibited variable expression among individuals with the same genotypes, masking statistically significant expression differences between alleles. Genome-wide information theory and expression analyses (RNAseq) in tumor exomes and genomes confirmed splicing effects for 7 of the HapMap SNP and 14 SNPs identified from tumor genomes. q-RT-PCR resolved rare splice isoforms with read abundance too low for statistical significance in ValidSpliceMut. Nevertheless, the web-beacon provides evidence of unanticipated splicing outcomes, for example, intron retention due to compromised recognition of constitutive splice sites. Thus, ValidSpliceMut and q-RT-PCR represent complementary resources for identification of allele-specific, alternative splicing.

Pan-cancer repository of validated natural and cryptic mRNA splicing mutations

We present a major public resource of mRNA splicing mutations validated according to multiple lines of evidence of abnormal gene expression. Likely mutations present in all tumor types reported in the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) were identified based on the comparative strengths of splice sites in tumor versus normal genomes, and then validated by respectively comparing counts of splice junction spanning and abundance of transcript reads in RNA-Seq data from matched tissues and tumors lacking these mutations. The comprehensive resource features 341,486 of these validated mutations, the majority of which (69.9%) are not present in the Single Nucleotide Polymorphism Database (dbSNP 150). There are 131,347 unique mutations which weaken or abolish natural splice sites, and 222,071 mutations which strengthen cryptic splice sites (11,932 affect both simultaneously). 28,812 novel or rare flagged variants (with <1% population frequency in dbSNP) were observed in multiple tumor tissue types. An algorithm was developed to classify variants into splicing molecular phenotypes that integrates germline heterozygosity, degree of information change and impact on expression. The classification thresholds were calibrated against the ClinVar clinical database phenotypic assignments. Variants are partitioned into allele-specific alternative splicing, likely aberrant and aberrant splicing phenotypes. Single variants or chromosome ranges can be queried using a Global Alliance for Genomics and Health (GA4GH)-compliant, web-based Beacon "Validated Splicing Mutations" either separately or in aggregate alongside other Beacons through the public Beacon Network, as well as through our website. The website provides additional information, such as a visual representation of supporting RNAseq results, gene expression in the corresponding normal tissues, and splicing molecular phenotypes.

Screening of BRCA1/2 deep intronic regions by targeted gene sequencing identifies the first germline BRCA1 variant causing pseudoexon activation in a patient with breast/ovarian cancer

BACKGROUND

Genetic analysis of BRCA1 and BRCA2 for the diagnosis of hereditary breast and ovarian cancer (HBOC) is commonly restricted to coding regions and exon-intron boundaries. Although germline pathogenic variants in these regions explain about ~20% of HBOC cases, there is still an important fraction that remains undiagnosed. We have screened BRCA1/2 deep intronic regions to identify potential spliceogenic variants that could explain part of the missing HBOC susceptibility.

METHODS

We analysed BRCA1/2 deep intronic regions by targeted gene sequencing in 192 high-risk HBOC families testing negative for BRCA1/2 during conventional analysis. Rare variants (MAF <0.005) predicted to create/activate splice sites were selected for further characterisation in patient RNA. The splicing outcome was analysed by RT-PCR and Sanger sequencing, and allelic imbalance was also determined when heterozygous exonic loci were present.

RESULTS

A novel transcript was detected in BRCA1 c.4185+4105C>T variant carrier. This variant promotes the inclusion of a pseudoexon in mature mRNA, generating an aberrant transcript predicted to encode for a non-functional protein. Quantitative and allele-specific assays determined haploinsufficiency in the variant carrier, supporting a pathogenic effect for this variant. Genotyping of 1030 HBOC cases and 327 controls did not identify additional carriers in Spanish population.

CONCLUSION

Screening of BRCA1/2 intronic regions has identified the first BRCA1 deep intronic variant associated with HBOC by pseudoexon activation. Although the frequency of deleterious variants in these regions appears to be low, our study highlights the importance of studying non-coding regions and performing comprehensive RNA assays to complement genetic diagnosis.

Pan-cancer repository of validated natural and cryptic mRNA splicing mutations

We present a major public resource of mRNA splicing mutations validated according to multiple lines of evidence of abnormal gene expression. Likely mutations present in all tumor types reported in the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) were identified based on the comparative strengths of splice sites in tumor versus normal genomes, and then validated by respectively comparing counts of splice junction spanning and abundance of transcript reads in RNA-Seq data from matched tissues and tumors lacking these mutations. The comprehensive resource features 341,486 of these validated mutations, the majority of which (69.9%) are not present in the Single Nucleotide Polymorphism Database (dbSNP 150). There are 131,347 unique mutations which weaken or abolish natural splice sites, and 222,071 mutations which strengthen cryptic splice sites (11,932 affect both simultaneously). 28,812 novel or rare flagged variants (with <1% population frequency in dbSNP) were observed in multiple tumor tissue types. An algorithm was developed to classify variants into splicing molecular phenotypes that integrates germline heterozygosity, degree of information change and impact on expression. The classification thresholds were calibrated against the ClinVar clinical database phenotypic assignments. Variants are partitioned into allele-specific alternative splicing, likely aberrant and aberrant splicing phenotypes. Single variants or chromosome ranges can be queried using a Global Alliance for Genomics and Health (GA4GH)-compliant, web-based Beacon "Validated Splicing Mutations" either separately or in aggregate alongside other Beacons through the public Beacon Network, as well as through our website. The website provides additional information, such as a visual representation of supporting RNAseq results, gene expression in the corresponding normal tissues, and splicing molecular phenotypes.

BRCA1 and BRCA2 5' noncoding region variants identified in breast cancer patients alter promoter activity and protein binding

The widespread use of next generation sequencing for clinical testing is detecting an escalating number of variants in noncoding regions of the genome. The clinical significance of the majority of these variants is currently unknown, which presents a significant clinical challenge. We have screened over 6,000 early-onset and/or familial breast cancer (BC) cases collected by the ENIGMA consortium for sequence variants in the 5' noncoding regions of BC susceptibility genes BRCA1 and BRCA2, and identified 141 rare variants with global minor allele frequency < 0.01, 76 of which have not been reported previously. Bioinformatic analysis identified a set of 21 variants most likely to impact transcriptional regulation, and luciferase reporter assays detected altered promoter activity for four of these variants. Electrophoretic mobility shift assays demonstrated that three of these altered the binding of proteins to the respective BRCA1 or BRCA2 promoter regions, including NFYA binding to BRCA1:c.-287C>T and PAX5 binding to BRCA2:c.-296C>T. Clinical classification of variants affecting promoter activity, using existing prediction models, found no evidence to suggest that these variants confer a high risk of disease. Further studies are required to determine if such variation may be associated with a moderate or low risk of BC.

Non-Coding Variants in BRCA1 and BRCA2 Genes: Potential Impact on Breast and Ovarian Cancer Predisposition

BRCA1 and BRCA2 are major breast cancer susceptibility genes whose pathogenic variants are associated with a significant increase in the risk of breast and ovarian cancers. Current genetic screening is generally limited to BRCA1/2 exons and intron/exon boundaries. Most identified pathogenic variants cause the partial or complete loss of function of the protein. However, it is becoming increasingly clear that variants in these regions only account for a small proportion of cancer risk. The role of variants in non-coding regions beyond splice donor and acceptor sites, including those that have no qualitative effect on the protein, has not been thoroughly investigated. The key transcriptional regulatory elements of BRCA1 and BRCA2 are housed in gene promoters, untranslated regions, introns, and long-range elements. Within these sequences, germline and somatic variants have been described, but the clinical significance of the majority is currently unknown and it remains a significant clinical challenge. This review summarizes the available data on the impact of variants on non-coding regions of BRCA1/2 genes and their role on breast and ovarian cancer predisposition.

Assessment of the functional impact of germline BRCA1/2 variants located in non-coding regions in families with breast and/or ovarian cancer predisposition

PURPOSE

The molecular mechanism of breast and/or ovarian cancer susceptibility remains unclear in the majority of patients. While germline mutations in the regulatory non-coding regions of BRCA1 and BRCA2 genes have been described, screening has generally been limited to coding regions. The aim of this study was to evaluate the contribution of BRCA1/2 non-coding variants.

METHODS

Four BRCA1/2 non-coding regions were screened using high-resolution melting analysis/Sanger sequencing or next-generation sequencing on DNA extracted from index cases with breast and ovarian cancer predisposition (3926 for BRCA1 and 3910 for BRCA2). The impact of a set of variants on BRCA1/2 gene regulation was evaluated by site-directed mutagenesis, transfection, followed by Luciferase gene reporter assay.

RESULTS

We identified a total of 117 variants and tested twelve BRCA1 and 8 BRCA2 variants mapping to promoter and intronic regions. We highlighted two neighboring BRCA1 promoter variants (c.-130del; c.-125C > T) and one BRCA2 promoter variants (c.-296C > T) inhibiting significantly the promoter activity. In the functional assays, a regulating region within the intron 12 was found with the same enhancing impact as within the intron 2. Furthermore, the variants c.81-3980A > G and c.4186-2022C > T suppress the positive effect of the introns 2 and 12, respectively, on the BRCA1 promoter activity. We also found some variants inducing the promoter activities.

CONCLUSION

In this study, we highlighted some variants among many, modulating negatively the promoter activity of BRCA1 or 2 and thus having a potential impact on the risk of developing cancer. This selection makes it possible to conduct future validation studies on a limited number of variants.

Discovery and validation of information theory-based transcription factor and cofactor binding site motifs

Data from ChIP-seq experiments can derive the genome-wide binding specificities of transcription factors (TFs) and other regulatory proteins. We analyzed 765 ENCODE ChIP-seq peak datasets of 207 human TFs with a novel motif discovery pipeline based on recursive, thresholded entropy minimization. This approach, while obviating the need to compensate for skewed nucleotide composition, distinguishes true binding motifs from noise, quantifies the strengths of individual binding sites based on computed affinity and detects adjacent cofactor binding sites that coordinate with the targets of primary, immunoprecipitated TFs. We obtained contiguous and bipartite information theory-based position weight matrices (iPWMs) for 93 sequence-specific TFs, discovered 23 cofactor motifs for 127 TFs and revealed six high-confidence novel motifs. The reliability and accuracy of these iPWMs were determined via four independent validation methods, including the detection of experimentally proven binding sites, explanation of effects of characterized SNPs, comparison with previously published motifs and statistical analyses. We also predict previously unreported TF coregulatory interactions (e.g. TF complexes). These iPWMs constitute a powerful tool for predicting the effects of sequence variants in known binding sites, performing mutation analysis on regulatory SNPs and predicting previously unrecognized binding sites and target genes.

Prevalence and spectrum of germline rare variants in BRCA1/2 and PALB2 among breast cancer cases in Sarawak, Malaysia

PURPOSE

To characterize the spectrum of germline mutations in BRCA1, BRCA2, and PALB2 in population-based unselected breast cancer cases in an Asian population.

METHODS

Germline DNA from 467 breast cancer patients in Sarawak General Hospital, Malaysia, where 93% of the breast cancer patients in Sarawak are treated, was sequenced for the entire coding region of BRCA1; BRCA2; PALB2; Exons 6, 7, and 8 of TP53; and Exons 7 and 8 of PTEN. Pathogenic variants included known pathogenic variants in ClinVar, loss of function variants, and variants that disrupt splice site.

RESULTS

We found 27 pathogenic variants (11 BRCA1, 10 BRCA2, 4 PALB2, and 2 TP53) in 34 patients, which gave a prevalence of germline mutations of 2.8, 3.23, and 0.86% for BRCA1, BRCA2, and PALB2, respectively. Compared to mutation non-carriers, BRCA1 mutation carriers were more likely to have an earlier age at onset, triple-negative subtype, and lower body mass index, whereas BRCA2 mutation carriers were more likely to have a positive family history. Mutation carrier cases had worse survival compared to non-carriers; however, the association was mostly driven by stage and tumor subtype. We also identified 19 variants of unknown significance, and some of them were predicted to alter splicing or transcription factor binding sites.

CONCLUSION

Our data provide insight into the genetics of breast cancer in this understudied group and suggest the need for modifying genetic testing guidelines for this population with a much younger age at diagnosis and more limited resources compared with Caucasian populations.

From General Aberrant Alternative Splicing in Cancers and Its Therapeutic Application to the Discovery of an Oncogenic DMTF1 Isoform

Alternative pre-mRNA splicing is a crucial process that allows the generation of diversified RNA and protein products from a multi-exon gene. In tumor cells, this mechanism can facilitate cancer development and progression through both creating oncogenic isoforms and reducing the expression of normal or controllable protein species. We recently demonstrated that an alternative cyclin D-binding myb-like transcription factor 1 (DMTF1) pre-mRNA splicing isoform, DMTF1β, is increasingly expressed in breast cancer and promotes mammary tumorigenesis in a transgenic mouse model. Aberrant pre-mRNA splicing is a typical event occurring for many cancer-related functional proteins. In this review, we introduce general aberrant pre-mRNA splicing in cancers and discuss its therapeutic application using our recent discovery of the oncogenic DMTF1 isoform as an example. We also summarize new insights in designing novel targeting strategies of cancer therapies based on the understanding of deregulated pre-mRNA splicing mechanisms.