Non-canonical Motif Research Articles

Design and Ribosomal Incorporation of Noncanonical Disulfide-Directing Motifs for the Development of Multicyclic Peptide Libraries.

The engineering of naturally occurring disulfide-rich peptides (DRPs) has been significantly hampered by the difficulty of manipulating disulfide pairing. New DRPs that take advantage of fold-directing motifs and noncanonical thiol-bearing amino acids are easy-to-fold with expected disulfide connectivities, representing a new class of scaffolds for the development of peptide ligands and therapeutics. However, the limited diversity of the scaffolds and particularly the use of noncanonical amino acids [e.g., penicillamine (Pen)] that are difficult to be translated by ribosomes greatly hamper the further development and application of these DRPs. Here, we designed and synthesized noncanonical bisthiol motifs bearing sterically obstructed thiol groups analogous to the Pen thiol to direct the folding of peptides into specific bicyclic and tricyclic structures. These bisthiol motifs can be ribosomally incorporated into peptides through a commercially available PURE system integrated with genetic code reprograming, which enables, for the first time, the in vitro expression of bicyclic peptides with two noncanonical and orthogonal disulfide bonds. We further constructed a bicyclic peptide library encoded by mRNA, with which new bicyclic peptide ligands with nanomolar affinity to proteins were successfully selected. Therefore, this study provides a new, general, and robust method for discovering de novo DRPs with new structures and functions not derived from natural peptides, which would greatly benefit the field of peptide drug discovery.

Comprehensive genetic diagnosis of tandem repeat expansion disorders with programmable targeted nanopore sequencing.

More than 50 neurological and neuromuscular diseases are caused by short tandem repeat (STR) expansions, with 37 different genes implicated to date. We describe the use of programmable targeted long-read sequencing with Oxford Nanopore’s ReadUntil function for parallel genotyping of all known neuropathogenic STRs in a single assay. Our approach enables accurate, haplotype-resolved assembly and DNA methylation profiling of STR sites, from a list of predetermined candidates. This correctly diagnoses all individuals in a small cohort (n = 37) including patients with various neurogenetic diseases (n = 25). Targeted long-read sequencing solves large and complex STR expansions that confound established molecular tests and short-read sequencing and identifies noncanonical STR motif conformations and internal sequence interruptions. We observe a diversity of STR alleles of known and unknown pathogenicity, suggesting that long-read sequencing will redefine the genetic landscape of repeat disorders. Last, we show how the inclusion of pharmacogenomic genes as secondary ReadUntil targets can further inform patient care.

Hidden Information Revealed Using the Orthogonal System of Nucleic Acids.

In this study, the organization of genetic information in nucleic acids is defined using a novel orthogonal representation. Clearly defined base pairing in DNA allows the linear base chain and sequence to be mathematically transformed into an orthogonal representation where the G–C and A–T pairs are displayed in different planes that are perpendicular to each other. This form of base allocation enables the evaluation of any nucleic acid and predicts the likelihood of a particular region to form non-canonical motifs. The G4Hunter algorithm is currently a popular method of identifying G-quadruplex forming sequences in nucleic acids, and offers promising scores despite its lack of a substantial rational basis. The orthogonal representation described here is an effort to address this incongruity. In addition, the orthogonal display facilitates the search for other sequences that are capable of adopting non-canonical motifs, such as direct and palindromic repeats. The technique can also be used for various RNAs, including any aptamers. This powerful tool based on an orthogonal system offers considerable potential for a wide range of applications.

GATA3 Truncation Mutants Alter EMT Related Gene Expression via Partial Motif Recognition in Luminal Breast Cancer Cells.

GATA3 is known to be one of the most frequently mutated genes in breast cancer. More than 10% of breast tumors carry mutations in this gene. However, the functional consequence of GATA3 mutations is still largely unknown. Clinical data suggest that different types of GATA3 mutations may have distinct roles in breast cancer characterization. In this study, we have established three luminal breast cancer cell lines that stably express different truncation mutants (X308 splice site deletion, C321 frameshift, and A333 frameshift mutants) found in breast cancer patients. Transcriptome analysis identified common and distinct gene expression patterns in these GATA3 mutant cell lines. In particular, the impacts on epithelial-to-mesenchymal transition (EMT) related genes are similar across these mutant cell lines. Chromatin localization of the mutants is highly overlapped and exhibits non-canonical motif enrichment. Interestingly, the A333 frameshift mutant expressed cells displayed the most significant impact on the GATA3 binding compared to X308 splice site deletion and C321fs mutants expressed cells. Our results suggest the common and different roles of GATA3 truncation mutations during luminal breast cancer development.

