shaPRS: Leveraging shared genetic effects across traits or ancestries improves accuracy of polygenic scores.

We present shaPRS, a method that leverages widespread pleiotropy between traits or shared genetic effects across ancestries, to improve the accuracy of polygenic scores. The method uses genome-wide summary statistics from two diseases or ancestries to improve the genetic effect estimate and standard error at SNPs where there is homogeneity of effect between the two datasets. When there is significant evidence of heterogeneity, the genetic effect from the disease or population closest to the target population is maintained. We show via simulation and a series of real-world examples that shaPRS substantially enhances the accuracy of polygenic risk scores (PRSs) for complex diseases and greatly improves PRS performance across ancestries. shaPRS is a PRS pre-processing method that is agnostic to the actual PRS generation method, and as a result, it can be integrated into existing PRS generation pipelines and continue to be applied as more performant PRS methods are developed over time.

Noncoding RNA and its associated proteins as regulatory elements of the immune system.

The rapid changes in gene expression that accompany developmental transitions, stress responses and proliferation are controlled by signal-mediated coordination of transcriptional and post-transcriptional mechanisms. In recent years, understanding of the mechanics of these processes and the contexts in which they are employed during hematopoiesis and immune challenge has increased. An important aspect of this progress is recognition of the importance of RNA-binding proteins and noncoding RNAs. These have roles in the development and function of the immune system and in pathogen life cycles, and they represent an important aspect of intracellular immunity.

EPISPOT: An epigenome-driven approach for detecting and interpreting hotspots in molecular QTL studies.

We present EPISPOT, a fully joint framework which exploits large panels of epigenetic annotations as variant-level information to enhance molecular quantitative trait locus (QTL) mapping. Thanks to a purpose-built Bayesian inferential algorithm, EPISPOT accommodates functional information for both cis and trans actions, including QTL hotspot effects. It effectively couples simultaneous QTL analysis of thousands of genetic variants and molecular traits with hypothesis-free selection of biologically interpretable annotations which directly contribute to the QTL effects. This unified, epigenome-aided learning boosts statistical power and sheds light on the regulatory basis of the uncovered hits; EPISPOT therefore marks an essential step toward improving the challenging detection and functional interpretation of trans-acting genetic variants and hotspots. We illustrate the advantages of EPISPOT in simulations emulating real-data conditions and in a monocyte expression QTL study, which confirms known hotspots and finds other signals, as well as plausible mechanisms of action. In particular, by highlighting the role of monocyte DNase-I sensitivity sites from >150 epigenetic annotations, we clarify the mediation effects and cell-type specificity of major hotspots close to the lysozyme gene. Our approach forgoes the daunting and underpowered task of one-annotation-at-a-time enrichment analyses for prioritizing cis and trans QTL hits and is tailored to any transcriptomic, proteomic, or metabolomic QTL problem. By enabling principled epigenome-driven QTL mapping transcriptome-wide, EPISPOT helps progress toward a better functional understanding of genetic regulation.

BBmix: a Bayesian beta-binomial mixture model for accurate genotyping from RNA-sequencing.

MOTIVATION: While many pipelines have been developed for calling genotypes using RNA-sequencing (RNA-Seq) data, they all have adapted DNA genotype callers that do not model biases specific to RNA-Seq such as allele-specific expression (ASE). RESULTS: Here, we present Bayesian beta-binomial mixture model (BBmix), a Bayesian beta-binomial mixture model that first learns the expected distribution of read counts for each genotype, and then deploys those learned parameters to call genotypes probabilistically. We benchmarked our model on a wide variety of datasets and showed that our method generally performed better than competitors, mainly due to an increase of up to 1.4% in the accuracy of heterozygous calls, which may have a big impact in reducing false positive rate in applications sensitive to genotyping error such as ASE. Moreover, BBmix can be easily incorporated into standard pipelines for calling genotypes. We further show that parameters are generally transferable within datasets, such that a single learning run of less than 1 h is sufficient to call genotypes in a large number of samples. AVAILABILITY AND IMPLEMENTATION: We implemented BBmix as an R package that is available for free under a GPL-2 licence at https://gitlab.com/evigorito/bbmix and https://cran.r-project.org/package=bbmix with accompanying pipeline at https://gitlab.com/evigorito/bbmix_pipeline.

