Significant impact of miRNA-target gene networks on genetics of human complex traits

MicroRNAs in atrial fibrillation target genes in structural remodelling

We aim to elucidate how miRNAs regulate the mRNA signature of atrial fibrillation (AF), to gain mechanistic insight and identify candidate targets for future therapies. We present combined miRNA-mRNA sequencing using atrial tissues of patient without AF (n = 22), with paroxysmal AF (n = 22) and with persistent AF (n = 20). mRNA sequencing previously uncovered upregulated epithelial to mesenchymal transition, endothelial cell proliferation and extracellular matrix remodelling involving glycoproteins and proteoglycans in AF. MiRNA co-sequencing discovered miRNAs regulating the mRNA expression changes. Key downregulated miRNAs included miR-135b-5p, miR-138-5p, miR-200a-3p, miR-200b-3p and miR-31-5p and key upregulated miRNAs were miR-144-3p, miR-15b-3p, miR-182-5p miR-18b-5p, miR-4306 and miR-206. MiRNA expression levels were negatively correlated with the expression levels of a multitude of predicted target genes. Downregulated miRNAs associated with increased gene expression are involved in upregulated epithelial and endothelial cell migration and glycosaminoglycan biosynthesis. In vitro inhibition of miR-135b-5p and miR-138-5p validated an effect of miRNAs on multiple predicted targets. Altogether, the discovered miRNAs may be explored in further functional studies as potential targets for anti-fibrotic therapies in AF.

The SYSCID map: a graphical and computational resource of molecular mechanisms across rheumatoid arthritis, systemic lupus erythematosus and inflammatory bowel disease

Chronic inflammatory diseases (CIDs), including inflammatory bowel disease (IBD), rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) are thought to emerge from an impaired complex network of inter- and intracellular biochemical interactions among several proteins and small chemical compounds under strong influence of genetic and environmental factors. CIDs are characterised by shared and disease-specific processes, which is reflected by partially overlapping genetic risk maps and pathogenic cells (e.g., T cells). Their pathogenesis involves a plethora of intracellular pathways. The translation of the research findings on CIDs molecular mechanisms into effective treatments is challenging and may explain the low remission rates despite modern targeted therapies. Modelling CID-related causal interactions as networks allows us to tackle the complexity at a systems level and improve our understanding of the interplay of key pathways. Here we report the construction, description, and initial applications of the SYSCID map (https://syscid.elixir-luxembourg.org/), a mechanistic causal interaction network covering the molecular crosstalk between IBD, RA and SLE. We demonstrate that the map serves as an interactive, graphical review of IBD, RA and SLE molecular mechanisms, and helps to understand the complexity of omics data. Examples of such application are illustrated using transcriptome data from time-series gene expression profiles following anti-TNF treatment and data from genome-wide associations studies that enable us to suggest potential effects to altered pathways and propose possible mechanistic biomarkers of treatment response.

Multi-ancestry genome-wide association analyses identify novel genetic mechanisms in rheumatoid arthritis

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10-8), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA.

Identification of functionally important miRNA targeted genes associated with child obesity trait in genome-wide association studies

BACKGROUND

Genome-wide association studies (GWAS) have uncovered thousands of genetic variants that are associated with complex human traits and diseases. miRNAs are single-stranded non-coding RNAs. In particular, genetic variants located in the 3'UTR region of mRNAs may play an important role in gene regulation through their interaction with miRNAs. Existing studies have not been thoroughly conducted to elucidate 3'UTR variants discovered through GWAS. The goal of this study is to analyze patterns of GWAS functional variants located in 3'UTRs about their relevance in the network between hosting genes and targeting miRNAs, and elucidate the association between the genes harboring these variants and genetic traits.

METHODS

We employed MIGWAS, ANNOVAR, MEME, and DAVID software packages to annotate the variants obtained from GWAS for 31 traits and elucidate the association between their harboring genes and their related traits. We identified variants that occurred in the motif regions that may be functionally important in affecting miRNA binding. We also conducted pathway analysis and functional annotation on miRNA targeted genes harboring 3'UTR variants for a trait with the highest percentage of 3'UTR variants occurring.

RESULTS

The Child Obesity trait has the highest percentage of 3'UTR variants (75%). Of the 16 genes related to the Child Obesity trait, 5 genes (ETV7, GMEB1, NFIX, ZNF566, ZBTB40) had a significant association with the term DNA-Binding (p < 0.05). EQTL analysis revealed 2 relevant tissues and 10 targeted genes associated with the Child Obesity trait. In addition, Red Blood Cells (RBC), Hemoglobin (HB), and Package Cell Volume (PCV) have overlapping variants. In particular, the PIM1 variant occurred inside the HB Motif region 37,174,641-37,174,660, and LUC7L3 variant occurred inside RBC Motif region 50,753,918-50,753,937.

CONCLUSION

