Long-read targeted sequencing uncovers clinicopathological associations for C9orf72-linked diseases.

To examine the length of a hexanucleotide expansion in C9orf72, which represents the most frequent genetic cause of frontotemporal lobar degeneration and motor neuron disease, we employed a targeted amplification-free long-read sequencing technology: No-Amp sequencing. In our cross-sectional study, we assessed cerebellar tissue from 28 well-characterized C9orf72 expansion carriers. We obtained 3507 on-target circular consensus sequencing reads, of which 814 bridged the C9orf72 repeat expansion (23%). Importantly, we observed a significant correlation between expansion sizes obtained using No-Amp sequencing and Southern blotting (P = 5.0 × 10-4). Interestingly, we also detected a significant survival advantage for individuals with smaller expansions (P = 0.004). Additionally, we uncovered that smaller expansions were significantly associated with higher levels of C9orf72 transcripts containing intron 1b (P = 0.003), poly(GP) proteins (P = 1.3 × 10- 5), and poly(GA) proteins (P = 0.005). Thorough examination of the composition of the expansion revealed that its GC content was extremely high (median: 100%) and that it was mainly composed of GGGGCC repeats (median: 96%), suggesting that expanded C9orf72 repeats are quite pure. Taken together, our findings demonstrate that No-Amp sequencing is a powerful tool that enables the discovery of relevant clinicopathological associations, highlighting the important role played by the cerebellar size of the expanded repeat in C9orf72-linked diseases.

Repetitive element transcripts are elevated in the brain of C9orf72 ALS/FTLD patients.

Significant transcriptome alterations are detected in the brain of patients with amyotrophic lateral sclerosis (ALS), including carriers of the C9orf72 repeat expansion and C9orf72-negative sporadic cases. Recently, the expression of repetitive element transcripts has been associated with toxicity and, while increased repetitive element expression has been observed in several neurodegenerative diseases, little is known about their contribution to ALS. To assess whether aberrant expression of repetitive element sequences are observed in ALS, we analysed RNA sequencing data from C9orf72-positive and sporadic ALS cases, as well as healthy controls. Transcripts from multiple classes and subclasses of repetitive elements (LINEs, endogenous retroviruses, DNA transposons, simple repeats, etc.) were significantly increased in the frontal cortex of C9orf72 ALS patients. A large collection of patient samples, representing both C9orf72 positive and negative ALS, ALS/FTLD, and FTLD cases, was used to validate the levels of several repetitive element transcripts. These analyses confirmed that repetitive element expression was significantly increased in C9orf72-positive compared to C9orf72-negative or control cases. While previous studies suggest an important link between TDP-43 and repetitive element biology, our data indicate that TDP-43 pathology alone is insufficient to account for the observed changes in repetitive elements in ALS/FTLD. Instead, we found that repetitive element expression positively correlated with RNA polymerase II activity in postmortem brain, and pharmacologic modulation of RNA polymerase II activity altered repetitive element expression in vitro. We conclude that increased RNA polymerase II activity in ALS/FTLD may lead to increased repetitive element transcript expression, a novel pathological feature of ALS/FTLD.

Corticobasal degeneration with TDP-43 pathology presenting with progressive supranuclear palsy syndrome: a distinct clinicopathologic subtype.

