The discovery of a novel eight-mRNA-lncRNA signature predicting survival of hepatocellular carcinoma patients

Increasing evidence indicates that the expressions of messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) undergo a frequent and aberrant change in carcinogenesis and cancer development. But some research was carried out on mRNA-lncRNA signatures for prediction of hepatocellular carcinoma (HCC) prognosis. We aimed to establish an mRNA-lncRNA signature to improve the ability to predict HCC patients’ survival. The subjects from the cancer genome atlas (TCGA) data set were randomly divided into two parts: training data set (n = 246) and testing data set (n = 124). Using computational methods, we selected eight gene signatures (five mRNAs and three lncRNAs) to generate the risk score model, which were significantly correlated with overall survival of patients with HCC in both training and testing data set. The signature had the ability to classify the patients in training data set into a high-risk group and low-risk group with significantly different overall survival (hazard ratio = 4.157, 95% confidence interval = 2.648-6.526, P < 0.001). The prognostic value was further validated in testing data set and the entire data set. Further analysis revealed that this signature was independent of tumor stage. In addition, Gene Set Enrichment Analysis suggested that high risk score group was associated with cell proliferation and division related pathways. Finally, we developed a well-performed nomogram integrating the prognostic signature and other clinical information to predict 3- and 5-year overall survival. In conclusion, the prognostic mRNAs and lncRNAs identified in our study indicate their potential role in HCC biogenesis. The risk score model based on the mRNA-lncRNA may be an efficient classification tool to evaluate the prognosis of patients’ with HCC.     

Genome‐wide methylomic analyses identify prognostic epigenetic signature in lower grade glioma

Glioma is the most malignant and aggressive type of brain tumour with high heterogeneity and mortality. Although some clinicopathological factors have been identified as prognostic biomarkers, the individual variants and risk stratification in patients with lower grade glioma (LGG) have not been fully elucidated. The primary aim of this study was to identify an efficient DNA methylation combination biomarker for risk stratification and prognosis in LGG. We conducted a retrospective cohort study by analysing whole genome DNA methylation data of 646 patients with LGG from the TCGA and GEO database. Cox proportional hazard analysis was carried out to screen and construct biomarker model that predicted overall survival (OS). The Kaplan‐Meier survival curves and time‐dependent ROC were constructed to prove the efficiency of the signature. Then, another independent cohort was used to further validate the finding. A two‐CpG site DNA methylation signature was identified by multivariate Cox proportional hazard analysis. Further analysis indicated that the signature was an independent survival predictor from other clinical factors and exhibited higher predictive accuracy compared with known biomarkers. This signature was significantly correlated with immune‐checkpoint blockade, immunotherapy‐related signatures and ferroptosis regulator genes. The expression pattern and functional analysis showed that these two genes corresponding with two methylation sites contained in the model were correlated with immune infiltration level, and involved in MAPK and Rap1 signalling pathway. The signature may contribute to improve the risk stratification of patients and provide a more accurate assessment for precision medicine in the clinic.     

采用数据挖掘技术对湖北省人类狂犬病开展生物信息学研究

狂犬病作为一种急性人畜共患疾病,其致死率接近100％.而湖北省作为我国中部狂犬病高发地之一,开展该省狂犬病的流行病学调查不仅有助于了解我国当前面临的狂犬病疫情风险,还能为当地及全国狂犬病防控工作提供有效的参考意见.首先运用描述性分析发现湖北省狂犬病发病人数随时间呈下降趋势,而狂犬病暴露人数则呈上升趋势,且近年来中西部地...     

Assessing vegetation recovery from energy development using a dynamic reference approach

Ecologically relevant references are useful for evaluating ecosystem recovery, but references that are temporally static may be less useful when environmental conditions and disturbances are spatially and temporally heterogeneous. This challenge is particularly acute for ecosystems dominated by sagebrush (Artemisia spp.), where communities may require decades to recover from disturbance. We demonstrated application of a dynamic reference approach to studying sagebrush recovery using three decades of sagebrush cover estimates from remote sensing (1985–2018). We modelled recovery on former oil and gas well pads (n = 1200) across southwestern Wyoming, USA, relative to paired references identified by the Disturbance Automated Reference Toolset. We also used quantile regression to account for unmodelled heterogeneity in recovery, and projected recovery from similar disturbance across the landscape. Responses to weather and site‐level factors often differed among quantiles, and sagebrush recovery on former well pads increased more when paired reference sites had greater sagebrush cover. Little (<5%) of the landscape was projected to recover within 100 years for low to mid quantiles, and recovery often occurred at higher elevations with cool and moist annual conditions. Conversely, 48%–78% of the landscape recovered quickly (within 25 years) for high quantiles of sagebrush cover. Our study demonstrates advantages of using dynamic reference sites when studying vegetation recovery, as well as how additional inferences obtained from quantile regression can inform management.     

