Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

The role of retrotransposable elements in ageing and age-associated diseases.

The genomes of virtually all organisms contain repetitive sequences that are generated by the activity of transposable elements (transposons). Transposons are mobile genetic elements that can move from one genomic location to another; in this process, they amplify and increase their presence in genomes, sometimes to very high copy numbers. In this Review we discuss new evidence and ideas that the activity of retrotransposons, a major subgroup of transposons overall, influences and even promotes the process of ageing and age-related diseases in complex metazoan organisms, including humans. Retrotransposons have been coevolving with their host genomes since the dawn of life. This relationship has been largely competitive, and transposons have earned epithets such as 'junk DNA' and 'molecular parasites'. Much of our knowledge of the evolution of retrotransposons reflects their activity in the germline and is evident from genome sequence data. Recent research has provided a wealth of information on the activity of retrotransposons in somatic tissues during an individual lifespan, the molecular mechanisms that underlie this activity, and the manner in which these processes intersect with our own physiology, health and well-being.

LINE-1 retrotransposon expression in cancerous, epithelial and neuronal cells revealed by 5' single-cell RNA-Seq.

LINE-1 retrotransposons are sequences capable of copying themselves to new genomic loci via an RNA intermediate. New studies implicate LINE-1 in a range of diseases, especially in the context of aging, but without an accurate understanding of where and when LINE-1 is expressed, a full accounting of its role in health and disease is not possible. We therefore developed a method-5' scL1seq-that makes use of a widely available library preparation method (10x Genomics 5' single cell RNA-seq) to measure LINE-1 expression in tens of thousands of single cells. We recapitulated the known pattern of LINE-1 expression in tumors-present in cancer cells, absent from immune cells-and identified hitherto undescribed LINE-1 expression in human epithelial cells and mouse hippocampal neurons. In both cases, we saw a modest increase with age, supporting recent research connecting LINE-1 to age related diseases.

LINE-1 expression in cancer correlates with p53 mutation, copy number alteration, and S phase checkpoint.

Retrotransposons are genomic DNA sequences that copy themselves to new genomic locations via RNA intermediates; LINE-1 is the only active and autonomous retrotransposon in the human genome. The mobility of LINE-1 is largely repressed in somatic tissues but is derepressed in many cancers, where LINE-1 retrotransposition is correlated with p53 mutation and copy number alteration (CNA). In cell lines, inducing LINE-1 expression can cause double-strand breaks (DSBs) and replication stress. Reanalyzing multiomic data from breast, ovarian, endometrial, and colon cancers, we confirmed correlations between LINE-1 expression, p53 mutation status, and CNA. We observed a consistent correlation between LINE-1 expression and the abundance of DNA replication complex components, indicating that LINE-1 may also induce replication stress in human tumors. In endometrial cancer, high-quality phosphoproteomic data allowed us to identify the DSB-induced ATM-MRN-SMC S phase checkpoint pathway as the primary DNA damage response (DDR) pathway associated with LINE-1 expression. Induction of LINE-1 expression in an in vitro model led to increased phosphorylation of MRN complex member RAD50, suggesting that LINE-1 directly activates this pathway.

Mapping of long interspersed element-1 (L1) insertions by TIPseq provides information about sub chromosomal genetic variation in human embryos.

