Identification of a gene expression signature of vascular invasion and recurrence in stage I lung adenocarcinoma via bulk and spatial transcriptomics.

Microscopic vascular invasion (VI) is predictive of recurrence and benefit from lobectomy in stage I lung adenocarcinoma (LUAD) but is difficult to assess in resection specimens and cannot be accurately predicted prior to surgery. Thus, new biomarkers are needed to identify this aggressive subset of stage I LUAD tumors. To assess molecular and microenvironment features associated with angioinvasive LUAD we profiled 162 resected stage I tumors with and without VI by RNA-seq and explored spatial patterns of gene expression in a subset of 15 samples by high-resolution spatial transcriptomics (stRNA-seq). Despite the small size of invaded blood vessels, we identified a gene expression signature of VI from the bulk RNA-seq discovery cohort (n=103) and found that it was associated with VI foci, desmoplastic stroma, and high-grade patterns in our stRNA-seq data. We observed a stronger association with high-grade patterns from VI+ compared with VI- tumors. Using the discovery cohort, we developed a transcriptomic predictor of VI, that in an independent validation cohort (n=60) was associated with VI (AUROC=0.86; p=5.42×10-6) and predictive of recurrence-free survival (HR=1.98; p=0.024), even in VI- LUAD (HR=2.76; p=0.003). To determine our VI predictor's robustness to intra-tumor heterogeneity we used RNA-seq data from multi-region sampling of stage I LUAD cases in TRACERx, where the predictor scores showed high correlation (R=0.87, p<2.2×10-16) between two randomly sampled regions of the same tumor. Our study suggests that VI-associated gene expression changes are detectable beyond the site of intravasation and can be used to predict the presence of VI. This may enable the prediction of angioinvasive LUAD from biopsy specimens, allowing for more tailored medical and surgical management of stage I LUAD.

Systematic dissection, preservation, and multiomics in whole human and bovine hearts.

OBJECTIVES: We sought to develop a rigorous, systematic protocol for the dissection and preservation of human hearts for biobanking that expands previous success in postmortem transcriptomics to multiomics from paired tissue. BACKGROUND: Existing cardiac biobanks consist largely of biopsy tissue or explanted hearts in select diseases and are insufficient for correlating whole organ phenotype with clinical data. METHODS: We demonstrate optimal conditions for multiomics interrogation (ribonucleic acid (RNA) sequencing, untargeted metabolomics) in hearts by evaluating the effect of technical variables (storage solution, temperature) and simulated postmortem interval (PMI) on RNA and metabolite stability. We used bovine (n=3) and human (n=2) hearts fixed in PAXgene or snap-frozen with liquid nitrogen. RESULTS: Using a paired Wald test, only two of the genes assessed were differentially expressed between left ventricular samples from bovine hearts stored in PAXgene at 0 and 12 hours PMI (FDR q<0.05). We obtained similar findings in human left ventricular samples, suggesting stability of RNA transcripts at PMIs up to 12 hours. Different library preparation methods (mRNA poly-A capture vs. rRNA depletion) resulted in similar quality metrics with both library preparations achieving >95% of reads properly aligning to the reference genomes across all PMIs for bovine and human hearts. PMI had no effect on RNA Integrity Number or quantity of RNA recovered at the time points evaluated. Of the metabolites identified (855 total) using untargeted metabolomics of human left ventricular tissue, 503 metabolites remained stable across PMIs (0, 4, 8, 12 hours). Most metabolic pathways retained several stable metabolites. CONCLUSIONS: Our data demonstrate a technically rigorous, reproducible protocol that will enhance cardiac biobanking practices and facilitate novel insights into human CVD. CONDENSED ABSTRACT: Cardiovascular disease (CVD) is the leading cause of mortality worldwide. Current biobanking practices insufficiently capture both the diverse array of phenotypes present in CVDs and the spatial heterogeneity across cardiac tissue sites. We have developed a rigorous and systematic protocol for the dissection and preservation of human cardiac biospecimens to enhance the availability of whole organ tissue for multiple applications. When combined with longitudinal clinical phenotyping, our protocol will enable multiomics in hearts to deepen our understanding of CVDs.

Clinical Study of Aspirin and Zileuton on Biomarkers of Tobacco-Related Carcinogenesis in Current Smokers.

