A model-based clustering via mixture of hierarchical models with covariate adjustment for detecting differentially expressed genes from paired design.

The causes of many complex human diseases are still largely unknown. Genetics plays an important role in uncovering the molecular mechanisms of complex human diseases. A key step to characterize the genetics of a complex human disease is to unbiasedly identify disease-associated gene transcripts on a whole-genome scale. Confounding factors could cause false positives. Paired design, such as measuring gene expression before and after treatment for the same subject, can reduce the effect of known confounding factors. However, not all known confounding factors can be controlled in a paired/match design. Model-based clustering, such as mixtures of hierarchical models, has been proposed to detect gene transcripts differentially expressed between paired samples. To the best of our knowledge, no model-based gene clustering methods have the capacity to adjust for the effects of covariates yet. In this article, we proposed a novel mixture of hierarchical models with covariate adjustment in identifying differentially expressed transcripts using high-throughput whole-genome data from paired design. Both simulation study and real data analysis show the good performance of the proposed method.

Increasing Enzyme Mannose-6-Phosphate Levels but Not Miglustat Coadministration Enhances the Efficacy of Enzyme Replacement Therapy in Pompe Mice.

Pompe disease is a rare glycogen storage disorder caused by a deficiency in the lysosomal enzyme acid α-glucosidase, which leads to muscle weakness, cardiac and respiratory failure, and early mortality. Alglucosidase alfa, a recombinant human acid α-glucosidase, was the first approved treatment of Pompe disease, but its uptake into skeletal muscle via the cation-independent mannose-6-phosphate (M6P) receptor (CIMPR) is limited. Avalglucosidase alfa has received marketing authorization in several countries for infantile-onset and/or late-onset Pompe disease. This recently approved enzyme replacement therapy (ERT) was glycoengineered to maximize CIMPR binding through high-affinity interactions with ∼7 bis-M6P moieties. Recently, small molecules like the glucosylceramide synthase inhibitor miglustat were reported to increase the stability of recombinant human acid α-glucosidase, and it was suggested that an increased serum half-life would result in better glycogen clearance. Here, the effects of miglustat on alglucosidase alfa and avalglucosidase alfa stability, activity, and efficacy in Pompe mice were evaluated. Although miglustat increased the stability of both enzymes in fluorescent protein thermal shift assays and when incubated in neutral pH buffer over time, it reduced their enzymatic activity by ∼50%. Improvement in tissue glycogen clearance and transcriptional dysregulation in Pompe mice correlated with M6P levels but not with miglustat coadministration. These results further substantiate the crucial role of CIMPR binding in lysosomal targeting of ERTs. SIGNIFICANCE STATEMENT: This work describes important new insights into the treatment of Pompe disease using currently approved enzyme replacement therapies (ERTs) coadministered with miglustat. Although miglustat increased the stability of ERTs in vitro, there was no positive impact to glycogen clearance and transcriptional correction in Pompe mice. However, increasing mannose-6-phosphate levels resulted in increased cell uptake in vitro and increased glycogen clearance and transcriptional correction in Pompe mice, further underscoring the crucial role of cation-independent mannose-6-phosphate receptor-mediated lysosomal targeting for ERTs.

An efficient mutagenesis system to improve the propamocarb tolerance in Lecanicillium lecanii (Zimmermann) Zare & Gams.

