Batch-effect correction in single-cell RNA sequencing data using JIVE.

Motivation: In single-cell RNA sequencing analysis, addressing batch effects-technical artifacts stemming from factors such as varying sequencing technologies, equipment, and capture times-is crucial. These factors can cause unwanted variation and obfuscate the underlying biological signal of interest. The joint and individual variation explained (JIVE) method can be used to extract shared biological patterns from multi-source sequencing data while adjusting for individual non-biological variations (i.e. batch effect). However, its current implementation is originally designed for bulk sequencing data, making it computationally infeasible for large-scale single-cell sequencing datasets. Results: In this study, we enhance JIVE for large-scale single-cell data by boosting its computational efficiency. Additionally, we introduce a novel application of JIVE for batch-effect correction on multiple single-cell sequencing datasets. Our enhanced method aims to decompose single-cell sequencing datasets into a joint structure capturing the true biological variability and individual structures, which capture technical variability within each batch. This joint structure is then suitable for use in downstream analyses. We benchmarked the results against four popular tools, Seurat v5, Harmony, LIGER, and Combat-seq, which were developed for this purpose. JIVE performed best in terms of preserving cell-type effects and in scenarios in which the batch sizes are balanced. Availability and implementation: The JIVE implementation used for this analysis can be found at https://github.com/oconnell-statistics-lab/scJIVE.

The Orbitofrontal Cortex Is Required for Learned Modulation of Innate Olfactory Behavior.

Animals have evolved innate responses to cues including social, food, and predator odors. In the natural environment, animals are faced with choices that involve balancing risk and reward where innate significance may be at odds with internal need. The ability to update the value of a cue through learning is essential for navigating changing and uncertain environments. However, the mechanisms involved in this modulation are not well defined in mammals. We have established a new olfactory assay that challenges a thirsty mouse to choose an aversive odor over an attractive odor in foraging for water, thus overriding their innate behavioral response to odor. Innately, mice prefer the attractive odor port over the aversive odor port. However, decreasing the probability of water at the attractive port leads mice to prefer the aversive port, reflecting a learned override of the innate response to the odors. The orbitofrontal cortex (OFC) is a fourth-order olfactory brain area, involved in flexible value association, with behaviorally relevant outputs throughout the limbic system. We performed optogenetic and chemogenetic silencing experiments that demonstrate the OFC is necessary for this learned modulation of innate aversion to odor. Further, we characterized odor evoked c-fos expression in learned and control mice and found significant suppression of activity in the bed nucleus of the stria terminalis, lateral septum, and central and medial amygdala. These findings reveal that the OFC is necessary for the learned override of innate behavior and may signal to limbic structures to modulate innate response to odor.

Transformation of valence signaling in a mouse striatopallidal circuit.

The ways in which sensory stimuli acquire motivational valence through association with other stimuli is one of the simplest forms of learning. Although we have identified many brain nuclei that play various roles in reward processing, a significant gap remains in understanding how valence encoding transforms through the layers of sensory processing. To address this gap, we carried out a comparative investigation of the mouse anteromedial olfactory tubercle (OT), and the ventral pallidum (VP) - 2 connected nuclei of the basal ganglia which have both been implicated in reward processing. First, using anterograde and retrograde tracing, we show that both D1 and D2 neurons of the anteromedial OT project primarily to the VP and minimally elsewhere. Using two-photon calcium imaging, we then investigated how the identity of the odor and reward contingency of the odor are differently encoded by neurons in either structure during a classical conditioning paradigm. We find that VP neurons robustly encode reward contingency, but not identity, in low-dimensional space. In contrast, the OT neurons primarily encode odor identity in high-dimensional space. Although D1 OT neurons showed larger responses to rewarded odors than other odors, consistent with prior findings, we interpret this as identity encoding with enhanced contrast. Finally, using a novel conditioning paradigm that decouples reward contingency and licking vigor, we show that both features are encoded by non-overlapping VP neurons. These results provide a novel framework for the striatopallidal circuit in which a high-dimensional encoding of stimulus identity is collapsed onto a low-dimensional encoding of motivational valence.

Genome-wide association study of major anxiety disorders in 122,341 European-ancestry cases identifies 58 loci and highlights GABAergic signaling.

