The role of the 3' untranslated region in post-transcriptional regulation of protein expression in mammalian cells

CTELS: A Cell-Free System for the Analysis of Translation Termination Rate

Translation termination is the final step in protein biosynthesis when the synthesized polypeptide is released from the ribosome. Understanding this complex process is important for treatment of many human disorders caused by nonsense mutations in important genes. Here, we present a new method for the analysis of translation termination rate in cell-free systems, CTELS (for C-terminally extended luciferase-based system). This approach was based on a continuously measured luciferase activity during in vitro translation reaction of two reporter mRNA, one of which encodes a C-terminally extended luciferase. This extension occupies a ribosomal polypeptide tunnel and lets the completely synthesized enzyme be active before translation termination occurs, i.e., when it is still on the ribosome. In contrast, luciferase molecule without the extension emits light only after its release. Comparing the translation dynamics of these two reporters allows visualization of a delay corresponding to the translation termination event. We demonstrated applicability of this approach for investigating the effects of cis- and trans-acting components, including small molecule inhibitors and read-through inducing sequences, on the translation termination rate. With CTELS, we systematically assessed negative effects of decreased 3' UTR length, specifically on termination. We also showed that blasticidin S implements its inhibitory effect on eukaryotic translation system, mostly by affecting elongation, and that an excess of eRF1 termination factor (both the wild-type and a non-catalytic AGQ mutant) can interfere with elongation. Analysis of read-through mechanics with CTELS revealed a transient stalling event at a "leaky" stop codon context, which likely defines the basis of nonsense suppression.

Long non-coding RNA XIST contributes to osteoarthritis progression via miR-149-5p/DNMT3A axis

Long non-coding RNAs (lncRNAs) are largely involved in the development of osteoarthritis (OA), a chronic and degenerative joint disease. The objective of this paper is to research the functional role and molecular mechanism of lncRNA X inactive specific transcript (XIST) in OA. The levels of XIST, microRNA-149-5p (miR-149-5p), and DNA methyltransferase 3A (DNMT3A) were measured. Cell viability and apoptosis rate were determined. Associated protein levels were examined through Western blot. Dual-luciferase reporter assay and RNA immunoprecipitation (RIP) assay were implemented for confirming the target relation. And the role of XIST on OA in vivo was investigated by a rat model. XIST was expressed at a high level in OA cartilage tissues and IL-1β-treated chondrocytes. XIST knockdown promoted cell viability but restrained cell apoptosis and extracellular matrix (ECM) protein degradation in IL-1β-treated chondrocytes. XIST directly targeted miR-149-5p and miR-149-5p down-regulation restored si-XIST-mediated pro-proliferative and anti-apoptotic or ECM degradative effects. DNMT3A was a target gene of miR-149-5p and DNMT3A overexpression ameliorated miR-149-5p-induced promotion of cell viability but repression of apoptosis and ECM degradation. Knockdown of XIST reduced DNMT3A level by motivating miR-149-5p expression. The inhibitory influence of XIST down-regulation on OA evolvement was also achieved by miR-149-5p/DNMT3A axis in vivo. In a word, knockdown of XIST can repress the development of OA by miR-149-5p/DNMT3A axis. This study discovers the XIST/miR-149-5p/DNMT3A axis in regulating OA evolution, which is beneficial for understanding the molecular pathomechanism and can lay a good foundation for targeted therapy of OA treatment.

Zebrafish rbm8a and magoh mutants reveal EJC developmental functions and new 3'UTR intron-containing NMD targets

