Natural selection on human microRNA binding sites inferred from SNP data

Pyrococcus furiosus Argonaute-Based Fluorometric Biosensor for One-Tube Detection of Cancer-Associated Single Nucleotide Polymorphisms in MicroRNAs

MicroRNA-related single nucleotide polymorphisms (miR-SNPs) are promising biomarkers for cancer diagnostics, yet accurate detection methods remain limited. Here, we introduce a ligation-triggered Pyrococcus furiosus Argonaute (PfAgo) cleavage (LTAC) strategy for the sensitive detection of miR-SNPs, demonstrated using the rs11614913 SNP in miR-196a2, which is associated with nonsmall cell lung cancer (NSCLC). The mutant miR-196a2T serves as a scaffold for the formation of guide DNA (gDNA) catalyzed by the SplintR ligase, leading to PfAgo activation and enhanced fluorescence. In contrast, wild-type miR-196a2C cannot facilitate gDNA formation and thus fails to activate PfAgo. This method exhibits a linear relationship with the logarithm of the miR-196a2T concentration over a range of 0.2 pM to 100 nM, achieving a low detection limit of 0.15 pM. Analysis of NSCLC patient samples using LTAC reveals elevated levels of the rs11614913 SNP in miR-196a2 compared to healthy controls, underscoring the diagnostic potential of LTAC.

Identification and Validation of Lipid Metabolism Gene FASN-Associated miRNA in Wilms Tumor

miRNA has been a research hotspot in recent years and its scope of action is very wide, involving the regulation of cell proliferation, differentiation, apoptosis, and other biological behaviors. This study intends to explore the role of miRNA in the lipid metabolism and development of Wilms tumor (WT) by detecting and analyzing the differences in the expression profiles of miRNAs between the tumor and adjacent normal tissue. Gene detection was performed in tumor tissues and adjacent normal tissues of three cases of WT to screen differentially expressed miRNAs (DEMs). According to our previous research, FASN, which participates in the lipid metabolism pathway, may be a target of WT. The starBase database was used to predict FASN-targeted miRNAs. The above two groups of miRNAs were intersected to obtain FASN-targeted DEMs and then GO Ontology (GO) functional enrichment analysis of FASN-targeted DEMs was performed. Finally, the FASN-targeted DEMs were compared and further verified by qRT‒PCR. Through gene sequencing and differential analysis, 287 DEMs were obtained, including 132 upregulated and 155 downregulated miRNAs. The top ten DEMs were all downregulated. Fourteen miRNAs targeted by the lipid metabolism-related gene FASN were predicted by starBase. After intersection with the DEMs, three miRNAs were finally obtained, namely, miR-107, miR-27a-3p, and miR-335-5p. GO enrichment analysis was mainly concentrated in the Parkin-FBXW7-Cul1 ubiquitin ligase complex and response to prostaglandin E. Further experimental verification showed that miR-27a-3p was significantly correlated with WT (P = 0.0018). Imbalanced expression of miRNAs may be involved in the occurrence and development of WT through lipid metabolism. The expression of miR-27a-3p is related to the malignant degree of WT, and it may become the target of diagnosis, prognosis, and treatment of WT in the later stage.

Argonaute protein assisted drug discovery for miRNA-181c-5p and target gene ATM translation repression: a computational approach

The miRNA binds to AGO's seed region, prompting the exploration of small molecules that can offset miRNA repression of target mRNA. This miRNA-181c-5p was found to be upregulated in the chronic traumatic encephalopathy, a prevalent neurodegenerative disease in contact sports and military personals. The research aimed to identify compounds that disrupt the AGO-assisted loop formation between miRNA-181c-5p and ATM, consequently repressing the translation of ATM. Target genes from commonly three databases (DIANA-microT-CDS, miRDB, RNA22 and TargetScan) were subjected to functional annotation and clustering analysis using DAVID bioinformatics tool. Haddock server were employed to make miRNA-181c-5p:ATM-AGO complex. A total of 2594 small molecules were screened using Glide XP based on their highest binding affinity towards the complex, through a three-phase docking approach. The top 5 compounds (DB00674-Galantamine, DB00371-Meprobamate, DB00694-Daunorubicin, DB00837-Progabide, and DB00851-Dacarbazine) were further analyzed for stability in the miRNA-181c-5p:ATM-AGO-ligand complex interaction using GROMACS (version 2023.2). Hence, these findings suggest that these molecules hold potential for facilitating AGO-assisted repression of ATM gene translation.

Genome-Wide Dissection of Selection on microRNA Target Genes Involved in Rice Flower Development

Although genome-wide studies have identified a number of candidate regions evolving under selection in domesticated animals and cultivated plants, few attempts have been made, from the point of a definite biological process, to assess sequence variation and characterize the regimes of the selection on miRNA-associated motifs. Here, we performed a genome-wide dissection of nucleotide variation and selection of miRNA targets associated with rice flower development. By sampling and resequencing 26 miRNA targets for globally diverse representative populations of Asian cultivated rice and wild relatives, we found that purifying selection has reduced genetic variation at the conserved miRNA binding sites on the whole, and highly conserved miRNA binding sequences were maintained in the studied rice populations. Conversely, non-neutral evolution of positive and/or artificial selection accelerates the elevated variations at nonconserved binding sites in a population-specific behavior which may have contributed to flower development-related phenotypic variation. Taken together, our results elucidate that miRNA targets involved in flower development are under distinctive selection regimes during rice evolution.

Genetic association between FADS and ELOVL polymorphisms and the circulating levels of EPA/DHA in humans: a scoping review

BACKGROUND

Docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) are two omega-3 fatty acids that can be synthesized out of their precursor alpha-linolenic acid (ALA). FADS and ELOVL genes encode the desaturase and elongase enzymes required for EPA and DHA synthesis from ALA; however, single nucleotide polymorphisms (SNPs) in FADS and ELOVL genes could modify the levels of EPA and DHA synthesized from ALA although there is no consensus in this area. This review aims to investigate EPA and DHA circulating levels in human blood and their association with FADS or ELOVL.

METHODS

PubMed, Cochrane, and Scopus databases were used to identify research articles. They were subsequently reviewed by two independent investigators.

RESULTS

Initially, 353 papers were identified. After removing duplicates and articles not meeting inclusion criteria, 98 full text papers were screened. Finally, this review included 40 studies investigating FADS and/or ELOVL polymorphisms. A total of 47 different SNPs in FADS genes were reported. FADS1 rs174537, rs174547, rs174556 and rs174561 were the most studied SNPs, with minor allele carriers having lower levels of EPA and DHA. SNPs in the FADS genes were in high linkage disequilibrium. SNPs in FADS were correlated with levels of EPA and DHA. No conclusion could be drawn with the ELOVL polymorphisms since the number of studies was too low.

CONCLUSION

Specific SNPs in FADS gene, such as rs174537, have strong associations with circulating levels of EPA and DHA. Continued investigation regarding the impact of genetic variants related to EPA and DHA synthesis is warranted.

Convergent and divergent evolution of microRNA-mediated regulation in metazoans

The evolution of microRNAs (miRNAs) has been studied extensively to understand their roles in gene regulation and evolutionary processes. This review focuses on how miRNA-mediated regulation has evolved in bilaterian animals, highlighting both convergent and divergent evolution. Since animals and plants display significant differences in miRNA biogenesis and target recognition, the 'independent origin' hypothesis proposes that miRNA pathways in these groups independently evolved from the RNA interference (RNAi) pathway, leading to modern miRNA repertoires through convergent evolution. However, recent evidence raises the alternative possibility that the miRNA pathway might have already existed in the last common ancestor of eukaryotes, and that the differences in miRNA pathway and miRNA repertoires among animal and plant lineages arise from lineage-specific innovations and losses of miRNA pathways, miRNA acquisition, and loss of miRNAs after eukaryotic divergence. The repertoire of miRNAs has considerably expanded during bilaterian evolution, primarily through de novo creation and duplication processes, generating new miRNAs. Although ancient functionally established miRNAs are rarely lost, many newly emerged miRNAs are transient and lineage specific, following a birth-death evolutionary pattern aligning with the 'out-of-the-testis' and 'transcriptional control' hypotheses. Our focus then shifts to the convergent molecular evolution of miRNAs. We summarize how miRNA clustering and seed mimicry contribute to this phenomenon, and we review how miRNAs from different sources converge to degrade maternal messenger RNAs (mRNAs) during animal development. Additionally, we describe how miRNAs evolve across species due to changes in sequence, seed shifting, arm switching, and spatiotemporal expression patterns, which can result in variations in target sites among orthologous miRNAs across distant strains or species. We also provide a summary of the current understanding regarding how the target sites of orthologous miRNAs can vary across strains or distantly related species. Although many paralogous miRNAs retain their seed or mature sequences after duplication, alterations can occur in the seed or mature sequences or expression patterns of paralogous miRNAs, leading to functional diversification. We discuss our current understanding of the functional divergence between duplicated miRNAs, and illustrate how the functional diversification of duplicated miRNAs impacts target site evolution. By investigating these topics, we aim to enhance our current understanding of the functions and evolutionary dynamics of miRNAs. Additionally, we shed light on the existing challenges in miRNA evolutionary studies, particularly the complexity of deciphering the role of miRNA-mediated regulatory network evolution in shaping gene expression divergence and phenotypic differences among species.