Corticobasal degeneration (CBD) is a clinically heterogeneous tauopathy, which has overlapping clinicopathologic and genetic characteristics with progressive supranuclear palsy (PSP). This study aimed to elucidate whether transactive response DNA-binding protein of 43 kDa (TDP-43) pathology contributes to clinicopathologic heterogeneity of CBD. Paraffin-embedded sections of the midbrain, pons, subthalamic nucleus, and basal forebrain from 187 autopsy-confirmed CBD cases were screened with immunohistochemistry for phospho-TDP-43. In cases with TDP-43 pathology, additional brain regions (i.e., precentral, cingulate, and superior frontal gyri, hippocampus, medulla, and cerebellum) were immunostained. Hierarchical clustering analysis was performed based on the topographical distribution and severity of TDP-43 pathology, and clinicopathologic and genetic features were compared between the clusters. TDP-43 pathology was observed in 45% of CBD cases, most frequently in midbrain tegmentum (80% of TDP-43-positive cases), followed by subthalamic nucleus (69%). TDP-43-positive CBD was divided into TDP-limited (52%) and TDP-severe (48%) by hierarchical clustering analysis. TDP-severe patients were more likely to have been diagnosed clinically as PSP compared to TDP-limited and TDP-negative patients (80 vs 32 vs 30%, P < 0.001). The presence of downward gaze palsy was the strongest factor for the antemortem diagnosis of PSP, and severe TDP-43 pathology in the midbrain tectum was strongly associated with downward gaze palsy. In addition, tau burden in the olivopontocerebellar system was significantly greater in TDP-positive than TDP-negative CBD. Genetic analyses revealed that MAPT H1/H1 genotype frequency was significantly lower in TDP-severe than in TDP-negative and TDP-limited CBD (65 vs 89 vs 91%, P < 0.001). The homozygous minor allele frequencies in GRN rs5848 and TMEM106B rs3173615 were not significantly different between the three groups. In conclusion, the present study indicates that CBD with severe TDP-43 pathology is a distinct clinicopathologic subtype of CBD, characterized by PSP-like clinical presentations, severe tau pathology in the olivopontocerebellar system, and low frequency of MAPT H1 haplotype.

Conserved DNA methylation combined with differential frontal cortex and cerebellar expression distinguishes C9orf72-associated and sporadic ALS, and implicates SERPINA1 in disease.

We previously found C9orf72-associated (c9ALS) and sporadic amyotrophic lateral sclerosis (sALS) brain transcriptomes comprise thousands of defects, among which, some are likely key contributors to ALS pathogenesis. We have now generated complementary methylome data and combine these two data sets to perform a comprehensive "multi-omic" analysis to clarify the molecular mechanisms initiating RNA misregulation in ALS. We found that c9ALS and sALS patients have generally distinct but overlapping methylome profiles, and that the c9ALS- and sALS-affected genes and pathways have similar biological functions, indicating conserved pathobiology in disease. Our results strongly implicate SERPINA1 in both C9orf72 repeat expansion carriers and non-carriers, where expression levels are greatly increased in both patient groups across the frontal cortex and cerebellum. SERPINA1 expression is particularly pronounced in C9orf72 repeat expansion carriers for both brain regions, where SERPINA1 levels are strictly down regulated across most human tissues, including the brain, except liver and blood, and are not measurable in E18 mouse brain. The altered biological networks we identified contain critical molecular players known to contribute to ALS pathology, which also interact with SERPINA1. Our comprehensive combined methylation and transcription study identifies new genes and highlights that direct genetic and epigenetic changes contribute to c9ALS and sALS pathogenesis.

A novel approach for multi-SNP GWAS and its application in Alzheimer's disease.

BACKGROUND: Genome-wide association studies (GWAS) have effectively identified genetic factors for many diseases. Many diseases, including Alzheimer's disease (AD), have epistatic causes, requiring more sophisticated analyses to identify groups of variants which together affect phenotype. RESULTS: Based on the GWAS statistical model, we developed a multi-SNP GWAS analysis to identify pairs of variants whose common occurrence signaled the Alzheimer's disease phenotype. CONCLUSIONS: Despite not having sufficient data to demonstrate significance, our preliminary experimentation identified a high correlation between GRIA3 and HLA-DRB5 (an AD gene). GRIA3 has not been previously reported in association with AD, but is known to play a role in learning and memory.

Evaluating the necessity of PCR duplicate removal from next-generation sequencing data and a comparison of approaches.

BACKGROUND: Analyzing next-generation sequencing data is difficult because datasets are large, second generation sequencing platforms have high error rates, and because each position in the target genome (exome, transcriptome, etc.) is sequenced multiple times. Given these challenges, numerous bioinformatic algorithms have been developed to analyze these data. These algorithms aim to find an appropriate balance between data loss, errors, analysis time, and memory footprint. Typical analysis pipelines require multiple steps. If one or more of these steps is unnecessary, it would significantly decrease compute time and data manipulation to remove the step. One step in many pipelines is PCR duplicate removal, where PCR duplicates arise from multiple PCR products from the same template molecule binding on the flowcell. These are often removed because there is concern they can lead to false positive variant calls. Picard (MarkDuplicates) and SAMTools (rmdup) are the two main softwares used for PCR duplicate removal. RESULTS: Approximately 92 % of the 17+ million variants called were called whether we removed duplicates with Picard or SAMTools, or left the PCR duplicates in the dataset. There were no significant differences between the unique variant sets when comparing the transition/transversion ratios (p = 1.0), percentage of novel variants (p = 0.99), average population frequencies (p = 0.99), and the percentage of protein-changing variants (p = 1.0). Results were similar for variants in the American College of Medical Genetics genes. Genotype concordance between NGS and SNP chips was above 99 % for all genotype groups (e.g., homozygous reference). CONCLUSIONS: Our results suggest that PCR duplicate removal has minimal effect on the accuracy of subsequent variant calls.