Variants located in 3'UTR can alter the binding affinity of miRNA and impact gene regulation, thus warranting further annotation and analysis. We have developed a bioinformatics bash pipeline to automatically annotate variants, determine the number of variants in different categories for each given trait, and check common variants across different traits. This is a valuable tool to annotate a large number of GWAS result files.

GWAS identifies candidate susceptibility loci and microRNA biomarkers for acute encephalopathy with biphasic seizures and late reduced diffusion

Acute encephalopathy with biphasic seizures and late reduced diffusion (AESD) is a severe encephalopathy preceded by viral infections with high fever. AESD is a multifactorial disease, however, few disease susceptibility genes have previously been identified. Here, we conducted a genome-wide association study (GWAS) and assessed functional variants in non-coding regions to study genetic susceptibility in AESD using 254 Japanese children with AESD and 799 adult healthy controls. We also performed a microRNA enrichment analysis using GWAS statistics to search for candidate biomarkers in AESD. The variant with the lowest p-value, rs1850440, was located in the intron of serine/threonine kinase 39 gene (STK39) on chromosome 2q24.3 (p = 2.44 × 10^-7, odds ratio = 1.71). The minor allele T of rs1850440 correlated with the stronger expression of STK39 in peripheral blood. This variant possessed enhancer histone modification marks in STK39, the encoded protein of which activates the p38 mitogen-activated protein kinase (MAPK) pathway. In the replication study, the odds ratios of three SNPs, including rs1850440, showed the same direction of association with that in the discovery stage GWAS. One of the candidate microRNAs identified by the microRNA enrichment analysis was associated with inflammatory responses regulated by the MAPK pathway. This study identified STK39 as a novel susceptibility locus of AESD, found microRNAs as potential biomarkers, and implicated immune responses and the MAPK cascade in its pathogenesis.

Genetic architecture of microRNA expression and its link to complex diseases in the Japanese population

Understanding the genetic effects on non-coding RNA (ncRNA) expression facilitates functional characterization of disease-associated genetic loci. Among several classes of ncRNAs, microRNAs (miRNAs) are key post-transcriptional gene regulators. Despite its biological importance, previous studies on the genetic architecture of miRNA expression focused mostly on the European individuals, underrepresented in other populations. Here, we mapped miRNA expression quantitative trait loci (miRNA-eQTL) for 343 miRNAs in 141 Japanese using small RNA sequencing (sRNA-seq) and whole-genome sequencing (WGS), identifying 1275 cis-miRNA-eQTL variants for 40 miRNAs (false discovery rate < 0.2). Of these, 25 miRNAs having eQTL were unreported in the European studies, including 5 miRNAs with their lead variant monomorphic in the European populations, which demonstrates the value of miRNA-eQTL analysis in diverse ancestral populations. MiRNAs with eQTL effect showed allele-specific expression (ASE) (e.g. miR-146a-3p), and ASE analysis further detected cis-regulatory variants not captured by the conventional miRNA-eQTL mapping (e.g. miR-933). We identified a copy number variation (CNV) associated with miRNA expression (e.g. miR-570-3p, P = 7.2 × 10-6), which contributes to a more comprehensive landscape of miRNA-eQTLs. To elucidate a post-transcriptional modification in miRNAs, we created a catalog of miRNA-editing sites, including ten canonical and six non-canonical sites. Finally, by integrating the miRNA-eQTLs and Japanese genome-wide association studies of 25 complex traits (mean n = 192 833), we conducted a transcriptome-wide association study (TWAS), identifying miR-1908-5p as a potential mediator for adult height, colorectal cancer, and type 2 diabetes (P < 9.1 × 10-5). Our study broadens the population diversity in ncRNA-eQTL studies and contributes to functional annotation of disease-associated loci found in non-European populations.

The nature of coxofemoral joint pathology across family Canidae

We evaluated coxofemoral joints from museum specimens of: Vulpes lagopus; Vulpes vulpes; Vulpes velox; Nyctereutes procyonoides; Urocyon cinereoargenteus; Aenocyon [Canis] dirus; Canis latrans; Canis lupus lupus; Canis lupus familiaris; C. l. familiaris × latrans; and Canis dingo. Acetabular components included: fossa; articular surface; medial and lateral articular margins; and periarticular surfaces. Acetabular components variably revealed: osteophyte-like features; varying appearance of articular margin rims (especially contour changes); rough bone surfaces (especially fossa and articular surface); and surface wear. Proximal femoral components included: articular surface; articular margin; periarticular surfaces; and joint capsule attachment. Femoral components variably revealed: rough bone surface; bone loss; articular margin osteophyte-like features; caudal post-developmental mineralized prominence; and enthesophytes along the joint capsule attachment. Non-metric multidimensional scaling was used to analyze right-left asymmetric relationships between observed traits, across taxa. Significantly different acetabular trait asymmetry involved only C. latrans-C. l. familiaris; V. vulpes-N. procyonoides, and U. cinereoargenteus-N. procyonoides. There were no significant lateralized differences in proximal femoral traits involving modern canids, ancient and modern C. l. familiaris, or modern vulpines. Thus, the observations were strongly bilateral. We hypothesized high similarity of traits across taxa. The data confirm the hypothesis and strongly suggest broad and deep morphological and mechanistic conservation that almost certainly pre-existed (at least) all modern canids. Further zoological studies are needed to evaluate phylogenic implications in greater detail.

