Search

Abstract Background The C-type lectin family 18 (CLEC18) with lipid and glycan binding capabilities is important to metabolic regulation and innate immune responses against viral infection. However, human CLEC18 comprises three paralogous genes with highly similar sequences, making it challenging to distinguish genetic variations, expression patterns, and biological functions of individual CLEC18 paralogs. Additionally, the evolutionary relationship between human CLEC18 and its counterparts in other species remains unclear. Methods To identify the sequence variation and evolutionary divergence of human CLEC18 paralogs, we conducted a comprehensive analysis using various resources, including human and non-human primate reference genome assemblies, human pangenome assemblies, and long-read-based whole-genome and -transcriptome sequencing datasets. Results We uncovered paralogous sequence variants (PSVs) and polymorphic variants (PVs) of human CLEC18 proteins, and identified distinct signatures specific to each CLEC18 paralog. Furthermore, we unveiled a novel segmental duplication for human CLEC18A gene. By comparing CLEC18 across human and non-human primates, our research showed that the CLEC18 paralogy probably occurred in the common ancestor of human and closely related non-human primates, and the lipid-binding CAP/SCP/TAPS domain of CLEC18 is more diverse than its glycan-binding CTLD. Moreover, we found that certain amino acids alterations at variant positions are exclusive to human CLEC18 paralogs. Conclusions Our findings offer a comprehensive profiling of the intricate variations and evolutionary characteristics of human CLEC18.

Abstract Background The human leukocyte antigen (HLA) genes, exhibiting significant genetic diversity, are associated with susceptibility to various clinical diseases and diverse in drug responses. High costs of HLA sequencing and the population-specific architecture of this genetic region necessitate the establishment of a population-specific HLA imputation reference panel. Moreover, there is a lack of understanding about the genetic and phenotypic landscape of HLA variations within the Taiwanese population. Methods We created models for a Taiwanese-specific HLA imputation reference panel. These models were trained with the array genotype data and HLA sequencing data from 845 Taiwanese subjects. HLA imputation was applied for 59,448 Taiwanese subjects to characterize the HLA allele and haplotype frequencies. Additionally, a phenome-wide association study (PheWAS) was conducted to identify the phenotypes associated with HLA variations. The association of the biallelic HLA variants with the binary and quantitative traits were evaluated with additive logistic and linear regression models, respectively. Furthermore, an omnibus test with likelihood-ratio test was applied for each HLA amino acid position in the multiallelic HLA amino acid polymorphisms to compare the difference between a fitted model and a null model following a χ2 distribution of n-1 degree of freedom at a position with n residues. Finally, we estimated the prevalence of adverse drug reactions (ADR)-related HLA alleles in the Taiwanese population. Results In this study, the reference panel models displayed remarkable accuracy, with averages of 99.3%, 98.9%, and 99.1% for 2-, 4-, 6-digit alleles of the eight classical HLA genes, respectively. For PheWAS, a total of 18,136 significant associations with HLA variants across 26 phenotypes are identified (p < 5×10-8), highlighting the pleiotropy feature of the HLA region. Among the independent signals, 15 are novel, including the association of HLA-B pos 138 variation with ankylosing spondylitis (AS), and rs9266290 and rs9266292 with allergy. Through an analysis spanning the entire HLA region, we identified clusters of phenotype correlations. Finally, the carriers of pharmacogenomic related HLA alleles, including HLA-C*01:02 (35.86%), HLA-B*58:01 (20.9%), and HLA-B*15:02 (8.38%), were characterized in the Taiwanese general population. Conclusions We successfully delivered the HLA imputation for 59,448 Taiwanese subjects and characterized the genetic and phenotypic landscapes of the HLA variations. In addition, we quantified the estimated prevalence of the ADR-related HLA alleles in the Taiwanese population. The developed HLA imputation reference panel could be used for estimation of population HLA allele frequencies, which can facilitate further studies in the role of HLA variants in a wider range of phenotypes in the population.

BackgroundRheumatoid arthritis (RA) is a chronic disease and may worsen over time. Today, nurse-led case management (NLCM) has been recommended to improve clinical outcomes for chronic disease patients, yet little is known regarding its impact on pain, fatigue, and C-reactive protein (CRP) among RA patients. We aimed to explore this issue among such groups via a two-group pre- and post-test approach.MethodsAll subjects were recruited from one hospital in Taiwan from January 2017 to June 2018 and assigned to either a 6-month NLCM program in addition to usual care or to a control group that received usual care only. All of them were followed for 2 years. Outcomes of interests were compared at four time points: baseline, the third day after NLCM completion, and at 6 and 24 months after NLCM. Effects between them were tested using the generalized estimating equations (GEE) model after adjusting for differences at baseline.ResultsA total of 50 patients in the NLCM group and 46 in the control group were recruited for data analysis. Results from the GEE model indicated that integrating NLCM into conventional care benefited patients in decreasing levels of pain and fatigue, as well as CRP value. These improvements were still observed for 2 years after NLCM.ConclusionNLCM was shown to be helpful in lowering pain, fatigue, and CRP, which implies that NLCM may be a reference in the provision of tailored care for those affected by rheumatism.

Abstract Immunoglobulin A nephropathy (IgAN) is an autoimmune disease characterized by abnormal IgA deposition in glomerulus. Current diagnosis of IgAN still depends on renal biopsy, an invasive method that might increase the risk of clinical outcomes. Therefore, we aimed to explore the characteristics of T cell repertoire in IgAN from peripheral blood samples for identifying innovative diagnostic biomarkers. Herein, we included 8 IgAN patients, 25 non-IgAN patients, and 10 healthy controls in the study. A high-throughput immune repertoire sequencing was conducted to investigate the T-cell receptor beta-chain (TCRβ) repertoire of peripheral blood. Characteristics of TCRβ repertoire were assessed for these three distinct groups. A reduced TCRβ repertoire diversity was observed in IgAN patients compared to non-IgAN and healthy individuals. A skewed distribution toward shorter TCRβ complementarity determining region (CDR3) length was found in non-IgAN relative to IgAN patients. In addition, the differences in usages of five TRBV genes (TRBV5-4, TRBV6-4, TRBV12-1, TRBV16, and TRBV21-1) were identified between IgAN, non-IgAN, and healthy subjects. Of note, the TRBV6-4 gene, which is associated with mucosal-associated invariant T (MAIT) cells, exhibited higher usage in IgAN patients, suggesting potential importance of MAIT cells in IgAN. In short, our findings supported TCR repertoire characteristics as potential biomarkers for IgAN diagnosis.

