Proteomic and genomic signatures of repeat instability in cancer and adjacent normal tissues

Variation within major internal repeats of KSHV in vivo

Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiologic agent of Kaposi's sarcoma (KS), yet the viral genetic factors that lead to the development of KS in KSHV-infected individuals have not been fully elucidated. Nearly, all previous analyses of KSHV genomic evolution and diversity have excluded the three major internal repeat regions: the two origins of lytic replication, internal repeats 1 and 2 (IR1 and IR2), and the latency-associated nuclear antigen (LANA) repeat domain (LANAr). These regions encode protein domains that are essential to the KSHV infection cycle but have been rarely sequenced due to their extended repetitive nature and high guanine and cytosine (GC) content. The limited data available suggest that their sequences and repeat lengths are more heterogeneous across individuals than in the remainder of the KSHV genome. To assess their diversity, the full-length IR1, IR2, and LANAr sequences, tagged with unique molecular identifiers (UMIs), were obtained by Pacific Biosciences' single-molecule real-time sequencing (SMRT-UMI) from twenty-four tumors and six matching oral swabs from sixteen adults in Uganda with advanced KS. Intra-host single-nucleotide variation involved an average of 0.16 per cent of base positions in the repeat regions compared to a nearly identical average of 0.17 per cent of base positions in the remainder of the genome. Tandem repeat unit (TRU) counts varied by only one from the intra-host consensus in a majority of individuals. Including the TRU indels, the average intra-host pairwise identity was 98.3 per cent for IR1, 99.6 per cent for IR2 and 98.9 per cent for LANAr. More individuals had mismatches and variable TRU counts in IR1 (twelve/sixteen) than in IR2 (two/sixteen). There were no open reading frames in the Kaposin coding sequence inside IR2 in at least fifty-five of ninety-six sequences. In summary, the KSHV major internal repeats, like the rest of the genome in individuals with KS, have low diversity. IR1 was the most variable among the repeats, and no intact Kaposin reading frames were present in IR2 of the majority of genomes sampled.

Recurrent repeat expansions in human cancer genomes

Expansion of a single repetitive DNA sequence, termed a tandem repeat (TR), is known to cause more than 50 diseases^1,2. However, repeat expansions are often not explored beyond neurological and neurodegenerative disorders. In some cancers, mutations accumulate in short tracts of TRs, a phenomenon termed microsatellite instability; however, larger repeat expansions have not been systematically analysed in cancer^3-8. Here we identified TR expansions in 2,622 cancer genomes spanning 29 cancer types. In seven cancer types, we found 160 recurrent repeat expansions (rREs), most of which (155/160) were subtype specific. We found that rREs were non-uniformly distributed in the genome with enrichment near candidate cis-regulatory elements, suggesting a potential role in gene regulation. One rRE, a GAAA-repeat expansion, located near a regulatory element in the first intron of UGT2B7 was detected in 34% of renal cell carcinoma samples and was validated by long-read DNA sequencing. Moreover, in preliminary experiments, treating cells that harbour this rRE with a GAAA-targeting molecule led to a dose-dependent decrease in cell proliferation. Overall, our results suggest that rREs may be an important but unexplored source of genetic variation in human cancer, and we provide a comprehensive catalogue for further study.

Genome-wide survey of D/E repeats in human proteins uncovers their instability and aids in identifying their role in the chromatin regulator ATAD2

D/E repeats are stretches of aspartic and/or glutamic acid residues found in over 150 human proteins. We examined genomic stability of D/E repeats and functional characteristics of D/E repeat-containing proteins vis-à-vis the proteins with poly-Q or poly-A repeats, which are known to undergo pathologic expansions. Mining of tumor sequencing data revealed that D/E repeat-coding regions are similar to those coding poly-Qs and poly-As in increased incidence of trinucleotide insertions/deletions but differ in types and incidence of substitutions. D/E repeat-containing proteins preferentially function in chromatin metabolism and are the more likely to be nuclear and interact with core histones, the longer their repeats are. One of the longest D/E repeats of unknown function is in ATAD2, a bromodomain family ATPase frequently overexpressed in tumors. We demonstrate that D/E repeat deletion in ATAD2 suppresses its binding to nascent and mature chromatin and to the constitutive pericentromeric heterochromatin, where ATAD2 represses satellite transcription.

•

Many human proteins contain runs of aspartic/glutamic acid residues (D/E repeats)

•

D/E repeats show increased incidence of in-frame insertions/deletions in tumors

•

Nuclear and histone-interacting proteins often have long D/E repeats

•

D/E repeat of the oncogene ATAD2 controls its binding to pericentric chromatin

Analysis of lineage-specific protein family variability in prokaryotes combined with evolutionary reconstructions

Background