Many post-transcriptional mechanisms operate via mRNA 3'UTRs to regulate protein expression, and such controls are crucial for development. We show that homozygous mutations in two zebrafish exon junction complex (EJC) core genes rbm8a and magoh leads to muscle disorganization, neural cell death, and motor neuron outgrowth defects, as well as dysregulation of mRNAs subjected to nonsense-mediated mRNA decay (NMD) due to translation termination ≥ 50 nts upstream of the last exon-exon junction. Intriguingly, we find that EJC-dependent NMD also regulates a subset of transcripts that contain 3'UTR introns (3'UI) < 50 nts downstream of a stop codon. Some transcripts containing such stop codon-proximal 3'UI are also NMD-sensitive in cultured human cells and mouse embryonic stem cells. We identify 167 genes that contain a conserved proximal 3'UI in zebrafish, mouse and humans. foxo3b is one such proximal 3'UI-containing gene that is upregulated in zebrafish EJC mutant embryos, at both mRNA and protein levels, and loss of foxo3b function in EJC mutant embryos significantly rescues motor axon growth defects. These data are consistent with EJC-dependent NMD regulating foxo3b mRNA to control protein expression during zebrafish development. Our work shows that the EJC is critical for normal zebrafish development and suggests that proximal 3'UIs may serve gene regulatory function in vertebrates.

CodAn: predictive models for precise identification of coding regions in eukaryotic transcripts

Motivation

Characterization of the coding sequences (CDSs) is an essential step in transcriptome annotation. Incorrect identification of CDSs can lead to the prediction of non-existent proteins that can eventually compromise knowledge if databases are populated with similar incorrect predictions made in different genomes. Also, the correct identification of CDSs is important for the characterization of the untranslated regions (UTRs), which are known to be important regulators of the mRNA translation process. Considering this, we present CodAn (Coding sequence Annotator), a new approach to predict confident CDS and UTR regions in full or partial transcriptome sequences in eukaryote species.

Results

Our analysis revealed that CodAn performs confident predictions on full-length and partial transcripts with the strand sense of the CDS known or unknown. The comparative analysis showed that CodAn presents better overall performance than other approaches, mainly when considering the correct identification of the full CDS (i.e. correct identification of the start and stop codons). In this sense, CodAn is the best tool to be used in projects involving transcriptomic data.

Availability

CodAn is freely available at https://github.com/pedronachtigall/CodAn.

Contact

aland@usp.br

Supplementary information

Supplementary data are available at Briefings in Bioinformatics online.

Functional Link between miR-200a and ELK3 Regulates the Metastatic Nature of Breast Cancer

Triple-negative breast cancer (TNBC) refers to breast cancer that does not have receptors for estrogen, progesterone, and HER2 protein. TNBC accounts for 10–20% of all cases of breast cancers and is characterized by its metastatic aggressiveness, poor prognosis, and limited treatment options. Here, we show that the metastatic nature of TNBC is critically regulated by a functional link between miR-200a and the transcription factor ELK3. We found that the expression levels of miR-200a and the ELK3 mRNA were negatively correlated in the luminal and TNBC subtypes of breast cancer cells. In vitro experiments revealed that miR-200a directly targets the 3’ untranslated region (UTR) of the ELK3 mRNA to destabilize the transcripts. Furthermore, ectopic expression of miR-200a impaired the migration and invasion of TNBC cells by reducing the expression level of the ELK3 mRNA. In in vivo studies, transfection of MDA-MB 231 cells (a claudin-low TNBC cell type) with exogenous miR-200a reduced their extravasation into the lung during 48 h after tail vein injection, and co-transfection of the cells with an expression plasmid harboring ELK3 that lacked an intact 3’UTR recovered their extravasation ability. Overall, our findings provide evidences that miR-200a and ELK3 is functionally linked to regulate invasive characteristics of breast cancers.

Variants in the Kisspeptin-GnRH Pathway Modulate the Hormonal Profile and Reproductive Outcomes

Kisspeptin has been identified as a key regulatory protein in the release of gonadotropin-releasing hormone (GnRH), which subsequently increases gonadotropin secretion during puberty to establish reproductive function and regulate the hypothalamic-pituitary-gonadal axis. The effects of variants in the KISS1, KISS1R, and GNRHR genes and their possible association with assisted reproduction outcomes remain to be elucidated. In this study, we used next-generation sequencing to investigate the associations of the genetic diversity at the candidate loci for KISS1, KISS1R, and GNRHR with the hormonal profiles and reproductive outcomes in 86 women who underwent in vitro fertilization treatments. Variants in the KISS1 and KISS1R genes were associated with luteinizing hormone (rs35431622:T>C), anti-Mullerian hormone (rs71745629delT), follicle-stimulating hormone (rs73507529:C>A), and estradiol (rs73507527:G>A, rs350130:A>G, and rs73507529:C>A) levels, as well as with reproductive outcomes such as the number of oocytes retrieved (s35431622:T>C), metaphasis II oocytes (rs35431622:T>C), and embryos (rs1132506:G>C). Additionally, variants in the GNRHR UTR3' (rs1038426:C>A, rs12508464:A>C, rs13150734:C>A, rs17635850:A>G, rs35683646:G>A, rs35610027:C>G, rs35845954:T>C, rs17635749:C>T, and rs7666201:C>T) were associated with low prolactin levels. A conjoint analysis of clinical, hormonal, and genetic variables using a generalized linear model identified two variants of the KISS1 gene (rs71745629delT and rs1132506:G>C) that were significantly associated with hormonal variations and reproductive outcomes. The findings suggest that variants in KISS1, KISS1R, and GNRHR genes can modulate hormone levels and reproductive outcomes.

Long Non-Coding RNA UCA1 Modulates Paclitaxel Resistance in Breast Cancer via miR-613/CDK12 Axis

Background

Paclitaxel (PTX) occupies a considerable status in the chemotherapies of breast cancer (BC), but the drug resistance keeps an obstructive factor of PTX treatment. This study was designed to explore the molecular mechanism of long non-coding RNA (lncRNA) urothelial carcinoma-associated 1 (UCA1) in PTX resistance of BC.

Methods

UCA1, microRNA-613 (miR-613) and cyclin-dependent kinase 12 (CDK12) expression was assayed through quantitative real-time polymerase chain reaction (qRT-PCR). Cell Counting Kit-8 (CCK-8) assay was implemented for evaluating the half inhibitory concentrations (IC₅₀) of PTX and cell viability. Cell apoptosis was examined by flow cytometry. The target relationship was explored using dual-luciferase reporter assay and RNA immunoprecipitation (RIP) assay. CDK12 protein level was detected through Western blot. Xenograft tumor assay was applied for assessing the influence of UCA1 on PTX resistance of BC in vivo.

Results

UCA1 expressed highly in PTX-resistant BC tissues and cells and regulated PTX resistance in BC cells by affecting cell viability and apoptosis in part. UCA1 negatively interacted with miR-613 and modulated PTX resistance via sponging miR-613. CDK12 was a downstream gene of miR-613 and miR-613 exerted the modulation of PTX resistance via targeting CDK12. Furthermore, UCA1 regulated CDK12 level through interacting with miR-613. The regulatory role of UCA1 in PTX resistance of BC was achieved by miR-613/CDK12 axis in vivo.

Conclusion

UCA1 mediated PTX resistance in BC through the miR-613/CDK12 axis, manifesting that UCA1 might improve the PTX treatment of BC as a significant therapeutic biomarker.

Dynamic Malignant Wave of Ribosome-Insulted Gut Niche via the Wnt-CTGF/CCN2 Circuit

Stress-driven ribosome dysfunction triggers an eIF2α-mediated integrated stress response to maintain cellular homeostasis. Among four key eIF2α kinases, protein kinase R (PKR) expression positively associates with poor prognoses for colorectal cancer (CRC) patients. We identified PKR-linked Wnt signaling networks that facilitate early inflammatory niche and epithelial-mesenchymal transitions of tumor tissues in response to ribosomal insults. However, the downstream Wnt signaling target fibrogenic connective tissue growth factor (CTGF/CCN2) regulates the nuclear translocation of β-catenin in a negative feedback manner. Moreover, dwindling expression of the Wnt/β-catenin pathway-regulator CTGF triggers noncanonical Wnt pathway-mediated exacerbation of intestinal cancer progression such as an increase in cancer stemness and acquisition of chemoresistance in the presence of ribosomal insults. The Wnt-CTGF-circuit-associated landscape of oncogenic signaling events was verified with clinical genomic profiling. This ribosome-associated wave of crosstalk between stress and oncogenes provides valuable insight into potential molecular interventions against intestinal malignancies.

•

PKR expression positively associates with poor prognoses for CRC patients

•

CTGF/CCN2 mediates tumor niche remodeling under PKR-activating ribosomal stress

•

CTGF/CCN2 antagonism of Wnt regulates cancer stemness and chemoresistance

Long non-coding RNA NEAT1/miR-338-3p axis impedes the progression of acute myeloid leukemia via regulating CREBRF

Background

Acute myeloid leukemia (AML) is a heterogeneous hematological disease. Our purpose of the research was to investigate the regulatory influence of long non-coding RNA (lncRNA) nuclear enriched abundant transcript 1 (NEAT1)/microRNA-338-3p (miR-338-3p)/CREB3 regulatory factor (CREBRF) in AML progression.

Methods

The associated RNA and protein levels were measured by quantitative real-time polymerase chain reaction (qRT-PCR) and Western blot, respectively. Cell growth was assessed through colony formation assay and 3-(4,5-dimethylthiazol-2-y1)-2, 5-diphenyl tetrazolium bromide (MTT) assay. Flow cytometry was exploited to determine the apoptosis rate. Cell migration and invasion were detected by transwell assay. The combination of miR-338-3p and NEAT1 or CREBRF was analyzed via the dual-luciferase reporter assay.

Results

NEAT1 and CREBRF were down-regulated in AML tissues and cells. NEAT1 up-regulation suppressed cell growth, migration and invasion but enhanced apoptosis of AML cells. Inhibition of CREBRF reverted the NEAT1-induced effects on AML cells. Moreover, NEAT1 directly targeted miR-338-3p and miR-338-3p targeted CREBRF. NEAT1/miR-338-3p could affect cellular behaviors of AML cells via the modulation of CREBRF.

Conclusion

NEAT1/miR-338-3p axis repressed the AML progression through regulating CREBRF, which might afford a favorable perspective for the AML treatment molecularly.

mRNA Metabolism in Cardiac Development and Disease: Life After Transcription

The central dogma of molecular biology illustrates the importance of mRNAs as critical mediators between genetic information encoded at the DNA level and proteomes/metabolomes that determine the diverse functional outcome at the cellular and organ levels. Although the total number of protein-producing (coding) genes in the mammalian genome is ~20,000, it is evident that the intricate processes of cardiac development and the highly regulated physiological regulation in the normal heart, as well as the complex manifestation of pathological remodeling in a diseased heart, would require a much higher degree of complexity at the transcriptome level and beyond. Indeed, in addition to an extensive regulatory scheme implemented at the level of transcription, the complexity of transcript processing following transcription is dramatically increased. RNA processing includes post-transcriptional modification, alternative splicing, editing and transportation, ribosomal loading, and degradation. While transcriptional control of cardiac genes has been a major focus of investigation in recent decades, a great deal of progress has recently been made in our understanding of how post-transcriptional regulation of mRNA contributes to transcriptome complexity. In this review, we highlight some of the key molecular processes and major players in RNA maturation and post-transcriptional regulation. In addition, we provide an update to the recent progress made in the discovery of RNA processing regulators implicated in cardiac development and disease. While post-transcriptional modulation is a complex and challenging problem to study, recent technological advancements are paving the way for a new era of exciting discoveries and potential clinical translation in the context of cardiac biology and heart disease.

RPSAP52 lncRNA Inhibits p21Waf1/CIP Expression by Interacting With the RNA Binding Protein HuR

Long noncoding RNAs have been recently demonstrated to have an important role in fundamental biological processes, and their deregulated expression has been found in several human neoplasias. Our group has recently reported a drastic overexpression of the long noncoding RNA (lncRNA) RPSAP52 (ribosomal protein SA pseudogene 52) in pituitary adenomas. We have shown that this lncRNA increased cell proliferation by upregulating the expression of the chromatinic proteins HMGA1 and HMGA2, functioning as a competing endogenous RNA (ceRNA) through competitively binding to microRNA-15a (miR-15a), miR-15b, and miR-16. The aim of this work was to identify further mechanisms by which RPSAP52 overexpression could contribute to the development of pituitary adenomas. We investigated the involvement of RPSAP52 in the modulation of the expression of cell cycle-related genes, such as p21Waf1/CIP, whose deregulation plays a critical role in pituitary cell transformation. We report that RPSAP52, interacting with the RNA binding protein HuR (human antigen R), favors the delocalization of miR-15a, miR-15b, and miR-16 on the cyclin-dependent kinase inhibitor p21Waf1/CIP1 that, accordingly, results in downregulation in pituitary adenomas. A RNA immunoprecipitation sequencing (RIPseq) analysis performed on cells overexpressing RPSAP52 identified 40 messenger RNAs (mRNAs) enriched in Argonaute 2 (AGO2) immunoprecipitated samples. Among them, we focused on GAS8 (growth arrest-specific protein 8) gene. Consistently, GAS8 expression was downregulated in all the analyzed pituitary adenomas with respect to normal pituitary and in RPSAP52-overepressing cells, supporting the role of RPSAP52 in addressing genes involved in growth inhibition and cell cycle arrest to miRNA-induced degradation. This study unveils another RPSAP52-mediated molecular mechanism in pituitary tumorigenesis.

Mitochondria in the Pulmonary Vasculature in Health and Disease: Oxygen-Sensing, Metabolism, and Dynamics

In lung vascular cells, mitochondria serve a canonical metabolic role, governing energy homeostasis. In addition, mitochondria exist in dynamic networks, which serve noncanonical functions, including regulation of redox signaling, cell cycle, apoptosis, and mitochondrial quality control. Mitochondria in pulmonary artery smooth muscle cells (PASMC) are oxygen sensors and initiate hypoxic pulmonary vasoconstriction. Acquired dysfunction of mitochondrial metabolism and dynamics contribute to a cancer-like phenotype in pulmonary arterial hypertension (PAH). Acquired mitochondrial abnormalities, such as increased pyruvate dehydrogenase kinase (PDK) and pyruvate kinase muscle isoform 2 (PKM2) expression, which increase uncoupled glycolysis (the Warburg phenomenon), are implicated in PAH. Warburg metabolism sustains energy homeostasis by the inhibition of oxidative metabolism that reduces mitochondrial apoptosis, allowing unchecked cell accumulation. Warburg metabolism is initiated by the induction of a pseudohypoxic state, in which DNA methyltransferase (DNMT)-mediated changes in redox signaling cause normoxic activation of HIF-1α and increase PDK expression. Furthermore, mitochondrial division is coordinated with nuclear division through a process called mitotic fission. Increased mitotic fission in PAH, driven by increased fission and reduced fusion favors rapid cell cycle progression and apoptosis resistance. Downregulation of the mitochondrial calcium uniporter complex (MCUC) occurs in PAH and is one potential unifying mechanism linking Warburg metabolism and mitochondrial fission. Mitochondrial metabolic and dynamic disorders combine to promote the hyperproliferative, apoptosis-resistant, phenotype in PAH PASMC, endothelial cells, and fibroblasts. Understanding the molecular mechanism regulating mitochondrial metabolism and dynamics has permitted identification of new biomarkers, nuclear and CT imaging modalities, and new therapeutic targets for PAH. © 2020 American Physiological Society. Compr Physiol 10:713-765, 2020.

Gene Polymorphisms in Boar Spermatozoa and Their Associations with Post-Thaw Semen Quality

Genetic markers have been used to assess the freezability of semen. With the advancement in molecular genetic techniques, it is possible to assess the relationships between sperm functions and gene polymorphisms. In this study, variant calling analysis of RNA-Seq datasets was used to identify single nucleotide polymorphisms (SNPs) in boar spermatozoa and to explore the associations between SNPs and post-thaw semen quality. Assessment of post-thaw sperm quality characteristics showed that 21 boars were considered as having good semen freezability (GSF), while 19 boars were classified as having poor semen freezability (PSF). Variant calling demonstrated that most of the polymorphisms (67%) detected in boar spermatozoa were at the 3’-untranslated regions (3’-UTRs). Analysis of SNP abundance in various functional gene categories showed that gene ontology (GO) terms were related to response to stress, motility, metabolism, reproduction, and embryo development. Genomic DNA was isolated from sperm samples of 40 boars. Forty SNPs were selected and genotyped, and several SNPs were significantly associated with motility and membrane integrity of frozen-thawed (FT) spermatozoa. Polymorphism in SCLT1 gene was associated with significantly higher motility and plasma membrane integrity of FT spermatozoa from boars of the GSF group compared with those of the PSF group. Likewise, polymorphisms in MAP3K20, MS4A2, and ROBO1 genes were significantly associated with reduced cryo-induced lipid peroxidation and DNA damage of FT spermatozoa from boars of the GSF group. Candidate genes with significant SNP associations, including APPL1, PLBD1, FBXO16, EML5, RAB3C, OXSR1,PRICKLE1, and MAP3K20 genes, represent potential markers for post-thaw semen quality, and they might be relevant for future improvement in the selection procedure of boars for cryopreservation. The findings of this study provide evidence indicating that polymorphisms in genes expressed in spermatozoa could be considered as factors associated with post-thaw semen quality.

Alternative polyadenylation drives oncogenic gene expression in pancreatic ductal adenocarcinoma

Alternative polyadenylation (APA) is a gene regulatory process that dictates mRNA 3'-UTR length, resulting in changes in mRNA stability and localization. APA is frequently disrupted in cancer and promotes tumorigenesis through altered expression of oncogenes and tumor suppressors. Pan-cancer analyses have revealed common APA events across the tumor landscape; however, little is known about tumor type-specific alterations that may uncover novel events and vulnerabilities. Here, we integrate RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project and The Cancer Genome Atlas (TCGA) to comprehensively analyze APA events in 148 pancreatic ductal adenocarcinomas (PDACs). We report widespread, recurrent, and functionally relevant 3'-UTR alterations associated with gene expression changes of known and newly identified PDAC growth-promoting genes and experimentally validate the effects of these APA events on protein expression. We find enrichment for APA events in genes associated with known PDAC pathways, loss of tumor-suppressive miRNA binding sites, and increased heterogeneity in 3'-UTR forms of metabolic genes. Survival analyses reveal a subset of 3'-UTR alterations that independently characterize a poor prognostic cohort among PDAC patients. Finally, we identify and validate the casein kinase CSNK1A1 (also known as CK1alpha or CK1a) as an APA-regulated therapeutic target in PDAC. Knockdown or pharmacological inhibition of CSNK1A1 attenuates PDAC cell proliferation and clonogenic growth. Our single-cancer analysis reveals APA as an underappreciated driver of protumorigenic gene expression in PDAC via the loss of miRNA regulation.

Nanomedicines to Deliver mRNA: State of the Art and Future Perspectives

The use of messenger RNA (mRNA) in gene therapy is increasing in recent years, due to its unique features compared to plasmid DNA: Transient expression, no need to enter into the nucleus and no risk of insertional mutagenesis. Nevertheless, the clinical application of mRNA as a therapeutic tool is limited by its instability and ability to activate immune responses; hence, mRNA chemical modifications together with the design of suitable vehicles result essential. This manuscript includes a revision of the strategies employed to enhance in vitro transcribed (IVT) mRNA functionality and efficacy, including the optimization of its stability and translational efficiency, as well as the regulation of its immunostimulatory properties. An overview of the nanosystems designed to protect the mRNA and to overcome the intra and extracellular barriers for successful delivery is also included. Finally, the present and future applications of mRNA nanomedicines for immunization against infectious diseases and cancer, protein replacement, gene editing, and regenerative medicine are highlighted.

Local adaptation fuels cryptic speciation in terrestrial annelids

Uncovering the genetic and evolutionary basis of cryptic speciation is a major focus of evolutionary biology. Next Generation Sequencing (NGS) allows the identification of genome-wide local adaptation signatures, but has rarely been applied to cryptic complexes - particularly in the soil milieu - as it is the case with integrative taxonomy. The earthworm genus Carpetania, comprising six previously suggested putative cryptic lineages, is a promising model to study the evolutionary phenomena shaping cryptic speciation in soil-dwelling lineages. Genotyping-By-Sequencing (GBS) was used to provide genome-wide information about genetic variability between 17 populations, and geometric morphometrics analyses of genital chaetae were performed to investigate unexplored cryptic morphological evolution. Genomic analyses revealed the existence of three cryptic species, with half of the previously-identified potential cryptic lineages clustering within them. Local adaptation was detected in more than 800 genes putatively involved in a plethora of biological functions (most notably reproduction, metabolism, immunological response and morphogenesis). Several genes with selection signatures showed shared mutations for each of the cryptic species, and genes under selection were enriched in functions related to regulation of transcription, including SNPs located in UTR regions. Finally, geometric morphometrics approaches partially confirmed the phylogenetic signal of relevant morphological characters such as genital chaetae. Our study therefore unveils that local adaptation and regulatory divergence are key evolutionary forces orchestrating genome evolution in soil fauna.

Inflammation-Induced Alternative Pre-mRNA Splicing in Mouse Alveolar Macrophages

Alveolar macrophages serve as central orchestrators of inflammatory responses in the lungs, both initiating their onset and promoting their resolution. However, the mechanisms that program macrophages for these dynamic responses are not fully understood. Over 95% of all mammalian genes undergo alternative pre-mRNA splicing. While alternative splicing has been shown to regulate inflammatory responses in macrophages in vitro, it has not been investigated on a genome-wide scale in vivo Here we used RNAseq to investigate alternative pre-mRNA splicing in alveolar macrophages isolated from lipopolysaccharide (LPS)-treated mice during the peak of inflammation and during its resolution. We found that lung inflammation induced substantial alternative pre-mRNA splicing in alveolar macrophages. The number of changes in isoform usage was greatest at the peak of inflammation and involved multiple classes of alternative pre-mRNA splicing events. Comparative pathway analysis of inflammation-induced changes in alternative pre-mRNA splicing and differential gene expression revealed overlap of pathways enriched for immune responses such as chemokine signaling and cellular metabolism. Moreover, alternative pre-mRNA splicing of genes in metabolic pathways differed in tissue resident vs. recruited (blood monocyte-derived) alveolar macrophages and corresponded to changes in core metabolism, including a switch to Warburg-like metabolism in recruited macrophages with increased glycolysis and decreased flux through the tricarboxylic acid cycle.

A Unique Regulation Region in the 3' UTR of HLA-G with a Promising Potential

Human leukocyte antigen G (HLA-G) is a non-classical human leukocyte antigen (HLA) class I protein that interacts with inhibitory receptors and is commonly overexpressed in various cancers, thereby establishing itself as an inhibitory checkpoint immune ligand. It is also expressed in trophoblast cells during pregnancy and protects the fetus from immune rejection. Despite its crucial role and its intriguing expression pattern, the regulation of HLA-G’s expression is only partially understood. HLA-G’s mRNA is expressed in many tissues but the protein expression is restricted only to the cells mentioned above. Therefore, we suggest that HLA-G is post-transcriptionally regulated. Here, we reveal a distinctive site present only in the 3′ Untranslated region (UTR) of HLA-G, which might explain its unique expression pattern. Consequently, we attempted to find binding factors such as RNA binding proteins (RBPs) and microRNAS (miRs) that regulate HLA-G expression by interacting with this distinct site present in its 3′ UTR. Our research indicates that this site should be further studied in order to reveal its significance.

A Cancer-Specific Ubiquitin Ligase Drives mRNA Alternative Polyadenylation by Ubiquitinating the mRNA 3' End Processing Complex

Alternative polyadenylation (APA) contributes to transcriptome complexity by generating mRNA isoforms with varying 3' UTR lengths. APA leading to 3' UTR shortening (3' US) is a common feature of most cancer cells; however, the molecular mechanisms are not understood. Here, we describe a widespread mechanism promoting 3' US in cancer through ubiquitination of the mRNA 3' end processing complex protein, PCF11, by the cancer-specific MAGE-A11-HUWE1 ubiquitin ligase. MAGE-A11 is normally expressed only in the male germline but is frequently re-activated in cancers. MAGE-A11 is necessary for cancer cell viability and is sufficient to drive tumorigenesis. Screening for targets of MAGE-A11 revealed that it ubiquitinates PCF11, resulting in loss of CFIm25 from the mRNA 3' end processing complex. This leads to APA of many transcripts affecting core oncogenic and tumor suppressors, including cyclin D2 and PTEN. These findings provide insights into the molecular mechanisms driving APA in cancer and suggest therapeutic strategies.

Long non-coding RNA MALAT1 facilitates the tumorigenesis, invasion and glycolysis of multiple myeloma via miR-1271-5p/SOX13 axis

Background

Long non-coding RNAs (lncRNAs) play crucial roles in the regulation and treatment of multiple myeloma (MM). The objective of this research was to study the functional mechanism of metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) in MM.

Methods

MALAT1, microRNA-1271-5p (miR-1271-5p), and SRY-Box 13 (SOX13) levels were examined by quantitative real-time polymerase chain reaction (qRT-PCR). Cell viability, apoptosis, and invasion were respectively assayed using 3-(4, 5-dimethylthiazol-2-y1)-2, 5-diphenyl tetrazolium bromide (MTT), flow cytometry, and transwell assay. Glycolysis was evaluated by glucose consumption, lactate production, ATP/ADP ratio, and the detection of related enzymes. Associated proteins were measured using Western blot. Target relation was verified via dual-luciferase reporter assay. Xenograft tumor assay was implemented to study the influence of MALAT1 on MM in vivo.

Results

The up-regulation of MALAT1 and the down-regulation of miR-1271-5p were found in MM serums and cells. MALAT1 knockdown suppressed cell viability, invasion, and glycolysis while expedited cell apoptosis in MM cells. MALAT1 directly targeted miR-1271-5p and miR-1271-5p depression reverted the effects of MALAT1 knockdown on MM cells. SOX13 was a target of miR-1271-5p and SOX13 overexpression weakened the effects of miR-1271-5p on MM. MALAT1 indirectly modulated SOX13 expression through targeting miR-1271-5p. MALAT1 down-regulation inhibited MM growth by miR-1271-5p/SOX13 axis in vivo.

Conclusion

LncRNA MALAT1 expedited MM tumorigenesis, invasion, and glycolysis via miR-1271-5p/SOX13 axis. MALAT1 might contribute to the therapy of MM as a promising indicator.

Association between single nucleotide polymorphisms in GPAT1 locus and pork quality in pigs

GPAT1 gene is considered to be a genetic marker for intramuscular fat content. The GPAT enzymes catalyze the first step in triacylglycerol synthesis. In the present study, the search for polymorphisms within the pig GPAT1 gene locus as well as association analyses between SNPs and meat quality traits were performed. The association analysis demonstrated that g.133513422C > T polymorphism affected IMF content in LL, SEMI post-mortem pH and shear force of cooked LL (p ≤ .05). While the association of g.133476803 T > C polymorphism was shown concerning IMF content, meat color (L*), shear force and energy of raw meat as well as some meat texture parameters (hardness, springiness, chewiness) (p ≤ .05). The g.133476733C > T SNP was significantly associated with LL 24 h post-mortem pH, raw meat toughness and TPA hardness and chewiness (p ≤ .05). The obtained results are promising but to confirm if the GPAT1 gene can be considered to be a genetic marker for improving the quality of pork, further research is required.

AR mRNA stability is increased with AR-antagonist resistance via 3'UTR variants

Advanced prostate cancer is often treated with AR antagonists which target the androgen receptor (AR) on which the growth of the tumour depends. Prostate cancer often develops AR-antagonist resistance via a plethora of mechanisms, many of which are as yet unknown, but it is thought that AR upregulation or AR ligand-binding site mutations, may be responsible. Here we describe the production of cell lines based on LNCaP and VCaP, with acquired resistance to the clinically relevant AR antagonists, bicalutamide and enzalutamide. In these resistant cells, we observed, via RNA-seq, that new variants in the 3'UTR of the AR mRNA were detectable and that the levels were increased both with AR-antagonist treatment and with hormonal starvation. Around 20% of AR transcripts showed a 3 kb deletion within the 6.7 kb 3'UTR sequence. Actinomycin D and luciferase fusion studies indicated that this shorter mRNA variant was inherently more stable in anti-androgen-resistant cell lines. Of additional interest was that the AR UTR variant could be detected in the sera of prostate cancer patients in a cohort of serum samples collected from patients of Gleason grades 6-10, with an increasing level correlated to increasing grade. We hypothesise that the shorter AR UTR variant is a survival adaptation to low hormone levels and/or AR-antagonist treatment in these cells, where a more stable mRNA may allow higher levels of AR expression under these conditions.

Identification of Novel Compounds Enhancing SR-BI mRNA Stability through High-Throughput Screening

Atherosclerosis is the pathological basis of most cardiovascular diseases. Reverse cholesterol transport (RCT) is a main mechanism of cholesterol homeostasis and involves the direct transport of high-density lipoprotein (HDL) cholesteryl ester by selective cholesterol uptake. Hepatic scavenger receptor class B member 1 (SR-BI) overexpression can effectively promote RCT and reduce atherosclerosis. SR-BI may be an important target for prevention or treatment of atherosclerotic disease. In our study, we inserted human SR-BI mRNA 3' untranslated region (3'UTR) downstream of the luciferase reporter gene, to establish a high-throughput screening model based on stably transfected HepG2 cells and to screen small-molecule compounds that can significantly enhance the mRNA stability of the SR-BI gene. Through multiple screenings of 25 755 compounds, the top five active compounds that have similar structures were obtained, with a positive rate of 0.19%. The five positive compounds could enhance the SR-BI expression and uptake of DiI-HDL in the hepatocyte HepG2. E238B-63 could also effectively extend the half-life of SR-BI mRNA and enhance the SR-BI mRNA and protein level and the uptake of DiI-HDL in hepatocytes in a time-dependent and dose-dependent manner. The structure-activity relationship analysis showed that the structure N-(3-hydroxy-2-pyridyl) carboxamide is possibly the key pharmacophore of the active compound, providing reference for acquiring candidate compounds with better activity. The positive small molecular compounds obtained in this study might become new drug candidates or lead compounds for the treatment of cardiovascular diseases and contribute to the further study of the posttranscriptional regulation mechanism of the SR-BI gene.