On the use of the variogram in checking for independence in spatial data

The variogram is a standard tool in the analysis of spatial data, and its shape provides useful information on the form of spatial correlation that may be present. However, it is also useful to be able to assess the evidence for the presence of any spatial correlation. A method of doing this, based on an assessment of whether the true function underlying the variogram is constant, is proposed. Nonparametric smoothing of the squared differences of the observed variables, on a suitably transformed scale, is used to estimate variogram shape. A statistic based on a ratio of quadratic forms is proposed and the test is constructed by investigating the distributional properties of this statistic under the assumption of an independent Gaussian process. The power of the test is investigated. Reference bands are proposed as a graphical follow-up. An example is discussed.     

Estimation of population attributable fractions from fitted incidence ratios and exposure survey data, with an application to electromagnetic fields and childhood leukemia

Standard presentations of epidemiological results focus on incidence-ratio estimates derived from regression models fit to specialized study data. These data are often highly nonrepresentative of populations for which public-health impacts must be evaluated. Basic methods are provided for interval estimation of attributable fractions from model-based incidence-ratio estimates combined with independent survey estimates of the exposure distribution in the target population of interest. These methods are illustrated in estimation of the potential impact of magnetic-field exposures on childhood leukemia in the United States, based on pooled data from 11 case-control studies and a U.S. sample survey of magnetic-field exposures.     

A semiparametric estimate of treatment effects with censored data

A semiparametric estimate of an average regression effect with right-censored failure time data has recently been proposed under the Cox-type model where the regression effect beta(t) is allowed to vary with time. In this article, we derive a simple algebraic relationship between this average regression effect and a measurement of group differences in k-sample transformation models when the random error belongs to the G(rho) family of Harrington and Fleming (1982, Biometrika 69, 553-566), the latter being equivalent to the conditional regression effect in a gamma frailty model. The models considered here are suitable for the attenuating hazard ratios that often arise in practice. The results reveal an interesting connection among the above three classes of models as alternatives to the proportional hazards assumption and add to our understanding of the behavior of the partial likelihood estimate under nonproportional hazards. The algebraic relationship provides a simple estimator under the transformation model. We develop a variance estimator based on the empirical influence function that is much easier to compute than the previously suggested resampling methods. When there is truncation in the right tail of the failure times, we propose a method of bias correction to improve the coverage properties of the confidence intervals. The estimate, its estimated variance, and the bias correction term can all be calculated with minor modifications to standard software for proportional hazards regression.     

Putting background information about relative risks into conjugate prior distributions

In Bayesian and empirical Bayes analyses of epidemiologic data, the most easily implemented prior specifications use a multivariate normal distribution for the log relative risks or a conjugate distribution for the discrete response vector. This article describes problems in translating background information about relative risks into conjugate priors and a solution. Traditionally, conjugate priors have been specified through flattening constants, an approach that leads to conflicts with the true prior covariance structure for the log relative risks. One can, however, derive a conjugate prior consistent with that structure by using a data-augmentation approximation to the true log relative-risk prior, although a rescaling step is needed to ensure the accuracy of the approximation. These points are illustrated with a logistic regression analysis of neonatal-death risk.     

Synthesis of evidence from epidemiological studies with interval-censored exposure due to grouping

We describe a method for assessing dose-response effects from a series of case-control and cohort studies in which the exposure information is interval censored. The interval censoring of the exposure variable is dealt with through the use of retrospective models in which the exposure is treated as a multinomial response and disease status as a binary covariate. Polychotomous logistic regression models are adopted in which the dose-response relationship between exposure and disease may be modeled in a discrete or continuous fashion. Partial conditioning is possible to eliminate some of the nuisance parameters. The methods are applied to the motivating study of the relationship between chorionic villus sampling and the occurrence of terminal transverse limb reduction.