In silico identification of human microRNAs pointing centrin genes in Leishmania donovani: Considering the RNAi-mediated gene control

Leishmaniasis, a parasitic disease caused by different species of the protozoa parasite Leishmania, is a neglected tropical human disease that is endemic in about a hundred countries worldwide. According to the World Health Organization (WHO), the annual incidence of cutaneous leishmaniasis (CL) is estimated to be 0.7-1.2 million cases globally, whereas the annual incidence of visceral leishmaniasis is estimated to be 0.2-0.4 million cases. In many eukaryotic organisms, including human beings and protozoan parasites, centrin genes encode proteins that play essential roles within the centrosome or basal body. Human microRNAs (miRNAs) have been linked to several infectious and non-infectious diseases associated with pathogen-host interactions, and they play the emphatic roles as gene expression regulators. In this study, we used the MirTarget bioinformatics tool, which is a machine learning-based approach implemented in miRDB, to predict the target of human miRNAs in Leishmania donovani centrin genes. For cross-validation, we utilized additional prediction algorithms, namely, RNA22 and RNAhybrid, targeting all five centrin isotypes. The centrin-3 (LDBPK_342160) and putative centrin-5 (NC_018236.1) genes in L. donovani were targeted by eight and twelve human miRNAs, respectively, among 2,635 known miRNAs (miRBase). hsa-miR-5193 consistently targeted both genes. Using TargetScan, TarBase, miRecords, and miRTarBase, we identified miRNA targets and off-targets in human homologs of centrin, inflammation, and immune-responsive genes. Significant targets were screened based on GO terminologies and KEGG pathway-enrichment analysis (Log10 p-value >0.0001). In silico tools that predict the biological roles of human miRNAs as primary gene regulators in pathogen-host interactions help unravel the regulatory patterns of these miRNAs, particularly in the early stages of inflammatory responses. It is also noted that these miRNAs played an important role in the late phase of adaptive immune response, inclusively their impacts on the immune system's response to L. donovani.

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.

Copper exposure induces ovarian granulosa cell apoptosis by activating the caspase-dependent apoptosis signaling pathway and corresponding changes in microRNA patterns

Environmental copper (Cu) contamination is a complex worldwide public health problem. However, information on the effects of Cu pollution on human reproduction is limited. Although our previous studies have indicated that Cu exposure disrupts ovarian folliculogenesis, the underlying mechanism needs to be further explored. In this study, human luteinized ovarian granulosa cells and a rat animal model were used to investigate whether Cu exposure affects ovarian follicle development by inducing apoptosis and to elucidate the possible mechanisms. The results showed that Cu exposure from weaning to sexual maturity significantly decreased the proportion of preantral follicles but increased the proportion of atretic follicles (P < 0.05). In addition, 6 mg/kg Cu increased the proportion of antral follicles, while 12 and 25 mg/kg Cu decreased it (P < 0.05). We also found that 6 mg/kg Cu exposure inhibited apoptosis of ovarian granulosa cells, while 12 and 25 mg/kg Cu promoted apoptosis (P < 0.05). Experiments on primary human luteinized ovarian granulosa cells suggested that higher levels of Cu exposure induced a significant increase in the mRNA levels of Bcl2 Bax , Fas, Caspase8, and Caspase3 (P < 0.05), and the protein levels of BAX, BCL2, CASPASE3, CASPASE8, CLE-CASPASE3, CLE-CASPASE8 and BAX/BCL2 were also increased (P < 0.05). miRNA chip analyses identified a total of 95 upregulated and 10 downregulated miRNAs in human luteinized granulosa cells exposed to Cu. Hsa-miR-19b-3p, hsa-miR-19a-3p, miR-548ar-3p, hsa-miR-652-5p, and hsa-miR-29b-5p were decreased after Cu exposure (P < 0.05). Additionally, the level of hsa-miR-144-5p was increased (P < 0.05). Together, our results reveal that Cu exposure induces abnormal ovarian folliculogenesis by inducing ovarian granulosa cell apoptosis, which is triggered by the caspase-dependent apoptosis signaling pathway, and that miRNAs may be involved in this process.

CTLA4 mRNA is downregulated by miR-155 in regulatory T cells, and reduced blood CTLA4 levels are associated with poor prognosis in metastatic melanoma patients

Cytotoxic T lymphocyte-associated antigen-4 (CTLA-4) is an immune checkpoint expressed in regulatory T (Treg) cells and activated T lymphocytes. Despite its potential as a treatment strategy for melanoma, CTLA-4 inhibition has limited efficacy. Using data from The Cancer Genome Atlas (TCGA) melanoma database and another dataset, we found that decreased CTLA4 mRNA was associated with a poorer prognosis in metastatic melanoma. To investigate further, we measured blood CTLA4 mRNA in 273 whole-blood samples from an Australian cohort and found that it was lower in metastatic melanoma than in healthy controls and associated with worse patient survival. We confirmed these findings using Cox proportional hazards model analysis and another cohort from the US. Fractionated blood analysis revealed that Treg cells were responsible for the downregulated CTLA4 in metastatic melanoma patients, which was confirmed by further analysis of published data showing downregulated CTLA-4 surface protein expression in Treg cells of metastatic melanoma compared to healthy donors. Mechanistically, we found that secretomes from human metastatic melanoma cells downregulate CTLA4 mRNA at the post-transcriptional level through miR-155 while upregulating FOXP3 expression in human Treg cells. Functionally, we demonstrated that CTLA4 expression inhibits the proliferation and suppressive function of human Treg cells. Finally, miR-155 was found to be upregulated in Treg cells from metastatic melanoma patients compared to healthy donors. Our study provides new insights into the underlying mechanisms of reduced CTLA4 expression observed in melanoma patients, demonstrating that post-transcriptional silencing of CTLA4 by miRNA-155 in Treg cells may play a critical role. Since CTLA-4 expression is downregulated in non-responder melanoma patients to anti-PD-1 immunotherapy, targeting miRNA-155 or other factors involved in regulating CTLA4 expression in Treg cells without affecting T cells could be a potential strategy to improve the efficacy of immunotherapy in melanoma. Further research is needed to understand the molecular mechanisms regulating CTLA4 expression in Treg cells and identify potential therapeutic targets for enhancing immune-based therapies.

Genetic analysis and allele-specific expression of SMAD7 3'UTR variants in human colorectal cancer reveal a novel somatic variant exhibiting allelic imbalance

BACKGROUND

Considering the impact of SMAD7 deregulation in colorectal cancer (CRC) progression and the significance of single nucleotide variant (SNV)-mediated disruptions of microRNA (miRNA)-dependent regulation for cancer susceptibility, our study aimed to analyze genetic variation in the SMAD7 3' untranslated region ( 3'UTR) in CRC, measure differences in allelic mRNA expression, and evaluate its interference with miRNA-mediated post-transcriptional regulation.

PATIENTS AND METHODS

This study included 80 patients with different CRC stages and six human colon cancer cell lines of various histological origins. SMAD7 3'UTR was analyzed by direct sequencing, followed by the relative quantification of differential allelic expression of detected variants by allele-specific qRT-PCR. In silico tools were employed for predictions of regulatory consequences of detected variants.

RESULTS

A total of four different SNVs in one cell line and nine patients were found, among which were a novel somatic point variant and three already known germline variants (rs16950113, rs1050799536, and rs1043778717). All evaluated SNVs exhibited variable extents of allelic imbalance in expression. In silico analysis predicted significant effects of SNVs on miRNA binding efficiency, with each SNV disrupting existing and creating new target sites for one or more miRNAs.