Lecanicillium lecanii (Zimmermann) Zare & Gams is used as an effective biopesticide for the control of sap-sucking insect pests on agricultural crops. However, low fungicide tolerance limits its large-scale field application. To improve the propamocarb tolerance in L. lecanii, a composite mutagenesis system was established by using UV-light (U), N-Methyl-N'-nitro-N-nitrosoguanidine (NTG) (N) and N+ ion-beam (I). The permutation type of three agents was a consecutive mutagenesis treatment (I/N/U) after an intermittent treatment (U + N + I). The "U" mutagenesis was performed at 254 nm for 60 s and at a distance of 45 cm under a 20 W germicidal lamp, the "N" mutagenesis was performed at a concentration of 1.0 mg/mL NTG for 60 min, and the "I" mutagenesis was performed by low energy N+ ion-beam using a dose of 10 × 1013 ions/cm2 at 30 keV. This composite mutagenesis system was recorded as the "U + N + I + I/N/U," and then the mutagenesis efficiency in improving propamocarb tolerance was assessed by analyzing changes of mutants in the propamocarb sensitivity, mitotic stability, mycelial growth speed on plates or in liquid, sporulation on plates or aphids, conidial germination, 50% lethal concentration (LC50) and 50% lethal time (LT50) to aphids, lipid constituent and cell membrane permeability and control against aphids in the presence or absence of propamocarb. Compared to the wild-type isolate with a 50% effective concentration (EC50) value of 503.6 µg/mL propamocarb, the Ll-IC-UNI produced by the "U + N + I + I/N/U" had the highest EC50 value of 3576.4 µg/mL and a tolerance ratio of 7.1. The mutant was mitotically stable in 20-passage cultivation and did not show any unfavorable changes in growth and virulence indicators. The mutant showed the highest ability to resist or avoid the damaging effects of propamocarb as reflected by the alternations of lipid constituents and membrane permeability. The interval time for applying fungal agent was significantly shortened in this mutant after spraying a field recommended dose of 550 µg/mL propamocarb. In conclude, the "U + N + I + I/N/U" composite mutagenesis mode was efficient and useful to improve the propamocarb-tolerance of L. lecanii and the obtained Ll-IC-UNI could have commercial potential for field application.

Evaluation of Neurofilament Light Chain as a Biomarker of Neurodegeneration in X-Linked Childhood Cerebral Adrenoleukodystrophy.

Cerebral adrenoleukodystrophy (CALD) is a devastating, demyelinating neuroinflammatory manifestation found in up to 40% of young males with an inherited mutation in ABCD1, the causative gene in adrenoleukodystrophy. The search for biomarkers which correlate to CALD disease burden and respond to intervention has long been sought after. We used the Olink Proximity Extension Assay (Uppsala, Sweden) to explore the cerebral spinal fluid (CSF) of young males with CALD followed by correlative analysis with plasma. Using the Target 96 Neuro Exploratory panel, we found that, of the five proteins significantly increased in CSF, only neurofilament light chain (NfL) showed a significant correlation between CSF and plasma levels. Young males with CALD had a 11.3-fold increase in plasma NfL compared with controls. Importantly, 9 of 11 young males with CALD who underwent HCT showed a mean decrease in plasma NfL of 50% at 1 year after HCT compared with pre-HCT levels. In conclusion, plasma NfL could be a great value in determining outcomes in CALD and should be scrutinized in future studies in patients prior to CALD development and after therapeutic intervention.

ETV2 regulates PARP-1 binding protein to induce ER stress-mediated death in tuberin-deficient cells.

Lymphangioleiomyomatosis (LAM) is a rare progressive disease, characterized by mutations in the tuberous sclerosis complex genes (TSC1 or TSC2) and hyperactivation of mechanistic target of rapamycin complex 1 (mTORC1). Here, we report that E26 transformation-specific (ETS) variant transcription factor 2 (ETV2) is a critical regulator of Tsc2-deficient cell survival. ETV2 nuclear localization in Tsc2-deficient cells is mTORC1-independent and is enhanced by spleen tyrosine kinase (Syk) inhibition. In the nucleus, ETV2 transcriptionally regulates poly(ADP-ribose) polymerase 1 binding protein (PARPBP) mRNA and protein expression, partially reversing the observed down-regulation of PARPBP expression induced by mTORC1 blockade during treatment with both Syk and mTORC1 inhibitors. In addition, silencing Etv2 or Parpbp in Tsc2-deficient cells induced ER stress and increased cell death in vitro and in vivo. We also found ETV2 expression in human cells with loss of heterozygosity for TSC2, lending support to the translational relevance of our findings. In conclusion, we report a novel ETV2 signaling axis unique to Syk inhibition that promotes a cytocidal response in Tsc2-deficient cells and therefore maybe a potential alternative therapeutic target in LAM.