Porphyrin induced structural destabilization of a parallel DNA G-quadruplex in human MRP1 gene promoter.

Porphyrins are among the first ligands that have been tested for their quadruplex binding and stabilization potential. We report the differential interaction of the positional cationic porphyrin isomers TMPyP3 and TMPyP4 with a parallel G-quadruplex (GQ) formed by 33-mer (TP) regulatory sequence present in the promoter region of the human multidrug resistance protein 1 (MRP1) transporter gene. This GQ element encompasses the three evolutionary conserved SP1 transcription factor binding sites. Taking into account that SP1 binds to a non-canonical GQ motif with higher affinity than to a canonical duplex DNA consensus motif, it is suggestive that GQ distortion by cationic porphyrin will have important implications in the regulation of MRP1 expression. Herein, we employed biophysical analysis using circular dichroism, visible absorption, UV-thermal melting and steady-state fluorescence spectroscopy, reporting destabilization of MRP1 GQ by cationic porphyrins. Results suggest that TMPyP4 and TMPyP3 interact with GQ with a binding affinity of 106 to 107 M-1 . Thermodynamic analysis indicated a significant decrease in melting temperature of GQ (ΔTm of 15.5°C-23.5°C), in the presence of 2 times excess of porphyrins. This study provides the biophysical evidence indicating the destabilisation of a parallel DNA G-quadruplex by cationic porphyrins.

Screening of Minimalist Noncanonical Sites in Duplex DNA and RNA Reveals Context and Motif-Selective Binding by Fluorogenic Base Probes.

We hypothesize that programmable hybridization to noncanonical nucleic acid motifs may be achieved by macromolecular display of binders to individual noncanonical pairs (NCPs). As each recognition element may individually have weak binding to an NCP, we developed a semi-rational approach to detect low affinity interactions between selected nitrogenous bases and noncanonical sites in duplex DNA and RNA. A set of fluorogenic probes was synthesized by coupling abiotic (triazines, pyrimidines) and native RNA bases to thiazole orange (TO) dye. This probe library was screened against duplex nucleic acid substrates bearing single abasic, single NCP, and tandem NCP sites. Probe engagement with NCP sites was reported by 100-1000× fluorescence enhancement over background. Binding is strongly context-dependent, reflective of both molecular recognition and stability: less stable motifs are more likely to bind a synthetic probe. Further, DNA and RNA substrates exhibit entirely different abasic and single NCP binding profiles. While probe binding in the abasic and single NCP screens was monotonous, much richer binding profiles were observed with the screen of tandem NCP sites in RNA, in part due to increased steric accessibility. In addition to known binding interactions between the triazine melamine (M) and T/U sites, the NCP screens identified new targeting elements for pyrimidine-rich motifs in single NCPs and 2×2 internal bulges. We anticipate that semi-rational approaches of this type will lead to programmable noncanonical hybridization strategies at the macromolecular level.

G-quadruplex RNA motifs influence gene expression in the malaria parasite Plasmodium falciparum.

G-quadruplexes are non-helical secondary structures that can fold in vivo in both DNA and RNA. In human cells, they can influence replication, transcription and telomere maintenance in DNA, or translation, transcript processing and stability of RNA. We have previously showed that G-quadruplexes are detectable in the DNA of the malaria parasite Plasmodium falciparum, despite a very highly A/T-biased genome with unusually few guanine-rich sequences. Here, we show that RNA G-quadruplexes can also form in P. falciparum RNA, using rG4-seq for transcriptome-wide structure-specific RNA probing. Many of the motifs, detected here via the rG4seeker pipeline, have non-canonical forms and would not be predicted by standard in silico algorithms. However, in vitro biophysical assays verified formation of non-canonical motifs. The G-quadruplexes in the P. falciparum transcriptome are frequently clustered in certain genes and associated with regions encoding low-complexity peptide repeats. They are overrepresented in particular classes of genes, notably those that encode PfEMP1 virulence factors, stress response genes and DNA binding proteins. In vitro translation experiments and in vivo measures of translation efficiency showed that G-quadruplexes can influence the translation of P. falciparum mRNAs. Thus, the G-quadruplex is a novel player in post-transcriptional regulation of gene expression in this major human pathogen.