Evolutionary rate is a key characteristic of gene families that is linked to the functional importance of the respective genes as well as specific biological functions of the proteins they encode. Accurate estimation of evolutionary rates is a challenging task that requires precise phylogenetic analysis. Here we present an easy to estimate protein family level measure of sequence variability based on alignment column homogeneity in multiple alignments of protein sequences from Clade-Specific Clusters of Orthologous Genes (csCOGs).

Results

We report genome-wide estimates of variability for 8 diverse groups of bacteria and archaea and investigate the connection between variability and various genomic and biological features. The variability estimates are based on homogeneity distributions across amino acid sequence alignments and can be obtained for multiple groups of genomes at minimal computational expense. About half of the variance in variability values can be explained by the analyzed features, with the greatest contribution coming from the extent of gene paralogy in the given csCOG. The correlation between variability and paralogy appears to originate, primarily, not from gene duplication, but from acquisition of distant paralogs and xenologs, introducing sequence variants that are more divergent than those that could have evolved in situ during the lifetime of the given group of organisms. Both high-variability and low-variability csCOGs were identified in all functional categories, but as expected, proteins encoded by integrated mobile elements as well as proteins involved in defense functions and cell motility are, on average, more variable than proteins with housekeeping functions. Additionally, using linear discriminant analysis, we found that variability and fraction of genomes carrying a given gene are the two variables that provide the best prediction of gene essentiality as compared to the results of transposon mutagenesis in Sulfolobus islandicus.

Conclusions

Variability, a measure of sequence diversity within an alignment relative to the overall diversity within a group of organisms, offers a convenient proxy for evolutionary rate estimates and is informative with respect to prediction of functional properties of proteins. In particular, variability is a strong predictor of gene essentiality for the respective organisms and indicative of sub- or neofunctionalization of paralogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13062-022-00337-7.

Kaiso Regulates DNA Methylation Homeostasis

Gain and loss of DNA methylation in cells is a dynamic process that tends to achieve an equilibrium. Many factors are involved in maintaining the balance between DNA methylation and demethylation. Previously, it was shown that methyl-DNA protein Kaiso may attract NCoR, SMRT repressive complexes affecting histone modifications. On the other hand, the deficiency of Kaiso resulted in reduced methylation of ICR in H19/Igf2 locus and Oct4 promoter in mouse embryonic fibroblasts. However, nothing is known about how Kaiso influences DNA methylation at the genome level. Here we show that deficiency of Kaiso led to whole-genome hypermethylation, using Kaiso deficient human renal cancer cell line obtained via CRISPR/CAS9 genome editing. However, Kaiso serves to protect genic regions, enhancers, and regions with a low level of histone modifications from demethylation. We detected hypomethylation of binding sites for Oct4 and Nanog in Kaiso deficient cells. Kaiso immunoprecipitated with de novo DNA methyltransferases DNMT3a/3b, but not with maintenance methyltransferase DNMT1. Thus, Kaiso may attract methyltransferases to surrounding regions and modulate genome methylation in renal cancer cells apart from being methyl DNA binding protein.

Mutation-selection balance and compensatory mechanisms in tumour evolution

Intratumour heterogeneity and phenotypic plasticity, sustained by a range of somatic aberrations, as well as epigenetic and metabolic adaptations, are the principal mechanisms that enable cancers to resist treatment and survive under environmental stress. A comprehensive picture of the interplay between different somatic aberrations, from point mutations to whole-genome duplications, in tumour initiation and progression is lacking. We posit that different genomic aberrations generally exhibit a temporal order, shaped by a balance between the levels of mutations and selective pressures. Repeat instability emerges first, followed by larger aberrations, with compensatory effects leading to robust tumour fitness maintained throughout the tumour progression. A better understanding of the interplay between genetic aberrations, the microenvironment, and epigenetic and metabolic cellular states is essential for early detection and prevention of cancer as well as development of efficient therapeutic strategies.

Amino acid homorepeats in proteins

Amino acid homorepeats, or homorepeats, are polypeptide segments found in proteins that contain stretches of identical amino acid residues. Although abnormal homorepeat expansions are linked to pathologies such as neurodegenerative diseases, homorepeats are prevalent in eukaryotic proteomes, suggesting that they are important for normal physiology. In this Review, we discuss recent advances in our understanding of the biological functions of homorepeats, which range from facilitating subcellular protein localization to mediating interactions between proteins across diverse cellular pathways. We explore how the functional diversity of homorepeat-containing proteins could be linked to the ability of homorepeats to adopt different structural conformations, an ability influenced by repeat composition, repeat length and the nature of flanking sequences. We conclude by highlighting how an understanding of homorepeats will help us better characterize and develop therapeutics against the human diseases to which they contribute.

FOLR1 was up-regulated in cervical squamous cell carcinoma and correlated with the patients’ progression free survival

Objective The aim of the present work was to evaluate the folate-receptor 1 (FOLR1) expression in cervical squamous cell carcinoma and its clinical significance.