Genetic perturbation of PU.1 binding and chromatin looping at neutrophil enhancers associates with autoimmune disease

Neutrophils play fundamental roles in innate immune response, shape adaptive immunity, and are a potentially causal cell type underpinning genetic associations with immune system traits and diseases. Here, we profile the binding of myeloid master regulator PU.1 in primary neutrophils across nearly a hundred volunteers. We show that variants associated with differential PU.1 binding underlie genetically-driven differences in cell count and susceptibility to autoimmune and inflammatory diseases. We integrate these results with other multi-individual genomic readouts, revealing coordinated effects of PU.1 binding variants on the local chromatin state, enhancer-promoter contacts and downstream gene expression, and providing a functional interpretation for 27 genes underlying immune traits. Collectively, these results demonstrate the functional role of PU.1 and its target enhancers in neutrophil transcriptional control and immune disease susceptibility.

Pinpointing miRNA and genes enrichment over trait-relevant tissue network in Genome-Wide Association Studies

BACKGROUND

Understanding gene regulation is important but difficult. Elucidating tissue-specific gene regulation mechanism is even more challenging and requires gene co-expression network assembled from protein-protein interaction, transcription factor and gene binding, and post-transcriptional regulation (e.g., miRNA targeting) information. The miRNA binding affinity could therefore be changed by SNP(s) located at the 3' untranslated regions (3'UTR) of the target messenger RNA (mRNA) which miRNA(s) interacts with. Genome-wide association study (GWAS) has reported significant numbers of loci hosting SNPs associated with many traits. The goal of this study is to pinpoint GWAS functional variants located in 3'UTRs and elucidate if the genes harboring these variants along with their targeting miRNAs are associated with genetic traits relevant to certain tissues.

METHODS

By applying MIGWAS, CoCoNet, ANNOVAR, and DAVID bioinformatics software and utilizing the gene expression database (e.g. GTEx data) to study GWAS summary statistics for 43 traits from 28 GWAS studies, we have identified a list of miRNAs and targeted genes harboring 3'UTR variants, which could contribute to trait-relevant tissue over miRNA-target gene network.

RESULTS

Our result demonstrated that strong association between traits and tissues exists, and in particular, the Primary Biliary Cirrhosis (PBC) trait has the most significant p-value for all 180 tissues among all 43 traits used for this study. We reported SNPs located in 3'UTR regions of genes (SFMBT2, ZC3HAV1, and UGT3A1) targeted by miRNAs for PBC trait and its tissue association network. After employing Gene Ontology (GO) analysis for PBC trait, we have also identified a very important miRNA targeted gene over miRNA-target gene network, PFKL, which encodes the liver subunit of an enzyme.

CONCLUSIONS

The non-coding variants identified from GWAS studies are casually assumed to be not critical to translated protein product. However, 3' untranslated regions (3'UTRs) of genes harbor variants can often change the binding affinity of targeting miRNAs playing important roles in protein translation degree. Our study has shown that GWAS variants could play important roles on miRNA-target gene networks by contributing the association between traits and tissues. Our analysis expands our knowledge on trait-relevant tissue network and paves way for future human disease studies.

Functional interaction between long non-coding RNA and microRNA in rheumatoid arthritis

MicroRNA (miRNA) has received widespread attention for its role in several key cellular processes such as cell differentiation, cell proliferation, apoptosis, and autoimmune diseases. Although we now have a good understanding of miRNA expression and function, our knowledge regarding the molecular mechanism of long non‐coding RNA (lncRNA) is still in its infancy. In this review, we will briefly introduce the definition and function of lncRNA and summarize the interactions between lncRNA and miRNA and their research progress in rheumatoid arthritis (RA). The expression of miR‐16, miR‐146a, miR‐155, and miR‐223 and the interactions between HOTAIR and miR138, ZFAS1 and miR‐27a, and GAPLINC and miR‐575 are representative examples that may augment the understanding of the pathogenesis of RA and help in the development of new biomarkers and target therapies.

Host transcriptome response to Borrelia burgdorferi sensu lato