PURPOSE: Retrotransposons play important roles during early development when they are transiently de-repressed during epigenetic reprogramming. Long interspersed element-1 (L1), the only autonomous retrotransposon in humans, comprises 17% of the human genome. We applied the Single Cell Transposon Insertion Profiling by Sequencing (scTIPseq) to characterize and map L1 insertions in human embryos. METHODS: Sixteen cryopreserved, genetically tested, human blastocysts, were accessed from consenting couples undergoing IVF at NYU Langone Fertility Center. Additionally, four trios (father, mother, and embryos) were also evaluated. scTIPseq was applied to map L1 insertions in all samples, using L1 locations reported in the 1000 Genomes as controls. RESULTS: Twenty-nine unknown and unique insertions were observed in the sixteen embryos. Most were intergenic; no insertions were located in exons or immediately upstream of genes. The location or number of unknown insertions did not differ between euploid and aneuploid embryos, suggesting they are not merely markers of aneuploidy. Rather, scTIPseq provides novel information about sub-chromosomal structural variation in human embryos. Trio analyses showed a parental origin of all L1 insertions in embryos. CONCLUSION: Several studies have measured L1 expression at different stages of development in mice, but this study for the first time reports unknown insertions in human embryos that were inherited from one parent, confirming no de novo L1 insertions occurred in parental germline or during embryogenesis. Since one-third of euploid embryo transfers fail, future studies would be useful for understanding whether these sub-chromosomal genetic variants or de novo L1 insertions affect embryo developmental potential.

Parasite manipulation of host phenotypes inferred from transcriptional analyses in a trematode-amphipod system.

Manipulation of host phenotypes by parasites is hypothesized to be an adaptive strategy enhancing parasite transmission across hosts and generations. Characterizing the molecular mechanisms of manipulation is important to advance our understanding of host-parasite coevolution. The trematode (Levinseniella byrdi) is known to alter the colour and behaviour of its amphipod host (Orchestia grillus) presumably increasing predation of amphipods which enhances trematode transmission through its life cycle. We sampled 24 infected and 24 uninfected amphipods from a salt marsh in Massachusetts to perform differential gene expression analysis. In addition, we constructed novel genomic tools for O. grillus including a de novo genome and transcriptome. We discovered that trematode infection results in upregulation of amphipod transcripts associated with pigmentation and detection of external stimuli, and downregulation of multiple amphipod transcripts implicated in invertebrate immune responses, such as vacuolar ATPase genes. We hypothesize that suppression of immune genes and the altered expression of genes associated with coloration and behaviour may allow the trematode to persist in the amphipod and engage in further biochemical manipulation that promotes transmission. The genomic tools and transcriptomic analyses reported provide new opportunities to discover how parasites alter diverse pathways underlying host phenotypic changes in natural populations.

Transposon insertion profiling by sequencing (TIPseq) identifies novel LINE-1 insertions in human sperm.

PURPOSE: Long interspersed nuclear element-1 (LINE-1 or L1) comprises 17% of the human genome. Retrotransposons may perturb gene integrity or alter gene expression by altering regulatory regions in the genome. The germline employs a number of mechanisms, including cytosine methylation, to repress retrotransposon transcription throughout most of life. Demethylation during germ cell and early embryo development de-represses retrotransposons. Intriguingly, de novo genetic variation appearing in sperm has been implicated in a number of disorders in offspring, including autism spectrum disorder, schizophrenia, and bipolar disorder. We hypothesize that human sperm exhibit de novo retrotransposition and employ a new sequencing method, single cell transposon insertion profiling by sequencing (scTIPseq) to map them in small amounts of human sperm. METHODS: Cross-sectional case-control study of sperm samples (n=10 men; ages 32-55 years old) from consenting men undergoing IVF at NYU Langone Fertility Center. scTIPseq identified novel LINE-1 insertions in individual sperm and TIPseqHunter, a custom bioinformatics pipeline, compared the architecture of sperm LINE-1 to known LINE-1 insertions from the European database of Human specific LINE-1 (L1Hs) retrotransposon insertions (euL1db). RESULTS: scTIPseq identified 17 novel insertions in sperm. New insertions were mainly intergenic or intronic. Only one sample did not exhibit new insertions. The location or number of novel insertions did not differ by paternal age. CONCLUSION: This study for the first time reports novel LINE-1 insertions in human sperm, demonstrating the feasibility of scTIPseq, and identifies new contributors to genetic diversity in the human germ line.