CONCLUSION

Imbalance observed in the expression of SNV alleles altering miRNA binding suggests that all investigated SNVs are potential contributing factors impacting SMAD7 expression regulation in CRC that further studies should investigate.

In silico analysis reveals hypoxia-induced miR-210-3p specifically targets SARS-CoV-2 RNA

Human coronaviruses (HCoVs) until the emergence of SARS in 2003 were associated with mild cold and upper respiratory tract infections. The ongoing pandemic caused by SARS-CoV-2 has enhanced the potential for infection and transmission as compared to other known members of this family. MicroRNAs (miRNA) are 21-25 nucleotides long non-coding RNA that bind to 3' UTR of genes and regulate almost every aspect of cellular function. Several human miRNAs have been known to target viral genomes, mostly to downregulate their expression and sometimes to upregulate also. In some cases, host miRNAs could be sequestered by the viral genome to create a condition for favourable virus existence. The ongoing SARS CoV-2 pandemic is unique based on its transmissibility and severity and we hypothesised that there could be a unique mechanism for its pathogenesis. In this study, we exploited in silico approach to identify human respiratory system-specific miRNAs targeting the viral genome of three highly pathogenic HCoVs (SARS-CoV-2 Wuhan strain, SARS-CoV, and MERS-CoV) and three low pathogenic HCoVs (OC43, NL63, and HKU1). We identified ten common microRNAs that target all HCoVs studied here. In addition, we identified unique miRNAs which targeted specifically one particular HCoV. miR-210-3p was the single unique lung-specific miRNA, which was found to target the NSP3, NSP4, and NSP13 genes of SARS-CoV-2. Further miR-210-NSP3, miR-210-NSP4, and miR-210-NSP13 SARS-CoV-2 duplexes were docked with the hAGO2 protein (PDB ID 4F3T) which showed Z-score values of -1.9, -1.7, and -1.6, respectively. The role of miR-210-3p as master hypoxia regulator and inflammation regulation may be important for SARS-CoV-2 pathogenesis. Overall, this analysis advocates that miR-210-3p be investigated experimentally in SARS-CoV-2 infection.Communicated by Ramaswamy H. Sarma.

MicroRNAs as clinical tools for diagnosis, prognosis, and therapy in prostate cancer

Prostate cancer (PCa) is one of the most commonly diagnosed cancers among men worldwide. Despite the presence of accumulated clinical strategies for PCa management, limited prognostic/sensitive biomarkers are available to follow up on disease occurrence and progression. MicroRNAs (miRNAs) are small non-coding RNAs that control gene expression through post-transcriptional regulation of their complementary target messenger RNA (mRNA). MiRNAs modulate fundamental biological processes and play crucial roles in the pathology of various diseases, including PCa. Multiple evidence proved an aberrant miRNA expression profile in PCa, which is actively involved in the carcinogenic process. The robust and pleiotropic impact of miRNAs on PCa suggests them as potential candidates to help more understand the molecular landscape of the disease, which is likely to provide tools for early diagnosis and prognosis as well as additional therapeutic strategies to manage prostate tumors. Here, we emphasize the most consistently reported dysregulated miRNAs and highlight the contribution of their altered downstream targets with PCa hallmarks. Also, we report the potential effectiveness of using miRNAs as diagnostic/prognostic biomarkers in PCa and the high-throughput profiling technologies that are being used in their detection. Another key aspect to be discussed in this review is the promising implication of miRNAs molecules as therapeutic tools and targets for fighting PCa.

Gene-environment interactions and their impact on human health

The molecular processes underlying human health and disease are highly complex. Often, genetic and environmental factors contribute to a given disease or phenotype in a non-additive manner, yielding a gene-environment (G × E) interaction. In this work, we broadly review current knowledge on the impact of gene-environment interactions on human health. We first explain the independent impact of genetic variation and the environment. We next detail well-established G × E interactions that impact human health involving environmental toxicants, pollution, viruses, and sex chromosome composition. We conclude with possibilities and challenges for studying G × E interactions.

Partial collapse of DNA tetrahedron for miRNA assay with duplex-specific nuclease-assisted amplification

We establish a facile electrochemical approach for detecting miRNA. Programmable DNA tetrahedron is designed using thiol groups for electrode modification, the amino group for the localization of electrochemical species and a hairpin structure that responds to target miRNA. In addition, duplex-specific nuclease-assisted amplification helps improve the sensitivity of this biosensor. The target-initiated partial collapse of the DNA tetrahedron event integrates recognition, electrode immobilization, signal recruitment and amplification. By measuring the sharp silver stripping peak, the highly sensitive detection of miRNA is achieved, which also performs satisfactorily challenging biological samples. This method is featured with simple operation, high sensitivity and practical utility, exhibiting great application potential in clinical diagnosis.

Identification of pathological-related and diagnostic potential circular RNAs in Stanford type A aortic dissection

Introduction

Stanford type A aortic dissection (TAAD) is one of the lethal macrovascular diseases caused by the invasion of blood into the media layer of ascending aortic wall. Inflammation, smooth muscle dysfunction, and extracellular matrix (ECM) degradation were regarded as the major pathology in affected tissue. However, the expression pattern and its regulation especially through circular RNAs (circRNAs) as an overall characteristic of TAAD molecular pathology remain unclear.

Methods

We employed CIRCexplorer2 to identify circRNAs based on the RNA sequencing (RNA-seq) data of human ascending aortic tissues to systematically assess the role of circRNA in the massive alterations of gene expression in TAAD aortas. The key circRNAs were determined by LASSO model and functionally annotated by competing endogenous RNAs (ceRNA) network and co-analysis with mRNA profile. The expression level and diagnostic capability of the 4 key circRNAs in peripheral serum were confirmed by real-time polymerase chain reaction (RT-PCR).

Results

The 4 key circRNAs, namely circPTGR1 (chr9:114341075-114348445[-]), circNOX4 (chr11:89069012-89106660[-]), circAMN1 (chr12:31854796-31862359[-]) and circUSP3 (chr15:63845913-63855207[+]), demonstrated a high power to discriminate between TAAD and control tissues, suggesting that these molecules stand for a major difference between the tissues at gene regulation level. Functionally, the ceRNA network of circRNA-miRNA-mRNA predicted by the online databases, combining gene set enrichment analysis (GSEA) and cell component prediction, revealed that the identified circRNAs covered all the aspects of primary TAAD pathology, centralized with increasing inflammatory factors and cells, and ECM destruction and loss of vascular inherent cells along with the circRNAs. Importantly, we validated the high concentration and diagnostic capability of the 4 key circRNAs in the peripheral serum in TAAD patients.

Discussion

This study reinforces the vital status of circRNAs in TAAD and the possibility of serving as promising diagnostic biomarkers.

Comprehensive analysis of roles of atrial-fibrillation-related genes in lung adenocarcinoma using bioinformatic methods

Atrial fibrillation (AF) is the most common tachyarrhythmia in the world. Lung cancer is the leading cause of cancer deaths in 93 countries. Previous studies demonstrated that the prevalence of AF was higher in patients with lung cancer. However, research on the associations between AF and lung cancer is still rare. In the present study, we first identified AF-related genes using weighted gene correlation network analysis. We then analyzed the expression profiles, prognosis, immune infiltration, and methylation characteristics of these genes in LUAD patients using bioinformatics analysis. We found several AF-related genes, including CBX3, BUB1, DSC2, P4HA1, and CYP4Z1, which differently expressed between tumor and normal tissues. Survival analysis demonstrated that CYP4Z1 was positively correlated with overall survival in LUAD patients, while CBX3, BUB1, DSC2, and P4HA1 were negatively correlated. Moreover, we found that the methylation level of DSC2 in normal lung tissues was significantly higher than that in tumor tissues, and six methylation sites in the DNA sequences of DSC2 were identified negatively correlated with its expression levels. Immune infiltration analysis suggested that levels of immune cell infiltration were related to gene expression levels in varying degrees. We identified AF-related genes and found these genes were correlated with prognosis, immune infiltration, and methylation levels in lung cancer patients. We also constructed a risk signature based on these genes in LUAD patients. We hoped that the current study could provide a novel insight into roles of AF-related genes in lung cancer patients.

Inhibition of O-GlcNAc transferase sensitizes prostate cancer cells to docetaxel