Polymorphism of the Myostatin (MSTN) Gene in Landes and Kielecka Geese Breeds

Myostatin, also known as growth differentiation factor 8 (GDF8), belongs to the TGF-β superfamily of proteins. MSTN is a highly conserved protein that acts as a negative regulator of skeletal muscle growth. Loss of myostatin functionality causes the phenotype to appear in the form of 'double musculature', among others in cattle, sheep, and house mice. The presented results of the research were carried out on two geese breeds-Landes and Kielecka. The aim of the study was to identify mutations in the MSTN gene and study their impact on body weight in both geese breeds in different periods of life. Analysis of the obtained results showed the existence of polymorphism in exon 3 (c.1231C>T) and suggested a possible association (p < 0.05) between BW and genotype in 12 weeks of life in male Kielecka geese breed. The identified polymorphism may be one of the factors important for improving body weight in the studied Kielecka breed, therefore, it is necessary to conduct further research on a larger population of geese breeds in order to more accurately estimate the effect of the identified SNP c.1231C>T on BW in geese.

TED-Seq Identifies the Dynamics of Poly(A) Length during ER Stress

Post-transcriptional RNA processing is a core mechanism of gene expression control in cell stress response. The poly(A) tail influences mRNA translation and stability, but it is unclear whether there are global roles of poly(A)-tail lengths in cell stress. To address this, we developed tail-end displacement sequencing (TED-seq) for an efficient transcriptome-wide profiling of poly(A) lengths and applied it to endoplasmic reticulum (ER) stress in human cells. ER stress induced increases in the poly(A) lengths of certain mRNAs, including known ER stress regulators, XBP1, DDIT3, and HSPA5. Importantly, the mRNAs with increased poly(A) lengths are both translationally de-repressed and stabilized. Furthermore, mRNAs in stress-induced RNA granules have shorter poly(A) tails than in the cytoplasm, supporting the view that RNA processing is compartmentalized. In conclusion, TED-seq reveals that poly(A) length is dynamically regulated upon ER stress, with potential consequences for both translation and mRNA turnover.

