Integrated analysis of single-cell RNA-seq, bulk RNA-seq, Mendelian randomization, and eQTL reveals T cell-related nomogram model and subtype classification in rheumatoid arthritis.

Objective: Rheumatoid arthritis (RA) is a systemic disease that attacks the joints and causes a heavy economic burden on humans worldwide. T cells regulate RA progression and are considered crucial targets for therapy. Therefore, we aimed to integrate multiple datasets to explore the mechanisms of RA. Moreover, we established a T cell-related diagnostic model to provide a new method for RA immunotherapy. Methods: scRNA-seq and bulk-seq datasets for RA were obtained from the Gene Expression Omnibus (GEO) database. Various methods were used to analyze and characterize the T cell heterogeneity of RA. Using Mendelian randomization (MR) and expression quantitative trait loci (eQTL), we screened for potential pathogenic T cell marker genes in RA. Subsequently, we selected an optimal machine learning approach by comparing the nine types of machine learning in predicting RA to identify T cell-related diagnostic features to construct a nomogram model. Patients with RA were divided into different T cell-related clusters using the consensus clustering method. Finally, we performed immune cell infiltration and clinical correlation analyses of T cell-related diagnostic features. Results: By analyzing the scRNA-seq dataset, we obtained 10,211 cells that were annotated into 7 different subtypes based on specific marker genes. By integrating the eQTL from blood and RA GWAS, combined with XGB machine learning, we identified a total of 8 T cell-related diagnostic features (MIER1, PPP1CB, ICOS, GADD45A, CD3D, SLFN5, PIP4K2A, and IL6ST). Consensus clustering analysis showed that RA could be classified into two different T-cell patterns (Cluster 1 and Cluster 2), with Cluster 2 having a higher T-cell score than Cluster 1. The two clusters involved different pathways and had different immune cell infiltration states. There was no difference in age or sex between the two different T cell patterns. In addition, ICOS and IL6ST were negatively correlated with age in RA patients. Conclusion: Our findings elucidate the heterogeneity of T cells in RA and the communication role of these cells in an RA immune microenvironment. The construction of T cell-related diagnostic models provides a resource for guiding RA immunotherapeutic strategies.

Association of cigarette smoking, smoking cessation with the risk of cardiometabolic multimorbidity in the UK Biobank.

BACKGROUND: To investigate the association between cigarette smoking, smoking cessation and the trajectory of cardiometabolic multimorbidity (CMM), and further to examine the association of age at smoking initiation and smoking cessation with CMM. METHODS: This study included 298,984 UK Biobank participants without cardiometabolic diseases (CMDs) (including type 2 diabetes, coronary heart diseases, stroke, and hypertension) at baseline. Smoking status was categorized into former, current, and never smokers, with age at smoking initiation and smoking cessation as a proxy for current and former smokers. The multi-state model was performed to evaluate the association between cigarette smoking, smoking cessation and CMM. RESULTS: During a median follow-up of 13.2 years, 59,193 participants developed first cardiometabolic disease (FCMD), 14,090 further developed CMM, and 16,487 died. Compared to former smokers, current smokers had higher risk at all transitions, with hazard ratio (95% confidence interval) = 1.59 (1.55 ∼ 1.63) vs. 1.18 (1.16 ∼ 1.21) (P = 1.48 × 10- 118) from health to FCMD, 1.40 (1.33 ∼ 1.47) vs. 1.09 (1.05 ∼ 1.14) (P = 1.50 × 10- 18) from FCMD to CMM, and 2.87 (2.72 ∼ 3.03) vs. 1.38 (1.32 ∼ 1.45) (P < 0.001) from health, 2.16 (1.98 ∼ 2.35) vs. 1.25 (1.16 ∼ 1.34) (P = 1.18 × 10- 46) from FCMD, 2.02 (1.79 ∼ 2.28) vs. 1.22 (1.09 ∼ 1.35) (P = 3.93 × 10- 17) from CMM to death; whereas quitting smoking reduced the risk attributed to cigarette smoking by approximately 76.5% across all transitions. Reduced risks of smoking cessation were also identified when age at quitting smoking was used as a proxy for former smokers. CONCLUSIONS: Cigarette smoking was associated with a higher risk of CMM across all transitions; however, smoking cessation, especially before the age of 35, was associated with a significant decrease in CMM risk attributed to cigarette smoking.