The host immune response to infection is a well-coordinated system of innate and adaptive immune cells working in concert to prevent the colonization and dissemination of a pathogen. While this typically leads to a beneficial outcome and the suppression of disease pathogenesis, the Lyme borreliosis bacterium, Borrelia burgdorferi sensu lato, can elicit an immune profile that leads to a deleterious state. As B. burgdorferi s.l. produces no known toxins, it is suggested that the immune and inflammatory response of the host are responsible for the manifestation of symptoms, including flu-like symptoms, musculoskeletal pain, and cognitive disorders. The past several years has seen a substantial increase in the use of microarray and sequencing technologies to investigate the transcriptome response induced by B. burgdorferi s.l., thus enabling researchers to identify key factors and pathways underlying the pathophysiology of Lyme borreliosis. In this review we present the major host transcriptional outcomes induced by the bacterium across several studies and discuss the overarching theme of the host inflammatory and immune response, and how it influences the pathology of Lyme borreliosis.

miR-30a-5p Inhibits Epithelial-to-Mesenchymal Transition by Targeting CDK6 in Nasal Polyps

BACKGROUND

Epithelial-to-Mesenchymal Transition (EMT) is considered as a crucial event in disease development and dysregulation of microRNAs (miRNAs) is involved in the regulation of EMT in various human diseases. Emerging evidences congregated over the years have demonstrated that miR-30a-5p was decreased in diseases and its overexpression inhibited the process of diseases via attenuating EMT. Although aberrant expression of miRNAs and occurrence of EMT were previously reported in Nasal Polyps (NPs), the role of miR-30a-5p in EMT of NPs is still remains unclear.

OBJECTIVE

The purpose of our present study was to explore the expression and potential function of miR-30a-5p in EMT of NPs.

METHODS

The expression of miR-30a-5p and mRNA expression level were detected by quantitative real-time PCR (qRT-PCR) in transforming growth factor β1 (TGF-β1) - induced EMT model and NPs patients. Western Blot (WB) and immunohistochemistry (IHC) were performed to evaluate the protein expression level of EMT markers. The cells mobility was assessed by Wound-Healing assay. Luciferase reporter assay was utilized to verify the relationship between Cyclin-dependent kinase 6 (CDK6) and miR-30a-5p.

RESULTS

Firstly, we observed that miR-30a-5p was down-regulated notably, accompanying with the alteration of EMT markers expression in NPs tissues and EMT model induced by TGF-β1 in primary Human Nasal Epithelial Cells (pHNECs) and A549 cells in vitro. Moreover, the functional assays demonstrated that overexpression of miR-30a-5p significantly inhibited EMT and cells mobility. Subsequently, CDK6 was validated as a direct target of miR-30a-5p. Finally, we performed the rescue experiments indicating that overexpression of CDK6 eliminated the suppressive effects of miR-30a-5p in TGF-β1-induced EMT in pHNECs and A549 cells.

CONCLUSION

Taken together, our results suggested that EMT was involved in NPs, and overexpression of miR-30a-5p could attenuate EMT via repressing the expression of the CDK6 in pHNECs and A549 cells.

Genetic influences on susceptibility to rheumatoid arthritis in African-Americans

Large meta-analyses of rheumatoid arthritis (RA) susceptibility in European (EUR) and East Asian (EAS) populations have identified >100 RA risk loci, but genome-wide studies of RA in African-Americans (AAs) are absent. To address this disparity, we performed an analysis of 916 AA RA patients and 1392 controls and aggregated our data with genotyping data from >100 000 EUR and Asian RA patients and controls. We identified two novel risk loci that appear to be specific to AAs: GPC5 and RBFOX1 (PAA < 5 × 10-9). Most RA risk loci are shared across different ethnicities, but among discordant loci, we observed strong enrichment of variants having large effect sizes. We found strong evidence of effect concordance for only 3 of the 21 largest effect index variants in EURs. We used the trans-ethnic fine-mapping algorithm PAINTOR3 to prioritize risk variants in >90 RA risk loci. Addition of AA data to those of EUR and EAS descent enabled identification of seven novel high-confidence candidate pathogenic variants (defined by posterior probability > 0.8). In summary, our trans-ethnic analyses are the first to include AAs, identified several new RA risk loci and point to candidate pathogenic variants that may underlie this common autoimmune disease. These findings may lead to better ways to diagnose or stratify treatment approaches in RA.

A joint effort: The interplay between the innate and the adaptive immune system in Lyme arthritis