RápidoPGS: a rapid polygenic score calculator for summary GWAS data without a test dataset.

MOTIVATION: Polygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RápidoPGS, a flexible and fast method to compute PGS requiring summary-level Genome-wide association studies (GWAS) datasets only, with little computational requirements and no test data required for parameter tuning. RESULTS: We show that RápidoPGS performs slightly less well than two out of three other widely used PGS methods (LDpred2, PRScs and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042 and 0.0064, respectively, but up to 17 000-fold faster with reduced computational requirements. RápidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog. AVAILABILITY AND IMPLEMENTATION: Our method is available with a GPL license as an R package from CRAN and GitHub. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

HVCN1 modulates BCR signal strength via regulation of BCR-dependent generation of reactive oxygen species.

Voltage-gated proton currents regulate generation of reactive oxygen species (ROS) in phagocytic cells. In B cells, stimulation of the B cell antigen receptor (BCR) results in the production of ROS that participate in B cell activation, but the involvement of proton channels is unknown. We report here that the voltage-gated proton channel HVCN1 associated with the BCR complex and was internalized together with the BCR after activation. BCR-induced generation of ROS was lower in HVCN1-deficient B cells, which resulted in attenuated BCR signaling via impaired BCR-dependent oxidation of the tyrosine phosphatase SHP-1. This resulted in less activation of the kinases Syk and Akt, impaired mitochondrial respiration and glycolysis and diminished antibody responses in vivo. Our findings identify unanticipated functions for proton channels in B cells and demonstrate the importance of ROS in BCR signaling and downstream metabolism.

Circadian control of innate immunity in macrophages by miR-155 targeting Bmal1.

The response to an innate immune challenge is conditioned by the time of day, but the molecular basis for this remains unclear. In myeloid cells, there is a temporal regulation to induction by lipopolysaccharide (LPS) of the proinflammatory microRNA miR-155 that correlates inversely with levels of BMAL1. BMAL1 in the myeloid lineage inhibits activation of NF-κB and miR-155 induction and protects mice from LPS-induced sepsis. Bmal1 has two miR-155-binding sites in its 3'-UTR, and, in response to LPS, miR-155 binds to these two target sites, leading to suppression of Bmal1 mRNA and protein in mice and humans. miR-155 deletion perturbs circadian function, gives rise to a shorter circadian day, and ablates the circadian effect on cytokine responses to LPS. Thus, the molecular clock controls miR-155 induction that can repress BMAL1 directly. This leads to an innate immune response that is variably responsive to challenges across the circadian day.

miR-155: an ancient regulator of the immune system.

MicroRNAs (miRNAs) are a newly recognized class of regulatory genes which repress the expression of protein-coding genes. Numerous studies have uncovered a complex role for miRNAs regulating many aspects of a variety of cellular processes including cell growth, differentiation, and lineage commitment. In the immune system, miR-155 is unique in its ability to shape the transcriptome of activated myeloid and lymphoid cells controlling diverse biological functions ranging from inflammation to immunological memory. Not surprisingly, a tight control of miR-155 expression is required to avoid malignant transformation, as evidenced by miR-155 overexpression in many cancers of B-cell origin. In this review, we discuss the potential of miR-155 as a molecular target for therapeutic intervention and discuss the function of miR-155 in the context of protective immunity. We first look back into the emergence of miR-155 in evolution, which is coincidental with the emergence of the ancestors of the antigen receptors. We then summarize what we have learned about the role of miR-155 in the regulation of lymphoid subsets at the cellular and molecular level in the context of recent progress in this field.

MiR-210 is induced by Oct-2, regulates B cells, and inhibits autoantibody production.

MicroRNAs (MiRs) are small, noncoding RNAs that regulate gene expression posttranscriptionally. In this study, we show that MiR-210 is induced by Oct-2, a key transcriptional mediator of B cell activation. Germline deletion of MiR-210 results in the development of autoantibodies from 5 mo of age. Overexpression of MiR-210 in vivo resulted in cell autonomous expansion of the B1 lineage and impaired fitness of B2 cells. Mice overexpressing MiR-210 exhibited impaired class-switched Ab responses, a finding confirmed in wild-type B cells transfected with a MiR-210 mimic. In vitro studies demonstrated defects in cellular proliferation and cell cycle entry, which were consistent with the transcriptomic analysis demonstrating downregulation of genes involved in cellular proliferation and B cell activation. These findings indicate that Oct-2 induction of MiR-210 provides a novel inhibitory mechanism for the control of B cells and autoantibody production.

Macrophage microRNA-155 promotes cardiac hypertrophy and failure.

BACKGROUND: Cardiac hypertrophy and subsequent heart failure triggered by chronic hypertension represent major challenges for cardiovascular research. Beyond neurohormonal and myocyte signaling pathways, growing evidence suggests inflammatory signaling pathways as therapeutically targetable contributors to this process. We recently reported that microRNA-155 is a key mediator of cardiac inflammation and injury in infectious myocarditis. Here, we investigated the impact of microRNA-155 manipulation in hypertensive heart disease. METHODS AND RESULTS: Genetic loss or pharmacological inhibition of the leukocyte-expressed microRNA-155 in mice markedly reduced cardiac inflammation, hypertrophy, and dysfunction on pressure overload. These alterations were macrophage dependent because in vivo cardiomyocyte-specific microRNA-155 manipulation did not affect cardiac hypertrophy or dysfunction, whereas bone marrow transplantation from wild-type mice into microRNA-155 knockout animals rescued the hypertrophic response of the cardiomyocytes and vice versa. In vitro, media from microRNA-155 knockout macrophages blocked the hypertrophic growth of stimulated cardiomyocytes, confirming that macrophages influence myocyte growth in a microRNA-155-dependent paracrine manner. These effects were at least partly mediated by the direct microRNA-155 target suppressor of cytokine signaling 1 (Socs1) because Socs1 knockdown in microRNA-155 knockout macrophages largely restored their hypertrophy-stimulating potency. CONCLUSIONS: Our findings reveal that microRNA-155 expression in macrophages promotes cardiac inflammation, hypertrophy, and failure in response to pressure overload. These data support the causative significance of inflammatory signaling in hypertrophic heart disease and demonstrate the feasibility of therapeutic microRNA targeting of inflammation in heart failure.

Chemical structure requirements and cellular targeting of microRNA-122 by peptide nucleic acids anti-miRs.

Anti-miRs are oligonucleotide inhibitors complementary to miRNAs that have been used extensively as tools to gain understanding of specific miRNA functions and as potential therapeutics. We showed previously that peptide nucleic acid (PNA) anti-miRs containing a few attached Lys residues were potent miRNA inhibitors. Using miR-122 as an example, we report here the PNA sequence and attached amino acid requirements for efficient miRNA targeting and show that anti-miR activity is enhanced substantially by the presence of a terminal-free thiol group, such as a Cys residue, primarily due to better cellular uptake. We show that anti-miR activity of a Cys-containing PNA is achieved by cell uptake through both clathrin-dependent and independent routes. With the aid of two PNA analogues having intrinsic fluorescence, thiazole orange (TO)-PNA and [bis-o-(aminoethoxy)phenyl]pyrrolocytosine (BoPhpC)-PNA, we explored the subcellular localization of PNA anti-miRs and our data suggest that anti-miR targeting of miR-122 may take place in or associated with endosomal compartments. Our findings are valuable for further design of PNAs and other oligonucleotides as potent anti-miR agents.

Enhanced susceptibility to Citrobacter rodentium infection in microRNA-155-deficient mice.

MicroRNAs (miRNAs) are small noncoding molecules that control gene expression posttranscriptionally, with microRNA-155 (miR-155) one of the first to be implicated in immune regulation. Here, we show that miR-155-deficient mice are less able to eradicate a mucosal Citrobacter rodentium infection than wild-type C57BL/6 mice. miR-155-deficient mice exhibited prolonged colonization associated with a higher C. rodentium burden in gastrointestinal tissue and spread into systemic tissues. Germinal center formation and humoral immune responses against C. rodentium were severely impaired in infected miR-155-deficient mice. A similarly susceptible phenotype was observed in µMT mice reconstituted with miR-155-deficient B cells, indicating that miR-155 is required intrinsically for mediating protection against this predominantly luminal bacterial pathogen.

MicroRNA fate upon targeting with anti-miRNA oligonucleotides as revealed by an improved Northern-blot-based method for miRNA detection.

MicroRNAs (miRNAs) are small non-coding RNAs involved in fine-tuning of gene regulation. Antisense oligonucleotides (ONs) are promising tools as anti-miRNA (anti-miR) agents toward therapeutic applications and to uncover miRNA function. Such anti-miR ONs include 2'-O-methyl (OMe), cationic peptide nucleic acids like K-PNA-K3, and locked nucleic acid (LNA)-based anti-miRs such as LNA/DNA or LNA/OMe. Northern blotting is a widely used and robust technique to detect miRNAs. However, miRNA quantification in the presence of anti-miR ONs has proved to be challenging, due to detection artifacts, which has led to poor understanding of miRNA fate upon anti-miR binding. Here we show that anti-miR ON bound to miR-122 can prevent the miRNA from being properly precipitated into the purified RNA fraction using the standard RNA extraction protocol (TRI-Reagent), yielding an RNA extract that does not reflect the real cellular levels of the miRNA. An increase in the numbers of equivalents of isopropanol during the precipitation step leads to full recovery of the targeted miRNA back into the purified RNA extract. Following our improved protocol, we demonstrate by Northern blotting, in conjunction with a PNA decoy strategy and use of high denaturing PAGE, that high-affinity anti-miRs (K-PNA-K3, LNA/DNA, and LNA/OMe) sequester miR-122 without causing miRNA degradation, while miR-122 targeting with a lower-affinity anti-miR (OMe) seems to promote degradation of the miRNA. The technical issues explored in this work will have relevance for other hybridization-based techniques for miRNA quantification in the presence of anti-miR ONs.

Efficient inhibition of miR-155 function in vivo by peptide nucleic acids.

MicroRNAs (miRNAs) play an important role in diverse physiological processes and are potential therapeutic agents. Synthetic oligonucleotides (ONs) of different chemistries have proven successful for blocking miRNA expression. However, their specificity and efficiency have not been fully evaluated. Here, we show that peptide nucleic acids (PNAs) efficiently block a key inducible miRNA expressed in the haematopoietic system, miR-155, in cultured B cells as well as in mice. Remarkably, miR-155 inhibition by PNA in primary B cells was achieved in the absence of any transfection agent. In mice, the high efficiency of the treatment was demonstrated by a strong overlap in global gene expression between B cells isolated from anti-miR-155 PNA-treated and miR-155-deficient mice. Interestingly, PNA also induced additional changes in gene expression. Our analysis provides a useful platform to aid the design of efficient and specific anti-miRNA ONs for in vivo use.

Cutting edge: the Foxp3 target miR-155 contributes to the development of regulatory T cells.

Foxp3 is a transcription factor that is essential for the normal development of regulatory T cells (Tregs). In the absence of microRNAs (miRNAs), Foxp3(+) Tregs develop but fail to maintain immune homeostasis, leading to a scurfy-like disease. Global analysis of the network of genes regulated by Foxp3 has identified the miRNA miR-155, which is highly expressed in Tregs, as a direct target of Foxp3. In this study we report that miR-155-deficient mice have reduced numbers of Tregs, both in the thymus and periphery, due to impaired development. However, we found no evidence for defective suppressor activity of miR-155-deficient Tregs, either in vitro or in vivo. Our results indicate that miR-155 contributes to Treg development, but that additional miRNAs control Treg function.

Detection of quantitative trait loci from RNA-seq data with or without genotypes using BaseQTL.

Detecting genetic variants associated with traits (quantitative trait loci, QTL) requires genotyped study individuals. Here we describe BaseQTL, a Bayesian method that exploits allele-specific expression to map molecular QTL from sequencing reads (eQTL for gene expression) even when no genotypes are available. When used with genotypes to map eQTL, BaseQTL has lower error rates and increased power compared with existing QTL mapping methods. Running without genotypes limits how many tests can be performed, but due to the proximity of QTL variants to gene bodies, the 2.8% of variants within a 100 kB window that could be tested contained 26% of eQTL detectable with genotypes. eQTL effect estimates were invariably consistent between analyses performed with and without genotypes. Often, sequencing data may be generated in the absence of genotypes on patients and controls in differential expression studies, and we identified an apparent psoriasis-specific eQTL for GSTP1 in one such dataset, providing new insights into disease-dependent gene regulation.

A crucial role for the p110delta subunit of phosphatidylinositol 3-kinase in B cell development and activation.

Mice lacking the p110delta catalytic subunit of phosphatidylinositol 3-kinase have reduced numbers of B1 and marginal zone B cells, reduced levels of serum immunoglobulins, respond poorly to immunization with type II thymus-independent antigen, and are defective in their primary and secondary responses to thymus-dependent antigen. p110delta(-/-) B cells proliferate poorly in response to B cell receptor (BCR) or CD40 signals in vitro, fail to activate protein kinase B, and are prone to apoptosis. p110delta function is required for BCR-mediated calcium flux, activation of phosphlipaseCgamma2, and Bruton's tyrosine kinase. Thus, p110delta plays a critical role in B cell homeostasis and function.

Vav GEFs regulate macrophage morphology and adhesion-induced Rac and Rho activation.

The Vav family of proteins have the potential to act as both signalling adapters and GEFs for Rho GTPases. They have therefore been proposed as regulators of the cytoskeleton in various cell types. We have used macrophages from mice deficient in all three Vav isoforms to determine how their function affects cell morphology and migration. Macrophages lacking Vav proteins adopt an elongated morphology and have enhanced migratory persistence in culture. To investigate the pathways through which Vav proteins exert their effects we analysed the responses of macrophages to the chemoattractant CSF-1 and to adhesion. We found that morphological and signalling responses of macrophages to CSF-1 did not require Vav proteins. In contrast, adhesion-induced cell spreading, RhoA and Rac1 activation and cell signalling were all dependent on Vav proteins. We propose that Vav proteins affect macrophage morphology and motile behaviour by coupling adhesion receptors to Rac1 and RhoA activity and regulating adhesion signalling events such as paxillin and ERK1/2 phosphorylation by acting as adapters.

MicroRNA-155 is essential for the optimal proliferation and survival of plasmablast B cells.

A fast antibody response can be critical to contain rapidly dividing pathogens. This can be achieved by the expansion of antigen-specific B cells in response to T-cell help followed by differentiation into plasmablasts. MicroRNA-155 (miR-155) is required for optimal T-cell-dependent extrafollicular responses via regulation of PU.1, although the cellular processes underlying this defect are largely unknown. Here, we show that miR-155 regulates the early expansion of B-blasts and later on the survival and proliferation of plasmablasts in a B-cell-intrinsic manner, by tracking antigen-specific B cells in vivo since the onset of antigen stimulation. In agreement, comparative analysis of the transcriptome of miR-155-sufficient and miR-155-deficient plasmablasts at the peak of the response showed that the main processes regulated by miR-155 were DNA metabolic process, DNA replication, and cell cycle. Thus, miR-155 controls the extent of the extrafollicular response by regulating the survival and proliferation of B-blasts, plasmablasts and, consequently, antibody production.

Vav family proteins are required for optimal regulation of PLCgamma2 by integrin alphaIIbbeta3.

Vav proteins belong to the family of guanine-nucleotide-exchange factors for the Rho/Rac family of small G-proteins. In addition, they serve as important adapter proteins for the activation of PLCgamma (phospholipase Cgamma) isoforms by ITAM (immunoreceptor tyrosine-based activation motif) receptors, including the platelet collagen receptor GPVI (glycoprotein VI). Vav proteins are also regulated downstream of integrins, including the major platelet integrin alphaIIbbeta3, which has recently been shown to regulate PLCgamma2. In the present study, we have investigated the role of Vav family proteins in filopodia and lamellipodia formation on fibrinogen using platelets deficient in Vav1 and Vav3. Wild-type mouse platelets undergo a limited degree of spreading on fibrinogen, characterized by the formation of numerous filopodia and limited lamellipodia structures. Platelets deficient in Vav1 and Vav3 exhibit reduced filopodia and lamellipodia formation during spreading on fibrinogen. This is accompanied by reduced alphaIIbbeta3-mediated PLCgamma2 tyrosine phosphorylation and reduced Ca(2+) mobilization. In contrast, the G-protein agonist thrombin stimulates full spreading of control and Vav1/3-deficient platelets. Consistent with this, stimulation of F-actin (filamentous actin) formation and Rac activation by thrombin is not altered in Vav-deficient cells. These results demonstrate that Vav1 and Vav3 are required for optimal spreading and regulation of PLCgamma2 by integrin alphaIIbbeta3, but that their requirement is by-passed upon G-protein receptor activation.