powerEQTL: an R package and shiny application for sample size and power calculation of bulk tissue and single-cell eQTL analysis.

SUMMARY: Genome-wide association studies (GWAS) have revealed thousands of genetic loci for common diseases. One of the main challenges in the post-GWAS era is to understand the causality of the genetic variants. Expression quantitative trait locus (eQTL) analysis is an effective way to address this question by examining the relationship between gene expression and genetic variation in a sufficiently powered cohort. However, it is frequently a challenge to determine the sample size at which a variant with a specific allele frequency will be detected to associate with gene expression with sufficient power. This is a particularly difficult task for single-cell RNAseq studies. Therefore, a user-friendly tool to estimate statistical power for eQTL analyses in both bulk tissue and single-cell data is needed. Here, we presented an R package called powerEQTL with flexible functions to estimate power, minimal sample size or detectable minor allele frequency for both bulk tissue and single-cell eQTL analysis. A user-friendly, program-free web application is also provided, allowing users to calculate and visualize the parameters interactively. AVAILABILITY AND IMPLEMENTATION: The powerEQTL R package source code and online tutorial are freely available at CRAN: https://cran.r-project.org/web/packages/powerEQTL/. The R shiny application is publicly hosted at https://bwhbioinfo.shinyapps.io/powerEQTL/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Drug Repurposing to Treat Glucocorticoid Resistance in Asthma.

Corticosteroid resistance causes significant morbidity in asthma, and drug repurposing may identify timely and cost-effective adjunctive treatments for corticosteroid resistance. In 95 subjects from the Childhood Asthma Management Program (CAMP) and 19 subjects from the Severe Asthma Research Program (SARP), corticosteroid response was measured by the change in percent predicted forced expiratory volume in one second (FEV1). In each cohort, differential gene expression analysis was performed comparing poor (resistant) responders, defined as those with zero to negative change in FEV1, to good responders, followed by Connectivity Map (CMap) analysis to identify inversely associated (i.e., negatively connected) drugs that reversed the gene expression profile of poor responders to resemble that of good responders. Mean connectivity scores weighted by sample size were calculated. The top five drug compound candidates underwent in vitro validation in NF-κB-based luciferase reporter A549 cells stimulated by IL-1ß ± dexamethasone. In CAMP and SARP, 134 and 178 respective genes were differentially expressed in poor responders. CMap analysis identified 46 compounds in common across both cohorts with connectivity scores < -50. γ-linolenic acid, ampicillin, exemestane, brinzolamide, and INCA-6 were selected for functional validation. γ-linolenic acid, brinzolamide, and INCA-6 significantly reduced IL-1ß induced luciferase activity and potentiated the anti-inflammatory effect of dexamethasone in A549/NF-κB-luc reporter cells. These results demonstrate how existing drugs, including γ-linolenic acid, brinzolamide, and INCA-6, may be repurposed to improve corticosteroid response in asthmatics.

MTOR suppresses autophagy-mediated production of IL25 in allergic airway inflammation.

INTRODUCTION: Airway epithelial cells are recognised as an essential controller for the initiation and perpetuation of asthmatic inflammation, yet the detailed mechanisms remain largely unknown. This study aims to investigate the roles and mechanisms of the mechanistic target of rapamycin (MTOR)-autophagy axis in airway epithelial injury in asthma. METHODS: We examined the MTOR-autophagy signalling in airway epithelium from asthmatic patients or allergic mice induced by ovalbumin or house dust mites, or in human bronchial epithelial (HBE) cells. Furthermore, mice with specific MTOR knockdown in airway epithelium and autophagy-related lc3b-/- mice were used for allergic models. RESULTS: MTOR activity was decreased, while autophagy was elevated, in airway epithelium from asthmatic patients or allergic mice, or in HBE cells treated with IL33 or IL13. These changes were associated with upstream tuberous sclerosis protein 2 signalling. Specific MTOR knockdown in mouse bronchial epithelium augmented, while LC3B deletion diminished allergen-induced airway inflammation and mucus hyperproduction. The worsened inflammation caused by MTOR deficiency was also ameliorated in lc3b-/- mice. Mechanistically, autophagy was induced later than the emergence of allergen-initiated inflammation, particularly IL33 expression. MTOR deficiency increased, while knocking out of LC3B abolished the production of IL25 and the eventual airway inflammation on allergen challenge. Blocking IL25 markedly attenuated the exacerbated airway inflammation in MTOR-deficiency mice. CONCLUSION: Collectively, these results demonstrate that allergen-initiated inflammation suppresses MTOR and induces autophagy in airway epithelial cells, which results in the production of certain proallergic cytokines such as IL25, further promoting the type 2 response and eventually perpetuating airway inflammation in asthma.