Articular joints are a major target of Borrelia burgdorferi, the causative agent of Lyme arthritis. Despite antibiotic treatment, recurrent or persistent Lyme arthritis is observed in a significant number of patients. The host immune response plays a crucial role in this chronic arthritic joint complication of Borrelia infections. During the early stages of B. burgdorferi infection, a major hinder in generating a proper host immune response is the lack of induction of a strong adaptive immune response. This may lead to a delayed hyperinflammatory reaction later in the disease. Several mechanisms have been suggested that might be pivotal for the development of Lyme arthritis and will be highlighted in this review, from molecular mimicry of matrix metallopeptidases and glycosaminoglycans, to autoimmune responses to live bacteria, or remnants of Borrelia spirochetes in joints. Murine studies have suggested that the inflammatory responses are initiated by innate immune cells, but this does not exclude the involvement of the adaptive immune system in this dysregulated immune profile. Genetic predisposition, via human leukocyte antigen-DR isotype and microRNA expression, has been associated with the development of antibiotic-refractory Lyme arthritis. Yet the ultimate cause for (antibiotic-refractory) Lyme arthritis remains unknown. Complex processes of different immune cells and signaling cascades are involved in the development of Lyme arthritis. When these various mechanisms are fully been unraveled, new treatment strategies can be developed to target (antibiotic-refractory) Lyme arthritis more effectively.

The Impact of Cigarette Smoking on Risk of Rheumatoid Arthritis: A Narrative Review

Rheumatoid arthritis (RA) is an autoimmune disease characterized by chronic inflammation and subsequent proliferation of synovial tissues, which eventually leads to cartilage and bone destruction without effective treatments. Anti-citrullinated cyclic peptide/protein antibody (ACPA) and rheumatoid factor (RF) are two main characteristic autoantibodies found in RA patients and are associated with unfavorable disease outcomes. Although etiologies and causes of the disease have not been fully clarified yet, it is likely that interactive contributions of genetic and environmental factors play a main role in RA pathology. Previous works have demonstrated several genetic and environmental factors as risks of RA development and/or autoantibody productions. Among these, cigarette smoking and HLA-DRB1 are the well-established environmental and genetic risks, respectively. In this narrative review, we provide a recent update on genetic contributions to RA and the environmental risks of RA with a special focus on cigarette smoking and its impacts on RA pathology. We also describe gene–environmental interaction in RA pathogenesis with an emphasis on cigarette smoking and HLA-DRB1.

Advances in genetics toward identifying pathogenic cell states of rheumatoid arthritis

Rheumatoid arthritis (RA) risk has a large genetic component (~60%) that is still not fully understood. This has hampered the design of effective treatments that could promise lifelong remission. RA is a polygenic disease with 106 known genome-wide significant associated loci and thousands of small effect causal variants. Our current understanding of RA risk has suggested cell-type-specific contexts for causal variants, implicating CD4 + effector memory T cells, as well as monocytes, B cells and stromal fibroblasts. While these cellular states and categories are still mechanistically broad, future studies may identify causal cell subpopulations. These efforts are propelled by advances in single cell profiling. Identification of causal cell subpopulations may accelerate therapeutic intervention to achieve lifelong remission.

Grimon: graphical interface to visualize multi-omics networks

Summary

Rapid advances in high-throughput sequencing technologies have enabled more efficient acquisition of massive amount of multi-omics data. However, interpretation of the underlying relationships across multi-omics networks has not been fully succeeded, partly due to the lack of effective methods in visualization. To aid interpretation of the results from such multi-omics data, we here present Grimon (Graphical interface to visualize multi-omics networks), an R package that visualizes high-dimensional multi-layered data sets in three-dimensional parallel coordinates. Grimon enables users to intuitively and interactively explore their analyzed data, helping their understanding of multiple inter-layer connections embedded in high-dimensional complex data.

Availability and implementation

Grimon is freely available at https://github.com/mkanai/grimon as an R package with example omics data sets.

Supplementary information

Supplementary data are available at bioinformatics online.

Complex Phenotypes: Mechanisms Underlying Variation in Human Stature

PURPOSE OF REVIEW

The goal of the review is to provide a comprehensive overview of the current understanding of the mechanisms underlying variation in human stature.

RECENT FINDINGS

Human height is an anthropometric trait that varies considerably within human populations as well as across the globe. Historically, much research focus was placed on understanding the biology of growth plate chondrocytes and how modifications to core chondrocyte proliferation and differentiation pathways potentially shaped height attainment in normal as well as pathological contexts. Recently, much progress has been made to improve our understanding regarding the mechanisms underlying the normal and pathological range of height variation within as well as between human populations, and today, it is understood to reflect complex interactions among a myriad of genetic, environmental, and evolutionary factors. Indeed, recent improvements in genetics (e.g., GWAS) and breakthroughs in functional genomics (e.g., whole exome sequencing, DNA methylation analysis, ATAC-sequencing, and CRISPR) have shed light on previously unknown pathways/mechanisms governing pathological and common height variation. Additionally, the use of an evolutionary perspective has also revealed important mechanisms that have shaped height variation across the planet. This review provides an overview of the current knowledge of the biological mechanisms underlying height variation by highlighting new research findings on skeletal growth control with an emphasis on previously unknown pathways/mechanisms influencing pathological and common height variation. In this context, this review also discusses how evolutionary forces likely shaped the genomic architecture of height across the globe.

A Polymer Physics Framework for the Entropy of Arbitrary Pseudoknots