L1EM: a tool for accurate locus specific LINE-1 RNA quantification.

MOTIVATION: LINE-1 elements are retrotransposons that are capable of copying their sequence to new genomic loci. LINE-1 derepression is associated with a number of disease states, and has the potential to cause significant cellular damage. Because LINE-1 elements are repetitive, it is difficult to quantify LINE-1 RNA at specific loci and to separate transcripts with protein coding capability from other sources of LINE-1 RNA. RESULTS: We provide a tool, L1EM that uses the expectation maximization algorithm to quantify LINE-1 RNA at each genomic locus, separating transcripts that are capable of generating retrotransposition from those that are not. We show the accuracy of L1EM on simulated data and against long read sequencing from HEK cells. AVAILABILITY AND IMPLEMENTATION: L1EM is written in python. The source code along with the necessary annotations are available at https://github.com/FenyoLab/L1EM and distributed under GPLv3. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization and visualization of RNA secondary structure Boltzmann ensemble via information theory.

BACKGROUND: The nearest neighbor model and associated dynamic programming algorithms allow for the efficient estimation of the RNA secondary structure Boltzmann ensemble. However because a given RNA secondary structure only contains a fraction of the possible helices that could form from a given sequence, the Boltzmann ensemble is multimodal. Several methods exist for clustering structures and finding those modes. However less focus is given to exploring the underlying reasons for this multimodality: the presence of conflicting basepairs. Information theory, or more specifically mutual information, provides a method to identify those basepairs that are key to the secondary structure. RESULTS: To this end we find most informative basepairs and visualize the effect of these basepairs on the secondary structure. Knowing whether a most informative basepair is present tells us not only the status of the particular pair but also provides a large amount of information about which other pairs are present or not present. We find that a few basepairs account for a large amount of the structural uncertainty. The identification of these pairs indicates small changes to sequence or stability that will have a large effect on structure. CONCLUSION: We provide a novel algorithm that uses mutual information to identify the key basepairs that lead to a multimodal Boltzmann distribution. We then visualize the effect of these pairs on the overall Boltzmann ensemble.

Predicting RNA hyper-editing with a novel tool when unambiguous alignment is impossible.

BACKGROUND: Repetitive elements are now known to have relevant cellular functions, including self-complementary sequences that form double stranded (ds) RNA. There are numerous pathways that determine the fate of endogenous dsRNA, and misregulation of endogenous dsRNA is a driver of autoimmune disease, particularly in the brain. Unfortunately, the alignment of high-throughput, short-read sequences to repeat elements poses a dilemma: Such sequences may align equally well to multiple genomic locations. In order to differentiate repeat elements, current alignment methods depend on sequence variation in the reference genome. Reads are discarded when no such variations are present. However, RNA hyper-editing, a possible fate for dsRNA, introduces enough variation to distinguish between repeats that are otherwise identical. RESULTS: To take advantage of this variation, we developed a new algorithm, RepProfile, that simultaneously aligns reads and predicts novel variations. RepProfile accurately aligns hyper-edited reads that other methods discard. In particular we predict hyper-editing of Drosophila melanogaster repeat elements in vivo at levels previously described only in vitro, and provide validation by Sanger sequencing sixty-two individual cloned sequences. We find that hyper-editing is concentrated in genes involved in cell-cell communication at the synapse, including some that are associated with neurodegeneration. We also find that hyper-editing tends to occur in short runs. CONCLUSIONS: Previous studies of RNA hyper-editing discarded ambiguously aligned reads, ignoring hyper-editing in long, perfect dsRNA - the perfect substrate for hyper-editing. We provide a method that simulation and Sanger validation show accurately predicts such RNA editing, yielding a superior picture of hyper-editing.

Quantification of LINE-1 RNA Expression from Bulk RNA-seq Using L1EM.

