A pan-cancer analysis reveals nonstop extension mutations causing SMAD4 tumour suppressor degradation

Underlying effect of SMAD4 gene polymorphism on risk prediction of papillary thyroid carcinoma in Chinese population

Objective

This research intends to explore how variations in the SMAD4 gene impact papillary thyroid carcinoma (PTC) among patients in China.

Methods

The rs10502913 and rs12968012 polymorphisms were genotyped in 405 subjects using SNP-scan high-throughput technology. Differential mRNA expression of SMAD4 was analyzed using data from TCGA and GSE33630, and protein level expression differences were analyzed using immunohistochemistry.

Results

The results showed that SMAD4 mRNA expression was lower in thyroid cancer (THCA) tissues than in normal tissues. Immunohistochemical results showed that the expression level of SMAD4 in normal tissue, thyroid papillary carcinoma tissue and poorly differentiated tissue was significantly different. We found that SMAD4 mismatch variants (rs10502913 and rs12968012) were associated with PTC susceptibility. Specifically, the SMAD4-rs10502913 genotypes (GA and AA) showed a notable correlation with a lower likelihood of PTC in comprehensive and segmented studies (genotype GA: OR (95% CI) = 0.270 (0.077-0.950), p = 0.041; genotype AA: OR (95% CI) = 0.103 (0.025-0.416), p = 0.001). We categorized the immunohistochemical results according to genotype and found that rs10502913-GG protein level was expressed at the lowest level, and both GA and AA were higher than GG (GG vs. AA, P < 0.05), and rs12968012-CG protein level was expressed at the lowest level, and both GG and CC were higher than CG (GG vs. CG, P < 0.01).

Conclusion

Two missense variants of SMAD4 (rs10502913 and rs12968012) are associated with reduced risk of papillary thyroid carcinoma, possibly by reducing protein expression leading to susceptibility to papillary thyroid carcinoma.

Constructing gene similarity networks using co-occurrence probabilities

Gene similarity networks play important role in unraveling the intricate associations within diverse cancer types. Conventionally, gauging the similarity between genes has been approached through experimental methodologies involving chemical and molecular analyses, or through the lens of mathematical techniques. However, in our work, we have pioneered a distinctive mathematical framework, one rooted in the co-occurrence of attribute values and single point mutations, thereby establishing a novel approach for quantifying the dissimilarity or similarity among genes. Central to our approach is the recognition of mutations as key players in the evolutionary trajectory of cancer. Anchored in this understanding, our methodology hinges on the consideration of two categorical attributes: mutation type and nucleotide change. These attributes are pivotal, as they encapsulate the critical variations that can precipitate substantial changes in gene behavior and ultimately influence disease progression. Our study takes on the challenge of formulating similarity measures that are intrinsic to genes' categorical data. Taking into account the co-occurrence probability of attribute values within single point mutations, our innovative mathematical approach surpasses the boundaries of conventional methods. We thereby provide a robust and comprehensive means to assess gene similarity and take a significant step forward in refining the tools available for uncovering the subtle yet impactful associations within the complex realm of gene interactions in cancer.

Target-dependent RNA polymerase as universal platform for gene expression control in response to intracellular molecules

Controlling gene expression in response to specific molecules is an essential technique for regulating cellular functions. However, current platforms with transcription and translation regulators have a limited number of detectable molecules to induce gene expression. Here to address these issues, we present a Target-dependent RNA polymerase (TdRNAP) that can induce RNA transcription in response to the intracellular target specifically recognized by single antibody. By substituting the fused antibody, we demonstrate that TdRNAPs respond to a wide variety of molecules, including peptides, proteins, RNA, and small molecules, and produce desired transcripts in human cells. Furthermore, we show that multiple TdRNAPs can construct orthogonal and multilayer genetic circuits. Finally, we apply TdRNAP to achieve cell-specific genome editing that is autonomously triggered by detecting the target gene product. TdRNAP can expand the molecular variety for controlling gene expression and provide the genetic toolbox for bioengineering and future therapeutic applications.

