Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis.

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

Identification of a candidate gene for the I locus determining the dominant white bulb color in onion (Allium cepa L.).

KEY MESSAGE: Through a map-based cloning approach, a gene coding for an R2R3-MYB transcription factor was identified as a causal gene for the I locus controlling the dominant white bulb color in onion. White bulb colors in onion (Allium cepa L.) are determined by either the C or I loci. The causal gene for the C locus was previously isolated, but the gene responsible for the I locus has not been identified yet. To identify candidate genes for the I locus, an approximately 7-Mb genomic DNA region harboring the I locus was obtained from onion and bunching onion (A. fistulosum) whole genome sequences using two tightly linked molecular markers. Within this interval, the AcMYB1 gene, known as a positive regulator of anthocyanin production, was identified. No polymorphic sequences were found between white and red AcMYB1 alleles in the 4,860-bp full-length genomic DNA sequences. However, a 4,838-bp LTR-retrotransposon was identified in the white allele, in the 79-bp upstream coding region from the stop codon. The insertion of this LTR-retrotransposon created a premature stop codon, resulting in the replacement of 26 amino acids with seven different residues. A molecular marker was developed based on the insertion of this LTR-retrotransposon to genotype the I locus. A perfect linkage between bulb color phenotypes and marker genotypes was observed among 5,303 individuals of segregating populations. The transcription of AcMYB1 appeared to be normal in both red and white onions, but the transcription of CHS-A, which encodes chalcone synthase and is involved in the first step of the anthocyanin biosynthesis pathway, was inactivated in the white onions. Taken together, an aberrant AcMYB1 protein produced from the mutant allele might be responsible for the dominant white bulb color in onions.

Specific DNMT3C flanking sequence preferences facilitate methylation of young murine retrotransposons.

The DNA methyltransferase DNMT3C appeared as a duplication of the DNMT3B gene in muroids and is required for silencing of young retrotransposons in the male germline. Using specialized assay systems, we investigate the flanking sequence preferences of DNMT3C and observe characteristic preferences for cytosine at the -2 and -1 flank that are unique among DNMT3 enzymes. We identify two amino acids in the catalytic domain of DNMT3C (C543 and V547) that are responsible for the DNMT3C-specific flanking sequence preferences and evolutionary conserved in muroids. Reanalysis of published data shows that DNMT3C flanking preferences are consistent with genome-wide methylation patterns in mouse ES cells only expressing DNMT3C. Strikingly, we show that CpG sites with the preferred flanking sequences of DNMT3C are enriched in murine retrotransposons that were previously identified as DNMT3C targets. Finally, we demonstrate experimentally that DNMT3C has elevated methylation activity on substrates derived from these biological targets. Our data show that DNMT3C flanking sequence preferences match the sequences of young murine retrotransposons which facilitates their methylation. By this, our data provide mechanistic insights into the molecular co-evolution of repeat elements and (epi)genetic defense systems dedicated to maintain genomic stability in mammals.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression.

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

Endogenous retroelement expression in the gut microenvironment of people living with HIV-1.

