Complex Selection on Human Polyadenylation Signals Revealed by Polymorphism and Divergence Data

Quantifying negative selection in human 3' UTRs uncovers constrained targets of RNA-binding proteins

Many non-coding variants associated with phenotypes occur in 3' untranslated regions (3' UTRs), and may affect interactions with RNA-binding proteins (RBPs) to regulate gene expression post-transcriptionally. However, identifying functional 3' UTR variants has proven difficult. We use allele frequencies from the Genome Aggregation Database (gnomAD) to identify classes of 3' UTR variants under strong negative selection in humans. We develop intergenic mutability-adjusted proportion singleton (iMAPS), a generalized measure related to MAPS, to quantify negative selection in non-coding regions. This approach, in conjunction with in vitro and in vivo binding data, identifies precise RBP binding sites, miRNA target sites, and polyadenylation signals (PASs) under strong selection. For each class of sites, we identify thousands of gnomAD variants under selection comparable to missense coding variants, and find that sites in core 3' UTR regions upstream of the most-used PAS are under strongest selection. Together, this work improves our understanding of selection on human genes and validates approaches for interpreting genetic variants in human 3' UTRs.

Regulation of closely juxtaposed proto-oncogene c-fms and HMGXB3 gene expression by mRNA 3' end polymorphism in breast cancer cells

Sense-antisense mRNA pairs generated by convergent transcription is a way of gene regulation. c-fms gene is closely juxtaposed to the HMGXB3 gene in the opposite orientation, in chromosome 5. The intergenic region (IR) between c-fms and HMGXB3 genes is 162 bp. We found that a small portion (∼4.18%) of HMGXB3 mRNA is transcribed further downstream, including the end of the c-fms gene generating antisense mRNA against c-fms mRNA. Similarly, a small portion (∼1.1%) of c-fms mRNA is transcribed further downstream, including the end of the HMGXB3 gene generating antisense mRNA against the HMGXB3 mRNA. Insertion of the strong poly(A) signal sequence in the IR results in decreased c-fms and HMGXB3 antisense mRNAs, resulting in up-regulation of both c-fms and HMGXB3 mRNA expression. miR-324-5p targets HMGXB3 mRNA 3' UTR, and as a result, regulates c-fms mRNA expression. HuR stabilizes c-fms mRNA, and as a result, down-regulates HMGXB3 mRNA expression. UALCAN analysis indicates that the expression pattern between c-fms and HMGXB3 proteins are opposite in vivo in breast cancer tissues. Together, our results indicate that the mRNA encoded by the HMGXB3 gene can influence the expression of adjacent c-fms mRNA, or vice versa.

Long-read sequencing uncovers a complex transcriptome topology in varicella zoster virus

BACKGROUND

Varicella zoster virus (VZV) is a human pathogenic alphaherpesvirus harboring a relatively large DNA molecule. The VZV transcriptome has already been analyzed by microarray and short-read sequencing analyses. However, both approaches have substantial limitations when used for structural characterization of transcript isoforms, even if supplemented with primer extension or other techniques. Among others, they are inefficient in distinguishing between embedded RNA molecules, transcript isoforms, including splice and length variants, as well as between alternative polycistronic transcripts. It has been demonstrated in several studies that long-read sequencing is able to circumvent these problems.

RESULTS

In this work, we report the analysis of the VZV lytic transcriptome using the Oxford Nanopore Technologies sequencing platform. These investigations have led to the identification of 114 novel transcripts, including mRNAs, non-coding RNAs, polycistronic RNAs and complex transcripts, as well as 10 novel spliced transcripts and 25 novel transcription start site isoforms and transcription end site isoforms. A novel class of transcripts, the nroRNAs are described in this study. These transcripts are encoded by the genomic region located in close vicinity to the viral replication origin. We also show that the ORF63 exhibits a complex structural variation encompassing the splice sites of VZV latency transcripts. Additionally, we have detected RNA editing in a novel non-coding RNA molecule.

CONCLUSIONS

Our investigations disclosed a composite transcriptomic architecture of VZV, including the discovery of novel RNA molecules and transcript isoforms, as well as a complex meshwork of transcriptional read-throughs and overlaps. The results represent a substantial advance in the annotation of the VZV transcriptome and in understanding the molecular biology of the herpesviruses in general.

Alternative Polyadenylation of Mammalian Transcripts Is Generally Deleterious, Not Adaptive

Alternative polyadenylation (APA) produces from the same gene multiple mature RNAs with varying 3' ends. Although APA is commonly believed to generate beneficial functional diversity and be adaptive, we hypothesize that most genes have one optimal polyadenylation site and that APA is caused largely by deleterious polyadenylation errors. The error hypothesis, but not the adaptive hypothesis, predicts that, as the expression level of a gene increases, its polyadenylation diversity declines, relative use of the major (presumably optimal) polyadenylation site increases, and that of each minor (presumably nonoptimal) site decreases. It further predicts that the number of polyadenylation signals per gene is smaller than the random expectation and that polyadenylation signals for major but not minor sites are under purifying selection. All of these predictions are confirmed in mammals, suggesting that numerous defective RNAs are produced in normal cells, many phenotypic variations at the molecular level are nonadaptive, and cellular life is noisier than is appreciated.

microRNAs in Cancer Susceptibility

microRNAs (miRNAs) are a small RNA species without protein-coding potential. However, they are key modulators of protein translation. Many studies have linked miRNAs with cancer initiation, progression, diagnosis, and prognosis, and recent studies have also linked them with cancer etiology and susceptibility, especially through single-nucleotide polymorphisms (SNPs). This review discusses some of the recent advances in miRNA-SNP literature-including SNPs in miRNA genes, miRNA target sites, and the processing machinery. In addition, we highlight some emerging areas of interest, including isomiRs and non-3'UTR focused miRNA-binding mechanisms that could provide further novel insight into the relationship between miR-SNPs and cancer. Finally, we note that additional epidemiological and experimental research is needed to close the gap in our understanding of the genotype-phenotype relationship between miRNA-SNPs and cancer.