The expression of O-GlcNAc transferase (OGT) and its catalytic product, O-GlcNAcylation (O-GlcNAc), are elevated in many types of cancers, including prostate cancer (PC). Inhibition of OGT serves as a potential strategy for PC treatment alone or combinational therapy. PC is the second common cancer type in male worldwide, for which chemotherapy is still the first-line treatment. However, the function of inhibition of OGT on chemotherapeutic response in PC cells is still unknown. In this study, we show that inhibition of OGT by genetic knockdown using shRNA or by chemical inhibition using OGT inhibitors sensitize PC cells to docetaxel, which is the most common chemotherapeutic agent in PC chemotherapy. Furthermore, we identified that microRNA-140 (miR-140) directly binds to OGT mRNA 3' untranslated region and inhibits OGT expression. Moreover, docetaxel treatment stimulates miR-140 expression, whereas represses OGT expression in PC cells. Overexpression of miR-140 enhanced the drug sensitivity of PC cells to docetaxel, which could be reversed by overexpression of OGT. Overall, this study demonstrates miR-140/OGT axis as therapeutic target in PC treatment and provides a promising adjuvant therapeutic strategy for PC therapy.

miRNAs as Biomarkers in Diabetes: Moving towards Precision Medicine

Diabetes mellitus (DM) is a complex metabolic disease with many specifically related complications. Early diagnosis of this disease could prevent the progression to overt disease and its related complications. There are several limitations to using existing biomarkers, and between 24% and 62% of people with diabetes remain undiagnosed and untreated, suggesting a large gap in current diagnostic practices. Early detection of the percentage of insulin-producing cells preceding loss of function would allow for effective therapeutic interventions that could delay or slow down the onset of diabetes. MicroRNAs (miRNAs) could be used for early diagnosis, as well as for following the progression and the severity of the disease, due to the fact of their pancreatic specific expression and stability in various body fluids. Thus, many studies have focused on the identification and validation of such groups or "signatures of miRNAs" that may prove useful in diagnosing or treating patients. Here, we summarize the findings on miRNAs as biomarkers in diabetes and those associated with direct cellular reprogramming strategies, as well as the relevance of miRNAs that act as a bidirectional switch for cell therapy of damaged pancreatic tissue and the studies that have measured and tracked miRNAs as biomarkers in insulin resistance are addressed.

Integrated analysis using ToppMiR uncovers altered miRNA- mRNA regulatory networks in pediatric hepatocellular carcinoma-A pilot study

BACKGROUND

Pediatric hepatocellular carcinoma (HCC) is a group of liver cancers whose mechanisms behind their pathogenesis and progression are poorly understood.

AIM

We aimed to identify alterations in the expression of miRNAs and their putative target mRNAs in not only tumor tissues of patients with pediatric HCC but also in corresponding non-tumorous background livers by using liver tissues without underlying liver disease as a control.

METHODS AND RESULTS

We performed a small-scale miRNA and mRNA profiling of pediatric HCC (consisting of fibrolamellar carcinoma [FLC] and non-FLC HCC) and paired liver tissues to identify miRNAs whose expression levels differed significantly from control livers without underlying liver disease. ToppMiR was used to prioritize both miRNAs and their putative target mRNAs in a gene-annotation network, and the mRNA profile was used to refine the prioritization. Our analysis generated prioritized lists of miRNAs and mRNAs from the following three sets of analyses: (a) pediatric HCC versus control; (b) FLC versus control; and (c) corresponding non-tumorous background liver tissues from the same patients with pediatric HCC versus control. No liver disease liver tissues were used as the control group for all analyses. Many miRNAs whose expressions were deregulated in pediatric HCC were consistent with their roles in adult HCC and/or other non-hepatic cancers. Our gene ontology analysis of target mRNAs revealed enrichment of biological processes related to the sustenance and propagation of cancer and significant downregulation of metabolic processes.

CONCLUSION

Our pilot study indicates that alterations in miRNA-mRNA networks were detected in not only tumor tissues but also corresponding non-tumorous liver tissues from patients with pediatric HCC, suggesting multi-faceted roles of miRNAs in disease progression. Our results may lead to novel hypotheses for future large-scale studies.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

DNRLCNN: A CNN Framework for Identifying MiRNA-Disease Associations Using Latent Feature Matrix Extraction with Positive Samples

Emerging evidence indicates that miRNAs have strong relationships with many human diseases. Investigating the associations will contribute to elucidating the activities of miRNAs and pathogenesis mechanisms, and providing new opportunities for disease diagnosis and drug discovery. Therefore, it is of significance to identify potential associations between miRNAs and diseases. The existing databases about the miRNA-disease associations (MDAs) only provide the known MDAs, which can be regarded as positive samples. However, the unknown MDAs are not sufficient to regard as reliable negative samples. To deal with this uncertainty, we proposed a convolutional neural network (CNN) framework, named DNRLCNN, based on a latent feature matrix extracted by only positive samples to predict MDAs. First, by only considering the positive samples into the calculation process, we captured the latent feature matrix for complex interactions between miRNAs and diseases in low-dimensional space. Then, we constructed a feature vector for each miRNA and disease pair based on the feature representation. Finally, we adopted a modified CNN for the feature vector to predict MDAs. As a result, our model achieves better performance than other state-of-the-art methods which based CNN in fivefold cross-validation on both miRNA-disease association prediction task (average AUC of 0.9030) and miRNA-phenotype association prediction task (average AUC of 0. 9442). In addition, we carried out case studies on two human diseases, and all the top-50 predicted miRNAs for lung neoplasms are confirmed by HMDD v3.2 and dbDEMC 2.0 databases, 98% of the top-50 predicted miRNAs for heart failure are confirmed. The experiment results show that our model has the capability of inferring potential disease-related miRNAs.

Bioinformatics Analysis of Differential Gene and MicroRNA Expression in Lung Adenocarcinoma: Genetic Effects on Patient Prognosis, as Indicated by the TCGA Database

Objective:

To investigate the differential expression of genes and microRNAs (miRNAs) in patients with lung adenocarcinoma and the relationship between such changes and patient prognosis.

Methods:

We analyzed the expression levels of genes and miRNAs in lung adenocarcinoma tissues and adjacent normal tissues using The Cancer Genome Atlas database (TCGA). We analyzed the function of the differentially expressed genes and miRNAs in a co-expression network. Finally, we performed survival analysis of differential genes and miRNAs in the co-expression network using clinical data from the TCGA database.

Results:

We successfully identified 6064 differentially expressed genes: 5324 upregulated genes and 740 downregulated genes. And we identified 161 differentially expressed miRNAs: 126 upregulated miRNAs and 35 downregulated miRNAs. We identified several genes that were related to each other in the co-expression network. Further analysis revealed that the high expression levels of G6PC, APOB, F2, PAQR9, and PAQR9-AS1 genes were associated with poor prognosis. However, there was no significant correlation between the expression of hsa-mir-122 with regards to patient prognosis.

Conclusions:

Our data showed that hsa-mir-122 and a number of related genes may affect the prognosis of patients with lung adenocarcinoma by regulating the cytoskeleton, thus promoting tumor angiogenesis and the metastasis of tumor cells. The high expression levels of some differentially expressed genes was associated with the low survival rate in patients with lung adenocarcinoma. However, the levels of hsa-mir-122 were not correlated with patient prognosis.

Novel insights into the mechanisms by which lncRNA HOTAIR regulates migration and invasion in HeLa cells

HOTAIR, as one of the few well-studied oncogenic lncRNAs, is involved in human tumorigenesis and is dys-regulated in most human cancers. The transcription co-activator factor YAP1 is broadly expressed in many tissues, and promotes cancer metastasis and progression. However, the precise biological roles of HOTAIR and YAP1 in cancer cells remain unclear. In this study, we showed that HOTAIR regulates H3K27 histone modification in the promoter of miR-200a to mediate miR-200a expression by recruiting EZH2. YAP1, as a potential target gene of miR-200a, aggravated the effects of miR-200a on the migration and invasion of HeLa cells. YAP1 activated the transcription of RPL23, which is a novel downstream transcriptional-regulator of YAP1. Agreement with this, the expression of YAP1 and RPL23 was dramatically decreased after injecting HeLa cells transfected with siHOTAIR in a xenograft mouse model. Accordingly, we propose a novel model of the molecular mechanism by which HOTAIR promotes the migration and invasion of cancer cells involving the miR-200a-3p/YAP1/RPL23 axis.