Abstract Background yqiC is required for colonizing the Salmonella enterica serovar Typhimurium (S. Typhimurium) in human cells; however, how yqiC regulates nontyphoidal Salmonella (NTS) genes to influence bacteria–host interactions remains unclear. Methods The global transcriptomes of S. Typhimurium yqiC-deleted mutant (ΔyqiC) and its wild-type strain SL1344 after 2 h of in vitro infection with Caco-2 cells were obtained through RNA sequencing to conduct comparisons and identify major yqiC-regulated genes, particularly those involved in Salmonella pathogenicity islands (SPIs), ubiquinone and menaquinone biosynthesis, electron transportation chains (ETCs), and carbohydrate/energy metabolism. A Seahorse XFp Analyzer and assays of NADH/NAD+ and H2O2 were used to compare oxygen consumption and extracellular acidification, glycolysis parameters, adenosine triphosphate (ATP) generation, NADH/NAD+ ratios, and H2O2 production between ΔyqiC and SL1344. Results After S. Typhimurium interacts with Caco-2 cells, yqiC represses gene upregulation in aspartate carbamoyl transferase, type 1 fimbriae, and iron–sulfur assembly, and it is required for expressing ilvB operon, flagellin, tdcABCD, and dmsAB. Furthermore, yqiC is required for expressing mainly SPI-1 genes and specific SPI-4, SPI-5, and SPI-6 genes; however, it diversely regulates SPI-2 and SPI-3 gene expression. yqiC significantly contributes to menD expression in menaquinone biosynthesis. A Kyoto Encyclopedia of Genes and Genomes analysis revealed the extensive association of yqiC with carbohydrate and energy metabolism. yqiC contributes to ATP generation, and the analyzer results demonstrate that yqiC is required for maintaining cellular respiration and metabolic potential under energy stress and for achieving glycolysis, glycolytic capacity, and glycolytic reserve. yqiC is also required for expressing ndh, cydA, nuoE, and sdhB but suppresses cyoC upregulation in the ETC of aerobically and anaerobically grown S. Typhimurium; priming with Caco-2 cells caused a reversed regulation of yiqC toward upregulation in these ETC complex genes. Furthermore, yqiC is required for maintaining NADH/NAD+ redox status and H2O2 production. Conclusions Specific unreported genes that were considerably regulated by the colonization-associated gene yqiC in NTS were identified, and the key role and tentative mechanisms of yqiC in the extensive modulation of virulence factors, SPIs, ubiquinone and menaquinone biosynthesis, ETCs, glycolysis, and oxidative stress were discovered.

Abstract Background Human traits, diseases susceptibility, and clinical outcomes vary hugely among individuals. Despite a fundamental understanding of genetic (or environmental) contributions, the detailed mechanisms of how genetic variation impacts molecular or cellular behaviours of a gene, and subsequently leads to such variability remain poorly understood. Methods Here, in addition to phenome-wide correlations, we leveraged multiomics to exploit mechanistic links, from genetic polymorphism to protein structural or functional changes and a cross-omics perturbation landscape of a germline variant. Results We identified a missense cis-acting expression quantitative trait locus in CLEC18A (rs75776403) in which the altered residue (T151→M151) disrupts the lipid-binding ability of the protein domain. The altered allele carriage led to a metabolic and proliferative shift, as well as immune deactivation, therefore determines human anthropometrics (body height), kidney, and hematological traits. Conclusions Collectively, we uncovered genetic pleiotropy in human complex traits and diseases via CLEC18A rs75776403-regulated pathways.

IntroductionCytoreductive surgery (CRS) and hyperthermic intraperitoneal chemotherapy (HIPEC) are considered for patients with peritoneal metastasis (PM). However, patients selection that relies on conventional prognostic factors is not yet optimal. In this study, we performed whole exome sequencing (WES) to establish tumor molecular characteristics and expect to identify prognosis profiles for PM management. MethodsIn this study, blood and tumor samples were collected from patients with PM before HIPEC. Tumor molecular signatures were determined using WES. Patient cohort was divided into responders and non-responders according to 12-month progression-free survival (PFS). Genomic characteristics between the two cohorts were compared to study potential targets. ResultsIn total, 15 patients with PM were enrolled in this study. Driver genes and enriched pathways were identified from WES results. AGAP5 mutation was found in all responders. This mutation was significantly associated with better OS (p = 0.00652). ConclusionsWe identified prognostic markers that might be useful to facilitate decision-making before CRS/HIPEC.

Dengue is an arboviral disease that has spread globally and become a major public health concern. A small proportion of patients may progress from symptomatic dengue fever (DF) to dengue hemorrhagic fever (DHF) or dengue shock syndrome (DSS). Findings from a previous genome-wide association study (GWAS) demonstrated that variations in the major histocompatibility complex (MHC) class I chain-related B (MICB) and the phospholipase C epsilon 1 (PLCE1) genes were related to DSS in a Vietnamese population. This study investigated associations of variations in MICB (rs3132468) and PLCE1 (rs3740360, rs3765524) with dengue severity and thrombocytopenia in both the Indonesian and Taiwanese populations. We sampled 160 patients from the Indonesian population and 273 patients from the Taiwanese population. None of the patients had DSS in the Taiwanese population. Based on age demographics, we found that dengue is more prevalent among younger individuals in the Indonesian population, whereas it has a greater impact on adults in the Taiwanese population. Our results showed the association between MICB rs3132468 and DSS. In addition, an association was identified between PLCE1 rs3740360 and DHF in secondary dengue in Indonesian patients. However, there is no association of MICB or PLCE1 variants with thrombocytopenia. This study highlights the value of genetic testing, which might be included in the clinical pathway for specific patients who can be protected from severe dengue.

Myocardial bridging (MB) is a congenital coronary artery anomaly and an important cause of angina. The genetic basis of MB is currently unknown. This study used a whole-exome sequencing technique and analyzed genotypic differences. Eight coronary angiography-confirmed cases of severe MB and eight age- and sex-matched control patients were investigated. In total, 139 rare variants that are potentially pathogenic for severe MB were identified in 132 genes. Genes with multiple rare variants or co-predicted by ClinVar and CADD/REVEL for severe MB were collected, from which heart-specific genes were selected under the guidance of tissue expression levels. Functional annotation indicated significant genetic associations with abnormal skeletal muscle mass, cardiomyopathies, and transmembrane ion channels. Candidate genes were reviewed regarding the functions and locations of each individual gene product. Among the gene candidates for severe MB, rare variants in DMD, SGCA, and TTN were determined to be the most crucial. The results suggest that altered anchoring proteins on the cell membrane and intracellular sarcomere unit of cardiomyocytes play a role in the development of the missed trajectory of coronary vessels. Additional studies are required to support the diagnostic application of cardiac sarcoglycan and dystroglycan complexes in patients with severe MB.

According to the National Comprehensive Cancer Network and the American Society of Clinical Oncology, the standard treatment for pancreatic cancer (PC) is gemcitabine and fluorouracil. Other chemotherapeutic agents have been widely combined. However, drug resistance remains a huge challenge, leading to the ineffectiveness of cancer therapy. Therefore, we are trying to discover new treatments for PC by utilizing genomic information to identify PC-associated genes as well as drug target genes for drug repurposing. Genomic information from a public database, the cBio Cancer Genomics Portal, was employed to retrieve the somatic mutation genes of PC. Five functional annotations were applied to prioritize the PC risk genes: Kyoto Encyclopedia of Genes and Genomes; biological process; knockout mouse; Gene List Automatically Derived For You; and Gene Expression Omnibus Dataset. DrugBank database was utilized to extract PC drug targets. To narrow down the most promising drugs for PC, CMap Touchstone analysis was applied. Finally, ClinicalTrials.gov and a literature review were used to screen the potential drugs under clinical and preclinical investigation. Here, we extracted 895 PC-associated genes according to the cBioPortal database and prioritized them by using five functional annotations; 318 genes were assigned as biological PC risk genes. Further, 216 genes were druggable according to the DrugBank database. CMap Touchstone analysis indicated 13 candidate drugs for PC. Among those 13 drugs, 8 drugs are in the clinical trials, 2 drugs were supported by the preclinical studies, and 3 drugs are with no evidence status for PC. Importantly, we found that midostaurin (targeted PRKA) and fulvestrant (targeted ESR1) are promising candidate drugs for PC treatment based on the genomic-driven drug repurposing pipelines. In short, integrated analysis using a genomic information database demonstrated the viability for drug repurposing. We proposed two drugs (midostaurin and fulvestrant) as promising drugs for PC.

