SpliceAPP: an interactive web server to predict splicing errors arising from human mutations.

BACKGROUND: Splicing variants are a major class of pathogenic mutations, with their severity equivalent to nonsense mutations. However, redundant and degenerate splicing signals hinder functional assessments of sequence variations within introns, particularly at branch sites. We have established a massively parallel splicing assay to assess the impact on splicing of 11,191 disease-relevant variants. Based on the experimental results, we then applied regression-based methods to identify factors determining splicing decisions and their respective weights. RESULTS: Our statistical modeling is highly sensitive, accurately annotating the splicing defects of near-exon intronic variants, outperforming state-of-the-art predictive tools. We have incorporated the algorithm and branchpoint information into a web-based tool, SpliceAPP, to provide an interactive application. This user-friendly website allows users to upload any genetic variants with genome coordinates (e.g., chr15 74,687,208 A G), and the tool will output predictions for splicing error scores and evaluate the impact on nearby splice sites. Additionally, users can query branch site information within the region of interest. CONCLUSIONS: In summary, SpliceAPP represents a pioneering approach to screening pathogenic intronic variants, contributing to the development of precision medicine. It also facilitates the annotation of splicing motifs. SpliceAPP is freely accessible using the link https://bc.imb.sinica.edu.tw/SpliceAPP . Source code can be downloaded at https://github.com/hsinnan75/SpliceAPP .

Immunometabolic features of natural killer cells are associated with infection outcomes in critical illness.

Immunosuppression increases the risk of nosocomial infection in patients with chronic critical illness. This exploratory study aimed to determine the immunometabolic signature associated with nosocomial infection during chronic critical illness. We prospectively recruited patients who were admitted to the respiratory care center and who had received mechanical ventilator support for more than 10 days in the intensive care unit. The study subjects were followed for the occurrence of nosocomial infection until 6 weeks after admission, hospital discharge, or death. The cytokine levels in the plasma samples were measured. Single-cell immunometabolic regulome profiling by mass cytometry, which analyzed 16 metabolic regulators in 21 immune subsets, was performed to identify immunometabolic features associated with the risk of nosocomial infection. During the study period, 37 patients were enrolled, and 16 patients (43.2%) developed nosocomial infection. Unsupervised immunologic clustering using multidimensional scaling and logistic regression analyses revealed that expression of nuclear respiratory factor 1 (NRF1) and carnitine palmitoyltransferase 1a (CPT1a), key regulators of mitochondrial biogenesis and fatty acid transport, respectively, in natural killer (NK) cells was significantly associated with nosocomial infection. Downregulated NRF1 and upregulated CPT1a were found in all subsets of NK cells from patients who developed a nosocomial infection. The risk of nosocomial infection is significantly correlated with the predictive score developed by selecting NK cell-specific features using an elastic net algorithm. Findings were further examined in an independent cohort of COVID-19-infected patients, and the results confirm that COVID-19-related mortality is significantly associated with mitochondria biogenesis and fatty acid oxidation pathways in NK cells. In conclusion, this study uncovers that NK cell-specific immunometabolic features are significantly associated with the occurrence and fatal outcomes of infection in critically ill population, and provides mechanistic insights into NK cell-specific immunity against microbial invasion in critical illness.

Upregulated PARP1 confers breast cancer resistance to CDK4/6 inhibitors via YB-1 phosphorylation.

