Joint analysis of single-cell and bulk tissue sequencing data to infer intratumor heterogeneity

Many computational methods have been developed to discern intratumor heterogeneity (ITH) using DNA sequence data from bulk tumor samples. These methods share an assumption that two mutations arise from the same subclone if they have similar mutant allele-frequencies (MAFs), and thus it is difficult or impossible to distinguish two subclones with similar MAFs. Single-cell DNA sequencing (scDNA-seq) data can be very informative for ITH inference. However, due to the difficulty of DNA amplification, scDNA-seq data are often very noisy. A promising new study design is to collect both bulk and single-cell DNA-seq data and jointly analyze them to mitigate the limitations of each data type. To address the analytic challenges of this new study design, we propose a computational method named BaSiC (B ulk tumor a nd Si ngle C ell), to discern ITH by jointly analyzing DNA-seq data from bulk tumor and single cells. We demonstrate that BaSiC has comparable or better performance than the methods using either data type. We further evaluate BaSiC using bulk tumor and single-cell DNA-seq data from a breast cancer patient and several leukemia patients.     

ROC-guided survival trees and ensembles

Tree-based methods are popular nonparametric tools in studying time-to-event outcomes. In this article, we introduce a novel framework for survival trees and ensembles, where the trees partition the dynamic survivor population and can handle time-dependent covariates. Using the idea of randomized tests, we develop generalized time-dependent receiver operating characteristic (ROC) curves for evaluating the performance of survival trees. The tree-building algorithm is guided by decision-theoretic criteria based on ROC, targeting specifically for prediction accuracy. To address the instability issue of a single tree, we propose a novel ensemble procedure based on averaging martingale estimating equations, which is different from existing methods that average the predicted survival or cumulative hazard functions from individual trees. Extensive simulation studies are conducted to examine the performance of the proposed methods. We apply the methods to a study on AIDS for illustration.     

Variable selection in joint frailty models of recurrent and terminal events

Recurrent event data are commonly encountered in biomedical studies. In many situations, they are subject to an informative terminal event, for example, death. Joint modeling of recurrent and terminal events has attracted substantial recent research interests. On the other hand, there may exist a large number of covariates in such data. How to conduct variable selection for joint frailty proportional hazards models has become a challenge in practical data analysis. We tackle this issue on the basis of the “minimum approximated information criterion” method. The proposed method can be conveniently implemented in SAS Proc NLMIXED for commonly used frailty distributions. Its finite-sample behavior is evaluated through simulation studies. We apply the proposed method to model recurrent opportunistic diseases in the presence of death in an AIDS study.     

Deregulation of UCA1 expression may be involved in the development of chemoresistance to cisplatin in the treatment of non-small-cell lung cancer via regulating the signaling pathway of microRNA-495/NRF2

Non-small-cell lung cancer (NSCLC) remains the leading cause of cancer death worldwide. As a platinum-based chemotherapeutic drug, cisplatin has been used for over 30 years in NSCLC treatment while its effects are diminished by drug resistance. Therefore, we aimed to study the potential role of UCA1 in the development of chemoresistance against cisplatin. Real-time polymerase chain reaction, western-blot analysis, and immunofluorescence were used to study the involvement of UCA1, miR-495, and NRF2 in chemoresistance against cisplatin. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay was performed to determine the effect of cisplatin on cell proliferation. Computational analysis and luciferase assay were carried out to explore the interaction among UCA1, miR-495, and NRF2. The cisplatin-R group exhibited lower levels of UCA1 and NRF2 expression but a higher level of miR-495 expression than the cisplatin-S group. The growth rate and half-maximal inhibitory concentration of cellular dipeptidyl peptidase (cisplatinum) of the cisplatin-R group were much higher than those in the cisplatin-S group. MiR-495 contained a complementary binding site of UCA1, and the luciferase activity of wild-type UCA1 was significantly reduced after the transfection of miR-495 mimics. MiR-495 directly targeted the 3′-untranslated region (3′-UTR) of NRF2, and the luciferase activity of wild-type NRF2 3′-UTR was evidently inhibited by miR-495 mimics. Finally, UCA1 and NRF2 expressions in the effective group were much lower than that in the ineffective group, along with a much higher level of miR-495 expression. We suggested for the first time that high expression of UCA1 contributed to the development of chemoresistance to cisplatin through the UCA1/miR-495/NRF2 signaling pathway.     

Cover Image

Synthetic microbial communities have become a focus of biotechnological research since they can overcome several of the limitations of single-specie cultures. A paradigmatic example is Clostridium cellulovorans DSM 743B, which can decompose lignocellulose but cannot produce butanol. Clostridium beijerinckii NCIMB 8052 however, is unable to use lignocellulose but can produce high amounts of butanol from simple sugars. In our previous studies, both organisms were cocultured to produce butanol by consolidated bioprocessing. However, such consolidated bioprocessing implementation strongly depends on pH regulation. Since low pH (pH 4.5–5.5) is required for butanol fermentation, C. cellulovorans cannot grow well and saccharify sufficient lignocellulose to feed both strains at a pH below 6.4. To overcome this bottleneck, this study engineered C. cellulovorans by adaptive laboratory evolution, inactivating cell wall lyases genes (Clocel_0798 and Clocel_2169), and overexpressing agmatine deiminase genes (augA, encoded by Cbei_1922) from C. beijerinckii NCIMB 8052. The generated strain WZQ36: 743B*6.0*3△lyt0798△lyt2169-(pXY1-P_thl-augA) can tolerate a pH of 5.5. Finally, the alcohol aldehyde dehydrogenase gene adhE1 from Clostridium acetobutylicum ATCC 824 was introduced into the strain to enable butanol production at low pH, in coordination with solvent fermentation of C. beijerinckii in consortium. The engineered consortium produced 3.94 g/L butanol without pH control within 83 hr, which is more than 5-fold of the level achieved by wild consortia under the same conditions. This exploration represents a proof of concept on how to combine metabolic and evolutionary engineering to coordinate coculture of a synthetic microbial community.     