The accurate prediction of RNA secondary structure from primary sequence has had enormous impact on research from the past 40 years. Although many algorithms are available to make these predictions, the inclusion of non-nested loops, termed pseudoknots, still poses challenges arising from two main factors: 1) no physical model exists to estimate the loop entropies of complex intramolecular pseudoknots, and 2) their NP-complete enumeration has impeded their study. Here, we address both challenges. First, we develop a polymer physics model that can address arbitrarily complex pseudoknots using only two parameters corresponding to concrete physical quantities-over an order of magnitude fewer than the sparsest state-of-the-art phenomenological methods. Second, by coupling this model to exhaustive enumeration of the set of possible structures, we compute the entire free energy landscape of secondary structures resulting from a primary RNA sequence. We demonstrate that for RNA structures of ∼80 nucleotides, with minimal heuristics, the complete enumeration of possible secondary structures can be accomplished quickly despite the NP-complete nature of the problem. We further show that despite our loop entropy model's parametric sparsity, it performs better than or on par with previously published methods in predicting both pseudoknotted and non-pseudoknotted structures on a benchmark data set of RNA structures of ≤80 nucleotides. We suggest ways in which the accuracy of the model can be further improved.

Integration of genetics and miRNA-target gene network identified disease biology implicated in tissue specificity

MicroRNAs (miRNAs) modulate the post-transcriptional regulation of target genes and are related to biology of complex human traits, but genetic landscape of miRNAs remains largely unknown. Given the strikingly tissue-specific miRNA expression profiles, we here expand a previous method to quantitatively evaluate enrichment of genome-wide association study (GWAS) signals on miRNA-target gene networks (MIGWAS) to further estimate tissue-specific enrichment. Our approach integrates tissue-specific expression profiles of miRNAs (∼1800 miRNAs in 179 cells) with GWAS to test whether polygenic signals enrich in miRNA-target gene networks and whether they fall within specific tissues. We applied MIGWAS to 49 GWASs (nTotal = 3 520 246), and successfully identified biologically relevant tissues. Further, MIGWAS could point miRNAs as candidate biomarkers of the trait. As an illustrative example, we performed differentially expressed miRNA analysis between rheumatoid arthritis (RA) patients and healthy controls (n = 63). We identified novel biomarker miRNAs (e.g. hsa-miR-762) by integrating differentially expressed miRNAs with MIGWAS results for RA, as well as novel associated loci with significant genetic risk (rs56656810 at MIR762 at 16q11; n = 91 482, P = 3.6 × 10-8). Our result highlighted that miRNA-target gene network contributes to human disease genetics in a cell type-specific manner, which could yield an efficient screening of miRNAs as promising biomarkers.

miRNA contributions to pediatric-onset multiple sclerosis inferred from GWAS

OBJECTIVE

Onset of multiple sclerosis (MS) occurs in childhood for approximately 5% of cases (pediatric MS, or ped-MS). Epigenetic influences are strongly implicated in MS pathogenesis in adults, including the contribution from microRNAs (miRNAs), small noncoding RNAs that affect gene expression by binding target gene mRNAs. Few studies have specifically examined miRNAs in ped-MS, but individuals developing MS at an early age may carry a relatively high burden of genetic risk factors, and miRNA dysregulation may therefore play a larger role in the development of ped-MS than in adult-onset MS. This study aimed to look for evidence of miRNA involvement in ped-MS pathogenesis.

METHODS

GWAS results from 486 ped-MS cases and 1362 controls from the U.S. Pediatric MS Network and Kaiser Permanente Northern California membership were investigated for miRNA-specific signals. First, enrichment of miRNA-target gene network signals was evaluated using MIGWAS software. Second, SNPs in miRNA genes and in target gene binding sites (miR-SNPs) were tested for association with ped-MS, and pathway analysis was performed on associated target genes.

RESULTS

MIGWAS analysis showed that miRNA-target gene signals were enriched in GWAS (P = 0.038) and identified 39 candidate biomarker miRNA-target gene pairs, including immune and neuronal signaling genes. The miR-SNP analysis implicated dysregulation of miRNA binding to target genes in five pathways, mainly involved in immune signaling.

INTERPRETATION

Evidence from GWAS suggests that miRNAs play a role in ped-MS pathogenesis by affecting immune signaling and other pathways. Candidate biomarker miRNA-target gene pairs should be further studied for diagnostic, prognostic, and/or therapeutic utility.

Genetics of rheumatoid arthritis: 2018 status

Study of the genetics of rheumatoid arthritis (RA) began about four decades ago with the discovery of HLA-DRB1 Since the beginning of this century, a number of non-HLA risk loci have been identified through genome-wide association studies (GWAS). We now know that over 100 loci are associated with RA risk. Because genetic information implies a clear causal relationship to the disease, research into the pathogenesis of RA should be promoted. However, only 20% of GWAS loci contain coding variants, with the remaining variants occurring in non-coding regions, and therefore, the majority of causal genes and causal variants remain to be identified. The use of epigenetic studies, high-resolution mapping of open chromatin, chromosomal conformation technologies and other approaches could identify many of the missing links between genetic risk variants and causal genetic components, thus expanding our understanding of RA genetics.