Resources and Methods for the Analysis of MicroRNA Function in Drosophila

Since the widespread discovery of microRNAs (miRNAs) 20 years ago, the Drosophila melanogaster model system has made important contributions to understanding the biology of this class of noncoding RNAs. These contributions are based on the amenability of this model system not only for biochemical analysis but molecular, genetic, and cell biological analyses as well. Nevertheless, while the Drosophila genome is now known to encode 258 miRNA precursors, the function of only a small minority of these have been well characterized. In this review, we summarize the current resources and methods that are available to study miRNA function in Drosophila with a particular focus on the large-scale resources that enable systematic analysis. Application of these methods will accelerate the discovery of ways that miRNAs are embedded into genetic networks that control basic features of metazoan cells.

The Interaction of Single Nucleotide Polymorphisms on Fibroblast Growth Factor 19 Superfamily Genes Is Associated With Alcohol Dependence-Related Aggression

Alcohol dependence (AD) is characterized by compulsive alcohol consumption, which involves behavioral impairments such as aggression. Members of fibroblast growth factor (FGF) 19 superfamily, including FGF19, FGF21, and FGF23, are major endocrine mediators that play an important role in alcohol metabolism and alcohol related disorders. The objective of the present study is to explore the possible associations among the interaction of single nucleotide polymorphisms (SNPs) of the FGF 19 superfamily, AD occurrence, and aggression in patients with AD. A total of 956 subjects were enrolled in this study, including 482 AD patients and 474 healthy controls (HCs). Michigan alcoholism screening test (MAST) was used to measure the level of AD, a Chinese version of the Buss-Perry Aggression Questionnaire was used to evaluate the aggressive behavior of subjects, and MassARRAY@ system was used to genotype rs948992 of FGF19, rs11665841 and rs11665896 of FGF21, rs7955866 and rs11063118 of FGF23. The results showed that AD patients presented a significantly higher level of aggression compared to HCs, and MAST scores were significantly positively associated Buss-Perry aggression scores (r = 0.402, p < 0.001) in AD patients. The interaction of FGF19 rs948992 TC × FGF21 rs11665896 GG presented the high-risk genotype combination predicting the high level of AD. In addition, the interaction of FGF19 rs948992 TC × FGF21 rs11665896 TG × FGF23 rs11063118 TT presented the high-risk genotype combination predicting the high level of aggression in AD patients. Our results added evidence linking the combination of rs948992 TC × rs11665896 TG × rs11063118 TT to aggressive behavior in AD patients and pointed out the potential usefulness of the SNPs of FGF19 superfamily as a predictor for the aggression in AD patients.

Cell-type-specific effects of genetic variation on chromatin accessibility during human neuronal differentiation

Common genetic risk for neuropsychiatric disorders is enriched in regulatory elements active during cortical neurogenesis. However, it remains poorly understood as to how these variants influence gene regulation. To model the functional impact of common genetic variation on the noncoding genome during human cortical development, we performed the assay for transposase accessible chromatin using sequencing (ATAC-seq) and analyzed chromatin accessibility quantitative trait loci (QTL) in cultured human neural progenitor cells and their differentiated neuronal progeny from 87 donors. We identified significant genetic effects on 988/1,839 neuron/progenitor regulatory elements, with highly cell-type and temporally specific effects. A subset (roughly 30%) of chromatin accessibility-QTL were also associated with changes in gene expression. Motif-disrupting alleles of transcriptional activators generally led to decreases in chromatin accessibility, whereas motif-disrupting alleles of repressors led to increases in chromatin accessibility. By integrating cell-type-specific chromatin accessibility-QTL and brain-relevant genome-wide association data, we were able to fine-map and identify regulatory mechanisms underlying noncoding neuropsychiatric disorder risk loci.

Association analysis of miRNA-related genetic polymorphisms in miR-143/145 and KRAS with colorectal cancer susceptibility and survival

Background: There is accumulating evidence of aberrant expression of miR-143 and miR-145 and their target gene KRAS in colorectal cancer (CRC). We hypothesize that single nucleotide polymorphisms (SNPs) within or near mRNA–microRNA (miRNA) binding sites may affect miRNA/target gene interaction, resulting in differential mRNA/protein expression and promoting the development and progression of CRC. Methods: We conducted a case–control study of 507 patients with CRC recruited from a tertiary hospital and 497 population-based controls to assess the association of genetic polymorphisms in miR-143/145 and the KRAS 3′ untranslated region (3′UTR) with susceptibility to CRC and patients’ survival. In addition, genetic variations of genomic regions located from 500 bp upstream to 500 bp downstream of the miR-143/miR-145 gene and the 3′UTR of KRAS were selected for analysis using the Haploview and HaploReg software. Results: Using publicly available expression profiling data, we found that miR-143/145 and KRAS expression were all reduced in rectal cancer tissue compared with adjacent non-neoplastic large intestinal mucosa. The rs74693964 C/T variant located 65 bp downstream of miR-145 genomic regions was observed to be associated with susceptibility to CRC (adjusted odds ratio (OR): 2.414, 95% CI: 1.385–4.206). Cumulative effects of miR-143 and miR-145 on CRC risk were observed (P_trend=0.03). Patients having CRC carrying variant genotype TT of KRAS rs712 had poorer survival (log-rank P=0.044, adjusted hazard ratio (HR): 4.328, 95% CI: 1.236–15.147). Conclusions: Our results indicate that miRNA-related polymorphisms in miR-143/145 and KRAS are likely to be deleterious and represent potential biomarkers for susceptibility to CRC and patients’ survival.

MicroRNAs in ascending thoracic aortic aneurysms

Thoracic Aortic Aneurysm (TAA) is characterized by the dilation of the aorta and is fatal if not diagnosed and treated appropriately. The underlying genetic mechanisms have not been completely delineated, so better knowledge of the physiopathology of TAAs is needed to improve detection and therapy. MicroRNAs (miRNAs) regulate gene expression post-transcriptionally and are known to be involved in cardiovascular diseases (CVDs). The current study aimed to identify miRNAs that can be used as possible biomarkers for the early diagnosis of patients with ascending TAAs (ATAAs). MiRNA expression was profiled by NanoString nCounter technology using 12 samples including tissue and pre- and post-surgical plasma from ATAA patients. Four miRNAs were selected and further validated by real time polymerase chain reaction (RT-PCR) in 22 plasma samples from which three miRNAs (hsa-miR140-5p, hsa-miR-191-5p and hsa-miR-214-3p) showed significant expression level differences between the two types of plasma samples. Further analyses of the corresponding predicted target genes by these miRNAs, revealed two genes (Myotubularin-related protein 4 (MTMR4) and Phosphatase 1 catalytic subunit β (PPP1CB)) whose expression was inversely correlated with the expression of their respective miRNAs. Overall, in this pilot study, we identified three miRNAs that might serve as potential biomarkers and therapeutic targets in ATAA.

Construction and Analysis of a ceRNA Network in Cardiac Fibroblast During Fibrosis Based on in vivo and in vitro Data

Aims

Activation of cardiac fibroblasts (CF) is crucial to cardiac fibrosis. We constructed a cardiac fibroblast-related competing endogenous RNA (ceRNA) network. Potential functions related to fibrosis of “hub genes” in this ceRNA network were explored.

Materials and Methods

The Gene Expression Omnibus database was searched for eligible datasets. Differentially expressed messenger (m)RNA (DE-mRNA) and long non-coding (lnc)RNA (DE-lncRNA) were identified. microRNA was predicted and validated. A predicted ceRNA network was constructed and visualized by Cytoscape, and ceRNA crosstalk was validated. A Single Gene Set Enrichment Analysis (SGSEA) was done, and the Comparative Toxicogenomics Database (CTD) was employed to analyze the most closely associated pathways and diseases of DE-mRNA in the ceRNA network. The functions of DE-mRNA and DE-lncRNA in the ceRNA network were validated by small interfering (si)RNA depletion.

Results

