A fast and scalable framework for large-scale and ultrahigh-dimensional sparse regression with application to the UK Biobank.

The UK Biobank is a very large, prospective population-based cohort study across the United Kingdom. It provides unprecedented opportunities for researchers to investigate the relationship between genotypic information and phenotypes of interest. Multiple regression methods, compared with genome-wide association studies (GWAS), have already been showed to greatly improve the prediction performance for a variety of phenotypes. In the high-dimensional settings, the lasso, since its first proposal in statistics, has been proved to be an effective method for simultaneous variable selection and estimation. However, the large-scale and ultrahigh dimension seen in the UK Biobank pose new challenges for applying the lasso method, as many existing algorithms and their implementations are not scalable to large applications. In this paper, we propose a computational framework called batch screening iterative lasso (BASIL) that can take advantage of any existing lasso solver and easily build a scalable solution for very large data, including those that are larger than the memory size. We introduce snpnet, an R package that implements the proposed algorithm on top of glmnet and optimizes for single nucleotide polymorphism (SNP) datasets. It currently supports ℓ1-penalized linear model, logistic regression, Cox model, and also extends to the elastic net with ℓ1/ℓ2 penalty. We demonstrate results on the UK Biobank dataset, where we achieve competitive predictive performance for all four phenotypes considered (height, body mass index, asthma, high cholesterol) using only a small fraction of the variants compared with other established polygenic risk score methods.

TIFA promotes colorectal cancer cell proliferation in an RSK- and PRAS40-dependent manner.

Previous studies have reported that TIFA plays different roles in various tumor types. However, the function of TIFA in colorectal cancer (CRC) remains unclear. Here, we showed that the expression of TIFA was markedly increased in CRC versus normal tissue, and positively correlated with CRC TNM stages. In agreement, we found that the CRC cell lines show increased TIFA expression levels versus normal control. The knockdown of TIFA inhibited cell proliferation but had no effect on cell apoptosis in vitro or in vivo. Moreover, the ectopic expression of TIFA enhanced cell proliferation ability in vitro and in vivo. In contrast, the expression of mutant TIFA (T9A, oligomerization site mutation; D6, TRAF6 binding site deletion) abolished TIFA-mediated cell proliferation enhancement. Exploration of the underlying mechanism revealed that the protein synthesis-associated kinase RSK and PRAS40 activation were responsible for TIFA-mediated CRC progression. In summary, these findings suggest that TIFA plays a role in mediating CRC progression. This could provide a promising target for CRC therapy.

Genome-wide analysis of cell-Free DNA methylation profiling with MeDIP-seq identified potential biomarkers for colorectal cancer.

BACKGROUND: Colorectal cancer is the most common malignancy and the third leading cause of cancer-related death worldwide. This study aimed to identify potential diagnostic biomarkers for colorectal cancer by genome-wide plasma cell-free DNA (cfDNA) methylation analysis. METHODS: Peripheral blood from colorectal cancer patients and healthy controls was collected for cfDNA extraction. Genome-wide cfDNA methylation profiling, especially differential methylation profiling between colorectal cancer patients and healthy controls, was performed by methylated DNA immunoprecipitation coupled with high-throughput sequencing (MeDIP-seq). Logistic regression models were established, and the accuracy of this diagnostic model for colorectal cancer was verified using tissue-sourced data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) due to the lack of cfDNA methylation data in public datasets. RESULTS: Compared with the control group, 939 differentially methylated regions (DMRs) located in promoter regions were found in colorectal cancer patients; 16 of these DMRs were hypermethylated, and the remaining 923 were hypomethylated. In addition, these hypermethylated genes, mainly PRDM14, RALYL, ELMOD1, and TMEM132E, were validated and confirmed in colorectal cancer by using publicly available DNA methylation data. CONCLUSIONS: MeDIP-seq can be used as an optimal approach for analyzing cfDNA methylomes, and 12 probes of four differentially methylated genes identified by MeDIP-seq (PRDM14, RALYL, ELMOD1, and TMEM132E) could serve as potential biomarkers for clinical application in patients with colorectal cancer.

2cChIP-seq and 2cMeDIP-seq: The Carrier-Assisted Methods for Epigenomic Profiling of Small Cell Numbers or Single Cells.

Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) can profile genome-wide epigenetic marks associated with regulatory genomic elements. However, conventional ChIP-seq is challenging when examining limited numbers of cells. Here, we developed a new technique by supplementing carrier materials of both chemically modified mimics with epigenetic marks and dUTP-containing DNA fragments during conventional ChIP procedures (hereafter referred to as 2cChIP-seq), thus dramatically improving immunoprecipitation efficiency and reducing DNA loss of low-input ChIP-seq samples. Using this strategy, we generated high-quality epigenomic profiles of histone modifications or DNA methylation in 10-1000 cells. By introducing Tn5 transposase-assisted fragmentation, 2cChIP-seq reliably captured genomic regions with histone modification at the single-cell level in about 100 cells. Moreover, we characterized the methylome of 100 differentiated female germline stem cells (FGSCs) and observed a particular DNA methylation signature potentially involved in the differentiation of mouse germline stem cells. Hence, we provided a reliable and robust epigenomic profiling approach for small cell numbers and single cells.