Invasive metastatic tumor-camouflaged ROS responsive nanosystem for targeting therapeutic brain injury after cardiac arrest.

Drug transmission through the blood-brain barrier (BBB) is considered an arduous challenge for brain injury treatment following the return of spontaneous circulation after cardiac arrest (CA-ROSC). Inspired by the propensity of melanoma metastasis to the brain, B16F10 cell membranes are camouflaged on 2-methoxyestradiol (2ME2)-loaded reactive oxygen species (ROS)-triggered "Padlock" nanoparticles that are constructed by phenylboronic acid pinacol esters conjugated D-a-tocopheryl polyethylene glycol succinate (TPGS-PBAP). The biomimetic nanoparticles (BM@TP/2ME2) can be internalized, mainly mediated by the mutual recognition and interaction between CD44v6 expressed on B16F10 cell membranes and hyaluronic acid on cerebral vascular endothelial cells, and they responsively release 2ME2 by the oxidative stress microenvironment. Notably, BM@TP/2ME2 can scavenge excessive ROS to reestablish redox balance, reverse neuroinflammation, and restore autophagic flux in damaged neurons, eventually exerting a remarkable neuroprotective effect after CA-ROSC in vitro and in vivo. This biomimetic drug delivery system is a novel and promising strategy for the treatment of cerebral ischemia-reperfusion injury after CA-ROSC.

Effect of probiotics on children with autism spectrum disorders: a meta-analysis.

BACKGROUND: Researches have found that alteration of intestinal flora may be closely related to the development of autism spectrum disorder (ASD). However, whether probiotics supplementation has a protective effect on ASD remains controversial. This meta-analysis aimed to analyze the outcome of probiotics in the treatment of ASD children. METHODS: The Pubmed, Cochrane Library, Web of Science and Embase were searched until Sep 2022. Randomized controlled trials (RCTs) relevant to the probiotics and placebo treatment on ASD children were screened. Quality assessment of the included RCTs was evaluated by the Cochrane collaboration's tool. The primary outcomes were ASD assessment scales, including ABC (aberrant behavior checklist) and CBCL (child behavior checklist) for evaluating the behavior improvement, SRS (social responsiveness scale) for social assessment, DQ (developmental quotient) for physical and mental development and CGI-I (clinical global impression improvement) for overall improvement. The secondary outcome was total 6-GSI (gastrointestinal severity index). RESULTS: In total, 6 RCTs from 6 studies with 302 children were included in the systemic review. Total 6-GSI (MD=-0.59, 95%CI [-1.02,-0.17], P < 0.05) decreased significantly after oral administration of probiotics. Whereas, there was no statistical difference in ABC, CBCL, SRS, DQ and CGI-I between probiotics and placebo groups in ASD children. CONCLUSION: Probiotics treatment could improve gastrointestinal symptoms, but there was no significant improvement in ASD.

A comprehensive comparison of microbial communities between aerobic granular sludge and flocculent sludge for nutrient removal in full-scale wastewater treatment plants.