Genome-Wide Association Study: Functional Variant rs2076295 Regulates Desmoplakin Expression in Airway Epithelial Cells.

Rationale: Genetic association studies have identified rs2076295 in association with idiopathic pulmonary fibrosis (IPF). We hypothesized that rs2076295 is the functional variant regulating DSP (desmoplakin) expression in human bronchial epithelial cells, and DSP regulates extracellular matrix-related gene expression and cell migration, which is relevant to IPF development.Objectives: To determine whether rs2076295 regulates DSP expression and the function of DSP in airway epithelial cells.Methods: Using CRISPR (clustered regularly interspaced short palindromic repeat)/Cas9 editing (including regional deletion, indel, CRISPR interference, and single-base editing), we modified rs2076295 and measured DSP expression in edited 16HBE14o- and primary airway epithelial cells. Cellular integrity, migration, and genome-wide gene expression changes were examined in 16HBE14o- single colonies with DSP knockout. The expression of DSP and its relevant matrix genes was measured by quantitative PCR and also analyzed in single-cell RNA-sequencing data from control and IPF lungs.Measurements and Main Results:DSP is expressed predominantly in bronchial and alveolar epithelial cells, with reduced expression in alveolar epithelial cells in IPF lungs. The deletion of the DNA region-spanning rs2076295 led to reduced expression of DSP, and the edited rs2076295GG 16HBE14o- line has lower expression of DSP than the rs2076295TT lines. Knockout of DSP in 16HBE14o- cells decreased transepithelial resistance but increased cell migration, with increased expression of extracellular matrix-related genes, including MMP7 and MMP9. Silencing of MMP7 and MMP9 abolished increased migration in DSP-knockout cells.Conclusions: rs2076295 regulates DSP expression in human airway epithelial cells. The loss of DSP enhances extracellular matrix-related gene expression and promotes cell migration, which may contribute to the pathogenesis of IPF.

The association of telomere length and telomerase activity with adverse outcomes in older patients with non-ST-elevation acute coronary syndrome.

BACKGROUND: Non-ST elevation acute coronary syndrome (NSTEACS) occurs more frequently in older patients with an increased occurrence of recurrent cardiac events following the index presentation. Telomeres are structures consisting of repeated DNA sequences as associated shelterin proteins at the ends of chromosomes. We aim to determine whether telomere length (TL) and telomerase activity (TA) predicted poor outcomes in older patients presenting with NSTEACS undergoing invasive care. METHOD: Older patients undergoing invasive management for NSTEACS were recruited to the ICON-1 biomarker study (NCT01933581). Peripheral blood mononuclear cells (PBMC) were recovered on 153 patients. DNA was isolated and mean TL was measured by quantitative PCR expressed as relative T (telomere repeat copy number) to S (single copy gene number) ratio (T/S ratio), and a telomere repeat amplification assay was used to assess TA during index presentation with NSTEACS. Primary clinical outcomes consisted of death, myocardial infarction (MI), unplanned revascularisation, stroke and significant bleeding recorded at 1 year. TL and TA were divided into tertile groups for analysis. Cox proportional hazards regression was performed. Ordinal regression was performed to evaluate the relationship between TL and TA and traditional cardiovascular risk factors at baseline. RESULTS: 298 patients were recruited in the ICON-1 study of which 153 had PBMC recovered. The mean age was 81.0 ± 4.0 years (64% male). Mean telomere length T/S ratio was 0.47 ± 0.25 and mean TA was 1.52 ± 0.61 units. The primary composite outcome occurred in 44 (28.8%) patients. There was no association between short TL or low TA and incidence of the primary composite outcome (Hazard Ratio [HR] 1.50, 95% Confidence Interval [CI] 0.68-3.34, p = 0.32 and HR 1.33, 95% CI 0.52-3.36, p = 0.51 respectively). CONCLUSION: TL and TA are not found to be associated with the incidence of adverse outcomes in older patients presenting with NSTEACS undergoing invasive care. CLINICAL TRIAL REGISTRATION: URL: https://www.clinicaltrials.gov Unique identifier: NCT01933581.