The major anxiety disorders (ANX; including generalized anxiety disorder, panic disorder, and phobias) are highly prevalent, often onset early, persist throughout life, and cause substantial global disability. Although distinct in their clinical presentations, they likely represent differential expressions of a dysregulated threat-response system. Here we present a genome-wide association meta-analysis comprising 122,341 European ancestry ANX cases and 729,881 controls. We identified 58 independent genome-wide significant ANX risk variants and 66 genes with robust biological support. In an independent sample of 1,175,012 self-report ANX cases and 1,956,379 controls, 51 of the 58 associated variants were replicated. As predicted by twin studies, we found substantial genetic correlation between ANX and depression, neuroticism, and other internalizing phenotypes. Follow-up analyses demonstrated enrichment in all major brain regions and highlighted GABAergic signaling as one potential mechanism underlying ANX genetic risk. These results advance our understanding of the genetic architecture of ANX and prioritize genes for functional follow-up studies.

Control of innate olfactory valence by segregated cortical amygdala circuits.

Animals perform innate behaviors that are stereotyped responses to specific evolutionarily relevant stimuli in the absence of prior learning or experience. These behaviors can be reduced to an axis of valence, whereby specific odors evoke approach or avoidance. The cortical amygdala (plCoA) mediates innate attraction and aversion to odor. However, little is known about how this brain area gives rise to behaviors of opposing motivational valence. Here, we sought to define the circuit features of plCoA that give rise to innate olfactory behaviors of valence. We characterized the physiology, gene expression, and projections of this structure, identifying a divergent, topographic organization that selectively controls innate attraction and avoidance to odor. First, we examined odor-evoked responses in these areas and found sparse encoding of odor identity, but not valence. We next considered a topographic organization and found that optogenetic stimulation of the anterior and posterior domains of plCoA elicits attraction and avoidance, respectively, suggesting a functional axis for valence. Using single cell and spatial RNA sequencing, we identified the molecular cell types in plCoA, revealing an anteroposterior gradient in cell types, whereby anterior glutamatergic neurons preferentially express Slc17a6 and posterior neurons express Slc17a7. Activation of these respective cell types recapitulates appetitive and aversive valence behaviors, and chemogenetic inhibition reveals partial necessity for valence responses to innate appetitive or aversive odors. Finally, we identified topographically organized circuits defined by projections, whereby anterior neurons preferentially project to medial amygdala, and posterior neurons preferentially project to nucleus accumbens, which are respectively sufficient and necessary for innate negative and positive olfactory valence. Together, these data advance our understanding of how the olfactory system generates stereotypic, hardwired attraction and avoidance, and supports a model whereby distinct, topographically distributed plCoA populations direct innate olfactory valence responses by signaling to divergent valence-specific targets, linking upstream olfactory identity to downstream valence behaviors, through a population code. This represents a novel circuit motif in which valence encoding is represented not by the firing properties of individual neurons, but by population level identity encoding that is routed through divergent targets to mediate distinct valence.

GAUSS: a summary-statistics-based R package for accurate estimation of linkage disequilibrium for variants, Gaussian imputation, and TWAS analysis of cosmopolitan cohorts.

MOTIVATION: As the availability of larger and more ethnically diverse reference panels grows, there is an increase in demand for ancestry-informed imputation of genome-wide association studies (GWAS), and other downstream analyses, e.g. fine-mapping. Performing such analyses at the genotype level is computationally challenging and necessitates, at best, a laborious process to access individual-level genotype and phenotype data. Summary-statistics-based tools, not requiring individual-level data, provide an efficient alternative that streamlines computational requirements and promotes open science by simplifying the re-analysis and downstream analysis of existing GWAS summary data. However, existing tools perform only disparate parts of needed analysis, have only command-line interfaces, and are difficult to extend/link by applied researchers. RESULTS: To address these challenges, we present Genome Analysis Using Summary Statistics (GAUSS)-a comprehensive and user-friendly R package designed to facilitate the re-analysis/downstream analysis of GWAS summary statistics. GAUSS offers an integrated toolkit for a range of functionalities, including (i) estimating ancestry proportion of study cohorts, (ii) calculating ancestry-informed linkage disequilibrium, (iii) imputing summary statistics of unobserved variants, (iv) conducting transcriptome-wide association studies, and (v) correcting for "Winner's Curse" biases. Notably, GAUSS utilizes an expansive, multi-ethnic reference panel consisting of 32 953 genomes from 29 ethnic groups. This panel enhances the range and accuracy of imputable variants, including the ability to impute summary statistics of rarer variants. As a result, GAUSS elevates the quality and applicability of existing GWAS analyses without requiring access to subject-level genotypic and phenotypic information. AVAILABILITY AND IMPLEMENTATION: The GAUSS R package, complete with its source code, is readily accessible to the public via our GitHub repository at https://github.com/statsleelab/gauss. To further assist users, we provided illustrative use-case scenarios that are conveniently found at https://statsleelab.github.io/gauss/, along with a comprehensive user guide detailed in Supplementary Text S1.