Found In Translation: a machine learning model for mouse-to-human inference.

Cross-species differences form barriers to translational research that ultimately hinder the success of clinical trials, yet knowledge of species differences has yet to be systematically incorporated in the interpretation of animal models. Here we present Found In Translation (FIT; http://www.mouse2man.org ), a statistical methodology that leverages public gene expression data to extrapolate the results of a new mouse experiment to expression changes in the equivalent human condition. We applied FIT to data from mouse models of 28 different human diseases and identified experimental conditions in which FIT predictions outperformed direct cross-species extrapolation from mouse results, increasing the overlap of differentially expressed genes by 20-50%. FIT predicted novel disease-associated genes, an example of which we validated experimentally. FIT highlights signals that may otherwise be missed and reduces false leads, with no experimental cost.

High-dimensional regression adjustments in randomized experiments.

We study the problem of treatment effect estimation in randomized experiments with high-dimensional covariate information and show that essentially any risk-consistent regression adjustment can be used to obtain efficient estimates of the average treatment effect. Our results considerably extend the range of settings where high-dimensional regression adjustments are guaranteed to provide valid inference about the population average treatment effect. We then propose cross-estimation, a simple method for obtaining finite-sample-unbiased treatment effect estimates that leverages high-dimensional regression adjustments. Our method can be used when the regression model is estimated using the lasso, the elastic net, subset selection, etc. Finally, we extend our analysis to allow for adaptive specification search via cross-validation and flexible nonparametric regression adjustments with machine-learning methods such as random forests or neural networks.

Copper Chaperone for Superoxide Dismutase Subtypes as a Prognostic Marker in Luminal B Breast Cancer.

Background: Copper chaperone for superoxide dismutase (CCS) is an essential component of the oxidation-reduction system. In breast cancer cells, CCS expression is highly up-regulated, which contributes to cellular proliferation and migration. Breast cancer is a multifaceted disease with different tumor prognoses and responses to clinical treatments, which may be associated with multiple molecular subtypes of CCS. Methods: The CCS expression patterns in breast cancer were investigated by TNMplot, cBioPortal, and HPA network database. The correlation of CCS expression with clinicopathological parameters was analyzed using the UALCAN database. The Cancer Genome Atlas (TCGA) data set was used to analyze the Clinical characteristics of CCS in luminal B patients. The bc-GenExMiner database was used to analyze the effects of BReast-CAncer susceptibility gene (BRCA)1/2, TP53 mutation status, estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor (HER) expression on CCS expression. The survival curves and prognostic value of CCS in luminal B breast cancer were performed through Kaplan-Meier curves, univariate and multivariate Cox regression using the PrognoScan, bc-GenExMiner, and Clinical bioinformatics analysis platform. Results: We found that CCS expression was associated with patient age, race, ER, and PR status. We also discovered that BRCA1/2 mutations had an effect on CCS expression. The luminal B subtype had the highest CCS expression, which was linked to poor survival compared with other subtypes. In addition, Kaplan-Meier curve analysis showed that luminal B patients with high CCS mRNA expression showed a poor survival and the CCS gene is an independent predictor of outcome in patients with luminal B breast cancer by univariate and multivariate Cox regression. Conclusions: Our findings emphasize the significant expression of CCS in luminal B breast cancer and its potential as an autonomous prognostic determinant for this specific molecular subtype. These findings suggest that CCS holds promise as a prospective marker for the treatment of luminal B breast cancer.

Integrative transcriptome analysis reveals alternative polyadenylation potentially contributes to GCRV early infection.

Introduction: Grass carp reovirus (GCRV), a member of the Aquareovirus genus in the Reoviridae family, is considered to be the most pathogenic aquareovirus. Productive viral infection requires extensive interactions between viruses and host cells. However, the molecular mechanisms underlying GCRV early infection remains elusive. Methods: In this study we performed transcriptome and DNA methylome analyses with Ctenopharyngodon idellus kidney (CIK) cells infected with GCRV at 0, 4, and 8 h post infection (hpi), respectively. Results: We found that at early infection stage the differentially expressed genes related to defense response and immune response in CIK cells are activated. Although DNA methylation pattern of CIK cells 8 hpi is similar to mock-infected cells, we identified a considerable number of genes that selectively utilize alternative polyadenylation sites. Particularly, we found that biological processes of cytoskeleton organization and regulation of microtubule polymerization are statistically enriched in the genes with altered 3'UTRs. Discussion: Our results suggest that alternative polyadenylation potentially contributes to GCRV early infection.

BMI1 fine-tunes gene repression and activation to safeguard undifferentiated spermatogonia fate.

Introduction: Spermatogenesis is sustained by the homeostasis of self-renewal and differentiation of undifferentiated spermatogonia throughout life, which is regulated by transcriptional and posttranscriptional mechanisms. B cell-specific Moloney murine leukemia virus integration site 1 (BMI1), one of spermatogonial stem cell markers, is a member of Polycomb repressive complex 1 (PRC1) and important to spermatogenesis. However, the mechanistic underpinnings of how BMI1 regulates spermatogonia fate remain elusive. Methods: We knocked down BMI1 by siRNA to investigate the role of BMI1 in undifferentiated spermatogonia. Differentially expressed genes were identified by RNA-seq and used for KEGG pathway analysis. We performed ChIP-seq analysis in wild type and BMI1 knockdown cells to explore the underlying molecular mechanisms exerted by BMI1. BMI1-associated alterations in repressive histone modifications were detected via Western blotting and ChIP-seq. Furthermore, we performed mass spectrometry and Co-immunoprecipitation assays to investigate BMI1 co-factors. Finally, we demonstrated the genomic regions occupied by both BMI1 and its co-factor. Results: BMI1 is required for undifferentiated spermatogonia maintenance by both repressing and activating target genes. BMI1 preserves PI3K-Akt signaling pathway for spermatogonia proliferation. Decrease of BMI1 affects the deposition of repressive histone modifications H2AK119ub1 and H3K27me3. BMI also positively regulates H3K27ac deposited genes which are associated with proliferation. Moreover, we demonstrate that BMI1 interacts with Sal-like 4 (SALL4), the transcription factor critical for spermatogonia function, to co-regulate gene expression. Discussion: Overall, our study reveals that BMI1 safeguards undifferentiated spermatogonia fate through multi-functional roles in regulating gene expression programs of undifferentiated spermatogonia.

TBX21 attenuates colorectal cancer progression via an ARHGAP29/RSK/GSK3ß dependent manner.

PURPOSE: Previous studies have shown that TBX21 (T-Box Transcription Factor 21) plays a vital role in coordinating multiple aspects of the immune response especially type 1 immune response as well as tumor progression. However, the function of TBX21 in colorectal cancer (CRC) remains unclear. METHODS: IHC to investigate TBX21 expression in CRC tissues. Cell proliferation and apoptosis assays to validate TBX21 function in vitro and in vivo. RNA-seq assay to explore target genes of TBX21. Human phospho-kinase array assay to explore down-stream signaling of TBX21. RESULTS: We disclosed that the expression of TBX21 was marked decreased in CRC versus normal tissue, and negatively correlated with CRC TNM stages. Surprisingly, we found that the CRC and normal cell lines show no TBX21 expression levels. Ectopic expression of TBX21 inhibited cell proliferation and promoted cell apoptosis in vitro. Moreover, RNA-sequence data first time showed that ARHGAP29 acts as the target gene of TBX21 to mediate down-stream signaling activation. Human phospho-kinase array data first time displayed that ectopic expression of TBX21 reduced kinase RSK and GSK3ß activation. In contrast, knocked down the expression of TBX21 or ARHGAP29 alternatively abolished TBX21 mediated cell proliferation suppression, cell apoptosis enhancement and RSK/GSK3ß activation. In addition, xenograft model studies demonstrated that TBX21 inhibits colorectal tumor progression via ARHGAP29/ RSK/ GSK3ß signaling in vivo. CONCLUSIONS: In summary, the aforementioned findings suggest a model of TBX21 in suppressing CRC progression. This may provide a promising target for CRC therapy.

E-cadherin maintains the undifferentiated state of mouse spermatogonial progenitor cells via ß-catenin.

BACKGROUND: Cadherins play a pivotal role in facilitating intercellular interactions between spermatogonial progenitor cells (SPCs) and their surrounding microenvironment. Specifically, E-cadherin serves as a cellular marker of SPCs in many species. Depletion of E-cadherin in mouse SPCs showed no obvious effect on SPCs homing and spermatogenesis. RESULTS: Here, we investigated the regulatory role of E-cadherin in regulating SPCs fate. Specific deletion of E-cadherin in germ cells was shown to promote SPCs differentiation, evidencing by reduced PLZF+ population and increased c-Kit+ population in mouse testes. E-cadherin loss down-regulated the expression level of ß-catenin, leading to the reduced ß-catenin in nuclear localization for transcriptional activity. Remarkably, increasing expression level of Cadherin-22 (CDH22) appeared specifically after E-cadherin deletion, indicating CDH22 played a synergistic effect with E-cadherin in SPCs. By searching for the binding partners of ß-catenin, Lymphoid enhancer-binding factor 1 (LEF1), T-cell factor (TCF3), histone deacetylase 4 (HDAC4) and signal transducer and activator 3 (STAT3) were identified as suppressors of SPCs differentiation by regulating acetylation of differentiation genes with PLZF. CONCLUSIONS: Two surface markers of SPCs, E-cadherin and Cadherin-22, synergically maintain the undifferentiation of SPCs via the pivotal intermediate molecule ß-catenin. LEF1, TCF3, STAT3 and HDAC4 were identified as co-regulatory factors of ß-catenin in regulation of SPC fate. These observations revealed a novel regulatory pattern of cadherins on SPCs fate.

Functional evaluation of various ICAM3 transcript variants in diffuse large B-Cell lymphoma.

Previous studies have identified several ICAM3 transcript variants and mainly investigated the function of the longest transcript of ICAM3 in various tumor progressions. However, the role of the other ICAM3 transcript variants remains unclear. Herein, we detected the expression of ICAM3 transcript variants 1-4 in DLBCL cells and tumor tissues, disclosed that variants 1, 3, and 4 were expressed in normal B cell lines and 3 DLBCL cell lines except SU-DHL-2 as well as tumor tissues, while variant 2 was not detected. Moreover, we found that ectopic expression of variants 1-4 enhanced cell proliferation by accelerating the cell cycle in SU-DHL2 cells in vitro. In addition, variants 1-4 overexpression showed no effects on SU-DHL2 cell apoptosis. Interestingly, the expression of variants 1, 3, and 4 promoted cell migration and EMT process while variant 2 had no effects. Collectively, the above results displayed the different roles of ICAM3 transcript variants in mediating DLBCL progression.

Platelet transcriptome identifies progressive markers and potential therapeutic targets in chronic myeloproliferative neoplasms.

Predicting disease progression remains a particularly challenging endeavor in chronic degenerative disorders and cancer, thus limiting early detection, risk stratification, and preventive interventions. Here, profiling the three chronic subtypes of myeloproliferative neoplasms (MPNs), we identify the blood platelet transcriptome as a proxy strategy for highly sensitive progression biomarkers that also enables prediction of advanced disease via machine-learning algorithms. The MPN platelet transcriptome reveals an incremental molecular reprogramming that is independent of patient driver mutation status or therapy. Subtype-specific markers offer mechanistic and therapeutic insights, and highlight impaired proteostasis and a persistent integrated stress response. Using a LASSO model with validation in two independent cohorts, we identify the advanced subtype MF at high accuracy and offer a robust progression signature toward clinical translation. Our platelet transcriptome snapshot of chronic MPNs demonstrates a proof-of-principle for disease risk stratification and progression beyond genetic data alone, with potential utility in other progressive disorders.

[Met-analysis on randomized controlled clinical trials of acupuncture and moxibustion on constipation].

OBJECTIVE: To assess the efficacy of acupuncture and moxibustion on constipation. METHODS: A retrieval on literatures concerning treatment of constipation with acupuncture was carried out in databases of VIP, CNKI, WANFANG and PubMed. And meta-analyses were conducted on randomized controlled trial (RCT) and controlled clinical trial (CCT) which met the enrolling requirements. RESULTS: A total number of 15 papers involving 1 052 patients were concluded. The result indicated that the curative rate of acupuncture and moxibustion on constipation is better than ordinary medication (RR = 1.92, 95% CI 1.61-2.30, Z = 7.18, P < 0.000 01). And statistical significance can be found between acupuncture-moxibustion treatment and the routine medicine treatment (RR = 1.26, 95% CI 1.18-1.34, Z = 7.26, P < 0.000 01). In the comparison of abdominal pain, defecation duration and general symptom scores, statistical significance can be found between the differences of acupuncture and moxibustion group and control group (abdominal pain: WMD = -0.22, 95% CI-0.32-0.12, Z = 4.28, P < 0.000 1; defecation duration: WMD = -0.47, 95% CI-0.79-0.15, Z = 2.85, P < 0.004; general symptom scores: WMD = -0.41, 95% CI-0.79-0.03, Z = 2.13, P = 0.03). CONCLUSION: Acupuncture and moxibustion is effective to treat constipation. It has certain advantage when compare with the routine medication treatment. However, since singleness still exists in the index of assessment on therapeutic effect of constipation, and the number of RCT and CCT literatures, especially high-quality, large samples and multi-center reports were still not enough, further studies are still necessary for approving the above conclusions.