Control of TGFβ signalling by ubiquitination independent function of E3 ubiquitin ligase TRIP12

Transforming growth factor β (TGFβ) pathway is a master regulator of cell proliferation, differentiation, and death. Deregulation of TGFβ signalling is well established in several human diseases including autoimmune disorders and cancer. Thus, understanding molecular pathways governing TGFβ signalling may help better understand the underlying causes of some of those conditions. Here, we show that a HECT domain E3 ubiquitin ligase TRIP12 controls TGFβ signalling in multiple models. Interestingly, TRIP12 control of TGFβ signalling is completely independent of its E3 ubiquitin ligase activity. Instead, TRIP12 recruits SMURF2 to SMAD4, which is most likely responsible for inhibitory monoubiquitination of SMAD4, since SMAD4 monoubiquitination and its interaction with SMURF2 were dramatically downregulated in TRIP12-/- cells. Additionally, genetic inhibition of TRIP12 in human and murine cells leads to robust activation of TGFβ signalling which was rescued by re-introducing wildtype TRIP12 or a catalytically inactive C1959A mutant. Importantly, TRIP12 control of TGFβ signalling is evolutionary conserved. Indeed, genetic inhibition of Drosophila TRIP12 orthologue, ctrip, in gut leads to a reduced number of intestinal stem cells which was compensated by the increase in differentiated enteroendocrine cells. These effects were completely normalised in Drosophila strain where ctrip was co-inhibited together with Drosophila SMAD4 orthologue, Medea. Similarly, in murine 3D intestinal organoids, CRISPR/Cas9 mediated genetic targeting of Trip12 enhances TGFβ mediated proliferation arrest and cell death. Finally, CRISPR/Cas9 mediated genetic targeting of TRIP12 in MDA-MB-231 breast cancer cells enhances the TGFβ induced migratory capacity of these cells which was rescued to the wildtype level by re-introducing wildtype TRIP12. Our work establishes TRIP12 as an evolutionary conserved modulator of TGFβ signalling in health and disease.

Mechanisms of readthrough mitigation reveal principles of GCN1-mediated translational quality control

Readthrough into the 3' untranslated region (3' UTR) of the mRNA results in the production of aberrant proteins. Metazoans efficiently clear readthrough proteins, but the underlying mechanisms remain unknown. Here, we show in Caenorhabditis elegans and mammalian cells that readthrough proteins are targeted by a coupled, two-level quality control pathway involving the BAG6 chaperone complex and the ribosome-collision-sensing protein GCN1. Readthrough proteins with hydrophobic C-terminal extensions (CTEs) are recognized by SGTA-BAG6 and ubiquitylated by RNF126 for proteasomal degradation. Additionally, cotranslational mRNA decay initiated by GCN1 and CCR4/NOT limits the accumulation of readthrough products. Unexpectedly, selective ribosome profiling uncovered a general role of GCN1 in regulating translation dynamics when ribosomes collide at nonoptimal codons, enriched in 3' UTRs, transmembrane proteins, and collagens. GCN1 dysfunction increasingly perturbs these protein classes during aging, resulting in mRNA and proteome imbalance. Our results define GCN1 as a key factor acting during translation in maintaining protein homeostasis.

Noncoding translation mitigation