Model-based clustering for identifying disease-associated SNPs in case-control genome-wide association studies.

Genome-wide association studies (GWASs) aim to detect genetic risk factors for complex human diseases by identifying disease-associated single-nucleotide polymorphisms (SNPs). The traditional SNP-wise approach along with multiple testing adjustment is over-conservative and lack of power in many GWASs. In this article, we proposed a model-based clustering method that transforms the challenging high-dimension-small-sample-size problem to low-dimension-large-sample-size problem and borrows information across SNPs by grouping SNPs into three clusters. We pre-specify the patterns of clusters by minor allele frequencies of SNPs between cases and controls, and enforce the patterns with prior distributions. In the simulation studies our proposed novel model outperforms traditional SNP-wise approach by showing better controls of false discovery rate (FDR) and higher sensitivity. We re-analyzed two real studies to identifying SNPs associated with severe bortezomib-induced peripheral neuropathy (BiPN) in patients with multiple myeloma (MM). The original analysis in the literature failed to identify SNPs after FDR adjustment. Our proposed method not only detected the reported SNPs after FDR adjustment but also discovered a novel BiPN-associated SNP rs4351714 that has been reported to be related to MM in another study.

DNA methylation is associated with inhaled corticosteroid response in persistent childhood asthmatics.

BACKGROUND: Response to inhaled corticosteroids is highly variable, and the association between DNA methylation and treatment response is not known. OBJECTIVE: To examine the association between peripheral blood DNA methylation and inhaled corticosteroid response in children with persistent asthma. METHODS: Epigenome-wide DNA methylation was analysed in individuals on inhaled corticosteroids in three independent and ethnically diverse cohorts-Childhood Asthma Management Program (CAMP); Children, Allergy, Milieu, Stockholm, Epidemiology (BAMSE); and Genetic Epidemiology of Asthma in Costa Rica Study (GACRS). Treatment response was evaluated using two definitions, the absence of emergency department visits and/or hospitalizations and the absence oral corticosteroid use while on inhaled corticosteroid therapy. CpG sites meeting nominal significance (P < 0.05) for each outcome were combined in a three-cohort meta-analysis with adjustment for multiple testing. DNA methylation was correlated with gene expression using Pearson and partial correlations. RESULTS: In 154 subjects from CAMP, 72 from BAMSE, and 168 from GACRS, relative hypomethylation of cg00066816 (171 bases upstream of IL12B) was associated with the absence of emergency department visits and/or hospitalizations (Q = 0.03) in all cohorts and lower IL12B expression (ρ = 0.34, P = 0.01) in BAMSE. Relative hypermethylation of cg04256470 (688 bases upstream of CORT) was associated with the absence of oral corticosteroid use (Q = 0.04) in all cohorts and higher CORT expression (ρ = 0.20, P = 0.045) in CAMP. CONCLUSION AND CLINICAL RELEVANCE: Differential DNA methylation of IL12B and CORT are associated with inhaled corticosteroid treatment response in persistent childhood asthmatics. Pharmaco-methylation can identify novel markers of treatment sensitivity in asthma.