Objective: Due to the shared pathogenesis of asthma and rheumatoid arthritis (RA), patients with asthma were found to have a higher risk of RA. While the benefits and safety of Chinese herbal medicine (CHM) for asthma have been reported, the scientific evidence regarding its effect on RA is limited. This longitudinal cohort study aimed to determine the relation between CHM use and RA risk in patients with asthma.Methods: Using the nationwide claims data, we enrolled 33,963 patients 20–80 years of age who were newly diagnosed with asthma and simultaneously free of RA between 2000 and 2007. From this sample, we utilized propensity score matching to create sets of participants as treatment and control groups, which comprised 13,440 CHM users and 13,440 non-CHM users. The incidence rate and hazard ratio (HR) for RA between the two groups were estimated at the end of 2013. A Cox proportional hazards model was constructed to examine the impact of the CHM use on the risk of RA.Results: The cumulative incidence of RA was substantially lower in the CHM user group. In the follow-up period, 214 patients in the CHM user group (1.92 per 1,000 person-years) and 359 patients in the non-CHM user group (2.92 per 1,000 person-years) developed RA (adjusted HR = 0.63, 95% confidence interval: 0.54–0.75). Of the commonly-prescribed formulae, nine CHM products were associated with a lower RA risk: Xiao-Qing-Long-Tang, Ma-Xing-Gan-Shi-Tang, Ding-Chuan-Tang, Xin-Yi-Qing-Fei-Tang, Bei Mu, Jie Geng, Xing Ren, Da Huang, and San Chi.Conclusion: This study found that patients with asthma who received CHM treatment, in addition to the conventional therapy, had a lower risk of RA. Use of CHM treatment may be integrated into conventional therapy to reduce subsequent RA risk among asthma patients.

BackgroundRecent evidence suggests that the use of orthopedic manual therapy (OMT) may lessen the subsequent risk of low back pain (LBP), but this association has not been examined among pregnant women who are at higher risk of LBP. This study aims to determine whether the addition of OMT to conventional LBP treatment before pregnancy could decrease the subsequent risk of LBP during pregnancy.MethodsFrom Taiwan's National Health Insurance Research Database, we identified 68,960 women, 20–55 years of age, with first pregnancy between 2001 and 2012. We then performed a nested case-control study in which 3,846 women with newly diagnosed LBP were matched to 3,846 controls according to age and cohort entry year. Multivariate conditional logistic regression was employed to estimate the association between OMT use before pregnancy and LBP during pregnancy.ResultsOMT users had a lower risk of LBP than did non-users, with an adjusted OR of 0.86 (95% CI, 0.78–0.93). Subgroup analysis showed that women with high intensity use of OMT treatment prior to pregnancy reported the lowest level of LBP during pregnancy by nearly 30%.ConclusionThe pre-pregnancy use of OMT treatment significantly decreased LBP risk during pregnancy, especially with high-intensity use. Thus, clinicians may consider recommending OMT for pregnant women to avoid possible obstetric complications during the pregnancy.

Cancer immunotherapies are treatments that use drugs or cells to activate patients’ own immune systems against cancer cells. Among them, cancer vaccines have recently been rapidly developed. Based on tumor-specific antigens referred to as neoantigens, these vaccines can be in various forms such as messenger (m)RNA and synthetic peptides to activate cytotoxic T cells and act with or without dendritic cells. Growing evidence suggests that neoantigen-based cancer vaccines possess a very promising future, yet the processes of immune recognition and activation to relay identification of a neoantigen through the histocompatibility complex (MHC) and T-cell receptor (TCR) remain unclear. Here, we describe features of neoantigens and the biological process of validating neoantigens, along with a discussion of recent progress in the scientific development and clinical applications of neoantigen-based cancer vaccines.

Aldosterone-producing adenoma (APA) is a benign adrenal tumor that results in persistent hyperaldosteronism. As one major subtype of primary aldosteronism, APA leads to secondary hypertension that is associated with immune dysregulation. However, how the adaptive immune system, particularly the T-cell population, is altered in APA patients remains largely unknown. Here, we performed TCR sequencing to characterize the TCR repertoire between two age-matched groups of patients: one with APA and the other one with essential hypertension (EH). Strikingly, we found a significant reduction of TCR repertoire diversity in the APA group. Analyses on TCR clustering and antigen annotation further showed that the APA group possessed lower diversity in TCR clonotypes with non-common antigen-specific features, compared with the EH group. In addition, our results indicated that the strength of correlation between generation probabilities and frequencies of TCR clonotypes was significantly higher in the APA group than that in the EH group. Finally, we observed that clinical features, including plasma aldosterone level, aldosterone–renin ratio, and blood sodium level, were positively associated with the strength of correlation between generation and abundance of TCR clonotypes in the APA group. Our findings unveiled the correlation between T-cell immune repertoire and APA, suggesting a critical role of such adrenal adenoma in the T-cell immunity of patients with hypertension.

Ya-Lang Huang et al. report a mechanism for TLR3-mediated signaling after immune simulation and influenza virus infection by way of the co-receptor CLEC18A. This study found that a single amino acid change in CLEC18A(S339R) can enhance the production of type I and type III interferons to suppress viral replication, and increase mice survival rate after flu infection infection.

Cancer is the leading cause of death in Taiwan, and the overall incidence rate has gradually increased. The four most common cancers in Taiwan are colorectal, lung, breast and liver cancers. With the rise in incidence, the clinical use and costs of all anticancer drugs have steadily increased. The costs of novel therapeutics, such as targeted therapies and immunotherapy were accounted almost two-third of all antineoplastic agents in Taiwan. Moving forward, it will be necessary to discuss the economic impacts to clinical use of new therapeutics, while continuing to monitor and improve the quality of cancer therapy. In this review, we describe the epidemiology, disease screening policies and medication treatment policies for colorectal, lung, breast and liver cancer. We focus on the recent developments in cancer therapeutics, discuss the use of biomarkers, and finally consider the costs and the recent advances of anticancer medications in Taiwan.

Adolescent idiopathic scoliosis (AIS) is a three-dimensional spinal curvature deformity that appears in the adolescent period. In this study, we performed whole-exome sequencing on 11 unrelated Taiwanese patients with a Cobb’s angle greater than 40 degrees. Our results identified more than 200 potential pathogenic rare variants, however, most of which were carried only by one individual. By in silico pathogenicity annotation studies, we found that TTN, CLCN1, and SOX8 were the most important genes, as multiple pathogenic variants were within these genes. Furthermore, biological functional annotation indicated critical roles of these AIS candidate genes in the skeletal muscle. Importantly, a pathogenic variant on SOX8 was shared by over 35% of the patients. These results highlighted TTN, CLCN1, and SOX8 as the most likely susceptibility genes for severe AIS.