Understanding the microbial community structure of sludge is crucial for improving the design, operation and optimisation of full-scale wastewater treatment plants (WWTPs). This study aimed to have a comprehensive comparison of microbial communities between aerobic granular sludge and flocculent sludge from two full-scale sequential batch reactors-based WWTPs with nutrient removal for the first time. To better understand key functional bacteria such as polyphosphate accumulating bacteria (PAOs), competitive bacteria such as glycogen accumulating bacteria (GAOs) and nitrifying bacteria for both nitrogen and phosphorus removal, another two full-scale WWTPs with only carbon (C) removal and C and nitrogen (N) removal were compared too. It was found that the richness and diversity of the microbial population in sludge increased with pollutant removal from only C, C and N, to C,N, P removal. For C, N P removal, granule structure led to a more diverse and rich microbial community structure than flocculent structure. Although more abundant nitrifying bacteria were enriched in granular sludge than flocculent sludge, the abundance of total putative PAOs was equivalent. However, the most typical putative PAOs such as Tetrasphaera and Candidatus Accumulibacter seemed to be more correlated with biological phosphorus removal performance, which might be more proper to be used as an indication for P removal potential. The higher abundance of GAOs in flocculent sludge with better phosphorus removal performance might suggest that further investigation is needed to understand the functions of GAOs. In addition, the equivalent abundances of PAOs in the WWTPs with only C removal and with C, N, and P removal, respectively, indicate that many newly reported putative PAOs might not contribute to P removal. This study provides insight into the microbial communities and functional bacteria in aerobic granular sludge and flocculent sludge in full-scale SBRs, which can provide microbes-informed optimisation of reactor operation for better nutrient removal.

Genetic variants in UNC93B1 predispose to childhood-onset systemic lupus erythematosus.

Rare genetic variants in toll-like receptor 7 (TLR7) are known to cause lupus in humans and mice. UNC93B1 is a transmembrane protein that regulates TLR7 localization into endosomes. In the present study, we identify two new variants in UNC93B1 (T314A, located proximally to the TLR7 transmembrane domain, and V117L) in a cohort of east Asian patients with childhood-onset systemic lupus erythematosus. The V117L variant was associated with increased expression of type I interferons and NF-κB-dependent cytokines in patient plasma and immortalized B cells. THP-1 cells expressing the variant UNC93B1 alleles exhibited exaggerated responses to stimulation of TLR7/-8, but not TLR3 or TLR9, which could be inhibited by targeting the downstream signaling molecules, IRAK1/-4. Heterozygous mice expressing the orthologous Unc93b1V117L variant developed a spontaneous lupus-like disease that was more severe in homozygotes and again hyperresponsive to TLR7 stimulation. Together, this work formally identifies genetic variants in UNC93B1 that can predispose to childhood-onset systemic lupus erythematosus.

Identification of potential shared gene signatures between gastric cancer and type 2 diabetes: a data-driven analysis.

Background: Gastric cancer (GC) and type 2 diabetes (T2D) contribute to each other, but the interaction mechanisms remain undiscovered. The goal of this research was to explore shared genes as well as crosstalk mechanisms between GC and T2D. Methods: The Gene Expression Omnibus (GEO) database served as the source of the GC and T2D datasets. The differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA) were utilized to identify representative genes. In addition, overlapping genes between the representative genes of the two diseases were used for functional enrichment analysis and protein-protein interaction (PPI) network. Next, hub genes were filtered through two machine learning algorithms. Finally, external validation was undertaken with data from the Cancer Genome Atlas (TCGA) database. Results: A total of 292 and 541 DEGs were obtained from the GC (GSE29272) and T2D (GSE164416) datasets, respectively. In addition, 2,704 and 336 module genes were identified in GC and T2D. Following their intersection, 104 crosstalk genes were identified. Enrichment analysis indicated that "ECM-receptor interaction," "AGE-RAGE signaling pathway in diabetic complications," "aging," and "cellular response to copper ion" were mutual pathways. Through the PPI network, 10 genes were identified as candidate hub genes. Machine learning further selected BGN, VCAN, FN1, FBLN1, COL4A5, COL1A1, and COL6A3 as hub genes. Conclusion: "ECM-receptor interaction," "AGE-RAGE signaling pathway in diabetic complications," "aging," and "cellular response to copper ion" were revealed as possible crosstalk mechanisms. BGN, VCAN, FN1, FBLN1, COL4A5, COL1A1, and COL6A3 were identified as shared genes and potential therapeutic targets for people suffering from GC and T2D.