The polarization dependence of 2D IR cross-peaks distinguishes parallel-stranded and antiparallel-stranded DNA G-quadruplexes

Guanine-rich nucleic acid sequences have a tendency to form four-stranded non-canonical motifs known as G-quadruplexes. These motifs may adopt a wide range of structures characterized by size, strand orientation, guanine base conformation, and fold topology. Using three K+-bound model systems, we show that vibrational coupling between guanine C6 = O and ring modes varies between parallel-stranded and antiparallel-stranded G-quadruplexes, and that such structures can be distinguished by comparison of the polarization dependences of cross-peaks in their two-dimensional infrared (2D IR) spectra. Combined with previously defined vibrational frequency trends, this analysis reveals key features of a 30-nucleotide unimolecular variant of the Bcl-2 proximal promoter that are consistent with its reported structure. This study shows that 2D IR spectroscopy is a convenient method for analyzing G-quadruplex structures that can be applied to complex sequences where traditional high-resolution methods are limited by solubility and disorder.

RiboDraw: semiautomated two-dimensional drawing of RNA tertiary structure diagrams

Publishing, discussing, envisioning, modeling, designing and experimentally determining RNA three-dimensional (3D) structures involve preparation of two-dimensional (2D) drawings that depict critical functional features of the subject molecules, such as noncanonical base pairs and protein contacts. Here, we describe RiboDraw, new software for crafting these drawings. We illustrate the features of RiboDraw by applying it to several RNAs, including the Escherichia coli tRNA-Phe, the P4–P6 domain of Tetrahymena ribozyme, a −1 ribosomal frameshift stimulation element from beet western yellows virus and the 5′ untranslated region of SARS-CoV-2. We show secondary structure diagrams of the 23S and 16S subunits of the E. coli ribosome that reflect noncanonical base pairs, ribosomal proteins and structural motifs, and that convey the relative positions of these critical features in 3D space. This software is a MATLAB package freely available at https://github.com/DasLab/RiboDraw.

Z-DNA as a Touchstone for Additive Empirical Force Fields and a Refinement of the Alpha/Gamma DNA Torsions for AMBER.

Although current AMBER force fields are relatively accurate for canonical B-DNA, many noncanonical structures are still described incorrectly. As noncanonical motifs are attracting increasing attention due to the role they play in living organisms, further improvement is desirable. Here, we have chosen the Z-DNA molecule, which can be considered a touchstone of the universality of empirical force fields, since the noncanonical α and γ backbone conformations native to Z-DNA are also found in protein-DNA complexes, i-motif DNA, and other noncanonical DNAs. We show that spurious α/γ conformations occurring in simulations with current AMBER force fields, OL15 and bsc1, are largely due to inaccurate α/γ parametrization. Moreover, stabilization of native Z-DNA substates involving γ = trans conformations appears to be in conflict with the correct description of the canonical B-DNA structure. Because the balance of the native and spurious conformations is influenced by nonadditive effects, this is a difficult case for an additive dihedral energy scheme such as AMBER. We propose new α/γ parameters, denoted OL21, and show that they improve the stability of native α/γ Z-DNA substates while keeping the canonical DNA description virtually unchanged, thus representing a reasonable compromise within the additive force field framework. Although further extensive testing is needed, the new modification appears to be a promising step toward a more reliable description of noncanonical DNA motifs and provides the best performance for Z-DNA molecules among current AMBER force fields.

Engineered DNase-inactive Cpf1 variants to improve targeting scope for base editing in E. coli

The development of base editing (BE) technology has opened a new avenue for research studies in bacteriology, particularly for bacterial species in which the DNA double-strand breaks (DSBs) introduced by CRISPR/Cas system would lead to cell death. However, a major limitation of BE-mediated gene editing is the restricted editable sites in the target bacterial genome due to highly diverse genomic compositions, such as GC content. Herein, we developed a broad-spectrum DNase-inactive Cpf1 (dCpf1) variant from Francisella novicida (bsdFnCpf1) through directed evolution. The resulting optimized mutant showed a substantially expanded targeting range, including previously non-canonical protospacer-adjacent motifs (PAMs), especially the GC-rich PAMs. Cytidine deaminase APOBEC1 and uracil DNA glycosylase inhibitor (UGI) were fused with bsdFnCpf1 to achieve specific C to T mutations at multiple target sites with canonical or non-canonical PAMs in the E. coli genome without compromising cell growth. We anticipate that bsdFnCpf1 could be applied for multiplex gene regulation and BE in species that have been reported to be suitable for Cpf1.