Transformation of valence signaling in a striatopallidal circuit.

The ways in which sensory stimuli acquire motivational valence through association with other stimuli is one of the simplest forms of learning. Though we have identified many brain nuclei that play various roles in reward processing, a significant gap remains in understanding how valence encoding transforms through the layers of sensory processing. To address this gap, we carried out a comparative investigation of the olfactory tubercle (OT), and the ventral pallidum (VP) - 2 connected nuclei of the basal ganglia which have both been implicated in reward processing. First, using anterograde and retrograde tracing, we show that both D1 and D2 neurons of the OT project primarily to the VP and minimally elsewhere. Using 2-photon calcium imaging, we then investigated how the identity of the odor and reward contingency of the odor are differently encoded by neurons in either structure during a classical conditioning paradigm. We find that VP neurons robustly encode reward contingency, but not identity, in low-dimensional space. In contrast, OT neurons primarily encode odor identity in high-dimensional space. Though D1 OT neurons showed larger response vectors to rewarded odors than other odors, we propose this is better interpreted as identity encoding with enhanced contrast rather than as valence encoding. Finally, using a novel conditioning paradigm that decouples reward contingency and licking vigor, we show that both features are encoded by non-overlapping VP neurons. These results provide a novel framework for the striatopallidal circuit in which a high-dimensional encoding of stimulus identity is collapsed onto a low-dimensional encoding of motivational valence.

KOMPUTE: imputing summary statistics of missing phenotypes in high-throughput model organism data.

Motivation: The International Mouse Phenotyping Consortium (IMPC) is striving to build a comprehensive functional catalog of mammalian protein-coding genes by systematically producing and phenotyping gene-knockout mice for almost every protein-coding gene in the mouse genome and by testing associations between gene loss-of-function and phenotype. To date, the IMPC has identified over 90 000 gene-phenotype associations, but many phenotypes have not yet been measured for each gene, resulting in largely incomplete data; ∼75.6% of association summary statistics are still missing in the latest IMPC summary statistics dataset (IMPC release version 16). Results: To overcome these challenges, we propose KOMPUTE, a novel method for imputing missing summary statistics in the IMPC dataset. Using conditional distribution properties of multivariate normal, KOMPUTE estimates the association Z-scores of unmeasured phenotypes for a particular gene as a conditional expectation given the Z-scores of measured phenotypes. Our evaluation of the method using simulated and real-world datasets demonstrates its superiority over the singular value decomposition matrix completion method in various scenarios. Availability and implementation: An R package for KOMPUTE is publicly available at https://github.com/statsleelab/kompute, along with usage examples and results for different phenotype domains at https://statsleelab.github.io/komputeExamples.

Hepatocyte-specific deletion of Bis causes senescence in the liver without deteriorating hepatic function.

Bcl-2-interacting cell death suppressor (BIS), also called as BAG3, regulates numerous physiological processes, such as apoptosis, protein quality control, and senescence. Whole-body Bis-knockout (KO) mice exhibit early lethality following cardiac and skeletal muscle dysfunction. The first attempt to generate organ-specific knockout mice resulted in constitutive or inducible heart-specific Bis-knockout mice, which exhibited cardiac dilation and underwent premature death. Here, we generated hepatocyte-specific Bis-knockout (Bis-HKO) mice and found no abnormalities in metabolic function and survival. However, depletion of HSPB8 and accumulation of p62 indicated impaired autophagy in Bis-HKO livers. Interestingly, the number of peroxisomes wrapped by phagophore membranes increased as evidenced by transmission electron microscopy analysis, indicating defects in the progression of pexophagy. In addition, increased dihydroethidine intensities and histone H3 K9me3-positive nuclei indicated increased oxidative stress and senescence induction in Bis-HKO livers. Mechanistically, p27 was upregulated in Bis-HKO livers. In SNU368 hepatocellular carcinoma cells, BIS depletion led to p27 upregulation, and increase in histone H3 K9me3 levels and senescence-associated ß-galactosidase staining; therefore, reproducing the in vivo senescence phenotype. Despite the observation of no metabolic abnormalities, BIS depletion led to defective autophagy, increased oxidative stress, and senescence in Bis-HKO livers. Collectively, our results suggest a role for BIS in maintaining liver regeneration potential under pathological conditions.