Investigating the Relationship between Birthweight and Breast Cancer from A Non-Linear and Mediation Perspective.

Background: Epidemiological studies have shown a positive relationship between birthweight and breast cancer; however, inconsistent, sometimes even controversial, observations emerged. We re-explored the association between them in the UK Biobank cohort. Methods: Relying on the UK Biobank cohort data of white British volunteers recruited between 2006 and 2010 (5,760 cases and 162,778 controls), we evaluated the causal mediation between birthweight and breast cancer, with age of menarche and age at menopause as two potential mediators under the traditional mediation analysis framework. The non-linear relationship between birthweight and breast cancer was also investigated by including the square of birthweight or discretized birthweight categories (<2.5, 2.5~4.0, or >4.0). Furthermore, we performed a stratification analysis in terms of the menopause status. Results: Birthweight can indirectly influence breast cancer risk in adulthood via the path of age of menarche or age at menopause, and found statistical evidence supporting the existence of suggestive non-linear association between birthweight and breast cancer (ß=0.062 and P=0.004 for the square of birthweight) although failing to discover a linear relationship (P=0.230). We also demonstrated such non-linear association seemed more pronounced and robust for premenopausal women compared with postmenopausal ones (27.5% vs. 19.5% increase in breast cancer risk). Conclusion: This study provided an in-depth insight into the observed relationship between birthweight and breast cancer and revealed that non-linear impact and causal mediation commonly drive the connection between the two traits.

Bioinformatics analysis of PANoptosis regulators in the diagnosis and subtyping of steroid-induced osteonecrosis of the femoral head.

In this study, we aimed to investigate the involvement of PANoptosis, a form of regulated cell death, in the development of steroid-induced osteonecrosis of the femoral head (SONFH). The underlying pathogenesis of PANoptosis in SONFH remains unclear. To address this, we employed bioinformatics approaches to analyze the key genes associated with PANoptosis. Our analysis was based on the GSE123568 dataset, allowing us to investigate both the expression profiles of PANoptosis-related genes (PRGs) and the immune profiles in SONFHallowing us to investigate the expression profiles of PRGs as well as the immune profiles in SONFH. We conducted cluster classification based on PRGs and assessed immune cell infiltration. Additionally, we used the weighted gene co-expression network analysis (WGCNA) algorithm to identify cluster-specific hub genes. Furthermore, we developed an optimal machine learning model to identify the key predictive genes responsible for SONFH progression. We also constructed a nomogram model with high predictive accuracy for assessing risk factors in SONFH patients, and validated the model using external data (area under the curve; AUC = 1.000). Furthermore, we identified potential drug targets for SONFH through the Coremine medical database. Using the optimal machine learning model, we found that 2 PRGs, CASP1 and MLKL, were significantly correlated with the key predictive genes and exhibited higher expression levels in SONFH. Our analysis revealed the existence of 2 distinct PANoptosis molecular subtypes (C1 and C2) within SONFH. Importantly, we observed significant variations in the distribution of immune cells across these subtypes, with C2 displaying higher levels of immune cell infiltration. Gene set variation analysis indicated that C2 was closely associated with multiple immune responses. In conclusion, our study sheds light on the intricate relationship between PANoptosis and SONFH. We successfully developed a risk predictive model for SONFH patients and different SONFH subtypes. These findings enhance our understanding of the pathogenesis of SONFH and offer potential insights into therapeutic strategies.

Association of Oxidative Stress on Cognitive Function: A Bidirectional Mendelian Randomisation Study.