The p53 mRNA: an integral part of the cellular stress response

A large number of signalling pathways converge on p53 to induce different cellular stress responses that aim to promote cell cycle arrest and repair or, if the damage is too severe, to induce irreversible senescence or apoptosis. The differentiation of p53 activity towards specific cellular outcomes is tightly regulated via a hierarchical order of post-translational modifications and regulated protein-protein interactions. The mechanisms governing these processes provide a model for how cells optimize the genetic information for maximal diversity. The p53 mRNA also plays a role in this process and this review aims to illustrate how protein and RNA interactions throughout the p53 mRNA in response to different signalling pathways control RNA stability, translation efficiency or alternative initiation of translation. We also describe how a p53 mRNA platform shows riboswitch-like features and controls the rate of p53 synthesis, protein stability and modifications of the nascent p53 protein. A single cancer-derived synonymous mutation disrupts the folding of this platform and prevents p53 activation following DNA damage. The role of the p53 mRNA as a target for signalling pathways illustrates how mRNA sequences have co-evolved with the function of the encoded protein and sheds new light on the information hidden within mRNAs.

Nuclear Tau, p53 and Pin1 Regulate PARN-Mediated Deadenylation and Gene Expression

While nuclear tau plays a role in DNA damage response (DDR) and chromosome relaxation, the mechanisms behind these functions are not fully understood. Here, we show that tau forms complex(es) with factors involved in nuclear mRNA processing such as tumor suppressor p53 and poly(A)-specific ribonuclease (PARN) deadenylase. Tau induces PARN activity in different cellular models during DDR, and this activation is further increased by p53 and inhibited by tau phosphorylation at residues implicated in neurological disorders. Tau's binding factor Pin1, a mitotic regulator overexpressed in cancer and depleted in Alzheimer's disease (AD), also plays a role in the activation of nuclear deadenylation. Tau, Pin1 and PARN target the expression of mRNAs deregulated in AD and/or cancer. Our findings identify novel biological roles of tau and toxic effects of hyperphosphorylated-tau. We propose a model in which factors involved in cancer and AD regulate gene expression by interactions with the mRNA processing machinery, affecting the transcriptome and suggesting insights into alternative mechanisms for the initiation and/or developments of these diseases.