BACKGROUND: Cyclic-dependent kinase (CDK) 4/6 kinases, as the critical drivers of the cell cycle, are involved in the tumor progression of various malignancies. Pharmacologic inhibitors of CDK4/6 have shown significant clinical prospects in treating hormone receptor-positive and human epidermal growth factor receptor-negative (HR + /HER2-) breast cancer (BC) patients. However, acquired resistance to CDK4/6 inhibitors (CDK4/6i), as a common issue, has developed rapidly. It is of great significance that the identification of novel therapeutic targets facilitates overcoming the CDK4/6i resistance. PARP1, an amplified gene for CDK4/6i-resistant patients, was found to be significantly upregulated during the construction of CDK4/6i-resistant strains. Whether PARP1 drives CDK4/6i resistance in breast cancer is worth further study. METHOD: PARP1 and p-YB-1 protein levels in breast cancer cells and tissues were quantified using Western blot (WB) analysis, immunohistochemical staining (IHC) and immunofluorescence (IF) assays. Bioinformatics analyses of Gene Expression Profiling Interactive Analysis (GEPIA), Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Cell Line Encyclopedia (CCLE) datasets were applied to explore the relationship between YB-1/PARP1 protein levels and CDK4/6i IC50. Cell Counting Kit-8 (CCK-8) and crystal violet staining assays were performed to evaluate cell proliferation rates and drug killing effects. Flow cytometry assays were conducted to assess apoptosis rates and the G1/S ratio in the cell cycle. An EdU proliferation assay was used to detect the DNA replication ratio after treatment with PARP1 and YB-1 inhibitors. A ChIP assay was performed to assess the interaction of the transcription factor YB-1 and associated DNA regions. A double fluorescein reporter gene assay was designed to assess the influence of WT/S102A/S102E YB-1 on the promoter region of PARP1. Subcutaneous implantation models were applied for in vivo tumor growth evaluations. RESULTS: Here, we reported that PARP1 was amplified in breast cancer cells and CDK4/6i-resistant patients, and knockdown or inhibition of PARP1 reversed drug resistance in cell experiments and animal models. In addition, upregulation of transcription factor YB-1 also occurred in CDK4/6i-resistant breast cancer, and YB-1 inhibition can regulate PARP1 expression. p-YB-1 and PARP1 were upregulated when treated with CDK4/6i based on the WB and IF results, and elevated PARP1 and p-YB-1 were almost simultaneously observed during the construction of MCF7AR-resistant strains. Inhibition of YB-1 or PAPR1 can cause decreased DNA replication, G1/S cycle arrest, and increased apoptosis. We initially confirmed that YB-1 can bind to the promoter region of PARP1 through a ChIP assay. Furthermore, we found that YB-1 phosphorylated at S102 was crucial for PARP1 transcription according to the double fluorescein reporter gene assay. The combination therapy of YB-1 inhibitors and CDK4/6i exerted a synergistic antitumor effect in vitro and in vivo. The clinical data suggested that HR + /HER2- patients with low expression of p-YB-1/PARP1 may be sensitive to CDK4/6i in breast cancer. CONCLUSION: These findings indicated that a ''YB-1/PARP1'' loop conferred resistance to CDK4/6 inhibitors. Furthermore, interrupting the loop can enhance tumor killing in the xenograft tumor model, which provides a promising strategy against drug resistance in breast cancer.

Metabolic regulation of Shewanella oneidensis for microbial electrosynthesis: From extracellular to intracellular.

Shewanella oneidensis MR-1 (S. oneidensis MR-1) has been shown to benefit from microbial electrosynthesis (MES) due to its exceptional electron transfer efficiency. In this study, genes involved in both extracellular electron uptake (EEU) and intracellular CO2 conversion processes were examined and regulated to enhance MES performance. The key genes identified for MES in the EEU process were mtrB, mtrC, mtrD, mtrE, omcA and cctA. Overexpression of these genes resulted in 1.5-2.1 times higher formate productivity than that of the wild-type strains (0.63 mmol/(L·µg protein)), as 0.94-1.61 mmol/(L·µg protein). In the intracellular CO2 conversion process, overexpression of the nadE, nadD, nadR, nadV, pncC and petC genes increased formate productivity 1.3-fold-3.4-fold. Moreover, overexpression of the formate dehydrogenase genes fdhA1, fdhB1 and fdhX1 in modified strains led to a 2.3-fold-3.1-fold increase in formate productivity compared to wild-type strains. The co-overexpression of cctA, fdhA1 and nadV in the mutant strain resulted in 5.59 times (3.50 mmol/(L·µg protein)) higher formate productivity than that of the wild-type strains. These findings revealed that electrons of MES derived from the electrode were utilized in the energy module for synthesizing ATP and NADH, followed by the synthesis of formate in formate dehydrogenase by the combinatorial effects of ATP, NADH, electrons and CO2. The results provide new insights into the mechanism of MES in S. oneidensis MR-1 and pave the way for genetic improvements that could facilitate the further application of MES.