The purpose of this study was to investigate the relationship between oxidative stress and cognitive function, encompassing cognitive performance, intelligence, memory, reaction time, speech and vision by a bidirectional Mendelian randomisation study. Independent genetic variants associated with glutathione S-transferase (GST), catalase (CAT), superoxide dismutase (SOD), glutathione peroxidase (GPX), peroxiredoxin (PRDX), sulfhydryl oxidase (SOX) and thyroid peroxidase (TPO) were explored using a genome-wide association study (GWAS). The inverse variance weighted (IVW) or Wald ratio method was employed to ascertain the relationship between antioxidant enzymes and cognitive function. The MR analyses indicated that the MR effect estimates of GST (ß = 0.0352, P = 0.0047, FDR = 0.0164) and TPO (ß = 0.0531, P = 0.0003, FDR = 0.0021) were significantly associated with cognitive performance elevation. Furthermore, genetically predicted GST (ß = 0.0334, P = 0.0043, FDR = 0.0151) and TPO (ß = 0.0496, P = 0.0031, FDR = 0.0151) were found to be associated with high intelligence. Additionally, there were also some associations of SOX (ß = 0.0243, P = 0.0283, FDR = 0.066) on high cognitive performance, TPO (ß = 0.1189, P = 0.0315, FDR = 0.2205) on larger maximum digits remembered correctly, and SOX (ß = - 0.2435, P = 0.0395, FDR = 0.1185) on reaction time. Nevertheless, the associations between antioxidant enzymes and speech and linguistic disorders, as well as visual disturbances, were not significant. We did not find reverse causation between antioxidant enzymes and cognitive function traits. This study provides evidence of potential causal relationships between oxidative stress and cognitive function.

Unearthing naturally-occurring cyclic antibacterial peptides and their structural optimization strategies.

Natural products with antibacterial activity are highly desired globally to combat against multidrug-resistant (MDR) bacteria. Antibacterial peptide (ABP), especially cyclic ABP (CABP), is one of the abundant classes. Most of them were isolated from microbes, demonstrating excellent bactericidal effects. With the improved proteolytic stability, CABPs are normally considered to have better druggability than linear peptides. However, most clinically-used CABP-based antibiotics, such as colistin, also face the challenges of drug resistance soon after they reached the market, urgently requiring the development of next-generation succedaneums. We present here a detail review on the novel naturally-occurring CABPs discovered in the past decade and some of them are under clinical trials, exhibiting anticipated application potential. According to their chemical structures, they were broadly classified into five groups, including (i) lactam/lactone-based CABPs, (ii) cyclic lipopeptides, (iii) glycopeptides, (iv) cyclic sulfur-rich peptides and (v) multiple-modified CABPs. Their chemical structures, antibacterial spectrums and proposed mechanisms are discussed. Moreover, engineered analogs of these novel CABPs are also summarized to preliminarily analyze their structure-activity relationship. This review aims to provide a global perspective on research and development of novel CABPs to highlight the effectiveness of derivatives design in identifying promising antibacterial agents. Further research efforts in this area are believed to play important roles in fighting against the multidrug-resistance crisis.

Pretransplant use of immune checkpoint inhibitors for hepatocellular carcinoma: A multicenter, retrospective cohort study.

Immune checkpoint inhibitors (ICIs) as a downstaging or bridging therapy for liver transplantation (LT) in hepatocellular carcinoma patients are rapidly increasing. However, the evidence about the feasibility and safety of pre-LT ICI therapy is limited and controversial. To this end, a multicenter, retrospective cohort study was conducted in 11 Chinese centers. The results showed that 83 recipients received pre-LT ICI therapy during the study period. The median post-LT follow-up was 8.1 (interquartile range 3.3-14.6) months. During the short follow-up, 23 (27.7%) recipients developed allograft rejection, and 7 of them (30.4%) were diagnosed by liver biopsy. Multivariate logistics regression analysis showed that the time interval between the last administration of ICI therapy and LT (TLAT) ≥ 30 days was an independent protective factor for allograft rejection (odds ratio = 0.096, 95% confidence interval 0.026-0.357; P < .001). Multivariate Cox analysis showed that allograft rejection was an independent risk factor for overall survival (hazard ratio = 9.960, 95% confidence interval 1.006-98.610; P = .043). We conclude that patients who receive a pre-LT ICI therapy with a TLAT shorter than 30 days have a much higher risk of allograft rejection than those with a TLAT longer than 30 days. The presence of rejection episodes might be associated with higher post-LT mortality.