BACKGROUND: Endogenous retroelements (EREs), including human endogenous retroviruses (HERVs) and long interspersed nuclear elements (LINEs), comprise almost half of the human genome. Our previous studies of the interferome in the gut suggest potential mechanisms regarding how IFNb may drive HIV-1 gut pathogenesis. As ERE activity is suggested to partake in type 1 immune responses and is incredibly sensitive to viral infections, we sought to elucidate underlying interactions between ERE expression and gut dynamics in people living with HIV-1 (PLWH). METHODS: ERE expression profiles from bulk RNA sequencing of colon biopsies and PBMC were compared between a cohort of PLWH not on antiretroviral therapy (ART) and uninfected controls. FINDINGS: 59 EREs were differentially expressed in the colon of PLWH when compared to uninfected controls (padj <0.05 and FC ≤ -1 or ≥ 1) [Wald's Test]. Of these 59, 12 EREs were downregulated in PLWH and 47 were upregulated. Colon expression of the ERE loci LTR19_12p13.31 and L1FLnI_1q23.1s showed significant correlations with certain gut immune cell subset frequencies in the colon. Furthermore L1FLnI_1q23.1s showed a significant upregulation in peripheral blood mononuclear cells (PBMCs) of PLWH when compared to uninfected controls suggesting a common mechanism of differential ERE expression in the colon and PBMC. INTERPRETATION: ERE activity has been largely understudied in genomic characterizations of human pathologies. We show that the activity of certain EREs in the colon of PLWH is deregulated, supporting our hypotheses that their underlying activity could function as (bio)markers and potential mediators of pathogenesis in HIV-1 reservoirs. FUNDING: US NIH grants NCI CA260691 (DFN) and NIAID UM1AI164559 (DFN).

Comprehensive profiling of L1 retrotransposons in mouse.

L1 elements are retrotransposons currently active in mammals. Although L1s are typically silenced in most normal tissues, elevated L1 expression is associated with a variety of conditions, including cancer, aging, infertility and neurological disease. These associations have raised interest in the mapping of human endogenous de novo L1 insertions, and a variety of methods have been developed for this purpose. Adapting these methods to mouse genomes would allow us to monitor endogenous in vivo L1 activity in controlled, experimental conditions using mouse disease models. Here, we use a modified version of transposon insertion profiling, called nanoTIPseq, to selectively enrich young mouse L1s. By linking this amplification step with nanopore sequencing, we identified >95% annotated L1s from C57BL/6 genomic DNA using only 200 000 sequencing reads. In the process, we discovered 82 unannotated L1 insertions from a single C57BL/6 genome. Most of these unannotated L1s were near repetitive sequence and were not found with short-read TIPseq. We used nanoTIPseq on individual mouse breast cancer cells and were able to identify the annotated and unannotated L1s, as well as new insertions specific to individual cells, providing proof of principle for using nanoTIPseq to interrogate retrotransposition activity at the single-cell level in vivo.

Regulatory SVA retrotransposons and classical HLA genotyped-transcripts associated with Parkinson's disease.

Introduction: Parkinson's disease (PD) is a neurodegenerative and polygenic disorder characterised by the progressive loss of neural dopamine and onset of movement disorders. We previously described eight SINE-VNTR-Alu (SVA) retrotransposon-insertion-polymorphisms (RIPs) located and expressed within the Human Leucocyte Antigen (HLA) genomic region of chromosome 6 that modulate the differential co-expression of 71 different genes including the HLA classical class I and class II genes in a Parkinson's Progression Markers Initiative (PPMI) cohort. Aims and methods: In the present study, we (1) reanalysed the PPMI genomic and transcriptomic sequencing data obtained from whole blood of 1521 individuals (867 cases and 654 controls) to infer the genotypes of the transcripts expressed by eight classical HLA class I and class II genes as well as DRA and the DRB3/4/5 haplotypes, and (2) examined the statistical differences between three different PD subgroups (cases) and healthy controls (HC) for the HLA and SVA transcribed genotypes and inferred haplotypes. Results: Significant differences for 57 expressed HLA alleles (21 HLA class I and 36 HLA class II alleles) up to the three-field resolution and four of eight expressed SVA were detected at p<0.05 by the Fisher's exact test within one or other of three different PD subgroups (750 individuals with PD, 57 prodromes, 60 individuals who had scans without evidence of dopamine deficits [SWEDD]), when compared against a group of 654 HCs within the PPMI cohort and when not corrected by the Bonferroni test for multiple comparisons. Fourteen of 20 significant alleles were unique to the PD-HC comparison, whereas 31 of the 57 alleles overlapped between two or more different subgroup comparisons. Only the expressed HLA-DRA*01:01:01 and -DQA1*03:01:01 protective alleles (PD v HC), the -DQA1*03:03:01 risk (HC v Prodrome) or protective allele (PD v Prodrome), the -DRA*01:01:02 and -DRB4*01:03:02 risk alleles (SWEDD v HC), and the NR_SVA_381 present genotype (PD v HC) at a 5% homozygous insertion frequency near HLA-DPA1, were significant (Pc<0.1) after Bonferroni corrections. The homologous NR_SVA_381 insertion significantly decreased the transcription levels of HLA-DPA1 and HLA-DPB1 in the PPMI cohort and its presence as a homozygous genotype is a risk factor (Pc=0.012) for PD. The most frequent NR_SVA_381 insertion haplotype in the PPMI cohort was NR_SVA_381/DPA1*02/DPB1*01 (3.7%). Although HLA C*07/B*07/DRB5*01/DRB1*15/DQB1*06 was the most frequent HLA 5-loci phased-haplotype (n, 76) in the PPMI cohort, the NR_SVA_381 insertion was present in only six of them (8%). Conclusions: These data suggest that expressed SVA and HLA gene alleles in circulating white blood cells are coordinated differentially in the regulation of immune responses and the long-term onset and progression of PD, the mechanisms of which have yet to be elucidated.

Identification and expression analysis of two steamer-like retrotransposons in the Chilean blue mussel (Mytilus chilensis).

BACKGROUND: Disseminated neoplasia (DN) is a proliferative cell disorder of the circulatory system of bivalve mollusks. The disease is transmitted between individuals and can also be induced by external chemical agents such as bromodeoxyuridine. In Mya arenaria, we have cloned and characterized an LTR-retrotransposon named Steamer. Steamer mRNA levels and gene copy number correlates with DN and can be used as a marker of the disease. So far, the only mollusk where a retrotransposon expression relates to DN is Mya arenaria. On the other hand, it has been reported that the Chilean blue mussel Mytilus chilensis can also suffers DN. Our aim was to identify retrotransposons in Mytilus chilensis and to study their expression levels in the context of disseminated neoplasia. RESULTS: Here we show that 7.1% of individuals collected in August 2018, from two farming areas, presents morphological characteristics described in DN. Using Steamer sequence to interrogate the transcriptome of M. chilensis we found two putative retrotransposons, named Steamer-like elements (MchSLEs). MchSLEs are present in the genome of M. chilensis and MchSLE1 is indeed an LTR-retrotransposon. Neither expression, nor copy number of the reported MchSLEs correlate with DN status but both are expressed at different levels among individual animals. We also report that in cultured M. chilensis haemocytes MchSLEs1 expression can be induced by bromodeoxyuridine. CONCLUSIONS: We conclude that SLEs present in Mytilus chilensis are differentially expressed among individuals and do not correlate with disseminated neoplasia. Treatment of haemocytes with a stressor like bromodeoxyuridine induces expression of MchSLE1 suggesting that in Mytilus chilensis environmental stressors can induce activation of LTR-retrotransposon.

Cytogenomic Characterization of Transposable Elements and Satellite DNA in Passiflora L. Species.

The species Passiflora alata, P. cincinnata, and P. edulis have great economic value due to the use of their fruits for human consumption. In this study, we compared the repetitive genome fractions of these three species. The compositions of the repetitive DNA of these three species' genomes were analyzed using clustering and identification of the repetitive sequences with RepeatExplorer. It was found that repetitive DNA content represents 74.70%, 66.86%, and 62.24% of the genome of P. alata, P. edulis, and P. cincinnata, respectively. LTR Ty3/Gypsy retrotransposons represent the highest genome proportions in P. alata and P. edulis, while Ty1/Copia comprises the largest proportion of P. cincinnata genome. Chromosomal mapping by Fluorescent In Situ Hybridization (FISH) showed that LTR retrotransposons have a dispersed distribution along chromosomes. The subtelomeric region of chromosomes is where 145 bp satellite DNA is located, suggesting that these elements may play important roles in genome structure and organization in these species. In this work, we obtained the first global characterization of the composition of repetitive DNA in Passiflora, showing that an increase in genome size is related to an increase in repetitive DNA, which represents an important evolutionary route for these species.

Genomes of historical specimens reveal multiple invasions of LTR retrotransposons in Drosophila melanogaster during the 19th century.

Transposable element invasions have a profound impact on the evolution of genomes and phenotypes. It is thus an important open question how often such TE invasions occur. To address this question, we utilize the genomes of historical specimens, sampled about 200 y ago. We found that the LTR retrotransposons Blood, Opus, and 412 spread in Drosophila melanogaster in the 19th century. These invasions constitute second waves, as degraded fragments were found for all three TEs. The composition of Opus and 412, but not of Blood, shows a pronounced geographic heterogeneity, likely due to founder effects during the invasions. Finally, we identified species from the Drosophila simulans complex as the likely origin of the TEs. We show that in total, seven TE families invaded D. melanogaster during the last 200y, thereby increasing the genome size by up to 1.2Mbp. We suggest that this high rate of TE invasions was likely triggered by human activity. Based on the analysis of strains and specimens sampled at different times, we provide a detailed timeline of TE invasions, making D. melanogaster the first organism where the invasion history of TEs during the last two centuries could be inferred.

The Mytilus chilensis Steamer-like Element-1 Retrotransposon Antisense mRNA Harbors an Internal Ribosome Entry Site That Is Modulated by hnRNPK.

LTR-retrotransposons are transposable elements characterized by the presence of long terminal repeats (LTRs) directly flanking an internal coding region. They share genome organization and replication strategies with retroviruses. Steamer-like Element-1 (MchSLE-1) is an LTR-retrotransposon identified in the genome of the Chilean blue mussel Mytilus chilensis. MchSLE-1 is transcribed; however, whether its RNA is also translated and the mechanism underlying such translation remain to be elucidated. Here, we characterize the MchSLE-1 translation mechanism. We found that the MchSLE-1 5' and 3'LTRs command transcription of sense and antisense RNAs, respectively. Using luciferase reporters commanded by the untranslated regions (UTRs) of MchSLE-1, we found that in vitro 5'UTR sense is unable to initiate translation, whereas the antisense 5'UTR initiates translation even when the eIF4E-eIF4G interaction was disrupted, suggesting the presence of an internal ribosomal entry site (IRES). The antisense 5'UTR IRES activity was tested using bicistronic reporters. The antisense 5'UTR has IRES activity only when the mRNA is transcribed in the nucleus, suggesting that nuclear RNA-binding proteins are required to modulate its activity. Indeed, heterogeneous nuclear ribonucleoprotein K (hnRNPK) was identified as an IRES trans-acting factor (ITAF) of the MchSLE-1 IRES. To our knowledge, this is the first report describing an IRES in an antisense mRNA derived from a mussel LTR-retrotransposon.

The structure, function, and evolution of plant centromeres.

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Systematic evaluation of retroviral LTRs as cis-regulatory elements in mouse embryos.

In mammals, many retrotransposons are de-repressed during zygotic genome activation (ZGA). However, their functions in early development remain elusive largely due to the challenge to simultaneously manipulate thousands of retrotransposon insertions in embryos. Here, we applied CRISPR interference (CRISPRi) to perturb the long terminal repeat (LTR) MT2_Mm, a well-known ZGA and totipotency marker that exists in ∼2,667 insertions throughout the mouse genome. CRISPRi robustly perturbed 2,485 (∼93%) MT2_Mm insertions and 1,090 (∼55%) insertions of the closely related MT2C_Mm in 2-cell embryos. Remarkably, such perturbation caused downregulation of hundreds of ZGA genes and embryonic arrest mostly at the morula stage. Mechanistically, MT2 LTRs are globally enriched for open chromatin and H3K27ac and function as promoters/enhancers downstream of OBOX/DUX proteins. Thus, we not only provide direct evidence to support the functional importance of MT2 activation in development but also systematically define cis-regulatory function of MT2 in embryos by integrating functional perturbation and multi-omic analyses.

RNA m6A modification regulates L1 retrotransposons in human spermatogonial stem cell differentiation in vitro and in vivo.

The maintenance of genome integrity in the germline is crucial for mammalian development. Long interspersed element type 1 (LINE-1, L1) is a mobile genetic element that makes up about 17% of the human genome and poses a threat to genome integrity. N6-methyl-adenosine (m6A) plays an essential role in regulating various biological processes. However, the function of m6A modification in L1 retrotransposons and human germline development remains largely unknown. Here we knocked out the m6A methyltransferase METTL3 or the m6A reader YTHDF2 in human embryonic stem cells (hESCs) and discovered that METTL3 and YTHDF2 are crucial for inducing human spermatogonial stem cells (hSSCs) from hESCs in vitro. The removal of METTL3 or YTHDF2 resulted in increased L1 retrotransposition and reduced the efficiency of SSC differentiation in vitro. Further analysis showed that YTHDF2 recognizes the METTL3-catalyzed m6A modification of L1 retrotransposons and degrades L1 mRNA through autophagy, thereby blocking L1 retrotransposition. Moreover, the study confirmed that m6A modification in human fetal germ cells promotes the degradation of L1 retrotransposon RNA, preventing the insertion of new L1 retrotransposons into the genome. Interestingly, L1 retrotransposon RNA was highly expressed while METTL3 was significantly downregulated in the seminal plasma of azoospermic patients with meiotic arrest compared to males with normal fertility. Additionally, we identified some potentially pathogenic variants in m6A-related genes in azoospermic men with meiotic arrest. In summary, our study suggests that m6A modification serves as a guardian of genome stability during human germline development and provides novel insights into the function and regulatory mechanisms of m6A modification in restricting L1 retrotransposition.

The Deficiency of Hypusinated eIF5A Decreases the Putrescine/Spermidine Ratio and Inhibits +1 Programmed Ribosomal Frameshifting during the Translation of Ty1 Retrotransposon in Saccharomyces cerevisiae.

Programmed ribosomal frameshifting (PRF) exists in all branches of life that regulate gene expression at the translational level. The eukaryotic translation initiation factor 5A (eIF5A) is a highly conserved protein essential in all eukaryotes. It is identified initially as an initiation factor and functions broadly in translation elongation and termination. The hypusination of eIF5A is specifically required for +1 PRF at the shifty site derived from the ornithine decarboxylase antizyme 1 (OAZ1) in Saccharomyces cerevisiae. However, whether the regulation of +1 PRF by yeast eIF5A is universal remains unknown. Here, we found that Sc-eIF5A depletion decreased the putrescine/spermidine ratio. The re-introduction of Sc-eIF5A in yeast eIF5A mutants recovered the putrescine/spermidine ratio. In addition, the Sc-eIF5A depletion decreases +1 PRF during the decoding of Ty1 retrotransposon mRNA, but has no effect on -1 PRF during the decoding of L-A virus mRNA. The re-introduction of Sc-eIF5A in yeast eIF5A mutants restored the +1 PRF rate of Ty1. The inhibition of the hypusine modification of yeast eIF5A by GC7 treatment or by mutating the hypusination site Lys to Arg caused decreases of +1 PRF rates in the Ty1 retrotransposon. Furthermore, mutational studies of the Ty1 frameshifting element support a model where the efficient removal of ribosomal subunits at the first Ty1 frame 0 stop codon is required for the frameshifting of trailing ribosomes. This dependency is likely due to the unique position of the frame 0 stop codon distance from the slippery sequence of Ty1. The results showed that eIF5A is a trans-regulator of +1 PRF for Ty1 retrotransposon and could function universally in yeast.

Locus-level L1 DNA methylation profiling reveals the epigenetic and transcriptional interplay between L1s and their integration sites.

Long interspersed element 1 (L1) retrotransposons are implicated in human disease and evolution. Their global activity is repressed by DNA methylation, but deciphering the regulation of individual copies has been challenging. Here, we combine short- and long-read sequencing to unveil L1 methylation heterogeneity across cell types, families, and individual loci and elucidate key principles involved. We find that the youngest primate L1 families are specifically hypomethylated in pluripotent stem cells and the placenta but not in most tumors. Locally, intronic L1 methylation is intimately associated with gene transcription. Conversely, the L1 methylation state can propagate to the proximal region up to 300 bp. This phenomenon is accompanied by the binding of specific transcription factors, which drive the expression of L1 and chimeric transcripts. Finally, L1 hypomethylation alone is typically insufficient to trigger L1 expression due to redundant silencing pathways. Our results illuminate the epigenetic and transcriptional interplay between retrotransposons and their host genome.

Every repeat is unique: Exploring the genomic impact of human L1 retrotransposons at locus-specific resolution.

Fully understanding the impact of the human retrotransposon L1 requires that each of ∼500,000 L1 copies be evaluated as a potentially unique genomic entity. In this issue of Cell Genomics, Lanciano et al.1 strive toward this goal, illuminating the reciprocal regulatory influence between individual L1s and their genomic integration sites.

Pharmacological inhibition of neddylation impairs long interspersed element 1 retrotransposition.

Aberrant long interspersed element 1 (LINE-1 or L1) activity can cause insertional mutagenesis and chromosomal rearrangements and has been detected in several types of cancers. Here, we show that neddylation, a post-translational modification process, is essential for L1 transposition. The antineoplastic drug MLN4924 is an L1 inhibitor that suppresses NEDD8-activating enzyme activity. Neddylation inhibition by MLN4924 selectively impairs ORF2p-mediated L1 reverse transcription and blocks the generation of L1 cDNA. Consistent with these results, MLN4924 treatment suppresses the retrotransposition activity of the non-autonomous retrotransposons short interspersed nuclear element R/variable number of tandem repeat/Alu and Alu, which rely on the reverse transcription activity of L1 ORF2p. The E2 enzyme UBE2M in the neddylation pathway, rather than UBE2F, is required for L1 ORF2p and retrotransposition. Interference with the functions of certain neddylation-dependent Cullin-really interesting new gene E3 ligases disrupts L1 reverse transcription and transposition activity. Our findings provide insights into the regulation of L1 retrotransposition and the identification of therapeutic targets for L1 dysfunctions.

Locus-specific LINE-1 expression in clinical ovarian cancer specimens at the single-cell level.

Long interspersed nuclear elements (LINE-1s/L1s) are a group of retrotransposons that can copy themselves within a genome. In humans, it is the most successful transposon in nucleotide content. L1 expression is generally mild in normal human tissues, but the activity has been shown to increase significantly in many cancers. Few studies have examined L1 expression at single-cell resolution, thus it is undetermined whether L1 reactivation occurs solely in malignant cells within tumors. One of the cancer types with frequent L1 activity is high-grade serous ovarian carcinoma (HGSOC). Here, we identified locus-specific L1 expression with 3' single-cell RNA sequencing in pre- and post-chemotherapy HGSOC sample pairs from 11 patients, and in fallopian tube samples from five healthy women. Although L1 expression quantification with the chosen technique was challenging due to the repetitive nature of the element, we found evidence of L1 expression primarily in cancer cells, but also in other cell types, e.g. cancer-associated fibroblasts. The expression levels were similar in samples taken before and after neoadjuvant chemotherapy, indicating that L1 transcriptional activity was unaffected by clinical platinum-taxane treatment. Furthermore, L1 activity was negatively associated with the expression of MYC target genes, a finding that supports earlier literature of MYC being an L1 suppressor.