Hypertension as a Novel Link for Shared Heritability in Age at Menarche and Cardiometabolic Traits.

CONTEXT: Extremely early age at menarche, also called precocious puberty, has been associated with various cardiometabolic traits, but their shared heritability remains unclear. OBJECTIVES: This work aimed to identify new shared genetic variants and their pathways for age at menarche and cardiometabolic traits and to investigate the influence of central precocious puberty on childhood cardiometabolic traits. METHODS: Using the conjunction false discovery rate method, this study analyzed genome-wide association study data from the menarche-cardiometabolic traits among 59 655 females of Taiwanese ancestry and systemically investigated pleiotropy between age at menarche and cardiometabolic traits. To support the novel hypertension link, we used the Taiwan Puberty Longitudinal Study (TPLS) to investigate the influence of precocious puberty on childhood cardiometabolic traits. RESULTS: We discovered 27 novel loci, with an overlap between age at menarche and cardiometabolic traits, including body fat and blood pressure. Among the novel genes discovered, SEC16B, CSK, CYP1A1, FTO, and USB1 are within a protein interaction network with known cardiometabolic genes, including traits for obesity and hypertension. These loci were confirmed through demonstration of significant changes in the methylation or expression levels of neighboring genes. Moreover, the TPLS provided evidence regarding a 2-fold higher risk of early-onset hypertension that occurred in girls with central precocious puberty. CONCLUSION: Our study highlights the usefulness of cross-trait analyses for identifying shared etiology between age at menarche and cardiometabolic traits, especially early-onset hypertension. The menarche-related loci may contribute to early-onset hypertension through endocrinological pathways.

Continuous fixation of dissolved arsenite from flooded soil by cooperating ferrihydrite with Geobacter sulfurreducens.

The fixed arsenic in soil is easy to be released into the aquatic environment in the form of arsenite (As(III)) with high toxicity and mobility due to the eutrophication of environment under anaerobic conditions. However, As(III) is difficult to be fixed in situ continuously by traditional methods, especially for the most efficient fixation method by iron ores. Based on that Fe(II) could promote the fixation of As(III), this study investigated the possibility that Geobacter sulfurreducens (G. sulfurreducens) cooperates with ferrihydrite to fix released As(III) from flooded soil in a glass column continuously under anaerobic conditions. During 42 days of operation of reactors that simulated the actual flooded soil environment, the concentration of released As(III) in the reactor with adding G. sulfurreducens and ferrihydrite is always lower than that in reactors with adding ferrihydrite or no treatment. Compared with reactors without treatment, the accumulated content of released As(III) (2455.0 ± 313.1 µg) decreased by 39.4% in the reactor with adding G. sulfurreducens and ferrihydrite on the last day, while that in reactors with adding ferrihydrite only decreased by 11.6%, respectively. These were caused by the cooperation of G. sulfurreducens and ferrihydrite, which increased the relative abundance of iron-reducing microorganisms to inhibit metabolisms of As-reducing microorganisms, inhibited the quick release of As(III) from solid soil, and promoted the release of iron to accelerate the formation of stable secondary ores with As. This study could provide an environmentally friendly method to fix dissolved As(III) pollutants from soil continuously.

Mechanism and modeling of human disease-associated near-exon intronic variants that perturb RNA splicing.

It is estimated that 10%-30% of disease-associated genetic variants affect splicing. Splicing variants may generate deleteriously altered gene product and are potential therapeutic targets. However, systematic diagnosis or prediction of splicing variants is yet to be established, especially for the near-exon intronic splice region. The major challenge lies in the redundant and ill-defined branch sites and other splicing motifs therein. Here, we carried out unbiased massively parallel splicing assays on 5,307 disease-associated variants that overlapped with branch sites and collected 5,884 variants across the 5' splice region. We found that strong splice sites and exonic features preserve splicing from intronic sequence variation. Whereas the splice-altering mechanism of the 3' intronic variants is complex, that of the 5' is mainly splice-site destruction. Statistical learning combined with these molecular features allows precise prediction of altered splicing from an intronic variant. This statistical model provides the identity and ranking of biological features that determine splicing, which serves as transferable knowledge and out-performs the benchmarking predictive tool. Moreover, we demonstrated that intronic splicing variants may associate with disease risks in the human population. Our study elucidates the mechanism of splicing response of intronic variants, which classify disease-associated splicing variants for the promise of precision medicine.