The association between vitamin D status and clinical events in high-risk older patients with non-ST elevation acute coronary syndrome undergoing invasive management.

There is a higher incidence of vitamin D deficiency in older adults. This may play a plausible mechanistic role in the occurrence of increased adverse events after non-ST elevation acute coronary syndrome (NSTEACS). This study investigated whether total vitamin D levels at the time of presentation predicted adverse outcomes in older adults undergoing invasive management of NSTEACS. Of the 629 patients screened, 300 high-risk older adults with NSTEACS managed by an invasive strategy were recruited. Serum total 25-hydroxyvitamin D was measured at index presentation. The primary outcome was defined as 1-year composite of all-cause mortality, acute coronary syndrome (ACS), unplanned repeat revascularisation, significant bleeding or stroke. Mean age was 80.5±4.8 years (61.9% male). Median vitamin D level was 29.5nmol/L [interquartile range IQR 16.0-53.0 nmol/L] and was split equally by the median for analysis forming two groups: high (median vitamin D 53.0 nmol/L [IQR 40.0-75.0]) and low (16.0 nmol/L [11.0-23.0]). The primary outcome occurred in 76 patients (25.9%); 32 (21.9%) in the low group and 44 (29.9%) in the high group, p = 0.12. Multivariable analyses showed no significant difference in the primary composite outcome at 1 year between the low and high group of baseline serum vitamin D (Hazard Ratio 1.20 [95% Confidence Interval 0.72-2.0], p = 0.48). Serum total vitamin D, measured at the time of angiography, was not associated with adverse outcomes at one year in this high-risk older cohort of patients with NSTEACS undergoing invasive management.

An Investigation about Gene Modules Associated with hDPSC Differentiation for Adolescents.

Dental pulp stem cells (DPSCs) have the property of self-renewal and multidirectional differentiation so that they have the potential for future regenerative therapy of various diseases. The latest breakthrough in the biology of stem cells and the development of regenerative biology provides an effective strategy for regenerative therapy. However, in the medium promoting differentiation during long-term passage, DPSCs would lose their differentiation capability. Some efforts have been made to find genes influencing human DPSC (hDPSC) differentiation based on hDPSCs isolated from adults. However, hDPSC differentiation is a very complex process, which involves multiple genes and multielement interactions. The purpose of this study is to detect sets of correlated genes (i.e., gene modules) that are associated to hDPSC differentiation at the crown-completed stage of the third molars, by using weighted gene coexpression network analysis (WGCNA). Based on the gene expression dataset GSE10444 from Gene Expression Omnibus (GEO), we identified two significant gene modules: yellow module (742 genes) and salmon module (9 genes). The WEB-based Gene SeT AnaLysis Toolkit showed that the 742 genes in the yellow module were enriched in 59 KEGG pathways (including Wnt signaling pathway), while the 9 genes in the salmon module were enriched in one KEGG pathway (neurotrophin signaling pathway). There were 660 (7) genes upregulated at P10 and 82 (2) genes downregulated at P10 in the yellow (salmon) module. Our results provide new insights into the differentiation capability of hDPSCs.

Pilot study of DNA methylation, molecular aging markers and measures of health and well-being in aging.