Characterisation of the metallothionein gene in the Sydney rock oyster and its expression upon metal exposure in oysters with different prior metal exposure histories

The metal-binding protein metallothionein (MT) is widely used as a biomarker of metal contamination. In this study, we cloned a MT gene (sgMT) from the Sydney rock oyster Saccostrea glomerata. The gene encodes a MT-I protein with a classical αβ domain structure and is expressed as two transcripts resulting from alternative polyadenylation. The gene promoter contains two putative metal-responsive elements (MREs) which are known to be required for metal-inducible transcription. A specific and efficient qPCR assay was developed to quantify sgMT mRNA expression. Further, we assessed whether prior metal exposure history influences sgMT mRNA expression upon subsequent metal exposure. Oysters with varying prior metal exposure histories (contaminated and reference) were exposed to Cu, Cd and Zn. Expression of sgMT generally increased with metal dose, and oysters with an elevated past metal exposure history exhibited higher sgMT expression under Cd and Zn stress, representing a potential acclimatory response to prior metal exposure.

Sequence determinants of polyadenylation-mediated regulation

The cleavage and polyadenylation reaction is a crucial step in transcription termination and pre-mRNA maturation in human cells. Despite extensive research, the encoding of polyadenylation-mediated regulation of gene expression within the DNA sequence is not well understood. Here, we utilized a massively parallel reporter assay to inspect the effect of over 12,000 rationally designed polyadenylation sequences (PASs) on reporter gene expression and cleavage efficiency. We find that the PAS sequence can modulate gene expression by over five orders of magnitude. By using a uniquely designed scanning mutagenesis data set, we gain mechanistic insight into various modes of action by which the cleavage efficiency affects the sensitivity or robustness of the PAS to mutation. Furthermore, we employ motif discovery to identify both known and novel sequence motifs associated with PAS-mediated regulation. By leveraging the large scale of our data, we train a deep learning model for the highly accurate prediction of RNA levels from DNA sequence alone (R = 0.83). Moreover, we devise unique approaches for predicting exact cleavage sites for our reporter constructs and for endogenous transcripts. Taken together, our results expand our understanding of PAS-mediated regulation, and provide an unprecedented resource for analyzing and predicting PAS for regulatory genomics applications.

Genome-wide profiling reveals alternative polyadenylation of mRNA in human non-small cell lung cancer

BACKGROUND

Lung cancer is the second most common cancer with an extremely poor overall survival rate. Post-transcriptional regulation of gene expression play many important roles in human cancer, and one of the potential mechanisms underlying this is alternative mRNA maturation at its 3' untranslated regions (3'-UTRs).

METHODS

Cancer tissues and paired adjacent normal lung tissues from 26 patients diagnosed with non-small cell lung cancer (NSCLC) were analyzed by in vitro transcription-sequencing alternative polyadenylation sites (IVT-SAPAS). 41,773,101 reads in average were obtained from each paired sample. A potential regulation of Cleavage Stimulation Factor Subunit 2 (CSTF2) on 3'UTR length of genes was tested in H460 cells.

RESULTS

1439 (10.26%) genes showed up-regulated expression and 1364 (9.72%) genes showed down-regulated expression in lung cancer tissue versus normal lung tissue, and shorten 3'UTR in cancer tissue was detected in cancer tissues collected from 96.2% (25/26) patients, indicating lung cancer tend to have shortened 3'UTRs of these identified genes. KEGG analysis showed 1855 genes with shorten 3'UTR were enriched in mTOR signaling, ubiquitin mediated proteolysis and RNA degradation. Knocking down CSTF2 expression in H460 cells results in 3'UTR elongation of genes that was identified to be with shortened length in cancer tissues.

CONCLUSION

Alternative polyadenylation (APA) site-switching of 3'UTRs is prevalent in NSCLC, and CSTF2 may serve as an oncogene regulates the 3'UTR length of cancer related genes in NSCLC.

Characterization of CCDC103 expression profiles: further insights in primary ciliary dyskinesia and in human reproduction

PROPOSE

To study CCDC103 expression profiles and understand how pathogenic variants in CCDC103 affect its expression profile at mRNA and protein level.

METHODS

To increase the knowledge about the CCDC103, we attempted genotype-phenotype correlations in two patients carrying novel homozygous (missense and frameshift) CCDC103 variants. Whole-exome sequencing, quantitative PCR, Western blot, electron microscopy, immunohistochemistry, immunocytochemistry, and immunogold labelling were performed to characterize CCDC103 expression profiles in reproductive and somatic cells.

RESULTS

Our data demonstrate that pathogenic variants in CCDC103 gene negatively affect gene and protein expression in both patients who presented absence of DA on their axonemes. Further, we firstly report that CCDC103 is expressed at different levels in reproductive tissues and somatic cells and described that CCDC103 protein forms oligomers with tissue-specific sizes, which suggests that CCDC103 possibly undergoes post-translational modifications. Moreover, we reported that CCDC103 was restricted to the midpiece of sperm and is present at the cytoplasm of the other cells.

CONCLUSIONS

