RNA-binding motif protein RBM47 promotes tumorigenesis in nasopharyngeal carcinoma through multiple pathways.

RNA binding motif proteins (RBMs) have been widely implicated in the tumorigenesis of multiple human cancers but scarcely studied in nasopharyngeal carcinoma (NPC). Here, we compare the mRNA levels of 29 RBMs between 87 NPC and 10 control samples. We find that RBM47 is frequently upregulated in NPC specimens, and its high expression is associated with the poor prognosis of patients with NPC. Biological experiments show that RBM47 plays an oncogenic role in NPC cells. Mechanically, RBM47 binds to the promoter and regulates the transcription of BCAT1, and its overexpression partially rescues the inhibitory effects of RBM47-knockdown on NPC cells. Moreover, transcriptome analysis reveals that RBM47 regulates alternative splicing of pre-mRNA, including those cancer-related, to a large extent in NPC cells. Furthermore, RBM47 binds to hnRNPM and cooperatively regulates multiple splicing events in NPC cells. In addition, we find that knockdown of hnRNPM inhibits proliferation and migration of NPC cells. Our study, taken together, shows that RBM47 promotes the progression of NPC through multiple pathways, acting as a transcriptional factor and a modulator of alternative splicing in cooperation with hnRNPM. Our study also highlights that RBM47 and hnRNPM could be prognostic factors and potential therapeutic targets for NPC.

RBFOX2/GOLIM4 Splicing Axis Activates Vesicular Transport Pathway to Promote Nasopharyngeal Carcinogenesis.

20-30% of patients with nasopharyngeal carcinoma (NPC) develop distant metastasis or recurrence leading to poor survival, of which the underlying key molecular events have yet to be addressed. Here alternative splicing events in 85 NPC samples are profiled using transcriptome analysis and it is revealed that the long isoform of GOLIM4 (-L) with exon-7 is highly expressed in NPC and associated with poor prognosis. Lines of evidence demonstrate the pro-tumorigenic function of GOLIM4-L in NPC cells. It is further revealed that RBFOX2 binds to a GGAA motif in exon-7 and promotes its inclusion forming GOLIM4-L. RBFOX2 knockdown suppresses the tumorigenesis of NPC cells, phenocopying GOLIM4-L knockdown, which is significantly rescued by GOLIM4-L overexpression. High expression of RBFOX2 is correlated with the exon-7 inclusion of GOLIM4 in NPC biopsies and associated with worse prognosis. It is observed that RBFOX2 and GOLIM4 can influence vesicle-mediated transport through maintaining the organization of Golgi apparatus. Finally, it is revealed that RAB26 interacts with GOLIM4 and mediates its tumorigenic potentials in NPC cells. Taken together, the findings provide insights into how alternative splicing contributes to NPC development, by highlighting a functional link between GOLIM4-L and its splicing regulator RBFOX2 activating vesicle-mediated transport involving RAB26.

Large-scale analysis of 2,152 Ig-seq datasets reveals key features of B cell biology and the antibody repertoire.

Antibody repertoire sequencing enables researchers to acquire millions of B cell receptors and investigate these molecules at the single-nucleotide level. This power and resolution in studying humoral responses have led to its wide applications. However, most of these studies were conducted with a limited number of samples. Given the extraordinary diversity, assessment of these key features with a large sample set is demanded. Thus, we collect and systematically analyze 2,152 high-quality heavy-chain antibody repertoires. Our study reveals that 52 core variable genes universally contribute to more than 99% of each individual's repertoire; a distal interspersed preferences characterize V gene recombination; the number of public clones between two repertoires follows a linear model, and the positive selection dominates at RGYW motif in somatic hypermutations. Thus, this population-level analysis resolves some critical features of the antibody repertoire and may have significant value to the large cadre of scientists.

Germline Polymorphisms and Length of Survival of Nasopharyngeal Carcinoma: An Exome-Wide Association Study in Multiple Cohorts.

Germline polymorphisms are linked with differential survival outcomes in cancers but are not well studied in nasopharyngeal carcinoma (NPC). Here, a two-phase association study is conducted to discover germline polymorphisms that are associated with the prognosis of NPC. The discovery phase includes two consecutive hospital cohorts of patients with NPC from Southern China. Exome-wide genotypes at 246 173 single nucleotide polymorphisms (SNPs) are determined, followed by survival analysis for each SNP under Cox proportional hazard regression model. Candidate SNP is replicated in another two independent cohorts from Southern China and Singapore. Meta-analysis of all samples (n = 5553) confirms that the presence of rs1131636-T, located in the 3'-UTR of RPA1, confers an inferior overall survival (HR = 1.33, 95% CI = 1.20-1.47, P = 6.31 × 10-8). Bioinformatics and biological assays show that rs1131636 has regulatory effects on upstream RPA1. Functional studies further demonstrate that RPA1 promotes the growth, invasion, migration, and radioresistance of NPC cells. Additionally, miR-1253 is identified as a suppressor for RPA1 expression, likely through regulation of its binding affinity to rs1131636 locus. Collectively, these findings provide a promising biomarker aiding in stratifying patients with poor survival, as well as a potential drug target for NPC.