Inhibitory Effect of Etravirine, a Non-Nucleoside Reverse Transcriptase Inhibitor, via Anterior Gradient Protein 2 Homolog Degradation against Ovarian Cancer Metastasis.

Anterior gradient protein 2 homolog (AGR2), an endoplasmic reticulum protein, is secreted in the tumor microenvironment. AGR2 is a member of the disulfide isomerase family, is highly expressed in multiple cancers, and promotes cancer metastasis. In this study, we found that etravirine, which is a non-nucleoside reverse transcriptase inhibitor, could induce AGR2 degradation via autophagy. Moreover, etravirine diminished proliferation, migration, and invasion in vitro. Moreover, in an orthotopic xenograft mouse model, the combination of etravirine and paclitaxel significantly suppressed cancer progression and metastasis. This drug may be a promising therapeutic agent for the treatment of ovarian cancer.

Protective Effects and Benefits of Olive Oil and Its Extracts on Women's Health.

Women and men share similar diseases; however, women have unique issues, including gynecologic diseases and diseases related to menstruation, menopause, and post menopause. In recent decades, scientists paid more attention to natural products and their derivatives because of their good tolerability and effectiveness in disease prevention and treatment. Olive oil is an essential component in the Mediterranean diet, a diet well known for its protective impact on human well-being. Investigation of the active components in olive oil, such as oleuropein and hydroxytyrosol, showed positive effects in various diseases. Their effects have been clarified in many suggested mechanisms and have shown promising results in animal and human studies, especially in breast cancer, ovarian cancer, postmenopausal osteoporosis, and other disorders. This review summarizes the current evidence of the role of olives and olive polyphenols in women's health issues and their potential implications in the treatment and prevention of health problems in women.

Increasing the resolution and precision of psychiatric genome-wide association studies by re-imputing summary statistics using a large, diverse reference panel.

Genotype imputation across populations of mixed ancestry is critical for optimal discovery in large-scale genome-wide association studies (GWAS). Methods for direct imputation of GWAS summary-statistics were previously shown to be practically as accurate as summary statistics produced after raw genotype imputation, while incurring orders of magnitude lower computational burden. Given that direct imputation needs a precise estimation of linkage-disequilibrium (LD) and that most of the methods using a small reference panel for example, ~2,500-subject coming from the 1000 Genome-Project, there is a great need for much larger and more diverse reference panels. To accurately estimate the LD needed for an exhaustive analysis of any cosmopolitan cohort, we developed DISTMIX2. DISTMIX2: (a) uses a much larger and more diverse reference panel compared to traditional reference panels, and (b) can estimate weights of ethnic-mixture based solely on Z-scores, when allele frequencies are not available. We applied DISTMIX2 to GWAS summary-statistics from the psychiatric genetic consortium (PGC). DISTMIX2 uncovered signals in numerous new regions, with most of these findings coming from the rarer variants. Rarer variants provide much sharper location for the signals compared with common variants, as the LD for rare variants extends over a lower distance than for common ones. For example, while the original PGC post-traumatic stress disorder GWAS found only 3 marginal signals for common variants, we now uncover a very strong signal for a rare variant in PKN2, a gene associated with neuronal and hippocampal development. Thus, DISTMIX2 provides a robust and fast (re)imputation approach for most psychiatric GWAS-studies.

TWAS pathway method greatly enhances the number of leads for uncovering the molecular underpinnings of psychiatric disorders.

Genetic signal detection in genome-wide association studies (GWAS) is enhanced by pooling small signals from multiple Single Nucleotide Polymorphism (SNP), for example, across genes and pathways. Because genes are believed to influence traits via gene expression, it is of interest to combine information from expression Quantitative Trait Loci (eQTLs) in a gene or genes in the same pathway. Such methods, widely referred to as transcriptomic wide association studies (TWAS), already exist for gene analysis. Due to the possibility of eliminating most of the confounding effects of linkage disequilibrium (LD) from TWAS gene statistics, pathway TWAS methods would be very useful in uncovering the true molecular basis of psychiatric disorders. However, such methods are not yet available for arbitrarily large pathways/gene sets. This is possibly due to the quadratic (as a function of the number of SNPs) computational burden for computing LD across large chromosomal regions. To overcome this obstacle, we propose JEPEGMIX2-P, a novel TWAS pathway method that (a) has a linear computational burden, (b) uses a large and diverse reference panel (33 K subjects), (c) is competitive (adjusts for background enrichment in gene TWAS statistics), and (d) is applicable as-is to ethnically mixed-cohorts. To underline its potential for increasing the power to uncover genetic signals over the commonly used nontranscriptomics methods, for example, MAGMA, we applied JEPEGMIX2-P to summary statistics of most large meta-analyses from Psychiatric Genetics Consortium (PGC). While our work is just the very first step toward clinical translation of psychiatric disorders, PGC anorexia results suggest a possible avenue for treatment.