LINE-1 retrotransposons have the potential to cause DNA damage, contribute to genome instability, and induce an interferon response. Thus, accurate measurements of their expression, especially in disease contexts where genome instability and the interferon response are relevant, are of particular importance. Illumina-based bulk RNA sequencing remains the most abundant datatype for measuring gene expression. However, "active" expression from its own internal promoter is only one source of LINE-1 aligning reads in an RNA-seq experiment. With about half a million LINE-1 sequences scattered throughout the genome, many are incorporated into other transcripts that have nothing to do with LINE-1 activity. We call this "passive" co-transcription. Here we will describe how to use L1EM, a computational method that separates active from passive LINE-1 expression at the locus-specific level.

Vitrification with Dimethyl Sulfoxide Induces Transcriptomic Alteration of Gene and Transposable Element Expression in Immature Human Oocytes.

Despite substantial advancements in the field of cryobiology, oocyte and embryo cryopreservation still compromise developmental competence. Furthermore, dimethyl sulfoxide (DMSO), one of the most commonly used cryoprotectants, has been found to exert potent effects on the epigenetic landscape of cultured human cells, as well as mouse oocytes and embryos. Little is known about its impact on human oocytes. Additionally, few studies investigate the effects of DMSO on transposable elements (TE), the control of which is essential for the maintenance of genomic instability. The objective of this study was to investigate the impact of vitrification with DMSO-containing cryoprotectant on the transcriptome, including on TEs, of human oocytes. Twenty-four oocytes at the GV stage were donated by four healthy women undergoing elective oocyte cryopreservation. Oocytes were paired such that half from each patient were vitrified with DMSO-containing cryoprotectant (Vitrified Cohort), while the other half were snap frozen in phosphate buffer, unexposed to DMSO (Non-Vitrified Cohort). All oocytes underwent RNA sequencing via a method with high fidelity for single cell analysis, and which allows for the analysis of TE expression through Switching Mechanism at the 5'-end of the RNA Transcript sequencing 2 (SMARTseq2), followed by functional enrichment analysis. Of the 27,837 genes identified by SMARTseq2, 7331 (26.3%) were differentially expressed (p < 0.05). There was a significant dysregulation of genes involved in chromatin and histone modification. Mitochondrial function, as well as the Wnt, insulin, mTOR, HIPPO, and MAPK signaling pathways were also altered. The expression of TEs was positively correlated with the expression of PIWIL2, DNMT3A, and DNMT3B, and negatively correlated with age. These findings suggest that the current standard process of oocyte vitrification, involving DMSO-containing cryoprotectant, induces significant transcriptome changes, including those involving TEs.

Affinity-Based Interactome Analysis of Endogenous LINE-1 Macromolecules.

During their proliferation and the host's concomitant attempts to suppress it, LINE-1 (L1) retrotransposons give rise to a collection of heterogeneous ribonucleoproteins (RNPs); their protein and RNA compositions remain poorly defined. The constituents of L1-associated macromolecules can differ depending on numerous factors, including, for example, position within the L1 life cycle, whether the macromolecule is productive or under suppression, and the cell type within which the proliferation is occurring. This chapter describes techniques that aid the capture and characterization of protein and RNA components of L1 macromolecules from tissues that natively express them. The protocols described have been applied to embryonal carcinoma cell lines that are popular model systems for L1 molecular biology (e.g., N2102Ep, NTERA-2, and PA-1 cells), as well as colorectal cancer tissues. N2102Ep cells are given as the use case for this chapter; the protocols should be applicable to essentially any tissue exhibiting endogenous L1 expression with minor modifications.

Proteogenomic data and resources for pan-cancer analysis.

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.

Proteogenomic insights suggest druggable pathways in endometrial carcinoma.