Quantitative analysis of chloramphenicol, thiamphenicol, florfenicol and florfenicol amine in eggs via liquid chromatography-electrospray ionization tandem mass spectrometry.

A widely applicable method involving liquid chromatography-electrospray ionization tandem triple quadrupole mass spectrometry (LC-ESI/MS/MS) was developed for the simultaneous determination of chloramphenicol, thiamphenicol, florfenicol, and florfenicol amine in eggs. Samples were extracted with ethyl acetate-acetonitrile-ammonium hydroxide (49:49:2, v/v) and defatted with hexane saturated with acetonitrile. The analysis was carried out on a mass spectrometer via an electrospray interface operated in the positive and negative ionization modes using deuterated chloramphenicol-d5 as the internal standard. The limits of detection and limits of quantification were 0.04-0.5 µg/kg and 0.1-1.5 µg/kg in eggs, respectively. The average recoveries of the four analytes from egg samples were 90.84-108.23%, with relative standard deviations of less than 9.61%. The corresponding intra-day and inter-day variations were found to be less than 8.11% and 11.30%, respectively. Finally, the new approach was successfully applied to the quantitative determination of these analytes in 50 commercial eggs from local supermarkets.

Selectively Constrained RNA Editing Regulation Crosstalks with piRNA Biogenesis in Primates.

Although millions of RNA editing events have been reported to modify hereditary information across the primate transcriptome, evidence for their functional significance remains largely elusive, particularly for the vast majority of editing sites in noncoding regions. Here, we report a new mechanism for the functionality of RNA editing-a crosstalk with PIWI-interacting RNA (piRNA) biogenesis. Exploiting rhesus macaque as an emerging model organism closely related to human, in combination with extensive genome and transcriptome sequencing in seven tissues of the same animal, we deciphered accurate RNA editome across both long transcripts and the piRNA species. Superimposing and comparing these two distinct RNA editome profiles revealed 4,170 editing-bearing piRNA variants, or epiRNAs, that primarily derived from edited long transcripts. These epiRNAs represent distinct entities that evidence an intersection between RNA editing regulations and piRNA biogenesis. Population genetics analyses in a macaque population of 31 independent animals further demonstrated that the epiRNA-associated RNA editing is maintained by purifying selection, lending support to the functional significance of this crosstalk in rhesus macaque. Correspondingly, these findings are consistent in human, supporting the conservation of this mechanism during the primate evolution. Overall, our study reports the earliest lines of evidence for a crosstalk between selectively constrained RNA editing regulation and piRNA biogenesis, and further illustrates that such an interaction may contribute substantially to the diversification of the piRNA repertoire in primates.

Emergence, Retention and Selection: A Trilogy of Origination for Functional De Novo Proteins from Ancestral LncRNAs in Primates.

While some human-specific protein-coding genes have been proposed to originate from ancestral lncRNAs, the transition process remains poorly understood. Here we identified 64 hominoid-specific de novo genes and report a mechanism for the origination of functional de novo proteins from ancestral lncRNAs with precise splicing structures and specific tissue expression profiles. Whole-genome sequencing of dozens of rhesus macaque animals revealed that these lncRNAs are generally not more selectively constrained than other lncRNA loci. The existence of these newly-originated de novo proteins is also not beyond anticipation under neutral expectation, as they generally have longer theoretical lifespan than their current age, due to their GC-rich sequence property enabling stable ORFs with lower chance of non-sense mutations. Interestingly, although the emergence and retention of these de novo genes are likely driven by neutral forces, population genetics study in 67 human individuals and 82 macaque animals revealed signatures of purifying selection on these genes specifically in human population, indicating a proportion of these newly-originated proteins are already functional in human. We thus propose a mechanism for creation of functional de novo proteins from ancestral lncRNAs during the primate evolution, which may contribute to human-specific genetic novelties by taking advantage of existed genomic contexts.

RNA editome in rhesus macaque shaped by purifying selection.

Understanding of the RNA editing process has been broadened considerably by the next generation sequencing technology; however, several issues regarding this regulatory step remain unresolved--the strategies to accurately delineate the editome, the mechanism by which its profile is maintained, and its evolutionary and functional relevance. Here we report an accurate and quantitative profile of the RNA editome for rhesus macaque, a close relative of human. By combining genome and transcriptome sequencing of multiple tissues from the same animal, we identified 31,250 editing sites, of which 99.8% are A-to-G transitions. We verified 96.6% of editing sites in coding regions and 97.5% of randomly selected sites in non-coding regions, as well as the corresponding levels of editing by multiple independent means, demonstrating the feasibility of our experimental paradigm. Several lines of evidence supported the notion that the adenosine deamination is associated with the macaque editome--A-to-G editing sites were flanked by sequences with the attributes of ADAR substrates, and both the sequence context and the expression profile of ADARs are relevant factors in determining the quantitative variance of RNA editing across different sites and tissue types. In support of the functional relevance of some of these editing sites, substitution valley of decreased divergence was detected around the editing site, suggesting the evolutionary constraint in maintaining some of these editing substrates with their double-stranded structure. These findings thus complement the "continuous probing" model that postulates tinkering-based origination of a small proportion of functional editing sites. In conclusion, the macaque editome reported here highlights RNA editing as a widespread functional regulation in primate evolution, and provides an informative framework for further understanding RNA editing in human.