The chemopreventive effect of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) on lung cancer risk is supported by epidemiologic and preclinical studies. Zileuton, a 5-lipoxygenase inhibitor, has additive activity with NSAIDs against tobacco carcinogenesis in preclinical models. We hypothesized that cyclooxygenase plus 5-lipoxygenase inhibition would be more effective than a placebo in modulating the nasal epithelium gene signatures of tobacco exposure and lung cancer. We conducted a randomized, double-blinded study of low-dose aspirin plus zileuton vs. double placebo in current smokers to compare the modulating effects on nasal gene expression and arachidonic acid metabolism. In total, 63 participants took aspirin 81 mg daily plus zileuton (Zyflo CR) 600 mg BID or the placebo for 12 weeks. Nasal brushes from the baseline, end-of-intervention, and one-week post intervention were profiled via microarray. Aspirin plus zilueton had minimal effects on the modulation of the nasal or bronchial gene expression signatures of smoking, lung cancer, and COPD but favorably modulated a bronchial gene expression signature of squamous dysplasia. Aspirin plus zileuton suppressed urinary leukotriene but not prostaglandin E2, suggesting shunting through the cyclooxygenase pathway when combined with 5-lipoxygenase inhibition. Continued investigation of leukotriene inhibitors is needed to confirm these findings, understand the long-term effects on the airway epithelium, and identify the safest, optimally dosed agents.

Multisite Evaluation of Next-Generation Methods for Small RNA Quantification.

Small RNAs (smRNAs) are important regulators of many biologic processes and are now most frequently characterized using Illumina sequencing. However, although standard RNA sequencing library preparation has become routine in most sequencing facilities, smRNA sequencing library preparation has historically been challenging because of high input requirements, laborious protocols involving gel purifications, inability to automate, and a lack of benchmarking standards. Additionally, studies have suggested that many of these methods are nonlinear and do not accurately reflect the amounts of smRNAs in vivo. Recently, a number of new kits have become available that permit lower input amounts and less laborious, gel-free protocol options. Several of these new kits claim to reduce RNA ligase-dependent sequence bias through novel adapter modifications and to lessen adapter-dimer contamination in the resulting libraries. With the increasing number of smRNA kits available, understanding the relative strengths of each method is crucial for appropriate experimental design. In this study, we systematically compared 9 commercially available smRNA library preparation kits as well as NanoString probe hybridization across multiple study sites. Although several of the new methodologies do reduce the amount of artificially over- and underrepresented microRNAs (miRNAs), we observed that none of the methods was able to remove all of the bias in the library preparation. Identical samples prepared with different methods show highly varied levels of different miRNAs. Even so, many methods excelled in ease of use, lower input requirement, fraction of usable reads, and reproducibility across sites. These differences may help users select the most appropriate methods for their specific question of interest.

Cross-site comparison of ribosomal depletion kits for Illumina RNAseq library construction.

BACKGROUND: Ribosomal RNA (rRNA) comprises at least 90% of total RNA extracted from mammalian tissue or cell line samples. Informative transcriptional profiling using massively parallel sequencing technologies requires either enrichment of mature poly-adenylated transcripts or targeted depletion of the rRNA fraction. The latter method is of particular interest because it is compatible with degraded samples such as those extracted from FFPE and also captures transcripts that are not poly-adenylated such as some non-coding RNAs. Here we provide a cross-site study that evaluates the performance of ribosomal RNA removal kits from Illumina, Takara/Clontech, Kapa Biosystems, Lexogen, New England Biolabs and Qiagen on intact and degraded RNA samples. RESULTS: We find that all of the kits are capable of performing significant ribosomal depletion, though there are differences in their ease of use. All kits were able to remove ribosomal RNA to below 20% with intact RNA and identify ~ 14,000 protein coding genes from the Universal Human Reference RNA sample at >1FPKM. Analysis of differentially detected genes between kits suggests that transcript length may be a key factor in library production efficiency. CONCLUSIONS: These results provide a roadmap for labs on the strengths of each of these methods and how best to utilize them.

Characterization of FUS mutations in amyotrophic lateral sclerosis using RNA-Seq.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease resulting in severe muscle weakness and eventual death by respiratory failure. Although little is known about its pathogenesis, mutations in fused in sarcoma/translated in liposarcoma (FUS) are causative for familial ALS. FUS is a multifunctional protein that is involved in many aspects of RNA processing. To elucidate the role of FUS in ALS, we overexpressed wild-type and two mutant forms of FUS in HEK-293T cells, as well as knocked-down FUS expression. This was followed by RNA-Seq to identify genes which displayed differential expression or altered splicing patterns. Pathway analysis revealed that overexpression of wild-type FUS regulates ribosomal genes, whereas knock-down of FUS additionally affects expression of spliceosome related genes. Furthermore, cells expressing mutant FUS displayed global transcription patterns more similar to cells overexpressing wild-type FUS than to the knock-down condition. This observation suggests that FUS mutants do not contribute to the pathogenesis of ALS through a loss-of-function. Finally, our results demonstrate that the R521G and R522G mutations display differences in their influence on transcription and splicing. Taken together, these results provide additional insights into the function of FUS and how mutations contribute to the development of ALS.

