Molecular variation between geographical populations and subtypes indicate potential genomic heterogeneity and novel genomic features within CCA. Here, we analyze exome-sequencing data of 87 perihilar cholangiocarcinoma (pCCA) and 261 intrahepatic cholangiocarcinoma (iCCA) cases from 3 Asian centers (including 43 pCCAs and 24 iCCAs from our center). iCCA tumours demonstrate a higher tumor mutation burden and copy number alteration burden (CNAB) than pCCA tumours, and high CNAB indicates a poorer pCCA prognosis. We identify 12 significantly mutated genes and 5 focal CNA regions, and demonstrate common mutations in post-transcriptional modification-related potential driver genes METTL14 and RBM10 in pCCA tumours. Finally we demonstrate the tumour-suppressive role of METTL14, a major RNA N6-adenosine methyltransferase (m6A), and illustrate that its loss-of-function mutation R298H may act through m6A modification on potential driver gene MACF1. Our results may be valuable for better understanding of how post-transcriptional modification can affect CCA development, and highlight both similarities and differences between pCCA and iCCA.

Circular RNAs (circRNAs), covalently closed continuous RNA loops, are generated from cognate linear RNAs through back splicing events, and alternative splicing events may generate different circRNA isoforms at the same locus. However, the challenges of reconstruction and quantification of alternatively spliced full-length circRNAs remain unresolved. On the basis of the internal structural characteristics of circRNAs, we developed CircAST, a tool to assemble alternatively spliced circRNA transcripts and estimate their expression by using multiple splice graphs. Simulation studies showed that CircAST correctly assembled the full sequences of circRNAs with a sensitivity of 85.63%–94.32% and a precision of 81.96%–87.55%. By assigning reads to specific isoforms, CircAST quantified the expression of circRNA isoforms with correlation coefficients of 0.85–0.99 between theoretical and estimated values. We evaluated CircAST on an in-house mouse testis RNA-seq dataset with RNase R treatment for enriching circRNAs and identified 380 circRNAs with full-length sequences different from those of their corresponding cognate linear RNAs. RT-PCR and Sanger sequencing analyses validated 32 out of 37 randomly selected isoforms, thus further indicating the good performance of CircAST, especially for isoforms with low abundance. We also applied CircAST to published experimental data and observed substantial diversity in circular transcripts across samples, thus suggesting that circRNA expression is highly regulated. CircAST can be accessed freely at https://github.com/xiaofengsong/CircAST.

Despite continuous expansion of the RNA binding protein (RBP) world, there is a lack of systematic understanding of RBPs in the mammalian testis, which harbors one of the most complex tissue transcriptomes. We adapted RNA interactome capture to mouse male germ cells, building an RBP atlas characterized by multiple layers of dynamics along spermatogenesis. Trapping of RNA–cross-linked peptides showed that the glutamic acid–arginine (ER) patch, a residue-coevolved polyampholytic element present in coiled coils, enhances RNA binding of its host RBPs. Deletion of this element in NONO (non-POU domain-containing octamer-binding protein) led to a defective mitosis-to-meiosis transition due to compromised NONO-RNA interactions. Whole-exome sequencing of over 1000 infertile men revealed a prominent role of RBPs in the human genetic architecture of male infertility and identified risk ER patch variants.

G-quadruplexes (G4) are noncanonical secondary genome structures. Aberrant formation of G4s can impair genome integrity. Investigation of the relationship between G4s and somatic structural variants (SV) in cancers could provide a better understanding of the role of G4 formation in cancer development and progression. In this study, we combined bioinformatic approaches and multiomics data to investigate the connection between G4s and the somatic SVs. Somatic SV breakpoints were significantly enriched in G4 regions, regardless of SV subtypes. This enrichment was only observed in regions demonstrated to form G4s in cells (“active quadruplexes”), rather than in regions with a sequence compatible with G4 formation but without confirmed G4 formation (“potential quadruplexes”). Several genomic features affected the connection between G4s and SVs, with the enrichment being notably strengthened at the boundary of topologically associated domains. Somatic breakpoints were also preferentially associated with G4 regions with earlier replication timing and open chromatin status. In patients with cancer with homologous recombination repair defects, G4s and somatic breakpoints were substantially more strongly associated. Machine learning models were constructed that showed that G4 propensity is a potent feature for predicting the density of SV breakpoints. Altogether, these findings suggest that the G4 structures play a critical role in modulating the production of somatic SVs in cancers. Significance: G-quadruplex structure formation constitutes a critical step in the production of somatic structural variants in cancers, suggesting G-quadruplex structures as potential targets for future cancer prevention and treatment strategies.

Cancer‐testis (CT) genes are a group of genes restrictedly expressed in testis and multiple cancers and can serve as candidate driver genes participating in the development of cancers. Our previous study identified a number of CT genes in nongerm cell tumors, but their expression pattern in testicular germ cell tumor (TGCT), a cancer type characterized by less genomic alterations, remained largely unknown. In this study, we systematically investigated the expression pattern of CT genes in TGCT samples and evaluated the transcriptome difference between TGCT and normal testis tissues, using datasets from the UCSC Xena platform, The Cancer Genome Atlas (TCGA) and the Genotype‐Tissue Expression (GTEx) project. Pathway enrichment analysis and survival analysis were conducted to evaluate the biological function and prognostic effect of expressed CT genes. We identified that 1036 testis‐specific expressed protein‐coding genes and 863 testis‐specific expressed long noncoding RNAs (lncRNAs) were expressed in TGCT samples, including 883 CT protein‐coding genes and 710 CT lncRNAs defined previously. The number of expressed CT genes was significantly higher in seminomas (P = 3.48 × 10⁻¹³) which were characterized by frequent mutations in driver genes (KIT, KRAS and NRAS). In contrast, the number of expressed CT genes showed a moderate negative correlation with the fraction of copy number altered genomes (cor = −0.28, P = 1.20 × 10⁻³). Unlike other cancers, our analysis revealed that 96.16% of the CT genes were down‐regulated in TGCT samples, while CT genes in stem cell maintenance related pathways were up‐regulated. Further survival analysis provided evidence that CT genes could also predict the prognosis of TGCT patients with both disease‐free interval and progression‐free interval as clinical endpoints. Taken together, our study provided a global view of CT genes in TGCT and provided evidence that CT genes played important roles in the progression and maintenance of TGCT.