Evolutionary interrogation of human biology in well-annotated genomic framework of rhesus macaque.

With genome sequence and composition highly analogous to human, rhesus macaque represents a unique reference for evolutionary studies of human biology. Here, we developed a comprehensive genomic framework of rhesus macaque, the RhesusBase2, for evolutionary interrogation of human genes and the associated regulations. A total of 1,667 next-generation sequencing (NGS) data sets were processed, integrated, and evaluated, generating 51.2 million new functional annotation records. With extensive NGS annotations, RhesusBase2 refined the fine-scale structures in 30% of the macaque Ensembl transcripts, reporting an accurate, up-to-date set of macaque gene models. On the basis of these annotations and accurate macaque gene models, we further developed an NGS-oriented Molecular Evolution Gateway to access and visualize macaque annotations in reference to human orthologous genes and associated regulations (www.rhesusbase.org/molEvo). We highlighted the application of this well-annotated genomic framework in generating hypothetical link of human-biased regulations to human-specific traits, by using mechanistic characterization of the DIEXF gene as an example that provides novel clues to the understanding of digestive system reduction in human evolution. On a global scale, we also identified a catalog of 9,295 human-biased regulatory events, which may represent novel elements that have a substantial impact on shaping human transcriptome and possibly underpin recent human phenotypic evolution. Taken together, we provide an NGS data-driven, information-rich framework that will broadly benefit genomics research in general and serves as an important resource for in-depth evolutionary studies of human biology.

RhesusBase: a knowledgebase for the monkey research community.

Although the rhesus macaque is a unique model for the translational study of human diseases, currently its use in biomedical research is still in its infant stage due to error-prone gene structures and limited annotations. Here, we present RhesusBase for the monkey research community (http://www.rhesusbase.org). We performed strand-specific RNA-Seq studies in 10 macaque tissues and generated 1.2 billion 90-bp paired-end reads, covering >97.4% of the putative exon in macaque transcripts annotated by Ensembl. We found that at least 28.7% of the macaque transcripts were previously mis-annotated, mainly due to incorrect exon-intron boundaries, incomplete untranslated regions (UTRs) and missed exons. Compared with the previous gene models, the revised transcripts show clearer sequence motifs near splicing junctions and the end of UTRs, as well as cleaner patterns of exon-intron distribution for expression tags and cross-species conservation scores. Strikingly, 1292 exon-intron boundary revisions between coding exons corrected the previously mis-annotated open reading frames. The revised gene models were experimentally verified in randomly selected cases. We further integrated functional genomics annotations from >60 categories of public and in-house resources and developed an online accessible database. User-friendly interfaces were developed to update, retrieve, visualize and download the RhesusBase meta-data, providing a 'one-stop' resource for the monkey research community.

Hominoid-specific de novo protein-coding genes originating from long non-coding RNAs.

Tinkering with pre-existing genes has long been known as a major way to create new genes. Recently, however, motherless protein-coding genes have been found to have emerged de novo from ancestral non-coding DNAs. How these genes originated is not well addressed to date. Here we identified 24 hominoid-specific de novo protein-coding genes with precise origination timing in vertebrate phylogeny. Strand-specific RNA-Seq analyses were performed in five rhesus macaque tissues (liver, prefrontal cortex, skeletal muscle, adipose, and testis), which were then integrated with public transcriptome data from human, chimpanzee, and rhesus macaque. On the basis of comparing the RNA expression profiles in the three species, we found that most of the hominoid-specific de novo protein-coding genes encoded polyadenylated non-coding RNAs in rhesus macaque or chimpanzee with a similar transcript structure and correlated tissue expression profile. According to the rule of parsimony, the majority of these hominoid-specific de novo protein-coding genes appear to have acquired a regulated transcript structure and expression profile before acquiring coding potential. Interestingly, although the expression profile was largely correlated, the coding genes in human often showed higher transcriptional abundance than their non-coding counterparts in rhesus macaque. The major findings we report in this manuscript are robust and insensitive to the parameters used in the identification and analysis of de novo genes. Our results suggest that at least a portion of long non-coding RNAs, especially those with active and regulated transcription, may serve as a birth pool for protein-coding genes, which are then further optimized at the transcriptional level.