Relations of DNA methylation markers to other biological aging markers and to psychosocial, behavioral, and health measures remain unclear. The sample included 23 participants (n = 11 cases with psychiatric diagnoses and n = 12 controls without current or lifetime psychiatric disorder), balanced by age and sex. Genomic DNA was extracted from blood samples; the following were performed: genome-wide DNA methylation assay using Illumina 850k methylationEPIC; PCR assays for relative telomere length (RTL) and mitochondrial DNA copy number (mtCN). Exposures were: case status; depression and anxiety symptoms; psychosocial support; subjective and objective cognition. Outcomes were: DNA methylation age (DNAm age); RTL; mtCN; extrinsic and intrinsic epigenetic age acceleration (EEAA and IEAA). Stronger correlation with chronological age was observed for DNAm age (ρ = 0.86; p < 0.0001) compared to RTL (ρ = -0.53; p < 0.01); mtCN was not correlated with age. DNAm age was more strongly correlated with behavioral and health variables than RTL or mtCN; e.g., correlations with DNAm age: body mass index (ρ = 0.36; p = 0.10); smoking pack-years (ρ = 0.37; p = 0.08); physical activity (ρ = -0.56; p = 0.01); alcohol intake (ρ = 0.56; p = 0.01). DNAm age was inversely correlated with psychosocial support (ρ = -0.42; p = 0.048) and Modified Mini-Mental State score (ρ = -0.44; p = 0.01). Anxiety, psychosocial support, and objective cognition were significantly related to accelerated aging; depression and subjective cognition were not. In conclusion, DNAm age correlated more strongly with chronological age and key psychosocial, behavioral, and health variables than RTL or mtCN. Signals for associations with epigenetic aging were observed for psychosocial and neurobehavioral variables.

Novel Data Transformations for RNA-seq Differential Expression Analysis.

We propose eight data transformations (r, r2, rv, rv2, l, l2, lv, and lv2) for RNA-seq data analysis aiming to make the transformed sample mean to be representative of the distribution center since it is not always possible to transform count data to satisfy the normality assumption. Simulation studies showed that for data sets with small (e.g., nCases = nControls = 3) or large sample size (e.g., nCases = nControls = 100) limma based on data from the l, l2, and r2 transformations performed better than limma based on data from the voom transformation in term of accuracy, FDR, and FNR. For datasets with moderate sample size (e.g., nCases = nControls = 30 or 50), limma with the rv and rv2 transformations performed similarly to limma with the voom transformation. Real data analysis results are consistent with simulation analysis results: limma with the r, l, r2, and l2 transformation performed better than limma with the voom transformation when sample sizes are small or large; limma with the rv and rv2 transformations performed similarly to limma with the voom transformation when sample sizes are moderate. We also observed from our data analyses that for datasets with large sample size, the gene-selection via the Wilcoxon rank sum test (a non-parametric two sample test method) based on the raw data outperformed limma based on the transformed data.

New Statistical Methods for Constructing Robust Differential Correlation Networks to characterize the interactions among microRNAs.

The interplay among microRNAs (miRNAs) plays an important role in the developments of complex human diseases. Co-expression networks can characterize the interactions among miRNAs. Differential correlation network is a powerful tool to investigate the differences of co-expression networks between cases and controls. To construct a differential correlation network, the Fisher's Z-transformation test is usually used. However, the Fisher's Z-transformation test requires the normality assumption, the violation of which would result in inflated Type I error rate. Several bootstrapping-based improvements for Fisher's Z test have been proposed. However, these methods are too computationally intensive to be used to construct differential correlation networks for high-throughput genomic data. In this article, we proposed six novel robust equal-correlation tests that are computationally efficient. The systematic simulation studies and a real microRNA data analysis showed that one of the six proposed tests (ST5) overall performed better than other methods.

Urine podoplanin heralds the onset of ischemia-reperfusion injury of the kidney.

Ischemia-reperfusion injury represents one of the most common causes of acute kidney injury, a serious and often deadly condition that affects up to 20% of all hospitalized patients in the United States. However, the current standard assay used universally for the diagnosis of acute kidney injury, serum creatinine, does not detect renal damage early in its course. Serendipitously, we found that the immunofluorescent signal of the constitutive podocyte marker podoplanin fades in the glomerulus and intensifies in the tubulointerstitial compartment of the kidney shortly after ischemia-reperfusion injury in 8- to 10-wk-old male C57Bl/6j mice. Therefore, we sought to define the appearance and course of the podoplanin-positive signal in the kidney after ischemia-reperfusion injury. The tubulointerstitial podoplanin-positive signal increased as early as 2 h but persisted for 7 days after ischemia-reperfusion injury. In addition, the strength of this tubulointerstitial signal was directly proportional to the severity of ischemia, and its location shifted from the tubules to interstitial cells over time. Finally, we detected podoplanin in the urine of mice after ischemia, and we observed that an increase in the urine podoplanin-to-creatinine ratio correlated strongly with the onset of renal ischemia-reperfusion injury. Our findings indicate that the measurement of urine podoplanin harbors promising potential for use as a novel biomarker for the early detection of ischemia-reperfusion injury of the kidney.