The GSE97358 and GSE116250 datasets (which described differentially expressed genes in human cardiac fibroblasts and failing ventricles, respectively) were used for analyses. Four-hundred-and-twenty DE-mRNA and 39 DE-lncRNA, and 369 DE-mRNA and 93 DE-lncRNA were identified, respectively, in the GSE97358 and GSE116250 datasets. Most of the genes were related to signal transduction, cytokine activity, and cell proliferation. Thirteen DE-mRNA with the same expression tendency were overlapped in the two datasets. Twenty-three candidate microRNAs were predicted and the expression of 11 were different. Only two DE-lncRNA were paired to any one of 11 microRNA. Finally, two mRNA [ADAM metallopeptidase domain 19, (ADAM19) and transforming growth factor beta induced, (TGFBI)], three microRNA (miR-9-5p, miR-124-3p, and miR-153-3p) and two lncRNA (LINC00511 and SNHG15) constituted our ceRNA network. siRNA against LINC00511 increased miR-124-3p and miR-9-5p expression, and decreased ADAM19 and TGFBI expression, whereas siRNA against SNHG15 increased miR-153-3p and decreased ADAM19 expression. ADAM19 and TGFBI were closely related to the TGF-β1 pathway and cardiac fibrosis, as shown by SGSEA and CTD, respectively. Depletion of two mRNA or two lncRNA could alleviate CF activation.

Conclusions

The CF-specific ceRNA network, including two lncRNA, three miRNA, and two mRNA, played a crucial role during cardiac fibrosis, which provided potential target genes in this field.

miR-875-5p exerts tumor-promoting function via down-regulation of CAPZA1 in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is one of the leading causes of cancer deaths worldwide. Currently, efficient genetic markers for diagnosis and treatment of ESCC are lacking. MicroRNAs (miRNAs) are global genetic regulators that control cancer gene expression by binding to the 3'untranslated regions (3'UTRs) of targeting mRNAs. In addition, miRNAs function as oncogenes or tumor suppressors in the progression of tumors. In the current study, we found that hsa-miR-875-5p (miR-875-5p) exhibited amplification in ESCC according to the TCGA database. Then, xCELLigence Real-Time Cell Analyzer (RTCA)-MP system and colony formation assays were employed to detect cell proliferationand colony formationability. The results showed that miR-875-5p promoted the proliferation ESCC cells. Subsequently, transwell results indicated that miR-875-5p promoted the invasion and migration of ESCC cells. Furthermore, we showed that miR-875-5p was able to bind to CAPZA13'UTR, which contains the single nucleotide polymorphism (SNP), rs373245753, as reported in our previous study involving WGS and WES on ESCC. Subsequently, mRNA affinity pull-down assays verifiedthat the SNP disrupts miR-875-5p binding to CAPZA1. The current study is the first demonstration that miR-875-5p may function as an oncogene via down-regulation of CAPZA1 expression in ESCC.

Inhibition of miR-188-5p Suppresses Progression of Experimental Abdominal Aortic Aneurysms

ABSTRACT

Abdominal aortic aneurysm (AAA) is an aging-related degenerative disease. miR-188-5p was reported to induce cell senescence and play a key role in aging-related disease. Therefore, in this study, we investigated miR-188-5p expression during progression in experimental AAAs. Furthermore, we investigated whether inhibition of miR-188-5p could suppress AAA progression. Experimental AAAs were created in 9-12-week-old male C57BL/6J mice by transient intra-aortic infusion of porcine pancreatic elastase. Expression of miR-188-5p levels were assessed in aneurysmal and control aortae during the progression of aneurysm. For inhibition experiment, miR-188 inhibiting group mice were injected with AAV2-miR188-5p sponge through tail vein and control group mice were injected with AAV2-CMV-GFP. Influences on experimental AAA progression were assessed by measurements of aortic diameter and histopathologic analysis at sacrifice. Meanwhile, immunohistochemistry and fluorescence in situ hybridization were used to determine the inflammatory cells infiltration and colocalization of miR-188-5p in aortic sections. Expression of miR-188-5p is upregulated during progression of AAA. Importantly, miR-188-5p inhibition treatment prevented enlargement of experimental aneurysms. Meanwhile, miR-188-5p inhibition regimens attenuated medial elastin degradation, smooth muscle cell depletion, and mural angiogenesis and the accumulation of macrophages, T cells, and angiogenesis. Furthermore, colocalization of miR188-5p with CD68 and CD3 was observed, which suggest miR-188-5p was expressed mainly in infiltrated macrophages and T cells. Expression of miR-188-5p is increased in experimental AAAs. Treatment with miR-188-5p inhibition limits experimental AAA progression, with histologic evidence of reduced neovessels and attenuated mural leukocyte infiltration. These findings underscore the potential significance of miR-188-5p in aneurysm pathogenesis and as a target for suppression of AAA disease.

Pervasive Selection against MicroRNA Target Sites in Human Populations

MicroRNA target sites are often conserved during evolution and purifying selection to maintain such sites is expected. On the other hand, comparative analyses identified a paucity of microRNA target sites in coexpressed transcripts, and novel target sites can potentially be deleterious. We proposed that selection against novel target sites pervasive. The analysis of derived allele frequencies revealed that, when the derived allele is a target site, the proportion of nontarget sites is higher than expected, particularly for highly expressed microRNAs. Thus, new alleles generating novel microRNA target sites can be deleterious and selected against. When we analyzed ancestral target sites, the derived (nontarget) allele frequency does not show statistical support for microRNA target allele conservation. We investigated the joint effects of microRNA conservation and expression and found that selection against microRNA target sites depends mostly on the expression level of the microRNA. We identified microRNA target sites with relatively high levels of population differentiation. However, when we analyze separately target sites in which the target allele is ancestral to the population, the proportion of single-nucleotide polymorphisms with high F_st significantly increases. These findings support that population differentiation is more likely in target sites that are lost than in the gain of new target sites. Our results indicate that selection against novel microRNA target sites is prevalent and, although individual sites may have a weak selective pressure, the overall effect across untranslated regions is not negligible and should be accounted when studying the evolution of genomic sequences.

miR‑7/SP1/TP53BP1 axis may play a pivotal role in NSCLC radiosensitivity

MicroRNA‑7 (miR‑7) has been identified as a tumor suppressor in non‑small cell lung cancer (NSCLC) and a radiosensitivity regulator. Numerous studies have revealed that specific protein 1 (SP1) plays a critical role in the tumorigenesis of various types of cancers and regulates radiosensitivity and tumor suppressor p53‑binding protein 1 (TP53BP1), which plays an essential role in DNA repair. However, it is not clear whether miR‑7 has a regulatory effect on SP1 and TP53BP1 in NSCLC. In the present study it was revealed that miR‑7 directly binds to the 3'UTR of SP1, thereby suppressing SP1 expression to regulate radiosensitivity. Overexpression of miR‑7 and SP1 and knockdown of miR‑7 and SP1 were performed using lentiviral transfection. Protein and mRNA abundance of SP1 and TP53BP1 were determined using western blotting and RT‑qPCR, respectively, while miR‑7 binding to SP1 was validated using a luciferase reporter assay. Biological function analysis indicated that miR‑7 negatively regulated SP1 and inhibited cell proliferation, migration, and invasion when combined with radiation. It was also revealed that the expression of TP53BP1 was positively regulated by SP1 or negatively regulated by miR‑7. In conclusion, SP1 was a target of miR‑7, and the decreased expression of SP1 resulting from miR‑7 overexpression in NSCLC was vital for improving radiosensitivity in NSCLC cells. Moreover, SP1 expression was detected in 95 paired NSCLC and adjacent normal tissues, and it was determined that SP1 was significantly upregulated in NSCLC tissues and that its upregulation was correlated with the degree of tissue differentiation. Thus, SP1 and/or miR‑7 may be potential molecular targets in NSCLC radiotherapy.

Impaired cell migration and structural defects in myeloid cells overexpressing miR-30b and miR-142-3p

Macrophages (MΦ) and dendritic cells (DC) play a fundamental role in shaping immune responses by sensing a plethora of Pathogen Associated Molecular Patterns (PAMPs), phagocytosis and antigen presentation to T lymphocytes. These important biological processes require efficient cell movement and an intact cellular morphology for dynamic interaction. The role of microRNAs (miRs) in this regard, however, is not well understood. In the present study, we show that miR-30b and miR-142-3p regulate migration and morphology of MΦ and DC. Transient overexpression of miR-30b and miR-142-3p attenuates migration and these cells display unique morphological deformities observed under electron microscopy. In addition, miR-142-3p overexpression in MΦ impaired phagocytosis of FITC-conjugated latex beads using live microscopy imaging. Interestingly, live cell imaging and F-actin staining revealed marked changes in the cell polarity and actin polymerization status, respectively. To identify miR-142-3p regulated pathways, we profiled global transcriptome changes in miR-142-3p or control mimic transfected DC. Expression of several genes were differentially altered by miR-142-3p and were associated with pathways related to cell movement, cell adhesion, and cytoskeletal rearrangement. Bioinformatics analysis identified a significant subset of downregulated genes with one or more predicted miR-142-3p binding sites in their 3'UTR strongly suggesting direct post-transcriptional impact of these miRNAs on multiple transcripts. Using dual luciferase assays, novel miR-142-3p binding sites were validated for three genes (Vinculin, Dab2 and Skap2) directly associated with cytoskeletal rearrangement and cell movement. In summary, our results show that miR-30b and miR-142-3p are regulators of myeloid cell cytoskeletal homeostasis and morphology.