We characterized a prospective endometrial carcinoma (EC) cohort containing 138 tumors and 20 enriched normal tissues using 10 different omics platforms. Targeted quantitation of two peptides can predict antigen processing and presentation machinery activity, and may inform patient selection for immunotherapy. Association analysis between MYC activity and metformin treatment in both patients and cell lines suggests a potential role for metformin treatment in non-diabetic patients with elevated MYC activity. PIK3R1 in-frame indels are associated with elevated AKT phosphorylation and increased sensitivity to AKT inhibitors. CTNNB1 hotspot mutations are concentrated near phosphorylation sites mediating pS45-induced degradation of ß-catenin, which may render Wnt-FZD antagonists ineffective. Deep learning accurately predicts EC subtypes and mutations from histopathology images, which may be useful for rapid diagnosis. Overall, this study identified molecular and imaging markers that can be further investigated to guide patient stratification for more precise treatment of EC.

RIP-seq reveals LINE-1 ORF1p association with p-body enriched mRNAs.

BACKGROUND: Long INterspersed Element-1 (LINE-1) is an autonomous retroelement able to "copy-and-paste" itself into new loci of the host genome through a process called retrotransposition. The LINE-1 bicistronic mRNA codes for two proteins, ORF1p, a nucleic acid chaperone, and ORF2p, a protein with endonuclease and reverse transcriptase activity. Both proteins bind LINE-1 mRNA in cis and are necessary for retrotransposition. While LINE-1 transcription is usually repressed in most healthy somatic cells through a plethora of mechanisms, ORF1p expression has been observed in nearly 50% of tumors, and new LINE-1 insertions have been documented in a similar fraction of tumors, including prostate cancer. RESULTS: Here, we utilized RNA ImmunoPrecipitation (RIP) and the L1EM analysis software to identify ORF1p bound RNA in prostate cancer cells. We identified LINE-1 loci that were expressed in parental androgen sensitive and androgen independent clonal derivatives. In all androgen independent cells, we found higher levels of LINE-1 RNA, as well as unique expression patterns of LINE-1 loci. Interestingly, we observed that ORF1p bound many non-LINE-1 mRNA in all prostate cancer cell lines evaluated, and polyA RNA, and RNA localized in p-bodies were especially enriched. Furthermore, the expression levels of RNAs identified in our ORF1p RIP correlated with RNAs expressed in LINE-1 positive tumors from The Cancer Genome Atlas (TCGA). CONCLUSION: Our results show a significant remodeling of LINE-1 loci expression in androgen independent cell lines when compared to parental androgen dependent cells. Additionally, we found that ORF1p bound a significant amount of non-LINE-1 mRNA, and that the enriched ORF1p bound mRNAs are also amplified in LINE-1 expressing TCGA prostate tumors, indicating the biological relevance of our findings to prostate cancer.

Human transposon insertion profiling by sequencing (TIPseq) to map LINE-1 insertions in single cells.

Long interspersed element-1 (LINE-1, L1) sequences, which comprise about 17% of human genome, are the product of one of the most active types of mobile DNAs in modern humans. LINE-1 insertion alleles can cause inherited and de novo genetic diseases, and LINE-1-encoded proteins are highly expressed in some cancers. Genome-wide LINE-1 mapping in single cells could be useful for defining somatic and germline retrotransposition rates, and for enabling studies to characterize tumour heterogeneity, relate insertions to transcriptional and epigenetic effects at the cellular level, or describe cellular phylogenies in development. Our laboratories have reported a genome-wide LINE-1 insertion site mapping method for bulk DNA, named transposon insertion profiling by sequencing (TIPseq). There have been significant barriers applying LINE-1 mapping to single cells, owing to the chimeric artefacts and features of repetitive sequences. Here, we optimize a modified TIPseq protocol and show its utility for LINE-1 mapping in single lymphoblastoid cells. Results from single-cell TIPseq experiments compare well to known LINE-1 insertions found by whole-genome sequencing and TIPseq on bulk DNA. Among the several approaches we tested, whole-genome amplification by multiple displacement amplification followed by restriction enzyme digestion, vectorette ligation and LINE-1-targeted PCR had the best assay performance. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