V-SVA: an R Shiny application for detecting and annotating hidden sources of variation in single-cell RNA-seq data.

SUMMARY: Single-cell RNA-sequencing (scRNA-seq) technology enables studying gene expression programs from individual cells. However, these data are subject to diverse sources of variation, including 'unwanted' variation that needs to be removed in downstream analyses (e.g. batch effects) and 'wanted' or biological sources of variation (e.g. variation associated with a cell type) that needs to be precisely described. Surrogate variable analysis (SVA)-based algorithms, are commonly used for batch correction and more recently for studying 'wanted' variation in scRNA-seq data. However, interpreting whether these variables are biologically meaningful or stemming from technical reasons remains a challenge. To facilitate the interpretation of surrogate variables detected by algorithms including IA-SVA, SVA or ZINB-WaVE, we developed an R Shiny application [Visual Surrogate Variable Analysis (V-SVA)] that provides a web-browser interface for the identification and annotation of hidden sources of variation in scRNA-seq data. This interactive framework includes tools for discovery of genes associated with detected sources of variation, gene annotation using publicly available databases and gene sets, and data visualization using dimension reduction methods. AVAILABILITY AND IMPLEMENTATION: The V-SVA Shiny application is publicly hosted at https://vsva.jax.org/ and the source code is freely available at https://github.com/nlawlor/V-SVA. CONTACT: leed13@miamioh.edu or duygu.ucar@jax.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Detection of correlated hidden factors from single cell transcriptomes using Iteratively Adjusted-SVA (IA-SVA).

Single cell RNA-sequencing (scRNA-seq) precisely characterizes gene expression levels and dissects variation in expression associated with the state (technical or biological) and the type of the cell, which is averaged out in bulk measurements. Multiple and correlated sources contribute to gene expression variation in single cells, which makes their estimation difficult with the existing methods developed for batch correction (e.g., surrogate variable analysis (SVA)) that estimate orthogonal transformations of these sources. We developed iteratively adjusted surrogate variable analysis (IA-SVA) that can estimate hidden factors even when they are correlated with other sources of variation by identifying a set of genes associated with each hidden factor in an iterative manner. Analysis of scRNA-seq data from human cells showed that IA-SVA could accurately capture hidden variation arising from technical (e.g., stacked doublet cells) or biological sources (e.g., cell type or cell-cycle stage). Furthermore, IA-SVA delivers a set of genes associated with the detected hidden source to be used in downstream data analyses. As a proof of concept, IA-SVA recapitulated known marker genes for islet cell subsets (e.g., alpha, beta), which improved the grouping of subsets into distinct clusters. Taken together, IA-SVA is an effective and novel method to dissect multiple and correlated sources of variation in scRNA-seq data.

Amniotic fluid index, single deepest pocket and transvaginal cervical length: Parameter of predictive delivery latency in preterm premature rupture of membranes.

OBJECTIVE: Prediction of delivery latency complicated with preterm premature rupture of membrane (PPROM) is crucial for reducing maternal and neonatal complications. Therefore, we investigated the correlations between latency period and cut-off values of ultrasonographic parameters, ultimately predicting delivery latency. MATERIALS AND METHODS: The retrospective study was performed on 121 PPROM patients enrolled between March 2010 and July 2015. Parameters including amniotic fluid index (AFI), single deepest pocket (SDP) and transvaginal cervical length (TVCL) were measured in 99 singleton pregnancies with PPROM. Latency was defined as the period from sonographic measurements to delivery day. The parameters were analyzed independently by Wilcoxon rank sum test and Fisher's exact test. Cut-off values were determined using a receiver operating characteristic (ROC) curve. RESULTS: In delivery latency within 3 days, AFI and SDP were decreased with significantly shorter TVCL. AFI and SDP had the highest sensitivity (82.2%) and SDP combined with TVCL showed the highest specificity (75.9%) in area under curve (AUC) value. The predicted median latency period was less than 2 days within the cutoff value of parameter (AFI ≤ 7.72, SDP ≤ 3.2 and TVCL ≤ 1.69). CONCLUSION: AFI and SDP combined with TVCL could be useful predictive parameters of the latency interval from PPROM to delivery.