Store-operated Ca2+ channel (SOC)-regulated Ca2+ entry is involved in inflammation and colorectal cancer (CRC) progression, but clinically applicable treatments targeting this mechanism are lacking. Recent studies have shown that nonsteroidal anti-inflammatory drugs (NSAIDs) not only inhibit inflammation but they also suppress Ca2+ entry via SOC (SOCE). Therefore, delineating the mechanisms of SOCE inhibition by NSAIDs may lead to new CRC treatments. In this study, we tested eight candidate NSAIDs in Ca2+ imaging experiments and found that Aspirin and Sulindac were the most effective at suppressing SOCE. Furthermore, time-lapse FRET imaging using TIRF microscopy and ground state depletion (GSD) super-resolution (SR) imaging revealed that SOC was inhibited by Aspirin and Sulindac via different mechanisms. Aspirin quickly interrupted the STIM1-Orai1 interaction, whereas Sulindac mainly suppressed STIM1 translocation. Additionally, Aspirin and Sulindac both inhibited metastasis-related endpoints in CRC cells. Both drugs were used throughout the study at doses that suppressed CRC cell migration and invasion without altering cell survival. This is the first study to reveal the differential inhibitory mechanisms of Aspirin and Sulindac on SOC activity. Thus, our results shed new light on the therapeutic potential of Aspirin for CRC and SOCE-related diseases.

Overweight and obese are risk factors for various diseases. In Taiwan, the combined prevalence of overweight and obesity has increased dramatically. Here, we conducted a genome-wide association study (GWAS) on four adiposity traits, including body-mass index (BMI), body fat percentage (BF%), waist circumference (WC), and waist-hip ratio (WHR), using the data for more than 21,000 subjects in Taiwan Biobank. Associations were evaluated between 6,546,460 single-nucleotide polymorphisms (SNPs) and adiposity traits, yielding 13 genome-wide significant (GWS) adiposity-associated trait-loci pairs. A known gene, FTO, as well as two BF%-associated loci (GNPDA2-GABRG1 [4p12] and RNU6-2-PIAS1 [15q23]) were identified as pleiotropic effects. Moreover, RALGAPA1 was found as a specific genetic predisposing factor to high BMI in a Taiwanese population. Compared to other populations, a slightly lower heritability of the four adiposity traits was found in our cohort. Surprisingly, we uncovered the importance of neural pathways that might influence BF%, WC and WHR in the Taiwanese (East Asian) population. Additionally, a moderate genetic correlation between the WHR and BMI (γg = 0.52; p = 2.37×10-9) was detected, suggesting different genetic determinants exist for abdominal adiposity and overall adiposity. In conclusion, the obesity-related genetic loci identified here provide new insights into the genetic underpinnings of adiposity in the Taiwanese population.

With obesity rate consistently increasing, a strong relationship between obesity and fatty liver disease has been discovered. More than 90% of bariatric surgery patients also have non-alcoholic fatty liver diseases (NAFLDs). NAFLD and non-alcoholic steatohepatitis (NASH), which are the hepatic manifestations of metabolic syndrome, can lead to liver carcinogenesis. Unfortunately, there is no effective medicine that can be used to treat NASH or liver cancer. Thus, it is critically important to understand the mechanism underlying the development of these diseases. Extensive evidence suggests that microRNA 22 (miR-22) can be a diagnostic marker for liver diseases as well as a treatment target. This review paper focuses on the roles of miR-22 in metabolism, steatosis, and liver carcinogenesis. Literature search is limited based on the publications included in the PubMed database in the recent 10 years.

Abstract Dramatic advances in immune therapy have emerged as a promising strategy in cancer therapeutics. In addition to chemotherapy and radiotherapy, inhibitors targeting immune-checkpoint molecules such as cytotoxic T-lymphocyte antigen-4 (CTLA-4), programmed cell death receptor-1 (PD-1) and its ligand (PD-L1) demonstrate impressive clinical benefits in clinical trials. In this review, we present background information about therapies involving PD-1/PD-L1 blockade and provide an overview of current clinical trials. Furthermore, we present recent advances involving predictive biomarkers associated with positive therapeutic outcomes in cancer immunotherapy.

Abstract Background Insulin-like growth factor 1 (IGF-1) is an important pleiotropic hormone that exerts neuroprotective and neuroreparative effects after a brain injury. However, the roles of IGF-1 variants in mild traumatic brain injury (mTBI) are not yet fully understood. This study attempted to elucidate the effects of IGF-1 variants on the risk and neuropsychiatric outcomes of mTBI. Methods Based on 176 recruited mTBI patients and 1517 control subjects from the Taiwan Biobank project, we first compared the genotypic distributions of IGF-1 variants between the two groups. Then, we analyzed associations of IGF-1 variants with neuropsychiatric symptoms after mTBI, including anxiety, depression, dizziness, and sleep disturbances. Functional annotation of IGF-1 variants was also performed through bioinformatics databases. Results The minor allele of rs7136446 was over-represented in mTBI patients compared to community-based control subjects. Patients carrying minor alleles of rs7136446 and rs972936 showed more dizziness and multiple neuropsychiatric symptoms after brain injury. Conclusions IGF-1 variants were associated with the risk and neuropsychiatric symptoms of mTBI. The findings highlight the important role of IGF-1 in the susceptibility and clinical outcomes of mTBI.

Abstract Background Genetic factors, dysregulation in the endocrine system, cytokine and paracrine factors are implicated in the pathogenesis of familial short stature (FSS). Nowadays, the treatment choice for FSS is limited, with only recombinant human growth hormone (rhGH) being available. Methods Herein, starting from the identification of 122 genetic loci related to FSS, we adopted a genetic-driven drug discovery bioinformatics pipeline based on functional annotation to prioritize crucial biological FSS-related genes. These genes were suggested to be potential targets for therapeutics. Results We discovered five druggable subnetworks, which contained seven FSS-related genes and 17 druggable targerts. Conclusions This study provides a valuable drug repositioning accompanied by corresponding targetable gene clusters for FSS therapy.

Atopic Dermatitis (AD) is a chronic and relapsing skin disease. The medications for treating AD are still limited, most of them are topical corticosteroid creams or antibiotics. The current study attempted to discover potential AD treatments by integrating a gene network and genomic analytic approaches. Herein, the Single Nucleotide Polymorphism (SNPs) associated with AD were extracted from the GWAS catalog. We identified 70 AD-associated loci, and then 94 AD risk genes were found by extending to proximal SNPs based on r2 > 0.8 in Asian populations using HaploReg v4.1. Next, we prioritized the AD risk genes using in silico pipelines of bioinformatic analysis based on six functional annotations to identify biological AD risk genes. Finally, we expanded them according to the molecular interactions using the STRING database to find the drug target genes. Our analysis showed 27 biological AD risk genes, and they were mapped to 76 drug target genes. According to DrugBank and Therapeutic Target Database, 25 drug target genes overlapping with 53 drugs were identified. Importantly, dupilumab, which is approved for AD, was successfully identified in this bioinformatic analysis. Furthermore, ten drugs were found to be potentially useful for AD with clinical or preclinical evidence. In particular, we identified filgotinub and fedratinib, targeting gene JAK1, as potential drugs for AD. Furthermore, four monoclonal antibody drugs (lebrikizumab, tralokinumab, tocilizumab, and canakinumab) were successfully identified as promising for AD repurposing. In sum, the results showed the feasibility of gene networking and genomic information as a potential drug discovery resource.