Effectiveness of entecavir vs tenofovir disoproxil fumarate for functional cure of chronic hepatitis B in an international cohort.

INTRODUCTION: Both entecavir (ETV) and tenofovir disoproxil fumarate (TDF) are first-line therapies for chronic hepatitis B (CHB), but their comparative effectiveness with regards to hepatitis B surface antigen (HBsAg) seroclearance remains unclear. METHODS: This international multicenter cohort study enrolled 7697 treatment-naïve CHB patients (median age 50 years; male 66.75%) initiated on either ETV (n = 5430) or TDF (n = 2267) without baseline malignancy or immunosuppression from 23 centers across 10 countries or regions. Patients were observed for HBsAg seroclearance until death, loss to follow-up, or treatment discontinuation or switching. The incidences of HBsAg seroclearance were adjusted for competing mortality and compared between ETV and TDF cohorts with inverse probability of treatment weighting (IPTW) and also by multivariable regression analysis. RESULTS: The study population was followed up for a median duration of 56.1 months with 36,929 11 person-years of observation. HBsAg seroclearance occurred in 70 ETV-treated and 21 TDF-treated patients, yielding 8-year cumulative incidence of 1.69% (95% confidence interval [CI] 1.32-2.17) for ETV and 1.34% (95% CI 0.85-2.10%), for TDF (p = 0.58). In the IPTW analysis with the two study cohorts more balanced in background covariates, the age-adjusted hazard ratio (HR) of TDF versus ETV for HBsAg seroclearance was 0.91 (95% CI 0.50-1.64; p = 0.75). Furthermore, there was no significant difference between the two medications in the multivariable competing risk regression model (adjusted sub-distributional HR 0.92 for TDF vs. ETV; 95% CI 0.56-1.53; p = 0.76). CONCLUSIONS: ETV and TDF did not differ significantly in the incidence of HBsAg seroclearance, which rarely occurred with either regimen.

Estimating the causal effect of embryo transfer day on clinical in vitro fertilization outcomes using propensity score matching.

BACKGROUND: For women undergoing in vitro fertilization (IVF), the clinical benefit of embryo transfer at the blastocyst stage (Day 5) versus cleavage stage (Day 3) remains controversial. The purpose of this study is to compare the implantation rate, clinical pregnancy rate and odds of live birth of Day 3 and Day 5 embryo transfer, and more importantly, to address the issue that patients were chosen to receive either transfer protocol due to their underlying clinical characteristics, i.e., confounding by indication. METHODS: We conducted a retrospective cohort study of 9,090 IVF cycles collected by Lee Women's Hospital in Taichung, Taiwan from 1998 to 2014. We utilized the method of propensity score matching to mimic a randomized controlled trial (RCT) where each patient with Day 5 transfer was matched by another patient with Day 3 transfer with respect to other clinical characteristics. Implantation rate, clinical pregnancy rate, and odds of live birth were compared for women underwent Day 5 transfer and Day 3 transfer to estimate the causal effects. We further investigated the causal effects in subgroups by stratifying age and anti-Mullerian hormone (AMH). RESULTS: Our analyses uncovered an evidence of a significant difference in implantation rate (p=0.04) favoring Day 5 transfer, and showed that Day 3 and Day 5 transfers made no difference in both odds of live birth (p=0.27) and clinical pregnancy rate (p=0.11). With the increase of gestational age, the trend toward non-significance of embryo transfer day in our result appeared to be consistent for subgroups stratified by age and AMH, while all analyses stratified by age and AMH were not statistically significant. CONCLUSIONS: We conclude that for women without strong indications for Day 3 or Day 5 transfer, there is a small significant difference in implantation rate in favor of Day 5 transfer. However, the two protocols have indistinguishable outcomes on odds of live birth and clinical pregnancy rate.

Age effect on in vitro fertilization pregnancy mediated by anti-Mullerian hormone (AMH) and modified by follicle stimulating hormone (FSH).