Translation is pervasive outside of canonical coding regions, occurring in long noncoding RNAs, canonical untranslated regions and introns1-4, especially in ageing4-6, neurodegeneration5,7 and cancer8-10. Notably, the majority of tumour-specific antigens are results of noncoding translation11-13. Although the resulting polypeptides are often nonfunctional, translation of noncoding regions is nonetheless necessary for the birth of new coding sequences14,15. The mechanisms underlying the surveillance of translation in diverse noncoding regions and how escaped polypeptides evolve new functions remain unclear10,16-19. Functional polypeptides derived from annotated noncoding sequences often localize to membranes20,21. Here we integrate massively parallel analyses of more than 10,000 human genomic sequences and millions of random sequences with genome-wide CRISPR screens, accompanied by in-depth genetic and biochemical characterizations. Our results show that the intrinsic nucleotide bias in the noncoding genome and in the genetic code frequently results in polypeptides with a hydrophobic C-terminal tail, which is captured by the ribosome-associated BAG6 membrane protein triage complex for either proteasomal degradation or membrane targeting. By contrast, canonical proteins have evolved to deplete C-terminal hydrophobic residues. Our results reveal a fail-safe mechanism for the surveillance of unwanted translation from diverse noncoding regions and suggest a possible biochemical route for the preferential membrane localization of newly evolved proteins.

Loss of function SMAD4 nonstop mutations in human cancer

BACKGROUND

SMAD4 is a tumour suppressor gene that is mutated in a variety of cancers. SMAD4 nonstop mutations, which convert stop codons to sense codons that extend transcription, have been identified in genomic databases but have not been characterised in human pathology samples. The frequency of SMAD4 nonstop mutations and the consequences of nonstop mutations on SMAD4 protein expression are unknown.

METHODS

We retrospectively analysed our cancer sequencing database of 38,002 tumour specimens and evaluated the spectrum of SMAD4 mutations. SMAD4 protein expression was evaluated by immunohistochemistry in tumours with SMAD4 nonstop mutations.

RESULTS

In total, 1956 SMAD4 mutations were identified in 1822 tumours. SMAD4 mutations were most common in tumours of the gastrointestinal tract and included nonsense variants (n = 344), frameshift indels (n = 258), splice site variants (n = 104), and missense variants at codon R361 (n = 245). In a subset of cases with immunohistochemistry, SMAD4 expression was lost in 23 of 25 tumours (92%) with protein truncating variants and in 7 of 27 tumours (26%) with missense variants. Four cases harboured SMAD4 nonstop mutations. SMAD4 nonstop mutations were identified in two pancreatic adenocarcinomas, one colonic adenocarcinoma, and one non-small cell lung carcinoma. Immunohistochemistry demonstrated loss of SMAD4 protein expression in each of the four tumours with SMAD4 nonstop mutations.

CONCLUSION

SMAD4 nonstop mutations are associated with loss of SMAD4 protein expression in multiple tumour types. SMAD4 nonstop mutations should be clinically interpreted as pathogenic loss of function alterations when identified in cancer sequencing panels.

Chimeric oligonucleotides combining guide RNA and single-stranded DNA repair template effectively induce precision gene editing

The ability to precisely alter the genome holds immense potential for molecular biology, medicine and biotechnology. The development of the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) into a genomic editing tool has vastly simplified genome engineering. Here, we explored the use of chemically synthesized chimeric oligonucleotides encoding a target-specific crRNA (CRISPR RNA) fused to a single-stranded DNA repair template for RNP-mediated precision genome editing. By generating three clinically relevant oncogenic driver mutations, two non-stop extension mutations, an FGFRi resistance mutation and a single nucleotide change, we demonstrate the ability of chimeric oligos to form RNPs and direct Cas9 to effectively induce genome editing. Further, we demonstrate that the polarity of the chimeric oligos is crucial: only chimeric oligos with the single-stranded DNA repair template fused to the 3'-end of the crRNA are functional for accurate editing, while templates fused to the 5'-end are ineffective. We also find that chimeras can perform editing with both symmetric and asymmetric single-stranded DNA repair templates. Depending on the target locus, the editing efficiency using chimeric RNPs is similar to or less than the efficiency of editing using the bipartite standard RNPs. Our results indicate that chimeric RNPs comprising RNA-DNA oligos formed from fusing the crRNA and DNA repair templates can successfully induce precise edits. While chimeric RNPs do not display an advantage over standard RNPs, they nonetheless represent a viable approach for one-molecule precision genome editing.

The HOXB13 variant X285K is associated with clinical significance and early age at diagnosis in African American prostate cancer patients