Overall, our data support the CCDC103 involvement in PCD and suggest that CCDC103 may have different assemblies and roles in cilia and sperm flagella biology that are still unexplored.

CETP (Cholesteryl Ester Transfer Protein) Concentration: A Genome-Wide Association Study Followed by Mendelian Randomization on Coronary Artery Disease

BACKGROUND

We aimed to identify independent genetic determinants of circulating CETP (cholesteryl ester transfer protein) to assess causal effects of variation in CETP concentration on circulating lipid concentrations and cardiovascular disease risk.

METHODS

A genome-wide association discovery and replication study on serum CETP concentration were embedded in the NEO study (Netherlands Epidemiology of Obesity). Based on the independent identified variants, Mendelian randomization was conducted on serum lipids (NEO study) and coronary artery disease (CAD; CARDIoGRAMplusC4D consortium).

RESULTS

In the discovery analysis (n=4248), we identified 3 independent variants (P<5×10^-8) that determine CETP concentration. These single-nucleotide polymorphisms were mapped to CETP and replicated in a separate subpopulation (n=1458). Per-allele increase (SE) in serum CETP was 0.32 (0.02) µg/mL for rs247616-C, 0.35 (0.02) µg/mL for rs12720922-A, and 0.12 (0.02) µg/mL for rs1968905-G. Combined, these 3 variants explained 16.4% of the total variation in CETP concentration. One microgram per milliliter increase in genetically determined CETP concentration strongly decreased high-density lipoprotein cholesterol (-0.23 mmol/L; 95% confidence interval, -0.26 to -0.20), moderately increased low-density lipoprotein cholesterol (0.08 mmol/L; 95% confidence interval, 0.00-0.16), and was associated with an odds ratio of 1.08 (95% confidence interval, 0.94-1.23) for CAD risk.

CONCLUSIONS

This is the first genome-wide association study identifying independent variants that largely determine CETP concentration. Although high-density lipoprotein cholesterol is not a causal risk factor for CAD, it has been unequivocally demonstrated that low-density lipoprotein cholesterol lowering is proportionally associated with a lower CAD risk. Therefore, the results of our study are fully consistent with the notion that CETP concentration is causally associated with CAD through low-density lipoprotein cholesterol.

ALG-1 Influences Accurate mRNA Splicing Patterns in the Caenorhabditis elegans Intestine and Body Muscle Tissues by Modulating Splicing Factor Activities

MicroRNAs (miRNAs) are known to modulate gene expression, but their activity at the tissue-specific level remains largely uncharacterized. To study their contribution to tissue-specific gene expression, we developed novel tools to profile putative miRNA targets in the Caenorhabditis elegans intestine and body muscle. We validated many previously described interactions and identified ∼3500 novel targets. Many of the candidate miRNA targets curated are known to modulate the functions of their respective tissues. Within our data sets we observed a disparity in the use of miRNA-based gene regulation between the intestine and body muscle. The intestine contained significantly more putative miRNA targets than the body muscle highlighting its transcriptional complexity. We detected an unexpected enrichment of RNA-binding proteins targeted by miRNA in both tissues, with a notable abundance of RNA splicing factors. We developed in vivo genetic tools to validate and further study three RNA splicing factors identified as putative miRNA targets in our study (asd-2, hrp-2, and smu-2), and show that these factors indeed contain functional miRNA regulatory elements in their 3'UTRs that are able to repress their expression in the intestine. In addition, the alternative splicing pattern of their respective downstream targets (unc-60, unc-52, lin-10, and ret-1) is dysregulated when the miRNA pathway is disrupted. A reannotation of the transcriptome data in C. elegans strains that are deficient in the miRNA pathway from past studies supports and expands on our results. This study highlights an unexpected role for miRNAs in modulating tissue-specific gene isoforms, where post-transcriptional regulation of RNA splicing factors associates with tissue-specific alternative splicing.

Regulation and related mechanism of GSN mRNA level by hnRNPK in lung adenocarcinoma cells

Gelsolin (GSN) is an actin filament-capping protein that plays a key role in cell migration. Here we show that heterogeneous nuclear ribonucleoprotein K (hnRNPK) regulates GSN expression level by binding to the 3'-untranslated region (3'UTR) of GSN mRNA in non-small cell lung cancers (NSCLC) H1299 cells which are highly metastatic and express high level of GSN. We found that hnRNPK overexpression increased the mRNA and protein level of GSN, whereas hnRNPK knockdown by siRNA decreased the mRNA and protein level of GSN in both H1299 and A549 cells, indicating a positive role of hnRNPK in the regulation of GSN expression. Furthermore, hnRNPK knockdown affected the migration ability of H1299 and A549 cells which could be rescued by ectopic expression of GSN in those cells. Conversely, GSN knockdown in hnRNPK-overexpressing cells could abort the stimulatory effect of hnRNPK on the cell migration. These results suggest that hnRNPK function in the regulation of cell migration is GSN-dependent. Taken together, these data unveiled a new mechanism of regulation of the GSN expression by hnRNPK and provides new clues for the discovery of new anti-metastatic therapy.

Comprehensive and quantitative analysis of G1 cyclins. A tool for studying the cell cycle

In eukaryotes, the cell cycle is driven by the actions of several cyclin dependent kinases (CDKs) and an array of regulatory proteins called cyclins, due to the cyclical expression patterns of the latter. In yeast, the accepted pattern of cyclin waves is based on qualitative studies performed by different laboratories using different strain backgrounds, different growing conditions and media, and different kinds of genetic manipulation. Additionally, only the subset of cyclins regulating Cdc28 was included, while the Pho85 cyclins were excluded. We describe a comprehensive, quantitative and accurate blueprint of G₁ cyclins in the yeast Saccharomyces cerevisiae that, in addition to validating previous conclusions, yields new findings and establishes an accurate G₁ cyclin blueprint. For the purposes of this research, we produced a collection of strains with all G₁ cyclins identically tagged using the same and most respectful procedure possible. We report the contribution of each G₁ cyclin for a broad array of growing and stress conditions, describe an unknown role for Pcl2 in heat-stress conditions and demonstrate the importance of maintaining the 3’UTR sequence of cyclins untouched during the tagging process.

A comprehensive analysis of core polyadenylation sequences and regulation by microRNAs in a set of cancer predisposition genes

Two core polyadenylation elements (CPE) located in the 3' untranslated region of eukaryotic pre-mRNAs play an essential role in their processing: the polyadenylation signal (PAS) AAUAAA and the cleavage site (CS), preferentially a CA dinucleotide. Herein, we characterized PAS and CS sequences in a set of cancer predisposition genes (CPGs) and performed an in silico investigation of microRNAs (miRNAs) regulation to identify potential tumor-suppressive and oncogenic miRNAs. NCBI and alternative polyadenylation databases were queried to characterize CPE sequences in 117 CPGs, including 81 and 17 known tumor suppressor genes and oncogenes, respectively. miRNA-mediated regulation analysis was performed using predicted and validated data sources. Based on NCBI analyses, we did not find an established PAS in 21 CPGs, and verified that the majority of PAS already described (74.4%) had the canonical sequence AAUAAA. Interestingly, "AA" dinucleotide was the most common CS (37.5%) associated with this set of genes. Approximately 90% of CPGs exhibited evidence of alternative polyadenylation (more than one functional PAS). Finally, the mir-192 family was significantly overrepresented as regulator of tumor suppressor genes (P < 0.01), which suggests a potential oncogenic function. Overall, this study provides a landscape of CPE in CPGs, which might be useful in development of future molecular analyses covering these frequently neglected regulatory sequences.

RNA-protein interactions: disorder, moonlighting and junk contribute to eukaryotic complexity

RNA-protein interactions are crucial for most biological processes in all organisms. However, it appears that the complexity of RNA-based regulation increases with the complexity of the organism, creating additional regulatory circuits, the scope of which is only now being revealed. It is becoming apparent that previously unappreciated features, such as disordered structural regions in proteins or non-coding regions in DNA leading to higher plasticity and pliability in RNA-protein complexes, are in fact essential for complex, precise and fine-tuned regulation. This review addresses the issue of the role of RNA-protein interactions in generating eukaryotic complexity, focusing on the newly characterized disordered RNA-binding motifs, moonlighting of metabolic enzymes, RNA-binding proteins interactions with different RNA species and their participation in regulatory networks of higher order.

TYR Gene in Llamas: Polymorphisms and Expression Study in Different Color Phenotypes