Genome-wide association study of prolactin levels in blood plasma and cerebrospinal fluid.

BACKGROUND: Prolactin is a polypeptide hormone secreted by the anterior pituitary gland that plays an essential role in lactation, tissue growth, and suppressing apoptosis to increase cell survival. Prolactin serves as a key player in many life-critical processes, including immune system and reproduction. Prolactin is also found in multiple fluids throughout the body, including plasma and cerebrospinal fluid (CSF). METHODS: In this study, we measured prolactin levels in both plasma and CSF, and performed a genome-wide association study. We then performed meta-analyses using METAL with a significance threshold of p < 5 × 10(-8) and removed SNPs where the direction of the effect was different between the two datasets. RESULTS: We identified 12 SNPs associated with increased prolactin levels in both biological fluids. CONCLUSIONS: Our efforts will help researchers understand how prolactin is regulated in both CSF and plasma, which could be beneficial in research for the immune system and reproduction.

Presenilin E318G variant and Alzheimer's disease risk: the Cache County study.

BACKGROUND: Alzheimer's disease is the leading cause of dementia in the elderly and the third most common cause of death in the United States. A vast number of genes regulate Alzheimer's disease, including Presenilin 1 (PSEN1). Multiple studies have attempted to locate novel variants in the PSEN1 gene that affect Alzheimer's disease status. A recent study suggested that one of these variants, PSEN1 E318G (rs17125721), significantly affects Alzheimer's disease status in a large case-control dataset, particularly in connection with the APOEε4 allele. METHODS: Our study looks at the same variant in the Cache County Study on Memory and Aging, a large population-based dataset. We tested for association between E318G genotype and Alzheimer's disease status by running a series of Fisher's exact tests. We also performed logistic regression to test for an additive effect of E318G genotype on Alzheimer's disease status and for the existence of an interaction between E318G and APOEε4. RESULTS: In our Fisher's exact test, it appeared that APOEε4 carriers with an E318G allele have slightly higher risk for AD than those without the allele (3.3 vs. 3.8); however, the 95 % confidence intervals of those estimates overlapped completely, indicating non-significance. Our logistic regression model found a positive but non-significant main effect for E318G (p = 0.895). The interaction term between E318G and APOEε4 was also non-significant (p = 0.689). CONCLUSIONS: Our findings do not provide significant support for E318G as a risk factor for AD in APOEε4 carriers. Our calculations indicated that the overall sample used in the logistic regression models was adequately powered to detect the sort of effect sizes observed previously. However, the power analyses of our Fisher's exact tests indicate that our partitioned data was underpowered, particularly in regards to the low number of E318G carriers, both AD cases and controls, in the Cache county dataset. Thus, the differences in types of datasets used may help to explain the difference in effect magnitudes seen. Analyses in additional case-control datasets will be required to understand fully the effect of E318G on Alzheimer's disease status.

Variants in ACPP are associated with cerebrospinal fluid Prostatic Acid Phosphatase levels.