Noncoding RNAs in multiple sclerosis

Multiple sclerosis (MS), a chronic inflammatory demyelinating disease of the central nervous system, is characterized by axonal degeneration and gliosis. Although the causes of MS remain unknown, gene dysregulation in the central nervous system has been associated with the disease pathogenesis. As such, the various regulators of gene expression may be contributing factors. The noncoding (nc) RNAs have piqued the interest of MS researchers due to their known functions in human physiology and various pathological processes, despite being generally characterized as transcripts without apparent protein-coding capacity. Accumulating evidence has indicated that ncRNAs participate in the regulation of MS by acting as epigenetic factors, especially the long (l) ncRNAs and the micro (mi) RNAs, and they are now recognized as key regulatory molecules in MS. In this review, we summarize the most current studies on the contribution of ncRNAs in MS pathogenic processes and discuss their potential applications in the diagnosis and treatment of MS.

The clinical significance of miR-335, miR-124, miR-218 and miR-484 downregulation in gastric cancer

Gastric cancer (GC) is one of the leading types of malignancy worldwide, particularly in Asian populations. Although the exact molecular mechanism of GC development remains unknown, microRNA (miRNA) has recently been shown to be involved. The current study aims to investigate the expression levels of bioinformatically ranked miRNAs in gastric tissues. Using bioinformatics tools, we prioritized miRNAs thought to be implicated in GC. Furthermore, polyA-qPCR was used to validate bioinformatics findings in 40 GC, 31 normal gastric tissue (NG) and 45 gastric dysplasia (GD) samples. As identified by bioinformatics analysis, miR-335 was shown to be the top-ranked miRNA implicated in GC. Moreover, a significant downregulation of miR-335, miR-124, miR-218 and miR-484 was found in GC and GD compared to NG samples. We found bioinformatics to be an efficient approach to finding candidate miRNAs relevant to GC development. Finally, the findings show that downregulation of miRNAs such as miR-124 and miR-218 in gastric tissue can be a significant indicator for neoplastic transformation.

Future Directions of Genomics Research in Rheumatic Diseases

Recent developments in human genome genotyping and sequencing technologies, such as genome-wide association studies and whole-genome sequencing analyses, have successfully identified several risk genes of rheumatic diseases. Fine-mapping studies using the HLA imputation method revealed that classical and non-classical HLA genes contribute to the risk of rheumatic diseases. Integration of human disease genomics with biological, medical, and clinical databases should contribute to the elucidation of disease pathogenicity and novel drug discovery. Disease risk genes identified by large-scale genetic studies are considered to be promising resources for novel drug discovery, including drug repositioning and biomarker microRNA screening for rheumatoid arthritis.

MicroRNA Expression Shows Inflammatory Dysregulation and Tumor-Like Proliferative Responses in Joints of Patients With Postinfectious Lyme Arthritis

OBJECTIVE

Lyme arthritis (LA) is caused by infection with Borrelia burgdorferi and usually resolves following spirochetal killing with antibiotics. However, in some patients, arthritis persists after antibiotic therapy. To provide insights into underlying pathogenic processes associated with antibiotic-refractory LA (postinfectious LA), we analyzed differences in microRNA (miRNA) expression between LA patients with active infection and those with postinfectious LA.

METHODS

MicroRNA expression was assayed in synovial fluid (SF) from LA patients before and after oral and intravenous antibiotic therapy, and in synovial tissue obtained months after antibiotic therapy from patients with postinfectious LA. SF and tissue from patients with other forms of arthritis, such as rheumatoid arthritis (RA) and osteoarthritis, were used for comparison.

RESULTS

SF from LA patients during active infection had marked elevations of white blood cells, particularly polymorphonuclear leukocytes, accompanied by elevated levels of microRNA-223 (miR-223). In contrast, SF from postantibiotic LA patients contained greater percentages of lymphocytes and mononuclear cells. SF from postantibiotic LA patients also exhibited marked inflammatory (miR-146a, miR-155), wound repair (miR-142), and proliferative (miR-17-92) miRNA signatures, and higher levels of these miRNAs correlated with longer arthritis duration. Levels of miR-146a, miR-155, miR-142, miR-223, and miR-17-92 were also elevated in synovial tissue in late postinfectious LA, and levels of let-7a were reduced, similar to RA.

CONCLUSION

During active infection, miRNA expression in SF reflected an immune response associated with bacterial killing, while in postinfectious LA, miRNA expression in SF and synovial tissue reflected chronic inflammation, synovial proliferation, and breakdown of wound repair processes, showing that the nature of the arthritis was altered after spirochetal killing.