Dynamic landscape and evolution of m6A methylation in human

m⁶A is a prevalent internal modification in mRNAs and has been linked to the diverse effects on mRNA fate. To explore the landscape and evolution of human m⁶A, we generated 27 m⁶A methylomes across major adult tissues. These data reveal dynamic m⁶A methylation across tissue types, uncover both broadly or tissue-specifically methylated sites, and identify an unexpected enrichment of m⁶A methylation at non-canonical cleavage sites. A comparison of fetal and adult m⁶A methylomes reveals that m⁶A preferentially occupies CDS regions in fetal tissues. Moreover, the m⁶A sub-motifs vary between fetal and adult tissues or across tissue types. From the evolutionary perspective, we uncover that the selection pressure on m⁶A sites varies and depends on their genic locations. Unexpectedly, we found that ∼40% of the 3′UTR m⁶A sites are under negative selection, which is higher than the evolutionary constraint on miRNA binding sites, and much higher than that on A-to-I RNA modification. Moreover, the recently gained m⁶A sites in human populations are clearly under positive selection and associated with traits or diseases. Our work provides a resource of human m⁶A profile for future studies of m⁶A functions, and suggests a role of m⁶A modification in human evolutionary adaptation and disease susceptibility.

Ambient Temperature is A Strong Selective Factor Influencing Human Development and Immunity

Solar energy, which is essential for the origin and evolution of all life forms on Earth, can be objectively recorded through attributes such as climatic ambient temperature (CAT), ultraviolet radiation (UVR), and sunlight duration (SD). These attributes have specific geographical variations and may cause different adaptation traits. However, the adaptation profile of each attribute and the selective role of solar energy as a whole during human evolution remain elusive. Here, we performed a genome-wide adaptation study with respect to CAT, UVR, and SD using the Human Genome Diversity Project-Centre Etude Polymorphism Humain (HGDP-CEPH) panel data. We singled out CAT as the most important driving force with the highest number of adaptive loci (6 SNPs at the genome-wide 1 × 10⁻⁷ level; 401 at the suggestive 1 × 10⁻⁵ level). Five of the six genome-wide significant adaptation SNPs were successfully replicated in an independent Chinese population (N = 1395). The corresponding 316 CAT adaptation genes were mostly involved in development and immunity. In addition, 265 (84%) genes were related to at least one genome-wide association study (GWAS)-mapped human trait, being significantly enriched in anthropometric loci such as those associated with body mass index (χ²; P < 0.005), immunity, metabolic syndrome, and cancer (χ²; P < 0.05). For these adaptive SNPs, balancing selection was evident in Euro-Asians, whereas obvious positive and/or purifying selection was observed in Africans. Taken together, our study indicates that CAT is the most important attribute of solar energy that has driven genetic adaptation in development and immunity among global human populations. It also supports the non-neutral hypothesis for the origin of disease-predisposition alleles in common diseases.

Strategies to Modulate MicroRNA Functions for the Treatment of Cancer or Organ Injury

Cancer and organ injury-such as that occurring in the perioperative period, including acute lung injury, myocardial infarction, and acute gut injury-are among the leading causes of death in the United States and impose a significant impact on quality of life. MicroRNAs (miRNAs) have been studied extensively during the last two decades for their role as regulators of gene expression, their translational application as diagnostic markers, and their potential as therapeutic targets for disease treatment. Despite promising preclinical outcomes implicating miRNA targets in disease treatment, only a few miRNAs have reached clinical trials. This likely relates to difficulties in the delivery of miRNA drugs to their targets to achieve efficient inhibition or overexpression. Therefore, understanding how to efficiently deliver miRNAs into diseased tissues and specific cell types in patients is critical. This review summarizes current knowledge on various approaches to deliver therapeutic miRNAs or miRNA inhibitors and highlights current progress in miRNA-based disease therapy that has reached clinical trials. Based on ongoing advances in miRNA delivery, we believe that additional therapeutic approaches to modulate miRNA function will soon enter routine medical treatment of human disease, particularly for cancer or perioperative organ injury. SIGNIFICANCE STATEMENT: MicroRNAs have been studied extensively during the last two decades in cancer and organ injury, including acute lung injury, myocardial infarction, and acute gut injury, for their regulation of gene expression, application as diagnostic markers, and therapeutic potentials. In this review, we specifically emphasize the pros and cons of different delivery approaches to modulate microRNAs, as well as the most recent exciting progress in the field of therapeutic targeting of microRNAs for disease treatment in patients.

Brain microRNAs among social and solitary bees

Evolutionary transitions to a social lifestyle in insects are associated with lineage-specific changes in gene expression, but the key nodes that drive these regulatory changes are unknown. We examined the relationship between social organization and lineage-specific microRNAs (miRNAs). Genome scans across 12 bee species showed that miRNA copy-number is mostly conserved and not associated with sociality. However, deep sequencing of small RNAs in six bee species revealed a substantial proportion (20-35%) of detected miRNAs had lineage-specific expression in the brain, 24-72% of which did not have homologues in other species. Lineage-specific miRNAs disproportionately target lineage-specific genes, and have lower expression levels than shared miRNAs. The predicted targets of lineage-specific miRNAs are not enriched for genes with caste-biased expression or genes under positive selection in social species. Together, these results suggest that novel miRNAs may coevolve with novel genes, and thus contribute to lineage-specific patterns of evolution in bees, but do not appear to have significant influence on social evolution. Our analyses also support the hypothesis that many new miRNAs are purged by selection due to deleterious effects on mRNA targets, and suggest genome structure is not as influential in regulating bee miRNA evolution as has been shown for mammalian miRNAs.

Genetic Variants in the 3'UTR of BRCA1 and BRCA2 Genes and their Putative Effects on the microRNA Mechanism in Hereditary Breast and Ovarian Cancer

Hereditary breast and ovarian cancer (HBOC) syndrome is mainly caused by mutations in the BRCA1 and BRCA2 genes. The 3’UTR region allows for the binding of microRNAs, which are involved in genetic tune regulation. We aimed to identify allelic variants on 3’UTR miRNA-binding sites in the BRCA1 and BRCA2 genes in HBOC patients. Blood samples were obtained from 50 patients with HBOC and from 50 controls. The 3’UTR regions of BRCA1 and BRCA2 were amplified by PCR and sequenced to identify genetic variants using bioinformatics tools. We detected nine polymorphisms in 3’UTR, namely: four in BRCA1 (rs3092995 (C/G), rs8176318 (C/T), rs111791349 (G/A), and rs12516 (C/T)) and five in BRCA2 (rs15869 (A/C), rs7334543 (A/G), rs1157836 (A/G), and rs75353978 (TT/del TT)). A new variant in position c.*457 (A/C) on 3’UTR of BRCA2 was also identified. The following three variants increased the risk of HBOC in the study population: rs111791349-A, rs15869-C, and c.*457-C (odds ratio (OR) range 3.7–15.4; p < 0.05). Genetic variants into the 3’UTR of BRCA1 and BRCA2 increased the risk of HBOC between 3.7–15.4 times in the study population. The presence/absence of these polymorphisms may influence the loss/creation of miRNA binding sites, such as hsa-miR-1248 in BRCA1 3′UTR or the hsa-miR-548 family binding site in BRCA2. Our results add new evidence of miRNA participation in the pathogenesis of HBOC.

miR-6086 inhibits ovarian cancer angiogenesis by downregulating the OC2/VEGFA/EGFL6 axis