Leptin, leptin receptors and hypoxia-induced factor-1α expression in the placental bed of patients with and without preeclampsia during pregnancy.

The mechanism underlying the pathogenesis of preeclampsia (PE) has been previously investigated but remains to be elucidated. Among numerous biomarkers that are associated with the pathogenesis of PE, leptin is most frequently investigated. Although studies concerning the association between PE and the expression of leptin in the serum and placenta have been conducted, the results are conflicting and inconsistent. Furthermore, the expression of leptin and its receptors in the placental bed and their association with PE, to the best of our knowledge, has not been previously reported. Therefore, to determine the association between the expression of leptin and its receptor, and pathogenesis and onset period of PE, placental bed tissues were obtained from cesarean section deliveries. The mRNA and protein expression levels of leptin and its receptor were investigated in normal pregnancies (n=18), pregnancies complicated with early­onset PE (n=9) and late­onset PE (n=9) by reverse transcription­quantitative polymerase chain reaction and western blotting, respectively. The results demonstrated that the mRNA and protein expression of leptin in the placental bed was significantly increased in the PE groups compared with normal controls and was associated with the onset period of PE. Furthermore, as evidenced by immunostaining, leptin was upregulated in endothelial cells of the placental bed in the PE groups, with a particularly strong upregulation in activated endothelial cells from patients with early­onset PE. The results of the present study indicate that altered expression of leptin in the placental bed may contribute to the pathogenesis of PE.

The Effect of Electrical Muscle Stimulation and In-bed Cycling on Muscle Strength and Mass of Mechanically Ventilated Patients: A Pilot Study.

BACKGROUND: Critically ill patients experience muscle weakness, which leads to functional disability. Both functional electrical stimulation (FES) and in-bed cycling can be an alternative measure for intensive care unit (ICU) patients who are not feasible for active exercise. The aim of this study was to examine whether FES and in-bed cycling have a positive effect on muscle mass in ICU patients. METHODS: Critically ill patients who received mechanical ventilation for at least 24 hours were included. After passive range of motion exercise, in-bed cycling was applied for 20 minutes, and FES was applied for 20 minutes on the left leg. The right leg received in-bed cycling and the left leg received both FES and in-bed cycling. Thigh circumferences and rectus femoris cross-sectional area (CSA) were assessed with ultrasonography before and after the intervention. Muscle strength was assessed by Medical Research Council scale. RESULTS: A total of 10 patients were enrolled in this study as a pilot study. Before and after the intervention, the CSA of right rectus femoris increased from 5.08 ± 1.51 cm2 to 6.01 ± 2.21 cm2 , which was statistically significant (P = 0.003). The thigh circumference was also increased and statistically significant (P = 0.006). There was no difference between left and right in regard to FES application. There is no significant change in muscle strength before and after the intervention (right and left, P = 0.317 and P = 0.368, respectively). CONCLUSIONS: In-bed cycling increased thigh circumferences rectus femoris CSA. Adding FES did not show differences.

JEPEGMIX2: improved gene-level joint analysis of eQTLs in cosmopolitan cohorts.

Motivation: To increase detection power, researchers use gene level analysis methods to aggregate weak marker signals. Due to gene expression controlling biological processes, researchers proposed aggregating signals for expression Quantitative Trait Loci (eQTL). Most gene-level eQTL methods make statistical inferences based on (i) summary statistics from genome-wide association studies (GWAS) and (ii) linkage disequilibrium patterns from a relevant reference panel. While most such tools assume homogeneous cohorts, our Gene-level Joint Analysis of functional SNPs in Cosmopolitan Cohorts (JEPEGMIX) method accommodates cosmopolitan cohorts by using heterogeneous panels. However, JEPGMIX relies on brain eQTLs from older gene expression studies and does not adjust for background enrichment in GWAS signals. Results: We propose JEPEGMIX2, an extension of JEPEGMIX. When compared to JPEGMIX, it uses (i) cis-eQTL SNPs from the latest expression studies and (ii) brains specific (sub)tissues and tissues other than brain. JEPEGMIX2 also (i) avoids accumulating averagely enriched polygenic information by adjusting for background enrichment and (ii) to avoid an increase in false positive rates for studies with numerous highly enriched (above the background) genes, it outputs gene q-values based on Holm adjustment of P-values. Availability and implementation: https://github.com/Chatzinakos/JEPEGMIX2. Supplementary information: Supplementary data are available at Bioinformatics online.