BackgroundLong noncoding RNA (lncRNA) mediates the pathogenesis of various diseases, including cancer and cardiovascular, infectious, and metabolic diseases. This study examined the role of lncRNA NTT in the development and progression of cancer.MethodsThe expression of NTT was determined using tissues containing complementary DNA (cDNA) from patients with liver, lung, kidney, oral, and colon cancers. The expression of cis-acting genes adjacent to the NTT locus (CTGF, STX7, MYB, BCLAF1, IFNGR1, TNFAIP3, and HIVEP2) was also assessed. We used knockdown and chromatin immunoprecipitation (ChIP) assays to identify the cis-acting genes that interact with NTT.ResultsNTT was most significantly downregulated in hepatocellular carcinoma (HCC), while a higher NTT level correlated with a shorter survival time of patients with HCC. Multivariate analysis indicated NTT was not an independent predictor for overall survival. MYB was significantly upregulated, and its increased expression was associated with dismal survival in HCC patients, similar to the results for NTT. NTT knockdown significantly decreased cellular migration. ChIP of HCC cell lines revealed that NTT is regulated by the transcription factor ATF3 and binds to the MYB promoter via the activated complex. Additionally, when NTT was knocked down, the expression of MYB target genes such as Bcl-xL, cyclinD1, and VEGF was also downregulated. NTT could play a positive or negative regulator for MYB with a context-dependent manner in both HCC tissues and animal model.ConclusionOur study suggests that NTT plays a key role in HCC progression via MYB-regulated target genes and may serve as a novel therapeutic target.

Immune checkpoint inhibitors (ICIs) have been approved to treat patients with various cancer types, including lung cancer, in many countries. This study aims to investigate the effectiveness and safety of ICIs under different treatment conditions of non-small cell lung cancer patients. A population-based retrospective cohort study was conducted using the electronic health records of three medical centers in Taiwan. From January 01, 2016, to November 30, 2018, a total of 91 ICIs and 300 traditional chemotherapy users who had undergone stage III and IV lung cancer treatment were included in the study. We performed the randomized matched pair design by selecting a Chemotherapy subject for each ICI patient in the sample population. All subjects were monitored from the date of taking ICIs or chemotherapy drugs until the event of death, loss to follow-up, or were occurred with any defined adverse events. Kaplan-Meier estimators and cox proportional hazard regression models were used to compute the overall survival, efficacy, and safety of the ICIs group. The median overall survival (OS) in the ICI and Chemo groups after matching was 11.2 months and 10.5 months, respectively. However, the results showed no significant OS differences between ICIs and chemo groups for both before and after matching (HR,1.30; 95%CI, 0.68-2.46; p=0.428 before matching and HR,0.96; 95CI%, 0.64-1.44; p=0.838 after matching). We observed that with the higher amount of PD-L1, the length of the patients’ overall survival was (positive vs. negative PD-L1, HR,0.21; 95%CI, 0.05-0.80; p=0.022). The incidences of serious adverse drug events above grade 3 in the ICIs and traditional chemo groups were 12.7% and 21.5%, respectively. We also found that the number of AEs was less in ICIs than in the Chemo group, and the AEs that occurred after treatments were observed earlier in the ICIs compared to the Chemo group. ICIs drugs were observed to be safer than traditional chemotherapy as they had a lower risk of serious adverse drug events. It is necessary to pay attention to immune-related side effects and provide appropriate treatment. Furthermore, the patient’s physical status and PD-L1 test can be used to evaluate the clinical effectiveness of ICIs.

Abstract Background Infection is the second most common cause of mortality for patients with end-stage renal disease (ESRD), accompanying with immune dysfunction. Endothelin (EDN) is known to be related to inflammation; however, it is unknown whether genetic variants of the EDN gene family are associated with increased risk of hospitalized infection events. Methods Nineteen tagging single-nucleotide polymorphisms (tSNPs) of the EDN gene family were selected for genotyping a cohort of 190 ESRD patients. Patient demographics were recorded, the subtypes of infection events were identified, and association analysis between the EDN genetic variants and hospitalized infection events was performed. Results In this study, 106 patients were hospitalized for infection events. The leading events were pneumonia, bacteremia, and cellulitis. The minor allele of rs260741, rs197173, and rs926632 SNPs of EDN3 were found to be associated with reduced risk of hospitalized bacteremia events. Conclusions The minor allele of rs260741, rs197173, and rs926632 in EDN3 were associated with reduced risk of hospitalized bacteremia events in ESRD patients.

Dengue is a viral infection caused by the dengue virus (DENV). Dengue infection is a self-limited acute febrile illness caused by four serotypes of DENV (DENV-1~4). Early recognition of high-risk patients would be helpful to reduce mortality rates and prevent severe dengue. Our study aimed to identify factors related to dengue hemorrhagic fever (DHF) based on admission-day data, and further to understand the distribution of biochemical laboratory data in dengue patients. This retrospective study was conducted in hospitals in Yogyakarta city, Indonesia, and involved febrile patients who were admitted to the hospital with a diagnosis of dengue during 2018 and 2020. Logistic regression models were used to identify variables related to DHF. In this study, 1087 patients were included as suspected dengue patients, among them 468 had dengue fever (DF) and 619 had DHF. Over half of the DHF patients were males (55.9%) with an average age of 17.9 years, and with a secondary infection (71.3%). By a multivariate analysis, on-admission laboratory data of thrombocytopenia and hemoglobin showed significant association with DHF. Furthermore, DHF patients had significantly prolonged hospitalizations compared to DF patients. In conclusion, on-admission platelet counts and hemoglobin laboratory data are useful as predictors of DHF especially for suspected dengue patients with the limitations of diagnostic tests.

Background Dual antiplatelet therapy (DAPT) is currently the standard treatment for the prevention of ischemic events after stent implantation. However, the optimal DAPT duration remains elusive for patients with chronic kidney disease (CKD). Therefore, we aimed to compare the effectiveness and safety between long-term and short-term DAPT after coronary stenting in patients with CKD. Methods This retrospective cohort study analyze data from the Taipei Medical University (TMU) Institutional and Clinical Database, which include anonymized electronic health data of 3 million patients that visited TMU Hospital, Wan Fang Hospital, and Shuang Ho Hospital. We enrolled patients with CKD after coronary stenting between 2008 and 2019. The patients were divided into the long-term (>6 months) and short-term DAPT group (≤ 6 months). The primary end point was major adverse cardiovascular events (MACE) from 6 months after the index date. The secondary outcomes were all-cause mortality and Thrombolysis in Myocardial Infarction (TIMI) bleeding. Results A total of 1899 patients were enrolled; of them, 1112 and 787 were assigned to the long-term and short-term DAPT groups, respectively. Long-term DAPT was associated with similar risk of MACE (HR: 1.05, 95% CI: 0.65–1.70, P = 0.83) compare with short-term DAPT. Different CKD risk did not modify the risk of MACE. There was also no significant difference in all-cause mortality (HR: 1.10, 95% CI: 0.75–1.61, P = 0.63) and TIMI bleeding (HR 1.19, 95% CI: 0.86–1.63, P = 0.30) between groups. Conclusions Among patients with CKD and coronary stenting, we found that long-term and short-term DAPT tied on the risk of MACE, all-cause mortality and TIMI bleeding.