Tyrosinase, encoded by TYR gene, is an enzyme that plays a major role in mammalian pigmentation. It catalyzes the oxidation of L-dihydroxy-phenylalanine (DOPA) to DOPA quinone, a precursor of both types of melanin: eumelanin and pheomelanin. TYR is commonly known as the albino locus since mutations in this gene result in albinism in several species. However, many other TYR mutations have been found to cause diluted phenotypes, like the Himalayan or chinchilla phenotypes in mice. The llama (Lama glama) presents a wide variety of coat colors ranging from non-diluted phenotypes (eumelanic and pheomelanic), through different degrees of dilution, to white. To investigate the possible contribution of TYR gene to coat color variation in llamas, we sequenced TYR exons and their flanking regions and genotyped animals with diluted, non-diluted, and white coat, including three blue-eyed white individuals. Moreover, we analyzed mRNA expression levels in skin biopsies by qPCR. TYR coding region presented nine SNPs, of which three were non-synonymous, c.428A > G, c.859G > T, and c.1490G > T. We also identified seven polymorphisms in non-coding regions, including two microsatellites, an homopolymeric repeat, and five SNPs: one in the promoter region (c.1-26C > T), two in the 3'-UTR, and two flanking the exons. Although no complete association was found between coat color and SNPs, c.1-26C > T was partially associated to diluted phenotypes. Additionally, the frequency of the G allele from c.428A > G was significantly higher in white compared to non-diluted. Results from qPCR showed that expression levels of TYR in white llamas were significantly lower (p < 0.05) than those in diluted and non-diluted phenotypes. Screening for variation in regulatory regions of TYR did not reveal polymorphisms that explain such differences. However, data from this study showed that TYR expression levels play a role in llama pigmentation.

Alternative polyadenylation coordinates embryonic development, sexual dimorphism and longitudinal growth in Xenopus tropicalis

RNA alternative polyadenylation contributes to the complexity of information transfer from genome to phenome, thus amplifying gene function. Here, we report the first X. tropicalis resource with 127,914 alternative polyadenylation (APA) sites derived from embryos and adults. Overall, APA networks play central roles in coordinating the maternal-zygotic transition (MZT) in embryos, sexual dimorphism in adults and longitudinal growth from embryos to adults. APA sites coordinate reprogramming in embryos before the MZT, but developmental events after the MZT due to zygotic genome activation. The APA transcriptomes of young adults are more variable than growing adults and male frog APA transcriptomes are more divergent than females. The APA profiles of young females were similar to embryos before the MZT. Enriched pathways in developing embryos were distinct across the MZT and noticeably segregated from adults. Briefly, our results suggest that the minimal functional units in genomes are alternative transcripts as opposed to genes.

CTRP13 attenuates vascular calcification by regulating Runx2

Vascular calcification is strongly associated with increased cardiovascular mortality and morbidity. C1q/TNF-related protein-13 (CTRP13) is a secreted adipokine that plays important roles in the cardiovascular system. However, the functional role of CTRP13 in the development of vascular calcification has yet to be explored. In this study, we collected blood samples from patients with chronic renal failure (CRF) and from rats with adenine-induced CRF. We found that the serum CTRP13 levels were decreased in patients and rats with CRF and were negatively associated with calcium deposition in the abdominal aorta. Compared to those of the controls, ectopic CTRP13 treatment significantly attenuated the calcium accumulation and alkaline phosphatase activity in the abdominal aorta of CRF rats, and β-glycerophosphate induced the formation of arterial rings and of vascular smooth muscle cells (VSMCs) and decreased the number of VSMCs that transitioned from a contractile to an osteogenic phenotype. The overexpression of Runx2 blocked CTRP13-reduced VSMC calcification. Mechanistically, CTRP13 repressed the phosphorylation of tristetraprolin (TTP), thereby activating TTP and increasing the TTP binding to the 3'untranslated region of the Runx2 mRNA, accelerating the Runx2 mRNA destabilization and degradation. In summary, these findings reveal that CTRP13 regulation is a novel method for the prevention of vascular calcification, representing a novel mechanism of the regulation of Runx2 expression in VSMCs.-Li, Y., Wang, W., Chao, Y., Zhang, F., Wang, C. CTRP13 attenuates vascular calcification by regulating Runx2.

Potential Impact of MicroRNA-423 Gene Variability in Coronary Artery Disease

AIM

Studies have evaluated the association of miRNA-423 C>A genotyping with the susceptibility to various diseases such cancers, atherosclerosis and inflammatory bowel disease but the results were contradictory. However, no studies have reported the association between miRNA-423 rs6505162 C>A polymorphism and susceptibility of coronary artery disease. MicroRNAs regulate expression of multiple genes involved in atherogenesis. Therefore, we investigated the association of microRNA-423C>T gene variations with susceptibility to coronary artery disease.

METHODOLOGY

This study was conducted on 100 coronary artery disease patients and 117 matched healthy controls. The genotyping of the microRNA-423 rs6505162C>A was performed by using Amplification refractory mutation system PCR method (ARMS-PCR).

RESULTS

A significant difference was observed in the genotype distribution among the coronary artery disease cases and sex-matched healthy controls (P=0.048). The frequencies of all three genotypes CC, CA, AA reported in the patient's samples were 55%, 41% and 4% and in the healthy controls samples were 55%, 41% and 4% respectively. Our findings showed that the microRNA-423 C>A variant was associated with an increased risk of coronary artery disease in codominant model (OR = 1.96, 95 % CI, 1.12-3.42; RR 1.35(1.05-1.75, p=0.017) of microRNA-423CA genotype and significant association in dominant model (OR 1.97, 95% CI (1.14-3.39), (CA+AA vs CC) and non-significant association for recessive model (OR=1.42, 95%CI=0.42-4.83, P=0.56, AA vs CC+CA).While, the A allele significantly increased the risk of coronary artery disease (OR =1.56, 95 % CI, 1.03-2.37; p=0.035) compared to C allele. Therefore, it was observed that more than 1.96, 1.97 and 1.56 fold increased risk of developing coronary artery disease.

CONCLUSION

Our findings indicated that microRNA-423 CA genotype and A allele are associated with an increased susceptibility to Coronary artery disease.

De novo emergence and potential function of human-specific tandem repeats in brain-related loci

Tandem repeats (TRs) are widespread in the genomes of all living organisms. In eukaryotes, they are found in both coding and noncoding regions and have potential roles in the regulation of cellular processes such as transcription, translation and in the modification of protein structure. Recent studies have highlighted TRs as a key regulator of gene expression and a potential contributor to human evolution. Thus, TRs are emerging as an important source of variation that can result in differential gene expression at intra- and inter-species levels. In this study, we performed a genome-wide survey to identify TRs that have emerged in the human lineage. We further examined these loci to explore their potential functional significance for human evolution. We identified 152 human-specific TR (HSTR) loci containing a repeat unit of more than ten bases, with most of them showing a repeat count of two. Gene set enrichment analysis showed that HSTR-associated genes were associated with biological functions in brain development and synapse function. In addition, we compared gene expression of human HSTR loci with orthologues from non-human primates (NHP) in seven different tissues. Strikingly, the expression level of HSTR-associated genes in brain tissues was significantly higher in human than in NHP. These results suggest the possibility that de novo emergence of TRs could have resulted in altered gene expression in humans within a short-time frame and contributed to the rapid evolution of human brain function.

Association between CLEC4E gene polymorphism of mincle and pulmonary tuberculosis infection in a northern Chinese population

BACKGROUND

Pulmonary tuberculosis caused by an intracellular pathogen, Mycobacterium tuberculosis continues to exist as a hazardous disease to human life globally. Genetic polymorphisms regulate resistance and susceptibility to tuberculosis. The C-type lectin receptor of family 4 member E (CLEC4E) confers protection against tuberculosis in laboratory animals but its function in influencing exposure or resistance to pulmonary tuberculosis (PTB) in humans remains obscure.

AIM

We conducted this research to analyze the effects or concomitance of CLEC4E gene variations with susceptibility to pulmonary tuberculosis in a northern Chinese population.

METHOD

In this study, 202 participants with pulmonary tuberculosis and 214 controls without PTB were enrolled. Two single nucleotide polymorphisms (SNPs) for CLEC4E on chromosome 12 were selected with a minor allele frequency of >0.05. All the SNPs were genotyped using high resolution melting analysis-PCR.

RESULTS

We estimated and compared two SNPs, rs10841845 and rs10841847. From our study findings, CLEC4E rs10841845 conferred protection against the development of pulmonary TB with a P value of <0.05 and odds ratio of <1 for all models of genetic inheritance. CLEC4E rs10841847 genotypes in co-dominant, Recessive, Dominant models and alleles had a significant statistical difference between patients and controls associated with resistance against the development of PTB (P<0.05 and OR<1).

CONCLUSION

Our findings suggest that variations at rs10841845 and rs10841847 of CLEC4E genes are associated with increased individual protection against PTB.

Hypoxia modulates protein phosphatase 2A through HIF-1α dependent and independent mechanisms in human aortic smooth muscle cells and ventricular cardiomyocytes

BACKGROUND AND PURPOSE