miRNAs have emerged as a pivotal component of gene regulatory networks, mediating cytokines secretion, cell cycle, and differentiation regulation. However, how miRNAs collaborate with transcription factors and downstream effector proteins that determine the fate of ovarian cancer cells remains to be understood, especially regarding to mechanism of tumor angiogenesis regulation. Based on the qRT-PCR and IHC analysis, we found that miR-6086 was maintained a very low level both in ovarian cancer cell lines and tissues. Further, we identified OC2 and EGFL6 as the direct targets of miR-6086 by luciferase assay and we observed an inverse relationship between the expression of miR-6086 and the OC2/VEGFA/EGFL6 axis. The Western blotting analysis suggested that OC2 could directly upregulate VEGFA and indirectly up-regulate EGFL6 through VEGFA. Moreover, miR-6086 could indirectly downregulate VEGFA through OC2. In addition, miR-6086, siOC2 and siEGFL6 could negatively regulate the tumor growth and angiogenesis of ovarian cancer (Skov3) in the animal studies, with the inhibition rates of 77.07%, 69.89%, and 73.62%, respectively (**p < 0.01). Moreover, the tumor cell proliferation, migration, and invasion of ovarian cancer cell lines (Caov3 and Skov3) and vascular formation (HUVECs) were significantly suppressed in vitro, by decreasing the AKT/MAPK pathways (*p < 0.05). Taken together, our results reveal that miR-6086 can suppress the angiogenesis networks in ovarian cancer by down-regulating the OC2/VEGFA/EGFL6 axis, directly or indirectly, which may provide potential targets for tumor therapeutics.

SG-LSTM-FRAME: a computational frame using sequence and geometrical information via LSTM to predict miRNA-gene associations

MOTIVATION

MircroRNAs (miRNAs) regulate target genes and are responsible for lethal diseases such as cancers. Accurately recognizing and identifying miRNA and gene pairs could be helpful in deciphering the mechanism by which miRNA affects and regulates the development of cancers. Embedding methods and deep learning methods have shown their excellent performance in traditional classification tasks in many scenarios. But not so many attempts have adapted and merged these two methods into miRNA-gene relationship prediction. Hence, we proposed a novel computational framework. We first generated representational features for miRNAs and genes using both sequence and geometrical information and then leveraged a deep learning method for the associations' prediction.

RESULTS

We used long short-term memory (LSTM) to predict potential relationships and proved that our method outperformed other state-of-the-art methods. Results showed that our framework SG-LSTM got an area under curve of 0.94 and was superior to other methods. In the case study, we predicted the top 10 miRNA-gene relationships and recommended the top 10 potential genes for hsa-miR-335-5p for SG-LSTM-core. We also tested our model using a larger dataset, from which 14 668 698 miRNA-gene pairs were predicted. The top 10 unknown pairs were also listed.

AVAILABILITY

Our work can be download in https://github.com/Xshelton/SG_LSTM.

CONTACT

luojiawei@hnu.edu.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Briefings in Bioinformatics online.

Integrative analysis reveals RNA G-quadruplexes in UTRs are selectively constrained and enriched for functional associations

G-quadruplex (G4) sequences are abundant in untranslated regions (UTRs) of human messenger RNAs, but their functional importance remains unclear. By integrating multiple sources of genetic and genomic data, we show that putative G-quadruplex forming sequences (pG4) in 5' and 3' UTRs are selectively constrained, and enriched for cis-eQTLs and RNA-binding protein (RBP) interactions. Using over 15,000 whole-genome sequences, we find that negative selection acting on central guanines of UTR pG4s is comparable to that of missense variation in protein-coding sequences. At multiple GWAS-implicated SNPs within pG4 UTR sequences, we find robust allelic imbalance in gene expression across diverse tissue contexts in GTEx, suggesting that variants affecting G-quadruplex formation within UTRs may also contribute to phenotypic variation. Our results establish UTR G4s as important cis-regulatory elements and point to a link between disruption of UTR pG4 and disease.

TP53/miR-34a-associated signaling targets SERPINE1 expression in human pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is a disease of aging. The TP53 gene product regulates cell growth, aging, and cancer. To determine the important targets of TP53 in PDAC, we examined the expression of 440 proteins on a reverse phase protein array (RPPA) in PDAC-derived MIA-PaCa-2 cells which either had WT-TP53 or lacked WT-TP53. MIA-PaCa-2 cells have a TP53 mutation as well as mutant KRAS and represent a good in vitro model to study PDAC. RPPA analysis demonstrated expression of tumor promoting proteins in cells that lacked WT-TP53; and this feature could be reversed significantly when the cells were transfected with vector encoding WT-TP53 or treated with berberine or a modified berberine (BBR). Expression of miR-34a-associated signaling was elevated in cells expressing WT-TP53 compared to cells expressing mTP53. Results from in vivo studies using human PDAC specimens confirmed the in vitro results as the expression of miR-34a and associated signaling was significantly decreased in PDAC specimens compared to non-cancerous tissues. This study determined SERPINE1 as a miR-34a target with relevance to the biology of PDAC. Thus, we have identified a key target (SERPINE1) of the TP53/miR-34a axis that may serve as a potential biomarker for early detection of pancreatic cancer.

Min3: Predict microRNA target gene using an improved binding-site representation method and support vector machine

MicroRNAs are single-stranded noncoding RNAs known to down-regulate target genes at the protein or mRNA level. Computational prediction of targets is essential for elucidating the detailed functions of microRNA. However, prediction specificity and sensitivity of the existing algorithms still need to be improved to generate useful hypotheses for subsequent experimental testing. A new microRNA binding-site representation method was developed, which uses four symbols "|", ":", "∼", and "∧" (indicating paired, unpaired, insertion, and bulge, respectively) to represent the status of each nucleotide base pair in the microRNA binding site. New features were established with the information of every two adjacent symbols. There are 12 possible combinations and the frequency of each defines a set of novel and useful features. A comprehensive training dataset is constructed for mammalian microRNAs with positive targets obtained from the microRNA target depository in the miRTarbase, while negative targets were derived from pseudo-microRNA bindings. An SVM model was established using the training dataset and a new software called Min3 was developed. Performance of Min3 was assessed with intensively studied examples of miR-155 and miR-92a. Prediction results showed that Min3 can discover 47% of experimental conformed targets on average. The overlapping is above 20% on average when compared with TargetScan and miRanda. Annotations of the public microRNA datasets showed that there is a negative effect (up-regulation) of the Min3 targets for the knock out/down of miR-155 and miR-92a. Six top ranked targets were selected for validation by wet-lab experiments, and five of them showed a regulation effect. The Min3 can be a good alternative to current microRNA target discovery software. This tool is available at https://sourceforge.net/projects/mirt3.

Mechanism of epithelial‑mesenchymal transition inhibited by miR‑203 in non‑small cell lung cancer

The aim of the present study was to investigate whether miR-203 can inhibit transforming growth factor-β (TGF-β)-induced epithelial-mesenchymal transition (EMT), and the migration and invasion ability of non-small cell lung cancer (NSCLC) cells by targeting SMAD3. In the present study, the expression levels of miR-203, SMAD3 mRNA and protein in NSCLC tissues were examined, as well as their corresponding paracancerous samples. The miR-203 mimics and miR-203 inhibitor were transfected into the H226 cell line. RT-qPCR was used to assess the expression levels of E-cadherin, Snail, N-cadherin and vimentin mRNA, and western blotting was performed to detect the expression levels of p-SMAD2, SMAD2, p-SMAD3, SMAD3 and SMAD4. The cell migration and invasion abilities were detected by Transwell assays. The target site of SMAD3 was predicted by the combined action between miR-203 and dual luciferase. The results revealed that the RNA levels of miR-203, compared with paracancerous tissues, were decreased in NSCLC tissues, while SMAD3 mRNA and protein levels were upregulated, and miR-203 inhibited SMAD3 expression. Induction of TGF-β led to decreased E-cadherin mRNA levels, upregulation of Snail, N-cadherin and vimentin mRNA levels (P<0.05), and significant increase in cell migration and invasion, whereas transfection of miR-203 mimics reversed the aforementioned results (P<0.05). Conversely, miR-203 inhibitor could further aggravate the aforementioned results (P<0.05). Western blot results revealed that transfection of miR-203 mimics significantly reduced the protein expression of SMAD3 and p-SMAD3 (P<0.05). Furthermore, the results of the Dual-Luciferase assay revealed that miR-203 inhibited SMAD3 expression by interacting with specific regions of its 3′-UTR. Overall, a novel mechanism is revealed, in which, miR-203 can inhibit SMAD3 by interacting with specific regions of the 3′-UTR of SMAD3, thereby restraining TGF-β-induced EMT progression and migration and invasion of NSCLC cells.