Marfan syndrome (MFS) is a rare disease that affects connective tissue, which causes abnormalities in several organ systems including the heart, eyes, bones, and joints. The autosomal dominant disorder was found to be strongly associated with FBN1, TGFBR1, and TGFBR2 mutations. Although multiple genetic mutations have been reported, data from Asian populations are still limited. As a result, we utilized the whole exome sequencing (WES) technique to identify potential pathogenic variants of MFS in a Taiwan cohort. In addition, a variety of annotation databases were applied to identify the biological functions as well as the potential mechanisms of candidate genes. In this study, we confirmed the pathogenicity of FBN1 to MFS. Our results indicated that TTN and POMT1 may be likely related to MFS phenotypes. Furthermore, we found nine unique variants highly shared in a MFS family cohort, of which eight are novel variants worthy of further investigation.

Asthma is a common and heterogeneous disease characterized by chronic airway inflammation. Currently, the two main types of asthma medicines are inhaled corticosteroids and long-acting β2-adrenoceptor agonists (LABAs). In addition, biological drugs provide another therapeutic option, especially for patients with severe asthma. However, these drugs were less effective in preventing severe asthma exacerbation, and other drug options are still limited. Herein, we extracted asthma-associated single nucleotide polymorphisms (SNPs) from the genome-wide association studies (GWAS) and phenome-wide association studies (PheWAS) catalog and prioritized candidate genes through five functional annotations. Genes enriched in more than two categories were defined as “biological asthma risk genes.” Then, DrugBank was used to match target genes with FDA-approved medications and identify candidate drugs for asthma. We discovered 139 biological asthma risk genes and identified 64 drugs targeting 22 of these genes. Seven of them were approved for asthma, including reslizumab, mepolizumab, theophylline, dyphylline, aminophylline, oxtriphylline, and enprofylline. We also found 17 drugs with clinical or preclinical evidence in treating asthma. In addition, eleven of the 40 candidate drugs were further identified as promising asthma therapy. Noteworthy, IL6R is considered a target for asthma drug repurposing based on its high target scores. Through in silico drug repurposing approach, we identified sarilumab and satralizumab as the most promising drug for asthma treatment.

Abstract Background Brain-derived neurotrophic factor (BDNF) is one of the most abundant neurotrophins in the adult brain, and it plays important roles in modulating synaptic plasticity and synaptogenesis. This study attempted to elucidate the role of the BDNF variant rs6265 in emotional symptoms following mild traumatic brain injury (mTBI). Methods To investigate the association between BDNF Val66Met polymorphism (rs6265) and emotional symptoms in mTBI patients, we recruited 192 mTBI patients and evaluated their Beck Anxiety Inventory (BAI) and Beck Depression Inventory (BDI) scores in the first and sixth week after mTBI. Results The patients carrying the T allele of rs6265 had significantly higher BAI scores in the first week following mTBI. In addition, the patients carrying the T allele also showed higher scores of BDI in the first week. In the gender-specific subgroup analysis, the male patients carrying the T allele of rs6265 had higher scores of both BAI and BDI in the first and sixth week. Meanwhile, female patients carrying the T allele also had significantly higher scores of BDI in the first week following mTBI. Conclusions This study provides evidence for the association between the BDNF variant rs6265 and emotional symptoms following mTBI.

This study analyzed the genetic diversity of ciprofloxacin (CIP) nonsusceptibility and the relationship between two major mechanisms and minimum inhibitory concentrations (MICs) of CIP in nontyphoidal Salmonella (NTS). Chromosomal mutations in quinolone resistance-determining regions (QRDRs) and plasmid-mediated quinolone resistance (PMQR) genes were searched from ResFinder, ARG-ANNOT, and PubMed for designing the sequencing regions in gyrA, gyrB, parC, and parE, and the 13 polymerase chain reactions for PMQR genes. We found that QRDR mutations were detected in gyrA (82.1%), parC (59.0%), and parE (20.5%) but not in gyrB among the 39 isolates. Five of the 13 PMQR genes were identified, including oqxA (28.2%), oqxB (28.2%), qnrS (18.0%), aac(6′)-Ib-cr (10.3%), and qnrB (5.1%), which correlated with the MICs of CIP within 0.25–2 μg/mL, and it was found that oxqAB contributed more than qnr genes to increase the MICs. All the isolates contained either QRDR mutations (53.8%), PMQR genes (15.4%), or both (30.8%). QRDR mutations (84.6%) were more commonly detected than PMQR genes (46.2%). QRDR mutation numbers were significantly associated with MICs (p < 0.001). Double mutations in gyrA and parC determined high CIP resistance (MICs ≥ 4 μg/mL). PMQR genes contributed to intermediate to low CIP resistance (MICs 0.25–2 μg/mL), thus providing insights into mechanisms underlying CIP resistance.

Thoracic aortic aneurysm (TAA) formation is a multifactorial process that results in diverse clinical manifestations and drug responses. Identifying the critical factors and their functions in Marfan syndrome (MFS) pathogenesis is important for exploring personalized medicine for MFS. Methylenetetrahydrofolate reductase (MTHFR), methionine synthase (MTR), and methionine synthase reductase (MTRR) polymorphisms have been correlated with TAA severity in MFS patients. However, the detailed relationship between the folate-methionine cycle and MFS pathogenesis remains unclear. Fbn1C1039G/+ mice were reported to be a disease model of MFS. To study the role of the folate-methionine cycle in MFS, Fbn1C1039G/+ mice were treated orally with methionine or vitamin B mixture (VITB), including vitamins B6, B9, and B12, for 20 weeks. VITB reduced the heart rate and circumference of the ascending aorta in Fbn1C1039G/+ mice. Our data showed that the Mtr and Smad4 genes were suppressed in Fbn1C1039G/+ mice, while VITB treatment restored the expression of these genes to normal levels. Additionally, VITB restored canonical transforming-growth factor β (TGF-β) signaling and promoted Loxl1-mediated collagen maturation in aortic media. This study provides a potential method to attenuate the pathogenesis of MFS that may have a synergistic effect with drug treatments for MFS patients.

COVID-19 has become a severe infectious disease and has caused high morbidity and mortality worldwide. Restriction rules such as quarantine and city lockdown have been implemented to mitigate the spread of infection, leading to significant economic impacts. Fortunately, development and inoculation of COVID-19 vaccines are being conducted on an unprecedented scale. The effectiveness of vaccines raises a hope that city lockdown might not be necessary in the presence of ongoing vaccination, thereby minimizing economic loss. The question, however, is how fast and what type of vaccines should be inoculated to control the disease without limiting economic activity. Here, we set up a simulation scenario of COVID-19 outbreak in a modest city with a population of 2.5 million. The basic reproduction number (R0) was ranging from 1.0 to 5.5. Vaccination rates at 1000/day, 10,000/day and 100,000/day with two types of vaccine (effectiveness v = 51% and 89%) were given. The results indicated that R0 was a critical factor. Neither high vaccination rate (10,000 persons/day) nor high-end vaccine (v = 89%) could control the disease when the scenario was at R0 = 5.5. Unless an extremely high vaccination rate was given (>4% of the entire population/per day), no significant difference was found between two types of vaccine. With the population scaled to 25 million, the required vaccination rate was >1,000,000/day, a quite unrealistic number. Nevertheless, with a slight reduction of R0 from 5 to 3.5, a significant impact of vaccine inoculation on disease control was observed. Thus, our study raised the importance of estimating transmission dynamics of COVID-19 in a city before determining the subsequent policy.