Although protein phosphatases regulate multiple cellular functions, their modulation under hypoxia remains unclear. We investigated expression of the protein phosphatase system under normoxic/hypoxic conditions and the mechanism by which hypoxia alters protein phosphatase 2A (PP2A) activity.

EXPERIMENTAL APPROACH

Human cardiovascular cells were cultured in cell type specific media under normoxic or hypoxic conditions (1% O₂ ). Effects on mRNA expression, phosphatase activity, post-translational modification, and involvement of hypoxia inducible factor 1α (HIF-1α) were assessed using RT-PCR, immunoblotting, an activity assay, and siRNA silencing.

KEY RESULTS

All components of the protein phosphatase system studied were expressed in each cell line. Hypoxia attenuated mRNA expression of the transcripts in a cell line- and time-dependent manner. In human aortic smooth muscle cells (HASMC) and AC16 cells, hypoxia decreased PP2Ac activity and mRNA expression without altering PP2Ac abundance. Hypoxia increased demethylated PP2Ac (DPP2Ac) and phosphatase methylesterase 1 (PME-1) abundance but decreased leucine carboxyl methyltransferase 1 (LCMT-1) abundance. HIF-1α siRNA prevented the hypoxia-mediated decrease in phosphatase activity and expression of the catalytic subunit of protein phosphatase 2A (PPP2CA), independently of altering pPP2Ac, DPP2Ac, LCMT-1, or PME-1 abundance.

CONCLUSION AND IMPLICATIONS

Cardiovascular cells express multiple components of the PP2A system. In HASMC and AC16 cells, hypoxia inhibits PP2A activity through HIF-1α-dependent and -independent mechanisms, with the latter being consistent with altered PP2A holoenzyme assembly. This indicates a complex inhibitory effect of hypoxia on the PP2A system, and highlights PP2A as a therapeutic target for diseases associated with dysregulated protein phosphorylation.

HLA-G 3' untranslated region variants +3187G/G, +3196G/G and +3035T define diametrical clinical status and disease outcome in epithelial ovarian cancer

Expression of the non-classical human leukocyte antigen-G (HLA-G) promotes cancer progression in various malignancies including epithelial ovarian cancer (EOC). As single nucleotide polymorphisms (SNPs) in the HLA-G 3' untranslated region (UTR) regulate HLA-G expression, we investigated HLA-G 3'UTR haplotypes arranged by SNPs in healthy controls (n = 75) and primary EOC patients (n = 79) and determined soluble HLA-G (sHLA-G) levels. Results were related to the clinical status and outcome. Although haplotype frequencies were similar in patients and controls, (i) sHLA-G levels were increased in EOC independent of the haplotype, (ii) homozygosity for UTR-1 or UTR-2 genotypes were significantly associated with metastases formation and presence of circulating tumor cells before therapy, whereas (iii) the UTR-5 and UTR-7 haplotypes were significantly associated with a beneficial clinical outcome regarding negative nodal status, early FIGO staging, and improved overall survival. Lastly, (iv) the ambivalent impact on clinical EOC aspects could be deduced to specific SNPs in the HLA-G 3'UTR: +3187G, +3196G and +3035T alleles. Our results give evidence that even if the genetic background of the HLA-G 3'UTR is identical between patients and controls, certain SNPs have the potential to contribute to diametrical clinical status/outcome in EOC.

Evolutionary and functional implications of 3' untranslated region length of mRNAs by comprehensive investigation among four taxonomically diverse metazoan species

Background

In eukaryotic organisms, it has been well acknowledged that 3′ untranslated regions (3′ UTRs) of mRNA are actively involved in post-transcriptional regulations of gene expression. Although both shortening and lengthening of 3′ UTRs of specific candidate genes were explicitly documented to have functional consequences, landscape of 3′ UTR lengths in relation to evolutionary dynamics and biological meanings remains to be elucidated when large-scale data become available.

Objectives

The primary objective of this study was to revealed different inter- and intra-species patterns on length distribution of 3′ UTRs in comparison with 5′ UTRs and coding regions.

Methods

In the present study, we investigated 3′ UTR lengths in a highly curated set of 57,135 mRNA sequences among four well-studied and taxonomically diverse metazoan species (fruit fly, zebrafish, mouse and human).

Results

The average length ratio of 3′–5′ UTRs considerably increased from fruit fly (twofold) to human (fivefold). Moreover, genes can be characterized by the obviously different patterns of evolutionary change on 3′ UTR lengths. By utilizing the Gene Ontology annotations, genes with differential lengths of 3′ UTRs were suggested to have the divergent functional implications. In human, we further revealed that ubiquitously transcribed genes had higher median lengths of 3′ UTRs than the genes of tissue-restricted expressions.

Conclusion

We conducted a comprehensive analysis and provided an overview regarding 3′ UTRs biology of mRNAs in animals, whereas the mechanistic explanations of 3′ UTRs length variation in correlation to regulation of gene expression still remain to be further studied.

Electronic supplementary material

The online version of this article (10.1007/s13258-019-00808-8) contains supplementary material, which is available to authorized users.

Cleavage factor 25 deregulation contributes to pulmonary fibrosis through alternative polyadenylation

Idiopathic pulmonary fibrosis (IPF) is a chronic and deadly disease with a poor prognosis and few treatment options. Pathological remodeling of the extracellular matrix (ECM) by myofibroblasts is a key factor that drives disease pathogenesis, although the underlying mechanisms remain unknown. Alternative polyadenylation (APA) has recently been shown to play a major role in cellular responses to stress by driving the expression of fibrotic factors and ECMs through altering microRNA sensitivity, but a connection to IPF has not been established. Here, we demonstrate that CFIm25, a global regulator of APA, is down-regulated in the lungs of patients with IPF and mice with pulmonary fibrosis, with its expression selectively reduced in alpha-smooth muscle actin (α-SMA) positive fibroblasts. Following the knockdown of CFIm25 in normal human lung fibroblasts, we identified 808 genes with shortened 3'UTRs, including those involved in the transforming growth factor-β signaling pathway, the Wnt signaling pathway, and cancer pathways. The expression of key pro-fibrotic factors can be suppressed by CFIm25 overexpression in IPF fibroblasts. Finally, we demonstrate that deletion of CFIm25 in fibroblasts or myofibroblast precursors using either the Col1a1 or the Foxd1 promoter enhances pulmonary fibrosis after bleomycin exposure in mice. Taken together, our results identified CFIm25 down-regulation as a novel mechanism to elevate pro-fibrotic gene expression in pulmonary fibrosis.

Identify gene expression pattern change at transcriptional and post-transcriptional levels

Gene transcription is regulated with distinct sets of regulatory factors at multiple levels. Transcriptional and post-transcriptional regulation constitute two major regulation modes of gene expression to either activate or repress the initiation of transcription and thereby control the number of proteins synthesized during translation. Disruptions of the proper regulation patterns at transcriptional and post-transcriptional levels are increasingly recognized as causes of human diseases. Consequently, identifying the differential gene expression at transcriptional and post-transcriptional levels respectively is vital to identify potential disease-associated and/or causal genes and understand their roles in the disease development. Here, we proposed a novel method with a linear mixed model that can identify a set of differentially expressed genes at transcriptional and post-transcriptional levels. The simulation and real data analysis showed our method could provide an accurate way to identify genes subject to aberrant transcriptional and post-transcriptional regulation and reveal the potential causal genes that contributed to the diseases.

HNRNP Q suppresses polyglutamine huntingtin aggregation by post-transcriptional regulation of vaccinia-related kinase 2

Misfolded proteins with abnormal polyglutamine (polyQ) expansion cause neurodegenerative disorders, including Huntington's disease. Recently, it was found that polyQ aggregates accumulate as a result of vaccinia-related kinase 2 (VRK2)-mediated degradation of TCP-1 ring complex (TRiC)/chaperonin-containing TCP-1 (CCT), which has an essential role in the prevention of polyQ protein aggregation and cytotoxicity. The levels of VRK2 are known to be much higher in actively proliferating cells but are maintained at a low level in the brain via an unknown mechanism. Here, we found that basal levels of neuronal cell-specific VRK2 mRNA are maintained by post-transcriptional, rather than transcriptional, regulation. Moreover, heterogeneous nuclear ribonucleoprotein Q (HNRNP Q) specifically binds to the 3'untranslated region of VRK2 mRNA in neuronal cells to reduce the mRNA stability. As a result, we found a dramatic decrease in CCT4 protein levels in response to a reduction in HNRNP Q levels, which was followed by an increase in polyQ aggregation in human neuroblastoma cells and mouse cortical neurons. Taken together, these results provide new insights into how neuronal HNRNP Q decreases VRK2 mRNA stability and contributes to the prevention of Huntington's disease, while also identifying new prognostic markers of HD.