Aminovinyl Cysteine Containing Peptides: A Unique Motif That Imparts Key Biological Activity.

Natural products that contain distinctive chemical functionality can serve as useful starting points to develop Nature’s compounds into viable therapeutics. Peptide natural products, an under-represented class of medicines, such as ribosomally synthesized and post-translationally modified peptides (RiPPs), often contain noncanonical amino acids and structural motifs that give rise to potent biological activity. However, these motifs can be difficult to obtain synthetically, thereby limiting the transition of RiPPs to the clinic. Aminovinyl cysteine containing peptides, which display potent antimicrobial or anticancer activity, possess an intricate C-terminal ring that is critical for bioactivity. To date, successful methods for the total chemical synthesis of such peptides are yet to be realized, although several advancements have been achieved. In this perspective, we review this burgeoning class of aminovinyl cysteine peptides and critically evaluate the chemical strategies to install the distinct aminovinyl cysteine motif.

Systematic identification of non-canonical transcription factor motifs

Sequence-specific transcription factors (TFs) recognize motifs of related nucleotide sequences at their DNA binding sites. Upon binding at these sites, TFs regulate critical molecular processes such as gene expression. It is widely assumed that a TF recognizes a single “canonical” motif, although recent studies have identified additional “non-canonical” motifs for some TFs. A comprehensive approach to identify non-canonical DNA binding motifs and the functional importance of those motifs’ matches in the human genome is necessary for fully understanding the mechanisms of TF-regulated molecular processes in human cells. To address this need, we developed a statistical pipeline for in vitro HT-SELEX data that identifies and characterizes the distributions of non-canonical TF motifs in a stringent manner. Analyzing ~170 human TFs’ HT-SELEX data, we found non-canonical motifs for 19 TFs (11%). These non-canonical motifs occur independently of the TFs’ canonical motifs. Non-canonical motif occurrences in the human genome show similar evolutionary conservation to canonical motif occurrences, explain TF binding in locations without canonical motifs, and occur within gene promoters and epigenetically marked regulatory sequences in human cell lines and tissues. Our approach and collection of non-canonical motifs expand current understanding of functionally relevant DNA binding sites for human TFs.

PUL21 is a viral phosphatase adaptor that promotes herpes simplex virus replication and spread

The herpes simplex virus (HSV)-1 protein pUL21 is essential for efficient virus replication and dissemination. While pUL21 has been shown to promote multiple steps of virus assembly and spread, the molecular basis of its function remained unclear. Here we identify that pUL21 is a virus-encoded adaptor of protein phosphatase 1 (PP1). pUL21 directs the dephosphorylation of cellular and virus proteins, including components of the viral nuclear egress complex, and we define a conserved non-canonical linear motif in pUL21 that is essential for PP1 recruitment. In vitro evolution experiments reveal that pUL21 antagonises the activity of the virus-encoded kinase pUS3, with growth and spread of pUL21 PP1-binding mutant viruses being restored in adapted strains where pUS3 activity is disrupted. This study shows that virus-directed phosphatase activity is essential for efficient herpesvirus assembly and spread, highlighting the fine balance between kinase and phosphatase activity required for optimal virus replication.

Machine Learning Aids Classification and Discrimination of Noncanonical DNA Folding Motifs by an Arrayed Host:Guest Sensing System.

An arrayed host:guest fluorescence sensor system can discriminate among and classify multiple different noncanonical DNA structures by exploiting selective molecular recognition. The sensor is highly selective and can discriminate between folds as similar as native G-quadruplexes and those with bulges or vacancies. The host and guest can form heteroternary complexes with DNA strands, with the host acting as mediator between the DNA and dye, modulating the emission. By applying machine learning algorithms to the sensing data, prediction of the folding state of unknown DNA strands is possible with high fidelity.

Rhodamine 6G-Ligand Influencing G-Quadruplex Stability and Topology.