A Manganese-Based Nanodriver Coordinates Tumor Prevention and Suppression through STING Activation in Glioblastoma.

Glioblastoma (GBM), the most prevalent and aggressive primary malignant brain tumor, exhibits profound immunosuppression and demonstrates a low response rate to current immunotherapy strategies. Manganese cations (Mn2+) directly activate the cGAS/STING pathway and induce the unique catalytic synthesis of 2'3'-cGAMP to facilitate type I IFN production, thereby enhancing innate immunity. Here, a telodendrimer and Mn2+-based nanodriver (PLHM) with a small size is developed, which effectively target lymph nodes through the blood circulation and exhibit tumor-preventive effects at low doses of Mn2+ (3.7 mg kg-1). On the other hand, the PLHM nanodriver also exhibits apparent antitumor effects in GBM-bearing mice via inducing in vivo innate immune responses. The combination of PLHM with doxorubicin nanoparticles (PLHM-DOX NPs) results in superior inhibition of tumor growth in GBM-bearing mice due to the synergistic potentiation of STING pathway functionality by Mn2+ and the presence of cytoplasmic DNA. These findings demonstrate that PLHM-DOX NPs effectively stimulate innate immunity, promote dendritic cell maturation, and orchestrate cascaded infiltration of CD8 cytotoxic T lymphocytes within glioblastomas characterized by low immunogenicity. These nanodivers chelated with Mn2+ show promising potential for tumor prevention and antitumor effects on glioblastoma by activating the STING pathway.

Multilocular thymic cysts can be easily misdiagnosed as malignant tumor on computer tomography: A case report.

BACKGROUND: Multilocular thymic cyst (MTC) is a rare mediastinal lesion which is considered to occur in the process of acquired inflammation. It is usually characterized by well-defined cystic density and is filled with transparent liquid. CASE SUMMARY: We report on a 39-year-old male with a cystic-solid mass in the anterior mediastinum. Computer tomography (CT) imaging showed that the mass was irregular with unclear boundaries. After injection of contrast agent, there was a slight enhancement of stripes and nodules. According to CT findings, it was diagnosed as thymic cancer. CONCLUSION: After surgery, MTC accompanied by bleeding and infection was confirmed by pathological examination. The main lesson of this case was that malignant thymic tumor and MTC of the anterior mediastinum sometimes exhibit similar CT findings. Caution is necessary in clinical work to avoid misdiagnosis.

Discovery of a Monoclonal Antibody That Targets Cell-Surface Pseudaminic Acid of Acinetobacter baumannii with Direct Bactericidal Effect.

The therapeutic effects of antibodies include neutralization of pathogens, activation of the host complement system, and facilitation of phagocytosis of pathogens. However, antibody alone has never been shown to exhibit bactericidal activity. In this study, we developed a monoclonal antibody that targets the bacterial cell surface component Pseudaminic acid (Pse). This monoclonal antibody, Pse-MAB1, exhibited direct bactericidal activity on Acinetobacter baumannii strains, even in the absence of the host complements or other immune factors, and was able to confer a protective effect against A. baumannii infections in mice. This study provides new insight into the potential of developing monoclonal antibody-based antimicrobial therapy of multidrug resistant bacterial infections, especially those which occurred among immunocompromised patients.

Joint Exposure to Multiple Air Pollutants, Genetic Susceptibility, and Incident Dementia: A Prospective Analysis in the UK Biobank Cohort.