First-in-class immune-modulating small molecule Icaritin in advanced hepatocellular carcinoma: preliminary results of safety, durable survival and immune biomarkers.

BACKGROUND: With poor prognosis and limited treatment options for advanced hepatocellular carcinoma (HCC), development of novel therapeutic agents is urgently needed. This single-arm phase I study sought to assess the safety and preliminary efficacy of icaritin in human as a potential oral immunotherapy in addition to the immune-checkpoint inhibitors. METHODS: Eligible advanced HCC patients with Child-Pugh Class A or B were administered with a fixed oral dose of icaritin at either 600 or 800 mg b.i.d. The primary endpoint was safety, and the secondary endpoints included time-to-progression (TTP), overall survival (OS) and the clinical benefit rate (CBR). Icaritin treatment induced immune biomarkers and immune-modulating activities in myeloid cells were also explored. RESULTS: No drug-related adverse events ≥ Grade 3 were observed in all 20 enrolled HCC patients. Among the 15 evaluable patients, 7 (46.7%) achieved clinical benefit, representing one partial response (PR, 6.7%) and 6 stable disease (SD, 40%). The median TTP was 141 days (range: 20-343 days), and the median OS was 192 days (range: 33-1036 days). Durable survival was observed in PR/SD patients with a median OS of 488 days (range: 72-773). TTP was significantly associated with the dynamic changes of peripheral neutrophils (p = 0.0067) and lymphocytes (p = 0.0337). Icaritin treatment induced changes in immune biomarkers-and immune-suppressive myeloid cells were observed. CONCLUSIONS: Icaritin demonstrated safety profiles and preliminary durable survival benefits in advanced HCC patients, which were correlated with its immune-modulation activities and immune biomarkers. These results suggested the potential of icaritin as a novel oral immunotherapy for advanced HCC in addition to antibody-based PD-1/PD-L1 blockade therapies. TRIAL REGISTRATION: Clinicaltrial.gov identifier. NCT02496949 (retrospectively registered, July 14, 2015).

Coronary artery lesion phenotype in frail older patients with non-ST-elevation acute coronary syndrome undergoing invasive care.

AIMS: The association of frailty with coronary plaque phenotype among older patients with non-ST-elevation acute coronary syndrome (NSTEACS) is not known. The aim of this study was to evaluate the association of frailty with coronary plaque phenotype among older patients with NSTEACS. METHODS AND RESULTS: Older patients with NSTEACS who underwent invasive angiography were recruited. Frailty was measured using the Fried frailty score. Following angiography, patients underwent greyscale and virtual histology intravascular ultrasound (VH-IVUS) imaging. Of the 90 patients, 26 (28.9%) were robust, 49 (54.4%) patients were pre-frail, and 15 (16.7%) were frail. Mean age was 80.9±3.8 years; 59 (65.6%) were male. Compared to robust patients, the pre-frail group had a significantly greater presence of high-risk lesions including VH thin-cap fibroatheroma (TCFA, p=0.011), minimum lumen area (MLA) ≤4 mm2 (p=0.016), TCFA+MLA ≤4 mm2 (p=0.005), TCFA+plaque burden (PB) ≥70% (p=0.005) and TCFA+PB ≥70%+MLA ≤4 mm2 (p=0.003). By age- and sex-adjusted logistic regression analysis, frailty was found to be strongly and independently associated with the presence of TCFA (odds ratio [OR] 2.81, 95% confidence interval [CI]:1.06-7.48, p=0.039). CONCLUSIONS: This is the first study to report the relationship between frailty phenotype and coronary plaque morphology among frail older NSTEACS patients. ClinicalTrials.gov Identifier: NCT01933581.