The involvement of G-quadruplex (G4) structures in nucleic acids in various molecular processes in cells such as replication, gene-pausing, the expression of crucial cancer-related genes and DNA damage repair is well known. The compounds targeting G4 usually bind directly to the G4 structure, but some ligands can also facilitate the G4 folding of unfolded G-rich sequences and stabilize them even without the presence of monovalent ions such as sodium or potassium. Interestingly, some G4-ligand complexes can show a clear induced CD signal, a feature which is indirect proof of the ligand interaction. Based on the dichroic spectral profile it is not only possible to confirm the presence of a G4 structure but also to determine its topology. In this study we examine the potential of the commercially available Rhodamine 6G (RhG) as a G4 ligand. RhG tends to convert antiparallel G4 structures to parallel forms in a manner similar to that of Thiazole Orange. Our results confirm the very high selectivity of this ligand to the G4 structure. Moreover, the parallel topology of G4 can be verified unambiguously based on the specific induced CD profile of the G4-RhG complex. This feature has been verified on more than 50 different DNA sequences forming various non-canonical structural motifs.

The Mutation Landscape of Follicular Lymphoma and Chronic Lymphocytic Leukemia Reveals a Distinctive Frequency and Composition of Aid-Related Mutation Signatures

Introduction: Somatic mutations associated with distinctive types of cancer can be globally classified by the trinucleotide context of single nucleotide variants into mutational signatures. Distinct signatures can be attributed to underlying causative mechanisms such as UV irradiation or smoking. B-cell lymphomas have only been analyzed globally and carry a signature characterized by APOBEC activity (Alexandrov, 2014).The two most common indolent B-cell neoplasms, chronic lymphocytic leukemia (CLL) and follicular lymphoma (FL), are characterized by distinctive pathogenic processes, but mutagenesis caused by the APOBEC enzyme activation-induced cytosine deaminase (AID) may be shared by both malignancies. AID deaminates deoxycytidine to deoxyuridine, which is processed by DNA repair enzymes into a larger spectrum of point mutations and also DNA double-strand breaks. AID activity and subsequent DNA repair create two different mutation signatures: the canonical signature characterized by C>T/G transitions at RCY motifs and the non-canonical defined by A>C transversions at WAN motifs. More recently a third AID signature characterized by C>T transitions at RCG motifs has been associated with AID-mediated CpG-methylation dependent mutagenesis.We here investigate the hypothesis that AID mutational signatures contribute to the pathogenesis of FL, a lymphoma characterized by constitutional expression of AID, but also of CLL since high AID expression levels may account for a more aggressive disease, and a subset of CLL B-cells with AID expression display dissociation between class switch recombination and somatic hypermutation.Methods: We analyzed the mutational landscape in whole exome sequencing data from 31 FL and 24 CLL cases. Somatic variants were called by comparison of tumor and germline DNA with an in-house developed pipeline. The relative contribution of mutations at canonical and non-canonical AID hotspot motifs was computed by direct quantification of variants occurring in known AID motifs. Student's t test was used for comparison between CLL and FL.Mutational signatures were defined according to the 96-base substitution model (Alexandrov, 2013) by an unsupervised machine learning approach, with a local implementation of the SomaticSignatures R package.Results: The frequency of C>T/G transitions at canonical AID hotspot was more prevalent in CLL (9.0%; 95% CI: 8.1-10.0%) than in FL (6.6%, 95% CI: 5.9-7.3%)(p<0.001). Likewise, A>C transversions at non-canonical AID sites were significantly more frequent in CLL (1.9%; 95% CI: 1.6-2.2%) as compared to FL (0.5%; 95% CI: 0.4-0.6%)(p<0.0001).In the unsupervised analysis of FL and CLL cases, 97.7% of the mutation spectrum variance was explained by only 5 signatures (Figure 1). In FL, the mutational landscape was dominated by a signature (S2) characterized by C>A transversions in the trinucleotide KCH context (40%, 95% CI: 33-36%). The second most frequent signature (S3; 22%, 95% CI: 15-30%) was characterized by C>T transitions in the context of a hybrid between the canonical AID and the CpG hotspot motifs (RCY/RCG).In contrast, the mutational landscape of CLL was dominated by a signature (S5; 42%, 95% CI: 34-49%) characterized by C>T transitions in the RCY motif and A>C transversions in the WAN context, consistent with an overlapping between the canonical and the non-canonical AID hotspots (Kasar, 2016).Conclusions: The mutational landscape of both CLL and FL seems to be strongly shaped by AID activity. In FL, AID-induced mutations are mainly restricted to canonical AID hotspots and CpG-methylation dependent mutagenesis sites. In strong contrast, both canonical and non-canonical AID motifs contributed to the mutational landscape of CLL. In addition, FL carries a high prevalence of C>A substitutions that are associated with exposure to polycyclic aromatic hydrocarbons that produce bulky DNA adducts followed by the nucleotide excision repair process. These compounds are present in several industrial processes (food products, beverages, tobacco) that have been suggested as occupational risk factors for follicular lymphomas (Mester, 2009).Overall, our data suggest a much stronger role of AID in CLL pathogenesis than hitherto appreciated. The differential contribution of canonical versus non-canonical AID activity points to different regulation of AID in CLL compared to FL and warrants further investigation. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Haplotype diversity and sequence heterogeneity of human telomeres.