Objectives: This study aimed to evaluate the joint effects of multiple air pollutants including PM2.5, PM10, NO2, and NOx with dementia and examined the modifying effects of genetic susceptibility. Methods: This study included 220,963 UK Biobank participants without dementia at baseline. Weighted air pollution score reflecting the joint exposure to multiple air pollutants were constructed by cross-validation analyses, and inverse-variance weighted meta-analyses were performed to create a pooled effect. The modifying effect of genetic susceptibility on air pollution score was assessed by genetic risk score and APOE ε4 genotype. Results: The HR (95% CI) of dementia for per interquartile range increase of air pollution score was 1.13 (1.07∼1.18). Compared with the lowest quartile (Q1) of air pollution score, the HR (95% CI) of Q4 was 1.26 (1.13∼1.40) (P trend = 2.17 × 10-5). Participants with high air pollution score and high genetic susceptibility had higher risk of dementia compared to those with low air pollution score and low genetic susceptibility. Conclusion: Our study provides evidence that joint exposure to multiple air pollutants substantially increases the risk of dementia, especially among individuals with high genetic susceptibility.

Tailoring raloxifene into single-component molecular crystals possessing multilevel stimuli-responsive room-temperature phosphorescence.

Simultaneously achieving room-temperature phosphorescence (RTP) and multiple-stimuli responsiveness in a single-component system is of significance but remains challenging. Crystallization has been recognized to be a workable strategy to fulfill the above task. However, how the molecular packing mode affects the intersystem crossing and RTP lifetime concurrently remains unclear so far. Herein, four economic small-molecular compounds, analogues of the famous drug raloxifene (RALO), are facilely synthesized and further explored as neat single-component and stimuli-responsive RTP emitters via crystallization engineering. Thanks to their simple structures and high ease to crystallize, these raloxifene analogues function as models to clarify the important role of molecular packing in the RTP and stimuli-responsiveness properties. Thorough combination of the single-crystal structure analysis and theoretical calculations clearly manifests that the tight antiparallel molecular packing mode is the key point to their RTP behaviors. Interestingly, harnessing the controllable and reversible phase transitions of the two polymorphs of RALO-OAc driven by mechanical force, solvent vapor, and heat, a single-component multilevel stimuli-responsive platform with tunable emission color is established and further exploited for optical information encryption. This work would shed light on the rational design of multi-stimuli responsive RTP systems based on single-component organics.

Incorporating genetic similarity of auxiliary samples into eGene identification under the transfer learning framework.

BACKGROUND: The term eGene has been applied to define a gene whose expression level is affected by at least one independent expression quantitative trait locus (eQTL). It is both theoretically and empirically important to identify eQTLs and eGenes in genomic studies. However, standard eGene detection methods generally focus on individual cis-variants and cannot efficiently leverage useful knowledge acquired from auxiliary samples into target studies. METHODS: We propose a multilocus-based eGene identification method called TLegene by integrating shared genetic similarity information available from auxiliary studies under the statistical framework of transfer learning. We apply TLegene to eGene identification in ten TCGA cancers which have an explicit relevant tissue in the GTEx project, and learn genetic effect of variant in TCGA from GTEx. We also adopt TLegene to the Geuvadis project to evaluate its usefulness in non-cancer studies. RESULTS: We observed substantial genetic effect correlation of cis-variants between TCGA and GTEx for a larger number of genes. Furthermore, consistent with the results of our simulations, we found that TLegene was more powerful than existing methods and thus identified 169 distinct candidate eGenes, which was much larger than the approach that did not consider knowledge transfer across target and auxiliary studies. Previous studies and functional enrichment analyses provided empirical evidence supporting the associations of discovered eGenes, and it also showed evidence of allelic heterogeneity of gene expression. Furthermore, TLegene identified more eGenes in Geuvadis and revealed that these eGenes were mainly enriched in cells EBV transformed lymphocytes tissue. CONCLUSION: Overall, TLegene represents a flexible and powerful statistical method for eGene identification through transfer learning of genetic similarity shared across auxiliary and target studies.