Abstract Background Cardiovascular (CV) complications are the main cause of death in end-stage renal disease (ESRD) patients. The high CV risks are attributable to the additive effects of multiple factors. Endothelin (EDN) is a potent vasoconstrictor and plays a role in regulating vascular homeostasis. However, whether variants of the EDN gene are associated with risks of CV events is not known. We conducted a study to investigate associations of variants of the EDN gene with CV events in ESRD patients. Methods A cohort of 190 ESRD patients was recruited, and 19 tagged single-nucleotide polymorphisms within the EDN gene family were selected for genotyping through a TaqMan assay. Data on clinical characteristics and hospitalizations for CV events were collected. Associations of genetic variants of the EDN gene with CV events were analyzed. Results In this cohort, 62% (n = 118) of patients were hospitalized for a CV event. The EDN1 rs4714384 (CC/TC vs. TT) polymorphism was associated with an increased risk of a CV event after multiple testing (p

Infertility is one of the important problems in the modern world. Male infertility is characterized by several clinical manifestations, including low sperm production (oligozoospermia), reduced sperm motility (asthenozoospermia), and abnormal sperm morphology (teratozoospermia). WDR4, known as Wuho, controls fertility in Drosophila. However, it is unclear whether WDR4 is associated with clinical manifestations of male fertility in human. Here, we attempted to determine the physiological functions of WDR4 gene. Two cohorts were applied to address this question. The first cohort was the general population from Taiwan Biobank. Genomic profiles from 68,948 individuals and 87 common physiological traits were applied for phenome-wide association studies (PheWAS). The second cohort comprised patients with male infertility from Wan Fang Hospital, Taipei Medical University. In total, 81 male participants were recruited for the genetic association study. Clinical records including gender, age, total testosterone, follicle-stimulating hormone (FSH), luteinizing hormone (LH), total sperm number, sperm motility, and sperm morphology were collected. In the first cohort, results from PheWAS exhibited no associations between WDR4 genetic variants and 87 common physiological traits. In the second cohort, a total of four tagging single-nucleotide polymorphisms (tSNPs) from WDR4 gene (rs2298666, rs465663, rs2248490, and rs3746939) were selected for genotyping. We found that SNP rs465663 solely associated with asthenozoospermia. Functional annotations through the GTEx portal revealed the correlation between TT or TC genotype and low expression of WDR4. Furthermore, we used mouse embryonic fibroblasts cells from mwdr4 heterozygous (+/‒) mice for functional validation by western blotting. Indeed, low expression of WDR4 contributed to ROS-induced DNA fragmentation. In conclusion, our results suggest a critical role of WDR4 gene variant as well as protein expression in asthenozoospermia.

Immunoglobulin G (IgG) N-glycosylation was discovered to have an association with inflammation status, which has the potential to be a novel biomarker for kidney diseases. In this study, we applied an ultra-high performance liquid chromatography–tandem mass spectrometry (UHPLC-MS/MS) method to plasma and urine samples from 57 individuals with different levels of kidney function. Natural abundances of total IgG, IgG1, IgG2, and IgG3 subclasses in plasma showed positive correlations to the estimated glomerular filtration rates (eGFRs). Eighteen IgG glycopeptides also showed positive correlations. In contrast, higher IgG amounts were found in urine samples from participants with lower eGFR values. After normalizing IgG glycopeptides from plasma to their respective protein amounts, H4N4F1S1-IgG1 (r = 0.37, p = 0.0047, significant) and H5N4F1S1-IgG1 (r = 0.25, p = 0.063, marginally significant) were the two glycopeptides that still had positive correlations with eGFRs. The results showed that the UHPLC-MS/MS method is capable of investigating IgG profiles, and monitoring IgG and glycosylation patterns is worthy of further clinical application for kidney disease.

Abstract Background Anemia is common among end-stage renal disease (ESRD) patients who undergone hemodialysis. The total reduction of red blood cell (RBC) count is associated with poor prognosis in these patients. Although erythropoietin (EPO) has been used as an effective treatment for ESRD patients with anemia, a large number of patients still present poor responses to EPO treatment. Methods We measured T-cell receptor sequencing profiles, including length of complementarity-deteremining region 3 (CDR3), intra- and inter-group (EPO resistant vs. responsive) clonotype diversity, V(D)J usage profiles and V-J combinations from ESRD patients and to investigate the correlation between these features and EPO treatment efficacy. Results Our results revealed statistical significance in the top 3 ~ 15 most abundant joint distributions of Vβ/Jβ among the two groups, suggesting the importance of V or J gene utilization in the EPO response of ESRD patients. Conclusions In summary, we provided evidence addressing the potential correlation between the immune repertoire and EPO response in ESRD patients. Trial registration TMU-JIRB 201309026. Registered 16 October 2013.

Abstract The insulin-like growth factor (IGF)-1 signaling is relevant in regulating cell growth and cytokine secretions by glioblastomas. MicroRNAs determine the cell fate in glioblastomas. However, relationships between IGF-1 signaling and miRNAs in glioblastoma pathogenesis are still unclear. Our aim was to validate the IGF-1-mediated mRNA/miRNA regulatory network in glioblastomas. Using in silico analyses of mRNA array and RNA sequencing data from The Cancer Genome Atlas (TCGA), we identified 32 core enrichment genes that were highly associated with IGF-1-promoted cytokine-cytokine receptor interactions. To investigate the IGF-1-downregulated miRNA signature, microarray-based approaches with IGF-1-treated U87-MG cells and array data in TCGA were used. Four miRNAs, including microRNA (miR)-9-5p, miR-9-3p, miR-181d, and miR-130b, exhibited an inverse correlation with IGF-1 levels. The miR-181d, that targeted the most IGF-1-related cytokine genes, was significantly reduced in IGF-1-treated glioma cells. Statistical models incorporating both high-IGF-1 and low-miR-181d statuses better predicted poor patient survival, and can be used as an independent prognostic factor in glioblastomas. The C-C chemokine receptor type 1 (CCR1) and interleukin (IL)-1b demonstrated inverse correlations with miR-181d levels and associations with patient survival. miR-181d significantly attenuated IGF-1-upregulated CCR1 and IL-1b gene expressions. These findings demonstrate a distinct role for IGF-1 signaling in glioma progression via miR-181d/cytokine networks.

Background Previous studies show that mTOR inhibitors decrease the risk of cancer development after kidney transplantation. However, the effect of cumulative doses of mTOR inhibitors on cancer after kidney transplantation is not well known. Methods In the current study, patients were registered into a national database in Taiwan. Between year 2000 and 2013, 4,563 patients received kidney transplantation. They were divided into two groups, according to mTOR inhibitors usage. The cumulative dose of mTOR inhibitors was recorded. Patients were followed-up until de novo cancer development, death, or the end of 2014. Results Patients were divided into two groups: mTOR inhibitors users (study group, n = 828) and mTOR inhibitors non-users (control group, n = 3,735). The median follow-up duration was 7.8 years. The risk of de novo cancer (hazards ratio (HR) 0.80, 95% CI [0.60–1.09], p = 0.16) and risk of death (HR 1.14, 95% CI [0.82–1.60], p = 0.43) was not different between mTOR inhibitor user and non-user groups. Neither high- nor low-dose exposure to mTOR inhibitors was associated with increased risk of cancer or mortality. Analysis of cancer subtypes showed no influence by mTOR inhibitors. In addition, the cause of mortality was not significantly different between the two groups. Discussion We could not find the association of mTOR inhibitors use and risk of de novo cancer development or mortality in patients with kidney transplantation in Chinese patients. Cumulative exposure to mTOR inhibitors did not change the results.

Intricate relationships among cytokines (including chemokines) shape the tumor microenvironment (TME) and reflect cell-cell interactions between malignant cells and other cells from the TME. Although our previous study indicated the transcriptional landscape of cytokines in 19 cancer types, the global pattern somatic copy number (SCN) alterations and the clinical relevance of cytokines have not been systematically investigated. Here, we reported a significant negative selection on cytokine genes. We also linked the SCN losses of cytokine genes to the abundance of immune infiltrates which affects cancer progression and patient prognoses. We also demonstrated and validated the correlations between SCN alterations of cytokine-containing loci and drug sensitivity. The results indicated the genomic loss of cytokines in malignant cells as a crucial theme for interrogating cancer progression, malignant cell-TME interactions, and therapeutics.

Among all cancers in women, breast cancer has the highest incidence. The mortality of breast cancer is highly associated with metastasis. Migration and malignant transformation of cancer cells have been reported to be modulated by store-operated calcium (SOC) channels, which control calcium signaling and cell proliferation pathways. Stromal interaction molecule 1 (STIM1) is a calcium sensor in the endoplasmic reticulum, triggering the activation of store-operated calcium signaling. However, the clinical relevance of STIM1 in breast cancer is still unclear. Here, we recruited 348 breast cancer patients and conducted a genetic association study to address this question. Four tagging germline single nucleotide variants (SNVs) in STIM1 were selected and RNA sequencing data of 525 breast cancer samples from The Cancer Genome Atlas (TCGA) database were evaluated. The results show that rs2304891 and rs3750996 were correlated with clinical stage of breast cancer. Expression quantitative trait loci (eQTL) analysis indicated that risk G allele of STIM1 contributed to the higher expression of STIM1. In addition, we found an increased risk of rs2304891 G allele and rs3750996 A allele in estrogen receptor (ER) positive and progesterone receptor (PR) positive patients. In conclusion, our results suggest that germline SNV, rs2304891 and rs3750996 as well as STIM1 expression are important biomarkers for the prediction of clinical outcomes in breast cancer patients.

Hepatocellular carcinoma (HCC) often develops from chronic hepatitis B (CHB) through replication of hepatitis B virus (HBV) infection. Calcium (Ca2+) signaling plays an essential role in HBV replication. Store-operated calcium (SOC) channels are a major pathway of Ca2+ entry into non-excitable cells such as immune cells and cancer cells. The basic components of SOC signaling include the STIM1 and ORAI1 genes. However, the roles of STIM1 and ORAI1 in HBV-mediated HCC are still unclear. Thus, long-term follow-up of HBV cohort was carried out in this study. This study recruited 3631 patients with chronic hepatitis (345 patients with HCC, 3286 patients without HCC) in a Taiwanese population. Genetic variants of the STIM1 and ORAI1 genes were detected using an Axiom CHB1 genome-wide array. Clinical associations of 40 polymorphisms were analyzed. Three of the STIM1 single-nucleotide polymorphisms (SNPs) (rs6578418, rs7116520, and rs11030472) and one SNP of ORAI1 (rs6486795) showed a trend of being associated with HCC disease (p < 0.05). However, after correction for multiple testing, none of the SNPs reached a significant level (q > 0.05); in contrast, neither STIM1 nor ORAI1 showed a significant association with HCC progression in CHB patients. Functional studies by both total internal reflection fluorescence images and transwell migration assay indicated the critical roles of SOC-mediated signaling in HCC migration. In conclusion, we reported a weak correlation between STIM1/ORAI1 polymorphisms and the risk of HCC progression in CHB patients.

Diabetic kidney disease (DKD) is the leading cause of morbidity and mortality in patients with diabetes mellitus (DM) and the most common variant of end-stage renal disease (ESRD) globally. The economic burden of ESRD treatment with dialysis is substantial. The incidence and prevalence of ESRD in Taiwan remain the highest worldwide. Therefore, identifying genetic factors affecting kidney function would have valuable clinical implications. We performed microarray experiments and identified that ubiquitin protein ligase E3C (UBE3C) is differentially expressed in two DKD patient groups with extreme (low and high) urine protein-to-creatinine ratios. A follow-up genotyping study was performed in a larger group to investigate any specific variants of UBE3C associated with DKD. A total of 263 patients were included in the study, comprising 172 patients with DKD and 91 control subjects (patients with DM without chronic kidney disease (CKD)). Two UBE3C variants (rs3802129(AA) and rs7807(CC)) were determined to be associated with reduced kidney function. The haplotype analysis revealed that rs3802129/rs3815217 (block 1) with A/G haplotype and rs8101/rs7807 (block 2) with T/C haplotype were associated with higher risks of CKD phenotypes. These findings suggest a clinical role of UBE3C variants in DKD risk.

Renal dysfunction is common in patients with diabetes mellitus (DM). Previous findings from a meta-analysis of GWAS indicated that the variation of RAB38/CTSC is highly associated with the urinary albumin-to-creatinine ratio (UACR) in European populations. In addition, RAB38 knockout rats showed an increase in urinary albumins. Although the prevalence of chronic kidney disease is high in Taiwan, the role of genetic variants in diabetic renal function is still unclear. In the current study, 275 diabetic nephropathy (DN) patients were recruited to perform a genetic association study. Our results indicated that rs1027027, rs302647, and rs302646 in RAB38 were significantly associated with urinary protein-to-creatinine ratio (UPCR) levels in DN patients. Importantly, after analysis stratified by gender, a significant genetic influence on UPCR levels was observed in the male population. The findings confirmed the roles of gender and variants of RAB38 in the risk of UPCR in Diabetic Nephropathy patients.

The effect of comorbidities and the immune profiles of the kidney cancer microenvironment play a major role in patients’ prognosis and survival. Using the National Health Insurance Research Database (Taiwan), we identified patients aged >20 years with a first diagnosis of kidney cancer between 2005 and 2014. Differences in demographic characteristics and comorbidities were examined using the Pearson chi-squared test or the t test. The Cox regression model was used to construct the nomogram. RNA-seq data were applied from The Cancer Genome Atlas database, and correlations between immune metagenes and clinical characteristics were determined using a linear regression model. In this nationwide cohort study, including 5090 patients with kidney cancer, predictors in our prediction models included age, sex, chronic kidney disease, dialysis requirements, renal stones, cerebrovascular disease, and metastasis tumor. In the tumor tissue profiles, significant positive correlations between immune metagenes and clinical stage or overall survival were observed among Natural Killer (NK) cells (CD56−), CD4+ T-helper 2 (Th2) cells, and activated Dendritic Cell (aDC). A negative correlation was observed between expression level of Dendritic Cell (DC) and overall survival. Patients with kidney cancer exhibit high prevalence of comorbid disease, especially in older patients. Comorbid disease types exert unique effects, and a particular comorbidity can affect cancer mortality. Moreover, the expression of immune metagenes can be utilized as potentialbiomarkers especially for further study of molecular mechanisms as well as microenvironments in kidney cancer.

In the original publication of this article [1] the name of the fifth author is uncorrect. The correct name of the fifth author should be Wei-Chiao Chang rather than Wei-Chao Chang. The original publication has been corrected.