Telomeres are regions of repetitive nucleotide sequences capping the ends of eukaryotic chromosomes that protect against deterioration, and whose lengths can be correlated with age and adverse health risk factors. Yet, given their length and repetitive nature, telomeric regions are not easily reconstructed from short-read sequencing, thus making telomere sequencing, mapping, and variant resolution challenging problems. Recently, long-read sequencing, with read lengths measuring in hundreds of kilobase pairs, has made it possible to routinely read into telomeric regions and inspect their sequence structure. Here, we describe a framework for extracting telomeric reads from whole-genome single-molecule sequencing experiments, including de novo identification of telomere repeat motifs and repeat types, and also describe their sequence variation. We find that long, complex telomeric stretches and repeats can be accurately captured with long-read sequencing, observe extensive sequence heterogeneity of human telomeres, discover and localize noncanonical telomere sequence motifs (both previously reported, as well as novel), and validate them in short-read sequence data. These data reveal extensive intra- and inter-population diversity of repeats in telomeric haplotypes, reveal higher paternal inheritance of telomeric variants, and represent the first motif composition maps of multi-kilobase-pair human telomeric haplotypes across three distinct ancestries (Ashkenazi, Chinese, and Utah), which can aid in future studies of genetic variation, aging, and genome biology.

Two distinct motifs for Zic-r.a drive specific gene expression in two cell lineages.

Zic-r.a, a maternal transcription factor, specifies posterior fate in ascidian embryos. However, its direct target, Tbx6-r.b, does not contain typical Zic-r.a-binding sites in its regulatory region. Using an in vitro selection assay, we found that Zic-r.a binds to sites dissimilar to the canonical motif, by which it activates Tbx6-r.b in a sub-lineage of muscle cells. These sites with non-canonical motifs have weak affinity for Zic-r.a; therefore, it activates Tbx6-r.b only in cells expressing Zic-r.a abundantly. Meanwhile, we found that Zic-r.a expressed zygotically in late embryos activates neural genes through canonical sites. Because different zinc-finger domains of Zic-r.a are important for driving reporters with canonical and non-canonical sites, it is likely that the non-canonical motif is not a divergent version of the canonical motif. In other words, our data indicate that the non-canonical motif represents a motif distinct from the canonical motif. Thus, Zic-r.a recognizes two distinct motifs to activate two sets of genes at two timepoints in development. This article has an associated 'The people behind the papers' interview.

Role of Regulatory Non-Coding RNAs in Aggressive Thyroid Cancer: Prospective Applications of Neural Network Analysis.

Thyroid cancer (TC) is the most common endocrine malignancy. Most TCs have a favorable prognosis, whereas anaplastic thyroid carcinoma (ATC) is a lethal form of cancer. Different genetic and epigenetic alterations have been identified in aggressive forms of TC such as ATC. Non-coding RNAs (ncRNAs) represent functional regulatory molecules that control chromatin reprogramming, including transcriptional and post-transcriptional mechanisms. Intriguingly, they also play an important role as coordinators of complex gene regulatory networks (GRNs) in cancer. GRN analysis can model molecular regulation in different species. Neural networks are robust computing systems for learning and modeling the dynamics or dependencies between genes, and are used for the reconstruction of large data sets. Canonical network motifs are coordinated by ncRNAs through gene production from each transcript as well as through the generation of a single transcript that gives rise to multiple functional products by post-transcriptional modifications. In non-canonical network motifs, ncRNAs interact through binding to proteins and/or protein complexes and regulate their functions. This article overviews the potential role of ncRNAs GRNs in TC. It also suggests prospective applications of deep neural network analysis to predict ncRNA molecular language for early detection and to determine the prognosis of TC. Validation of these analyses may help in the design of more effective and precise targeted therapies against aggressive TC.