LINE-1 ORF2p expression is nearly imperceptible in human cancers.

BACKGROUND: Long interspersed element-1 (LINE-1, L1) is the major driver of mobile DNA activity in modern humans. When expressed, LINE-1 loci produce bicistronic transcripts encoding two proteins essential for retrotransposition, ORF1p and ORF2p. Many types of human cancers are characterized by L1 promoter hypomethylation, L1 transcription, L1 ORF1p protein expression, and somatic L1 retrotransposition. ORF2p encodes the endonuclease and reverse transcriptase activities required for L1 retrotransposition. Its expression is poorly characterized in human tissues and cell lines. RESULTS: We report mass spectrometry-based tumor proteome profiling studies wherein ORF2p eludes detection. To test whether ORF2p could be detected with specific reagents, we developed and validated five rabbit monoclonal antibodies with immunoreactivity for specific epitopes on the protein. These reagents readily detect ectopic ORF2p expressed from bicistronic L1 constructs. However, endogenous ORF2p is not detected in human tumor samples or cell lines by western blot, immunoprecipitation, or immunohistochemistry despite high levels of ORF1p expression. Moreover, we report endogenous ORF1p-associated interactomes, affinity isolated from colorectal cancers, wherein we similarly fail to detect ORF2p. These samples include primary tumors harboring hundreds of somatically acquired L1 insertions. The new data are available via ProteomeXchange with identifier PXD013743. CONCLUSIONS: Although somatic retrotransposition provides unequivocal genetic evidence for the expression of ORF2p in human cancers, we are unable to directly measure its presence using several standard methods. Experimental systems have previously indicated an unequal stoichiometry between ORF1p and ORF2p, but in vivo, the expression of these two proteins may be more strikingly uncoupled. These findings are consistent with observations that ORF2p is not tolerable for cell growth.

SmCL3, a gastrodermal cysteine protease of the human blood fluke Schistosoma mansoni.

BACKGROUND: Blood flukes of the genus Schistosoma are platyhelminth parasites that infect 200 million people worldwide. Digestion of nutrients from the host bloodstream is essential for parasite development and reproduction. A network of proteolytic enzymes (proteases) facilitates hydrolysis of host hemoglobin and serum proteins. METHODOLOGY/PRINCIPAL FINDINGS: We identified a new cathepsin L termed SmCL3 using PCR strategies based on S. mansoni EST sequence data. An ortholog is present in Schistosoma japonicum. SmCL3 was heterologously expressed as an active enzyme in the yeast, Pichia pastoris. Recombinant SmCL3 has a broad pH activity range against peptidyl substrates and is inhibited by Clan CA protease inhibitors. Consistent with a function in degrading host proteins, SmCL3 hydrolyzes serum albumin and hemoglobin, is localized to the adult gastrodermis, and is expressed mainly in those life stages infecting the mammalian host. The predominant form of SmCL3 in the parasite exists as a zymogen, which is unusual for proteases. This zymogen includes an unusually long prodomain with alpha helical secondary structure motifs. The striking specificity of SmCL3 for amino acids with large aromatic side chains (Trp and Tyr) at the P2 substrate position, as determined with positional scanning-synthetic combinatorial library, is consistent with a molecular model that shows a large and deep S2 pocket. A sequence similarity network (SSN) view clusters SmCL3 and other cathepsins L in accordance with previous large-scale phylogenetic analyses that identify six super kingdoms. CONCLUSIONS/SIGNIFICANCE: SmCL3 is a gut-associated cathepsin L that may contribute to the network of proteases involved in degrading host blood proteins as nutrients. Furthermore, this enzyme exhibits some unusual sequence and biophysical features that may result in additional functions. The visualization of network inter-relationships among cathepsins L suggests that these enzymes are suitable 'marker sequences' for inclusion in future phylogenetic analyses.