Genome-scale modeling for Bacillus coagulans to understand the metabolic characteristics

Lactic acid is widely used in many industries, especially in the production of poly-lactic acid. Bacillus coagulans is a promising lactic acid producer in industrial fermentation due to its thermophilic property. In this study, we developed the first genome-scale metabolic model (GEM) of B. coagulans iBag597, together with an enzyme-constrained model ec-iBag597. We measured strain-specific biomass composition and integrated the data into a biomass equation. Then, we validated iBag597 against experimental data generated in this study, including amino acid requirements and carbon source utilization, showing that simulations were generally consistent with the experimental results. Subsequently, we carried out chemostats to investigate the effects of specific growth rate and culture pH on metabolism of B. coagulans. Meanwhile, we used iBag597 to estimate the intracellular metabolic fluxes for those conditions. The results showed that B. coagulans was capable of generating ATP via multiple pathways, and switched among them in response to various conditions. With ec-iBag597, we estimated the protein cost and protein efficiency for each ATP-producing pathway to investigate the switches. Our models pave the way for systems biology of B. coagulans, and our findings suggest that maintaining a proper growth rate and selecting an optimal pH are beneficial for lactate fermentation.     

MiR-34a inhibits the proliferation,migration, and invasion of oral squamous cell carcinoma by directly targeting SATB2

In various kinds of carcinomas, the special AT-rich sequence-binding protein 2 (SATB2) with its atypical expression promotes the metastasis and progression of the tumor, though in the oral squamous cell carcinoma (OSCC) its inherent mechanism and the status of SATB2 remain unclear. The role played by the SATB2 expression in the OSCC cell lines and tissue samples in the target of miR-34a downstream is the intended endeavor of this study. In te OSCCs the miR-34a expression was determined by quantitative real-time polymerase chain reaction (q-PCR), while the SATB2 expression in the cell lines and tissue samples in OSCC was analyzed with the q-PCR and the western blot. Studies in both in vitro and in vivo of the effects of miR-34a on the initiation of OSCC were conducted. As a direct target of the miR-34a the SATB2 was verified with the luciferase reporter assay. In cases where the miR-34a levels were low, the SATB2 in OSCCs seemed to be overexpressed. Besides, both in the in vitro and in vivo a suppression of migration, invasion, and cell growth was caused by miR-34a by down regulating the SATB2 expression. The SATB2 being a direct target of miR-34a was confirmed by the cotransfection of miR-34a mimics specifically the decrease in the expression of luciferase of SATB2–3′UTR-wt reporter. As a whole, our study confirmed the inhibition of miR-34a in the invasion, proliferation, and migration of the OSCCs, playing a potential tumor suppressor role with SATB2 as its downstream target.     

miR-181a/b-5p ameliorates inflammatory response in monocrotaline-induced pulmonary arterial hypertension by targeting endocan

Pulmonary arterial hypertension (PAH) is a progressive disorder characterized by vascular remodeling, endothelial cell (EC) dysfunction, and inflammation. The roles of microRNAs have received much critical attention. Thus, this study was attempted to show the biological function of miR-181a/b-5p (miR-181a/b) in monocrotaline (MCT)-induced PAH. Here, rats injected with MCT were used as PAH models. The expression of miR-181a/b and its effect on PAH pathologies were examined using miR-181a/b overexpression lentivirus. A luciferase reporter analysis was performed to measure the relationships between miR-181a/b and endocan. Additionally, primary rat pulmonary arterial endothelial cells (rPAECs) treated with tumor necrosis factor-α (TNF-α) were employed to further validate the regulatory mechanism of miR-181a/b in vitro. Our results showed that miR-181a/b expression was reduced in PAH, and its upregulation significantly attenuated the short survival period, right ventricular systolic pressure and mean pulmonary artery pressure increments, right ventricular remodeling, and lung injury. Furthermore, the increase of intercellular cell adhesion molecule-1 (ICAM1) and vascular cell adhesion molecule-1 (VCAM1) in PAH rats was inhibited by miR-181a/b overexpression. Similarly, our in vitro results showed that inducing miR-181a/b suppressed TNF-α-stimulated increase of ICAM1 and VCAM1 in rPAECs. Importantly, the increased expression of endocan in PAH model or TNF-α-treated rPAECs was restored by miR-181a/b upregulation. Further analysis validated the direct targeting relationships between miR-181a/b and endocan. Collectively, this study suggests that miR-181a/b targets endocan to ameliorate PAH symptoms by inhibiting inflammatory states, shedding new lights on the prevention and treatment of PAH.