Mutations in the profilin 1 gene cause familial amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a late-onset neurodegenerative disorder resulting from motor neuron death. Approximately 10% of cases are familial (FALS), typically with a dominant inheritance mode. Despite numerous advances in recent years, nearly 50% of FALS cases have unknown genetic aetiology. Here we show that mutations within the profilin 1 (PFN1) gene can cause FALS. PFN1 is crucial for the conversion of monomeric (G)-actin to filamentous (F)-actin. Exome sequencing of two large ALS families showed different mutations within the PFN1 gene. Further sequence analysis identified 4 mutations in 7 out of 274 FALS cases. Cells expressing PFN1 mutants contain ubiquitinated, insoluble aggregates that in many cases contain the ALS-associated protein TDP-43. PFN1 mutants also display decreased bound actin levels and can inhibit axon outgrowth. Furthermore, primary motor neurons expressing mutant PFN1 display smaller growth cones with a reduced F/G-actin ratio. These observations further document that cytoskeletal pathway alterations contribute to ALS pathogenesis.

Prevalence of Huntington's disease gene CAG repeat alleles in sporadic amyotrophic lateral sclerosis patients.

A higher prevalence of intermediate ataxin-2 CAG repeats in amyotrophic lateral sclerosis (ALS) patients has raised the possibility that CAG expansions in other polyglutamine disease genes could contribute to ALS neurodegeneration. We sought to determine whether expansions of the CAG repeat of the HTT gene that causes Huntington's disease, are associated with ALS. We compared the HTT CAG repeat length on a total of 3144 chromosomes from 1572 sporadic ALS patients and 4007 control chromosomes, and also tested its possible effects on ALS-specific parameters, such as age and site of onset and survival rate. Our results show that the CAG repeat in the HTT gene is not a risk factor for ALS nor modifies its clinical presentation. These findings suggest that distinct neuronal degeneration processes are involved in these two different neurodegenerative disorders.

Angiogenin variants in Parkinson disease and amyotrophic lateral sclerosis.

OBJECTIVE: Several studies have suggested an increased frequency of variants in the gene encoding angiogenin (ANG) in patients with amyotrophic lateral sclerosis (ALS). Interestingly, a few ALS patients carrying ANG variants also showed signs of Parkinson disease (PD). Furthermore, relatives of ALS patients have an increased risk to develop PD, and the prevalence of concomitant motor neuron disease in PD is higher than expected based on chance occurrence. We therefore investigated whether ANG variants could predispose to both ALS and PD. METHODS: We reviewed all previous studies on ANG in ALS and performed sequence experiments on additional samples, which allowed us to analyze data from 6,471 ALS patients and 7,668 controls from 15 centers (13 from Europe and 2 from the USA). We sequenced DNA samples from 3,146 PD patients from 6 centers (5 from Europe and 1 from the USA). Statistical analysis was performed using the variable threshold test, and the Mantel-Haenszel procedure was used to estimate odds ratios. RESULTS: Analysis of sequence data from 17,258 individuals demonstrated a significantly higher frequency of ANG variants in both ALS and PD patients compared to control subjects (p = 9.3 × 10(-6) for ALS and p = 4.3 × 10(-5) for PD). The odds ratio for any ANG variant in patients versus controls was 9.2 for ALS and 6.7 for PD. INTERPRETATION: The data from this multicenter study demonstrate that there is a strong association between PD, ALS, and ANG variants. ANG is a genetic link between ALS and PD.

Mutational analysis of TARDBP in neurodegenerative diseases.

Neurodegenerative diseases are often characterized by the presence of aggregates of misfolded proteins. TDP-43 is a major component of these aggregates in amyotrophic lateral sclerosis (ALS), but has also been observed in Alzheimer's (AD) and Parkinson's Diseases (PD). In addition, mutations in the TARDBP gene, encoding TDP-43, have been found to be a significant cause of familial ALS (FALS). All mutations, except for one, have been found in exon 6. To confirm this observation in ALS and to investigate whether TARDBP may play a role in the pathogenesis of AD and PD, we screened for mutations in exon 6 of the TARDBP gene in three cohorts composed of 376 AD, 463 PD (18% familial PD) and 376 ALS patients (50% FALS). We found mutations in ∼ 7% of FALS and ∼0.5% of sporadic ALS (SALS) patients, including two novel mutations, p.N352T and p.G384R. In contrast, we did not find TARDBP mutations in our cohort of AD and PD patients. These results suggest that mutations in TARDBP are not a significant cause of AD and PD.

Paraoxonase gene mutations in amyotrophic lateral sclerosis.

Three clustered, homologous paraoxonase genes (PON1, PON2, and PON3) have roles in preventing lipid oxidation and detoxifying organophosphates. Recent reports describe a genetic association between the PON genes and sporadic amyotrophic lateral sclerosis (ALS). We now report that in genomic DNA from individuals with familial and sporadic ALS, we have identified at least 7 PON gene mutations that are predicted to alter PON function.