Cellular network perturbations by disease-associated variants

Genetic and genome-wide association studies (GWAS) have identified a myriad of human disease-associated genomic variants. However, these studies do not reveal the mechanisms by which these variants perturb cellular networks, a necessary step to intervene and improve disease outcomes. This has been challenging because multiple variants are present in haplotype blocks, thereby confounding the identification of causal variants, and because most reside in noncoding regions. Here, we review recent advances in the identification of functional variants and gene-variant associations. In addition, we examine approaches used to study perturbations in protein-protein and protein-DNA interactions associated with disease, and discuss how these perturbations affect cellular networks.

Update on the genetic architecture of rheumatoid arthritis

Human genetic studies into rheumatoid arthritis (RA) have uncovered more than 100 genetic loci associated with susceptibility to RA and have refined the RA-association model for HLA variants. The majority of RA-risk variants are highly shared across multiple ancestral populations and are located in noncoding elements that might have allele-specific regulatory effects in relevant tissues. Emerging multi-omics data, high-density genotype data and bioinformatic approaches are enabling researchers to use RA-risk variants to identify functionally relevant cell types and biological pathways that are involved in impaired immune processes and disease phenotypes. This Review summarizes reported RA-risk loci and the latest insights from human genetic studies into RA pathogenesis, including how genetic data has helped to identify currently available drugs that could be repurposed for patients with RA and the role of genetics in guiding the development of new drugs.

miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients

PURPOSE

The proper validation of prognostic biomarkers is an important clinical issue in breast cancer research. MicroRNAs (miRNAs) have emerged as a new class of promising breast cancer biomarkers. In the present work, we developed an integrated online bioinformatic tool to validate the prognostic relevance of miRNAs in breast cancer.

METHODS

A database was set up by searching the GEO, EGA, TCGA, and PubMed repositories to identify datasets with published miRNA expression and clinical data. Kaplan-Meier survival analysis was performed to validate the prognostic value of a set of 41 previously published survival-associated miRNAs.

RESULTS

All together 2178 samples from four independent datasets were integrated into the system including the expression of 1052 distinct human miRNAs. In addition, the web-tool allows for the selection of patients, which can be filtered by receptors status, lymph node involvement, histological grade, and treatments. The complete analysis tool can be accessed online at: www.kmplot.com/mirpower . We used this tool to analyze a large number of deregulated miRNAs associated with breast cancer features and outcome, and confirmed the prognostic value of 26 miRNAs. A significant correlation in three out of four datasets was validated only for miR-29c and miR-101.

CONCLUSIONS

In summary, we established an integrated platform capable to mine all available miRNA data to perform a survival analysis for the identification and validation of prognostic miRNA markers in breast cancer.

A DNA Circuit for IsomiR Detection

A synthetic DNA oligonucleotide has been programmed to function as a biological circuit to detect 5'-IsomiRs. The circuit consists of two integrated DNA switches. The first is "activated" when a DNA probe is enzymatically modified by a reverse transcriptase that incorporates nucleotides complementary to the 5'-region of a microRNA (miRNA). Addition of the correct number and sequence of nucleotides enables the probe to assemble into an asymmetric DNA hairpin. The reconfigured hairpin probe then primes an internal polymerisation reaction, resulting in the synthesis of a symmetrical DNA hairpin. This activates the second switch, which then initiates the amplification of reverse-transcribed miRNA through a continuous cycle of DNA nicking and polymerisation. The DNA circuit enables sensitive and rapid detection of femtomoles of a miRNA transcript under isothermal conditions.

MiR-30a-5p/UBE3C axis regulates breast cancer cell proliferation and migration

Aberrant expression of ubiquitin Protein ligase E3C (UBE3C) has been documented in breast cancer (BC). MicroRNAs (miRNAs) were shown to play an important role in the regulation of tumor properties in BC. However, whether miRNAs contributes to UBE3C expression in BC cells remains poorly understood. In this study, we report that UBE3C was a direct target of miR-30a-5p. Expression of miR-30a-5p in BC cells reduced UBE3C expression. MCF-7 and MDA-MB-453 cells were transfected miR-30a-5p-overexpression, and found that cell proliferation and migration were inhibited. In contrast, when miR-30a-5p inhibitor were transfected into MCF-7 and MDA-MB-453 cells, cell proliferation and migration were promoted. We study demonstrated that upregulation of miR-30a-5p was significantly suppressed levels of cyclin B1, cyclin D1 and c-myc. Moreover, Correlation analysis indicated that expression of miR-30a-5p was highly negatively correlated with UBE3C, which was upregulated in BC specimens. These data highlight the important role of miR-30a-5p/UBE3C axis in BC development and progression. Therefore, miR-30a-5p activation or UBE3C inhibition may be provide a novel strategy for the treatment of BC.