BACKGROUND: Both follicle stimulating hormone (FSH) and anti-mullerian hormone (AMH) are widely used to assess the ovarian reserve in women for in vitro fertilization (IVF). However, studies also showed that both AMH and FSH are significantly associated with age: as age increases, AMH decreases and FSH increases. This study aims to investigate the mechanism of age effect on IVF live birth rate, particularly through mediation and interaction by AMH and FSH. METHODS: We conducted a retrospective cohort study of 13970 IVF cycles collected by eIVF from 2010 to 2016. A series of logistic mixed models were used to estimate the association of live birth and AMH (or FSH). The mediation effects and proportion mediated, were quantified by our previously proposed mediation analyses. We further investigated the FSH-modified mediation effects on live birth rate through AMH, accounting for the nonlinear age effect. RESULTS: Our analyses showed that age effect on live birth was mediated more by AMH than by FSH (18 vs. 6%). The mediation effect through AMH can be further modified by FSH level regardless of women's age. CONCLUSIONS: In summary, mediation model provides a new perspective elucidating the mechanism of IVF successful rate by age. The majority of the age effect on live birth rate remained unexplained by AMH and FSH, suggesting its importance and independent role in IVF.

Cigarette smoking increases the risk of nasopharyngeal carcinoma through the elevated level of IgA antibody against Epstein-Barr virus capsid antigen: A mediation analysis.

BACKGROUND: The study aims are to evaluate the associations between nasopharyngeal carcinoma (NPC) risk and cigarette smoking and to explore the effects of cigarette smoking on Epstein-Barr virus (EBV) infection for NPC risk. METHODS: 1235 male NPC cases and 1262 hospital-based male controls matched to cases were recruited across six collaborative hospitals between 2010 and 2014. Using a standardized questionnaire, information on cigarette smoking and other potential risk factors for NPC was obtained. Blood was collected and used for anti-EBV VCA IgA and anti-EBV EA-EBNA1 IgA testing using standard methods. Unconditional logistic regression analysis was used to estimate odds ratio (OR) with 95% confidence interval (CI) for each risk factor after adjusting for confounders. RESULTS: 63.6% of cases and 44.0% of controls reported ever smoking cigarettes. After full adjustment, current smokers had a significant 1.60-fold (95% CI = 1.30-1.97) and former smokers a borderline significant 1.27-fold (95% CI = 1.00-1.60) increased NPC risk compared to never smokers. NPC risk increased with increasing duration, intensity, and pack-years of cigarette smoking but not with age at smoking initiation. Among controls, anti-EBV VCA IgA seropositivity rate was higher in current smokers than never smokers (14.0% vs 8.4%; OR = 1.82; 95% CI = 1.19-2.79). Mediation analyses showed that more than 90% of the cigarette smoking effect on NPC risk is mediated through anti-EBV VCA IgA. CONCLUSION: This study confirms the association between long-term cigarette smoking and NPC and demonstrates that current smoking is associated with seropositivity of anti-EBV VCA IgA antibodies.

Cellular Uptake of Gold Nanoparticles and Their Movement in 3D Multicellular Tumor Spheroids: Effect of Molecular Weight and Grafting Density of Poly(2-hydroxyl ethyl acrylate).

It is known that the size of gold nanoparticles (GNPs) is not the only determining factor in the uptake by cells such as cancer cells. The surface functionalization plays a crucial role, in particular the nature of the ligand as well as the molecular weight and the grafting density. Here, poly(2-hydroxy ethyl) acrylate (pHEA) with molecular weights ranging from 10, 20 to 39 g mol-1 via reversible addition-fragmentation chain transfer polymerization is synthesized. These polymers are used directly to coat GNPs with sizes of 20, 40, and 70 nm as the trithiocarbonate functionality can strongly bind to the gold surface. The library of nine GNP is found to be nontoxic against lung carcinoma cells A549 and has negligible albumin protein absorption as determined by quartz crystal microbalance. Laser scanning confocal microscopy and flow cytometry reveal that GNP coated with medium length pHEA displays the highest cellular uptake while the effect of the size is not statistically significant. In contrast, multicellular tumor spheroids, which is a 3D model that simulates the tissue, enable the penetration of GNP coated with the longest pHEA chain while it also appears that smaller GNPs have now a clear advantage.