Evaluation of relative quantification of alternatively spliced transcripts using droplet digital PCR

Introduction: For the relative quantification of isoform expression, RT-qPCR has been the gold standard for over a decade. More recently, digital PCR is becoming widely implemented, as it is promised to be more accurate, sensitive and less affected by inhibitors, without the need for standard curves. In this study we evaluated RT-qPCR versus RT-droplet digital PCR (ddPCR) for the relative quantification of isoforms in controls and carriers of the splice site mutation BRCA1 c.212+3A>G, associated with increased expression of several isoforms. Materials and methods: RNA was extracted from EBV cell lines of controls and heterozygous BRCA1 c.212+3A>G carriers. Transcript-specific plasmids were available to determine the efficiency, precision, reproducibility and accuracy of each method. Results: Both ddPCR and RT-qPCR were able to accurately quantify all targets and showed the same LOB, LOD and LOQ; also precision and reproducibility were similar. Both techniques have the same dynamic range and linearity at biologically relevant template concentrations. However, a significantly higher cost and workload was required for ddPCR experiments. Conclusions: Our study recognizes the potential and validity of digital PCR but shows the value of a highly optimized qPCR for the relative quantification of isoforms. Cost efficiency and simplicity turned out to be better for RT-qPCR. Keywords: Reverse transcriptase polymerase chain reaction, Alternative splicing, Droplet digital PC

Long-read sequencing to unravel complex structural variants of CEP78 leading to cone-rod dystrophy and hearing loss

Inactivating variants as well as a missense variant in the centrosomal CEP78 gene have been identified in autosomal recessive cone-rod dystrophy with hearing loss (CRDHL), a rare syndromic inherited retinal disease distinct from Usher syndrome. Apart from this, a complex structural variant (SV) implicating CEP78 has been reported in CRDHL. Here we aimed to expand the genetic architecture of typical CRDHL by the identification of complex SVs of the CEP78 region and characterization of their underlying mechanisms. Approaches used for the identification of the SVs are shallow whole-genome sequencing (sWGS) combined with quantitative polymerase chain reaction (PCR) and long-range PCR, or ExomeDepth analysis on whole-exome sequencing (WES) data. Targeted or whole-genome nanopore long-read sequencing (LRS) was used to delineate breakpoint junctions at the nucleotide level. For all SVs cases, the effect of the SVs on CEP78 expression was assessed using quantitative PCR on patient-derived RNA. Apart from two novel canonical CEP78 splice variants and a frameshifting single-nucleotide variant (SNV), two SVs affecting CEP78 were identified in three unrelated individuals with CRDHL: a heterozygous total gene deletion of 235 kb and a partial gene deletion of 15 kb in a heterozygous and homozygous state, respectively. Assessment of the molecular consequences of the SVs on patient’s materials displayed a loss-of-function effect. Delineation and characterization of the 15-kb deletion using targeted LRS revealed the previously described complex CEP78 SV, suggestive of a recurrent genomic rearrangement. A founder haplotype was demonstrated for the latter SV in cases of Belgian and British origin, respectively. The novel 235-kb deletion was delineated using whole-genome LRS. Breakpoint analysis showed microhomology and pointed to a replication-based underlying mechanism. Moreover, data mining of bulk and single-cell human and mouse transcriptional datasets, together with CEP78 immunostaining on human retina, linked the CEP78 expression domain with its phenotypic manifestations. Overall, this study supports that the CEP78 locus is prone to distinct SVs and that SV analysis should be considered in a genetic workup of CRDHL. Finally, it demonstrated the power of sWGS and both targeted and whole-genome LRS in identifying and characterizing complex SVs in patients with ocular diseases

Mutations in SAMD7 cause autosomal-recessive macular dystrophy with or without cone dysfunction

Sterile alpha motif domain containing 7 (SAMD7) is a component of the Polycomb repressive complex 1, which inhibits transcription of many genes, including those activated by the transcription factor Cone-Rod Homeobox (CRX). Here we report bi-allelic mutations in SAMD7 as a cause of autosomal-recessive macular dystrophy with or without cone dysfunction. Four of these mutations affect splicing, while another mutation is a missense variant that alters the repressive effect of SAMD7 on CRX-dependent promoter activity, as shown by in vitro assays. Immunostaining of human retinal sections revealed that SAMD7 is localized in the nuclei of both rods and cones, as well as in those of cells belonging to the inner nuclear layer. These results place SAMD7 as a gene crucial for human retinal function and demonstrate a significant difference in the role of SAMD7 between the human and the mouse retina

Mutational spectrum in a worldwide study of 29,700 families with BRCA1 or BRCA2 mutations.

The prevalence and spectrum of germline mutations in BRCA1 and BRCA2 have been reported in single populations, with the majority of reports focused on White in Europe and North America. The Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) has assembled data on 18,435 families with BRCA1 mutations and 11,351 families with BRCA2 mutations ascertained from 69 centers in 49 countries on six continents. This study comprehensively describes the characteristics of the 1,650 unique BRCA1 and 1,731 unique BRCA2 deleterious (disease-associated) mutations identified in the CIMBA database. We observed substantial variation in mutation type and frequency by geographical region and race/ethnicity. In addition to known founder mutations, mutations of relatively high frequency were identified in specific racial/ethnic or geographic groups that may reflect founder mutations and which could be used in targeted (panel) first pass genotyping for specific populations. Knowledge of the population-specific mutational spectrum in BRCA1 and BRCA2 could inform efficient strategies for genetic testing and may justify a more broad-based oncogenetic testing in some populations

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer.

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC

Juggling offsets unlocks RNA-seq tools for fast scalable differential usage, aberrant splicing and expression analyses

AbstractRNA-sequencing (RNA-seq) is increasingly used to diagnose patients with rare diseases by prioritising genes with aberrant expression and/or splicing. State-of-the-art methods for detecting aberrant expression and splicing, however, are extremely slow. The latter, also discard much information because they only use junction reads to infer aberrant splicing. In this contribution, we show that replacing the offset for library size unlocks conventional bulk RNA-seq workflows for fast and scalable differential usage, aberrant splicing and expression analyses. Our method, saseR, is several orders of magnitude faster than the state-of-the-art methods and dramatically outperforms these in terms of sensitivity and specificity for aberrant splicing, while being on par with these inferring differential usage and aberrant expression. Finally, our framework is also very flexible and can be used for all applications that involve the analysis of proportions of short- or long RNA-seq read counts.</jats:p