BACKGROUND

Recently, a novel HOXB13 variant (X285K) was observed in men of African descent with prostate cancer (PCa) in Martinique. Little is known about this or other variants in HOXB13 which may play a role in PCa susceptibility in African-American (AA) men.

METHODS

We sequenced HOXB13 in an AA population of 1048 men undergoing surgical treatment for PCa at Johns Hopkins Hospital.

RESULTS

Seven non-synonymous germline variants were observed in the patient population. While six of these variants were seen only once, X285K was found in eight patients. In a case-case analysis, we find that carriers of this latter variant are at increased risk of clinically significant PCa (1.2% carrier rate in Gleason Score ≥7 PCa vs. 0% in Gleason Score <7 PCa, odds ratio, OR = inf; 95% Confidence Interval, 95%CI:1.05-inf, P = 0.028), as well as PCa with early age at diagnosis (2.4% carrier rate in patients <50 year vs. 0.5% carrier rate in patients ≥50 year, OR = 5.25, 95% CI:1.00-28.52, P = 0.03).

CONCLUSIONS

While this variant is rare in the AA population (~0.2% MAF), its ancestry-specific occurrence and apparent preferential association with risk for the more aggressive disease at an early age emphasizes its translational potential as an important, novel PCa susceptibility marker in the high-risk AA population.

Systematic analysis of migration factors by MigExpress identifies essential cell migration control genes in non-small cell lung cancer

Cell migration is an essential process in health and in disease, including cancer metastasis. A comprehensive inventory of migration factors is nonetheless lacking-in part due to the difficulty in assessing migration using high-throughput technologies. Hence, there are currently very few screens that systematically reveal factors controlling cell migration. Here, we introduce MigExpress as a platform for the 'identification of Migration control genes by differential Expression'. MigExpress exploits the combination of in-depth molecular profiling and the robust quantitative analysis of migration capacity in a broad panel of samples and identifies migration-associated genes by their differential expression in slow- versus fast-migrating cells. We applied MigExpress to investigate non-small cell lung cancer (NSCLC), which is the most frequent cause of cancer mortality mainly due to metastasis. In 54 NSCLC cell lines, we comprehensively determined mRNA and protein expression. Correlating the transcriptome and proteome profiles with the quantified migration properties led to the discovery and validation of FLNC, DSE, CPA4, TUBB6, and BICC1 as migration control factors in NSCLC cells, which were also negatively correlated with patient survival. Notably, FLNC was the least expressed filamin in NSCLC, but the only one controlling cell migration and correlating with patient survival and metastatic disease stage. In our study, we present MigExpress as a new method for the systematic analysis of migration factors and provide a comprehensive resource of transcriptomic and proteomic data of NSCLC cell lines related to cell migration.

Identification of novel protein-coding sequences in Eucalyptus grandis plants by high-resolution mass spectrometry

Eucalyptus species are widely used in the forestry industry, and a significant increase in the number of sequences available in database repositories has been observed for these species. In proteomics, a protein is identified by correlating the theoretical fragmentation spectrum derived from genomic/transcriptomic data against the experimental fragmentation mass spectrum acquired from large-scale analysis of protein mixtures. Proteogenomics is an alternative approach that can identify novel proteins encoded by regions previously considered as non-coding. This study aimed to confidently identify and confirm the existence of previously unknown protein-coding sequences in the Eucalyptus grandis genome. To this end, we used a modified spectral correlation strategy and a dedicated de novo peptide sequencing pipeline. Upon the strategy used here, we confidently identified 41 novel peptide forms and six peptides containing at least one single amino acid substitution. The most representative genomic class of novel peptides was identified as originating from alternative reading frames. In contrast, no clear single amino acid substitution pattern was identified. Validation of the identifications was carried out using a parallel reaction monitoring approach that provided further mass spectrometry support for the existence of the novel peptide sequences. Data are available via ProteomeXchange with identifier PXD022110.