BACKGROUND: Prostatic Acid Phosphatase (PAP) is an enzyme that is produced primarily in the prostate and functions as a cell growth regulator and potential tumor suppressor. Understanding the genetic regulation of this enzyme is important because PAP plays an important role in prostate cancer and is expressed in other tissues such as the brain. METHODS: We tested association between 5.8 M SNPs and PAP levels in cerebrospinal fluid across 543 individuals in two datasets using linear regression. We then performed meta-analyses using METAL =with a significance threshold of p < 5 × 10(-8) and removed SNPs where the direction of the effect was different between the two datasets, identifying 289 candidate SNPs that affect PAP cerebrospinal fluid levels. We analyzed each of these SNPs individually and prioritized SNPs that had biologically meaningful functional annotations in wANNOVAR (e.g. non-synonymous, stop gain, 3' UTR, etc.) or had a RegulomeDB score less than 3. RESULTS: Thirteen SNPs met our criteria, suggesting they are candidate causal alleles that underlie ACPP regulation and expression. CONCLUSIONS: Given PAP's expression in the brain and its role as a cell-growth regulator and tumor suppressor, our results have important implications in brain health such as cancer and other brain diseases including neurodegenerative diseases (e.g., Alzheimer's disease and Parkinson's disease) and mental health (e.g., anxiety, depression, and schizophrenia).

Variants in CCL16 are associated with blood plasma and cerebrospinal fluid CCL16 protein levels.

BACKGROUND: CCL16 is a chemokine predominantly expressed in the liver, but is also found in the blood and brain, and is known to play important roles in immune response and angiogenesis. Little is known about the gene's regulation. METHODS: Here, we test for potential causal SNPs that affect CCL16 protein levels in both blood plasma and cerebrospinal fluid in a genome-wide association study across two datasets. We then use METAL to performed meta-analyses with a significance threshold of p < 5x10(-8). We removed SNPs where the direction of the effect was different between the two datasets. RESULTS: We identify 10 SNPs associated with increased CCL16 protein levels in both biological fluids. CONCLUSIONS: Our results will help understand CCL16's regulation, allowing researchers to better understand the gene's effects on human health.

Interaction between variants in CLU and MS4A4E modulates Alzheimer's disease risk.

INTRODUCTION: Ebbert et al. reported gene-gene interactions between rs11136000-rs670139 (CLU-MS4A4E) and rs3865444-rs670139 (CD33-MS4A4E). We evaluate these interactions in the largest data set for an epistasis study. METHODS: We tested interactions using 3837 cases and 4145 controls from Alzheimer's Disease Genetics Consortium using meta-analyses and permutation analyses. We repeated meta-analyses stratified by apolipoprotein E (APOE) ε4 status, estimated combined odds ratio (OR) and population attributable fraction (cPAF), and explored causal variants. RESULTS: Results support the CLU-MS4A4E interaction and a dominant effect. An association between CLU-MS4A4E and APOE ε4 negative status exists. The estimated synergy factor, OR, and cPAF for rs11136000-rs670139 are 2.23, 2.45, and 8.0, respectively. We identified potential causal variants. DISCUSSION: We replicated the CLU-MS4A4E interaction in a large case-control series and observed APOE ε4 and possible dominant effect. The CLU-MS4A4E OR is higher than any Alzheimer's disease locus except APOE ε4, APP, and TREM2. We estimated an 8% decrease in Alzheimer's disease incidence without CLU-MS4A4E risk alleles and identified potential causal variants.

Variant Tool Chest: an improved tool to analyze and manipulate variant call format (VCF) files.

BACKGROUND: Since the advent of next-generation sequencing many previously untestable hypotheses have been realized. Next-generation sequencing has been used for a wide range of studies in diverse fields such as population and medical genetics, phylogenetics, microbiology, and others. However, this novel technology has created unanticipated challenges such as the large numbers of genetic variants. Each caucasian genome has more than four million single nucleotide variants, insertions and deletions, copy number variants, and structural variants. Several formats have been suggested for storing these variants; however, the variant call format (VCF) has become the community standard. RESULTS: We developed new software called the Variant Tool Chest (VTC) to provide much needed tools to work with VCF files. VTC provides a variety of tools for manipulating, comparing, and analyzing VCF files beyond the functionality of existing tools. In addition, VTC was written to be easily extended with new tools. CONCLUSIONS: Variant Tool Chest brings new and important functionality that complements and integrates well with existing software. VTC is available at https://github.com/mebbert/VariantToolChest.

Surveying the landscape of RNA isoform diversity and expression across 9 GTEx tissues using long-read sequencing data.

Even though alternative RNA splicing was discovered nearly 50 years ago (1977), we still understand very little about most isoforms arising from a single gene, including in which tissues they are expressed and if their functions differ. Human gene annotations suggest remarkable transcriptional complexity, with approximately 252,798 distinct RNA isoform annotations from 62,710 gene bodies (Ensembl v109; 2023), emphasizing the need to understand their biological effects. For example, 256 gene bodies have ≥50 annotated isoforms and 30 have ≥100, where one protein-coding gene (MAPK10) even has 192 distinct RNA isoform annotations. Whether such isoform diversity results from biological redundancy or spurious alternative splicing (i.e., noise), or whether individual isoforms have specialized functions (even if subtle) remains a mystery for most genes. Recent studies by Aguzzoli-Heberle et al., Leung et al., and Glinos et al. demonstrated long-read RNAseq enables improved RNA isoform quantification for essentially any tissue, cell type, or biological condition (e.g., disease, development, aging, etc.), making it possible to better assess individual isoform expression and function. While each study provided important discoveries related to RNA isoform diversity, deeper exploration is needed. We sought to quantify and characterize real isoform usage across tissues (compared to annotations). We used long-read RNAseq data from 58 GTEx samples across nine tissues (three brain, two heart, muscle, lung, liver, and cultured fibroblasts) generated by Glinos et al. and found considerable isoform diversity within and across tissues. Cerebellar hemisphere was the most transcriptionally complex tissue (22,522 distinct isoforms; 3,726 unique); liver was least diverse (12,435 distinct isoforms; 1,039 unique). We highlight gene clusters exhibiting high tissue-specific isoform diversity per tissue (e.g., TPM1 expresses 19 in heart's atrial appendage). We also validated 447 of the 700 new isoforms discovered by Aguzzoli-Heberle et al. and found that 88 were expressed in all nine tissues, while 58 were specific to a single tissue. This study represents a broad survey of the RNA isoform landscape, demonstrating isoform diversity across nine tissues and emphasizes the need to better understand how individual isoforms from a single gene body contribute to human health and disease.

Mapping medically relevant RNA isoform diversity in the aged human frontal cortex with deep long-read RNA-seq.

Determining whether the RNA isoforms from medically relevant genes have distinct functions could facilitate direct targeting of RNA isoforms for disease treatment. Here, as a step toward this goal for neurological diseases, we sequenced 12 postmortem, aged human frontal cortices (6 Alzheimer disease cases and 6 controls; 50% female) using one Oxford Nanopore PromethION flow cell per sample. We identified 1,917 medically relevant genes expressing multiple isoforms in the frontal cortex where 1,018 had multiple isoforms with different protein-coding sequences. Of these 1,018 genes, 57 are implicated in brain-related diseases including major depression, schizophrenia, Parkinson's disease and Alzheimer disease. Our study also uncovered 53 new RNA isoforms in medically relevant genes, including several where the new isoform was one of the most highly expressed for that gene. We also reported on five mitochondrially encoded, spliced RNA isoforms. We found 99 differentially expressed RNA isoforms between cases with Alzheimer disease and controls.

A phase II study of UCN-01 in combination with irinotecan in patients with metastatic triple negative breast cancer.

Mutations in TP53 lead to a defective G1 checkpoint and the dependence on checkpoint kinase 1 (Chk1) for G2 or S phase arrest in response to DNA damage. In preclinical studies, Chk1 inhibition resulted in enhanced cytotoxicity of several chemotherapeutic agents. The high frequency of TP53 mutations in triple negative breast cancer (TNBC: negative for estrogen receptor, progesterone receptor, and HER2) make Chk1 an attractive therapeutic target. UCN-01, a non-selective Chk1 inhibitor, combined with irinotecan demonstrated activity in advanced TNBC in our Phase I study. The goal of this trial was to further evaluate this treatment in women with TNBC. Patients with metastatic TNBC previously treated with anthracyclines and taxanes received irinotecan (100-125 mg/m(2) IV days 1, 8, 15, 22) and UCN-01 (70 mg/m(2) IV day 2, 35 mg/m(2) day 23 and subsequent doses) every 42-day cycle. Peripheral blood mononuclear cells (PBMC) and tumor specimens were collected. Twenty five patients were enrolled. The overall response (complete response (CR) + partial response (PR)) rate was 4 %. The clinical benefit rate (CR + PR + stable disease ≥6 months) was 12 %. Since UCN-01 inhibits PDK1, phosphorylated ribosomal protein S6 (pS6) in PBMC was assessed. Although reduced 24 h post UCN-01, pS6 levels rose to baseline by day 8, indicating loss of UCN-01 bioavailability. Immunostains of γH2AX and pChk1(S296) on serial tumor biopsies from four patients demonstrated an induction of DNA damage and Chk1 activation following irinotecan. However, Chk1 inhibition by UCN-01 was not observed in all tumors. Most tumors were basal-like (69 %), and carried mutations in TP53 (53 %). Median overall survival in patients with TP53 mutant tumors was poor compared to wild type (5.5 vs. 20.3 months, p = 0.004). This regimen had limited activity in TNBC. Inconsistent Chk1 inhibition was likely due to the pharmacokinetics of UCN-01. TP53 mutations were associated with a poor prognosis in metastatic TNBC.

PAM50 proliferation score as a predictor of weekly paclitaxel benefit in breast cancer.

To identify a group of patients who might benefit from the addition of weekly paclitaxel to conventional anthracycline-containing chemotherapy as adjuvant therapy of node-positive operable breast cancer. The predictive value of PAM50 subtypes and the 11-gene proliferation score contained within the PAM50 assay were evaluated in 820 patients from the GEICAM/9906 randomized phase III trial comparing adjuvant FEC to FEC followed by weekly paclitaxel (FEC-P). Multivariable Cox regression analyses of the secondary endpoint of overall survival (OS) were performed to determine the significance of the interaction between treatment and the (1) PAM50 subtypes, (2) PAM50 proliferation score, and (3) clinical and pathological variables. Similar OS analyses were performed in 222 patients treated with weekly paclitaxel versus paclitaxel every 3 weeks in the CALGB/9342 and 9840 metastatic clinical trials. In GEICAM/9906, with a median follow up of 8.7 years, OS of the FEC-P arm was significantly superior compared to the FEC arm (unadjusted HR = 0.693, p = 0.013). A benefit from paclitaxel was only observed in the group of patients with a low PAM50 proliferation score (unadjusted HR = 0.23, p < 0.001; and interaction test, p = 0.006). No significant interactions between treatment and the PAM50 subtypes or the various clinical-pathological variables, including Ki-67 and histologic grade, were identified. Finally, similar OS results were obtained in the CALGB data set, although the interaction test did not reach statistical significance (p = 0.109). The PAM50 proliferation score identifies a subset of patients with a low proliferation status that may derive a larger benefit from weekly paclitaxel.

TRIM29 functions as a tumor suppressor in nontumorigenic breast cells and invasive ER+ breast cancer.

Tripartite motif-containing 29 (TRIM29) is a member of the TRIM protein family that has been implicated in hematologic and solid tumor cancers. We found that TRIM29 functions as a tumor suppressor in both the nontumorigenic MCF10A [estrogen receptor (ER)-/TRIM29+] breast cell line and the invasive MCF7 (ER+/TRIM29-) breast cell line. Silencing TRIM29 in MCF10A cells resulted in preneoplastic changes that included loss of polarity in three-dimensional culture, increased proliferation, anchorage-independent growth, and increased migration and invasion. Conversely, the introduction of TRIM29 into MCF7 cells caused reversion to a less aggressive phenotype by antagonizing the growth effect of 17ß-estradiol. The interaction between TRIM29 and ER signaling in MCF7 cells was supported by a reduction in ERE binding in the presence of TRIM29 and suppression of ER-dependent gene expression of TFF1, FOS, and GREB1. By microarray analyses, we showed that younger women (<55 years of age) with early-stage, ER+ breast cancer who were given no adjuvant systemic therapy had a significantly lower risk of relapse when their tumor had high TRIM29 expression (P = 0.02). This effect was not observed in older women (>55 years of age) and thus may be due to menopause and loss of circulating estrogens. Our results suggest that loss of TRIM29 expression in normal breast luminal cells can contribute to malignant transformation and lead to progression of ER+ breast cancer in premenopausal women.

Using deep long-read RNAseq in Alzheimer's disease brain to assess medical relevance of RNA isoform diversity.

Due to alternative splicing, human protein-coding genes average over eight RNA isoforms, resulting in nearly four distinct protein coding sequences per gene. Long-read RNAseq (IsoSeq) enables more accurate quantification of isoforms, shedding light on their specific roles. To assess the medical relevance of measuring RNA isoform expression, we sequenced 12 aged human frontal cortices (6 Alzheimer's disease cases and 6 controls; 50% female) using one Oxford Nanopore PromethION flow cell per sample. Our study uncovered 53 new high-confidence RNA isoforms in medically relevant genes, including several where the new isoform was one of the most highly expressed for that gene. Specific examples include WDR4 (61%; microcephaly), MYL3 (44%; hypertrophic cardiomyopathy), and MTHFS (25%; major depression, schizophrenia, bipolar disorder). Other notable genes with new high-confidence isoforms include CPLX2 (10%; schizophrenia, epilepsy) and MAOB (9%; targeted for Parkinson's disease treatment). We identified 1,917 medically relevant genes expressing multiple isoforms in human frontal cortex, where 1,018 had multiple isoforms with different protein coding sequences, demonstrating the need to better understand how individual isoforms from a single gene body are involved in human health and disease, if at all. Exactly 98 of the 1,917 genes are implicated in brain-related diseases, including Alzheimer's disease genes such as APP (Aß precursor protein; five), MAPT (tau protein; four), and BIN1 (eight). As proof of concept, we also found 99 differentially expressed RNA isoforms between Alzheimer's cases and controls, despite the genes themselves not exhibiting differential expression. Our findings highlight the significant knowledge gaps in RNA isoform diversity and their medical relevance. Deep long-read RNA sequencing will be necessary going forward to fully comprehend the medical relevance of individual isoforms for a "single" gene.

Agreement in risk prediction between the 21-gene recurrence score assay (Oncotype DX®) and the PAM50 breast cancer intrinsic Classifier™ in early-stage estrogen receptor-positive breast cancer.

PURPOSE: To compare risk assignment by PAM50 Breast Cancer Intrinsic Classifier™ and Oncotype DX_Recurrence Score (RS) in the same population. METHODS: RNA was extracted from 151 estrogen receptor (ER)+ stage I-II breast cancers and gene expression profiled using PAM50 "intrinsic" subtyping test. RESULTS: One hundred eight cases had complete molecular information; 103 (95%) were classified as luminal A (n = 76) or luminal B (n = 27). Ninety-two percent (n = 98) had a low (n = 59) or intermediate (n = 39) RS. Among luminal A cancers, 70% had low (n = 53) and the remainder (n = 23) had an intermediate RS. Among luminal B cancers, nine were high (33%) and 13 were intermediate (48%) by the RS. Almost all cancers with a high RS were classified as luminal B (90%, n = 9). One high RS cancer was identified as basal-like and had low ER/ESR1 and low human epidermal growth factor receptor 2 (HER2) expression by quantitative polymerase chain reaction in both assays. The majority of low RS cases were luminal A (83%, n = 53). Importantly, half of the intermediate RS cancers were re-categorized as low risk luminal A subtype by PAM50. CONCLUSION: There is good agreement between the two assays for high (i.e., luminal B or RS > 31) and low (i.e., luminal B or RS < 18) prognostic risk assignment but PAM50 assigns more patients to the low risk category. About half of the intermediate RS group was reclassified as luminal A by PAM50.

Transcriptional response to VZV infection is modulated by RNA polymerase III in lung epithelial cell lines.

Ancestral RNA polymerase III (Pol III) is a multi-subunit polymerase responsible for transcription of short non-coding RNA, such as double-stranded short interspersed nuclear elements (SINEs). Although SINE ncRNAs are generally transcriptionally repressed, they can be induced in response to viral infections and can stimulate immune signaling pathways. Indeed, mutations in RNA Pol III have been associated with poor antiviral interferon response following infection with varicella zoster virus (VZV). In this study, we probed the role of Pol III transcripts in the detection and initial immune response to VZV by characterizing the transcriptional response following VZV infection of wild type A549 lung epithelial cells as well as A549 cells lacking specific RNA sensors MAVS and TLR3, or interferon-stimulated genes RNase L and PKR in presence or absence of functional RNA Pol III. Multiple components of the antiviral sensing and interferon signaling pathways were involved in restricting VZV replication in lung epithelial cells thus suggesting an innate defense system with built-in redundancy. In addition, RNA Pol III silencing altered the antiviral transcriptional program indicating that it plays an essential role in the sensing of VZV infection.