Methods FOLR1 mRNA expression level was detected in the cancer genome atlas (TCGA) database for multiple carcinomas. The FOLR1 mRNA relative expression between tumor tissue and normal cervix tissue of the cervical squamous cell cancer patients was compared by the online data analysis tool of GEPIA. The overall survival (OS) and progression free survival (PFS) between the FOLR1 high and low expression groups were compared by the log-rank test. Thirty one cervical squamous cancer patients and 20 healthy controls were included in and tested for serum FOLR1 protein level detection. Eighty one cervical squamous cell cancer patients who received surgery were included for FOLR1 protein expression detected by immunohistochemistry assay (IHC). The correlation between FOLR1 protein expression and patients’ clinical features was analyzed.

Results FOLR1 mRNA was up-regulated in tumor tissue compared to corresponding normal cervical tissue of cervical squamous cell carcinoma. Top 20 genes interacted with FOLR1 was identified through the network with the edges of 146. UBXN10 (r=0.668, P<0.01) and GBP6 (r-=0.606, P<0.01) were the top 2 genes that most correlated with FOLR1. The serum level of FR-α (FOLR1 coding protein) were 275.50±83.79 and 161.70±66.62 (ng/L) for the cervical cancer and healthy control subjects respectively with significant statistical difference (P<0.05). Using the serum FR-α as serological marker for cervical cancer detection, the diagnostic sensitivity, specificity and AUC were 80.0% (58.40% to 91.93%), 80.65% (63.72% to 90.81%) and 0.85(95%CI:0.74-0.96), respectively. Immunohistochemical assay indicated that of the 81 cancer tissue samples, 45 (55.6%) was FOLR1 protein positive. FOLR1 protein positive expression rate in FIGO stage Ⅲ/Ⅳ was significant higher than in the stage Ⅰ/Ⅱ with statistical difference (P<0.05). The progression free survival (PFS) was significant different between FOLR1 high and low expression group (HR=2.48, 95%CI:1.1-5.58, P=0.023). However, the overall survival (OS) was not statistical different between the two groups (HR=1.34, 95%CI:0.84-2.15, P=0.22).

Conclusion: FOLR1 was up-regulated in both serum and cancer tissue of cervical squamous cell carcinoma which may act as diagnostic and prognostic maker for cervical squamous cell cancer.

On the wrong DNA track: Molecular mechanisms of repeat-mediated genome instability

Expansions of simple tandem repeats are responsible for almost 50 human diseases, the majority of which are severe, degenerative, and not currently treatable or preventable. In this review, we first describe the molecular mechanisms of repeat-induced toxicity, which is the connecting link between repeat expansions and pathology. We then survey alternative DNA structures that are formed by expandable repeats and review the evidence that formation of these structures is at the core of repeat instability. Next, we describe the consequences of the presence of long structure-forming repeats at the molecular level: somatic and intergenerational instability, fragility, and repeat-induced mutagenesis. We discuss the reasons for gender bias in intergenerational repeat instability and the tissue specificity of somatic repeat instability. We also review the known pathways in which DNA replication, transcription, DNA repair, and chromatin state interact and thereby promote repeat instability. We then discuss possible reasons for the persistence of disease-causing DNA repeats in the genome. We describe evidence suggesting that these repeats are a payoff for the advantages of having abundant simple-sequence repeats for eukaryotic genome function and evolvability. Finally, we discuss two unresolved fundamental questions: (i) why does repeat behavior differ between model systems and human pedigrees, and (ii) can we use current knowledge on repeat instability mechanisms to cure repeat expansion diseases?

Role of telomerase in the tumour microenvironment

Telomeres are specialized genomic structures that protect chromosomal ends to maintain genomic stability. Telomeric length is primarily regulated by the telomerase complex, essentially consisting of an RNA template (TERC), an enzymatic subunit (telomerase reverse transcriptase, TERT). In humans, telomerase activity is repressed during embryonic differentiation and is absent in most somatic cells. However, it is upregulated or reactivated in 80%-90% of the primary tumours in humans. The human TERT (hTERT) plays a pivotal role in cellular immortality and tumourigenesis. However, the molecular mechanisms of telomerase functioning in cancer have not been fully understood beyond the telomere maintenance. Several research groups, including ours, have demonstrated that hTERT possesses vital functions independent of its telomere maintenance, including angiogenesis, inflammation, cancer cell stemness, and epithelial-mesenchymal transformation (EMT). All these telomere-independent activities of hTERT may contribute to the regulation of the dynamics and homeostasis of the tumour microenvironment (TME), thereby promoting tumour growth and development. Cancer progression and metastasis largely depend upon the interactions between cancer cells and their microenvironment. In this review, the involvement of TERT in the tumour microenvironment and the underlying implications in cancer therapeutics have been summarized.