Cubic regression-based degree of correction predicts the performance of whole bisulfitome amplified DNA methylation analysis

Epigenetic mechanisms, including DNA methylation, are important determinants in development and disease. There is a need for technologies capable of detecting small variations in methylation levels in an accurate and reproducible manner, even if only limited amounts of DNA are available (which is the case in many studies in humans). Quantitative methylation analysis of minute DNA amounts after whole bisulfitome amplification (qMAMBA) has been proposed as an alternative, but this technique has not been adequately standardized and no comparative study against conventional methods has been performed, that includes a wide range of methylation percentages and different target assays. We designed an experiment to compare the performance of qMAMBA and bisulfite-treated genomic (non-amplified) DNA pyrosequencing. Reactions were performed in duplicate for each technique in eight different target genes, using nine artificially constructed DNA samples with methylation levels ranging between 0% and 100% with intervals of 12.5%. Cubic polynomial curves were plotted from the experimental results and the real methylation values and the resulting equation was used to estimate new corrected data points. The use of the cubic regression-based correction benefits the accuracy and the power of discrimination in methylation studies. Additionally, dispersion of the new estimated data around a y = x line (R²) served to fix a cutoff that can discriminate, with a single 9-point curve experiment, whether whole bisulfitome amplification and subsequent qMAMBA can produce accurate methylation results. Finally, even with an optimized reagent kit, DNA samples subjected to whole bisulfitome amplification enhance the preferential amplification of unmethylated alleles, and subtle changes in methylation levels cannot be detected confidently.     

Identification of novel prognostic alternative splicing signature in papillary renal cell carcinoma

Papillary renal cell carcinoma (pRCC) is a heterogeneous disease containing multifocal or solitary tumors with an aggressive phenotype. Increasing evidence has indicated the involvement of aberrant splicing variants in renal cell cancer, while systematic profiling of aberrant alternative splicing (AS) in pRCC was lacking and largely unknown. In the current study, comprehensive profiling of AS events were performed based on the integration of pRCC cohort from the Cancer Genome Atlas database and SpliceSeq software. With rigorous screening and univariate Cox analysis, a total of 2077 prognoses AS events from 1642 parent genes were identified. Then, stepwise least absolute shrinkage and selection operator method-penalized Cox regression analyses with 10-fold cross-validation followed by multivariate Cox regression were used to construct the prognostic AS signatures within each AS type. And a final 21 AS event-based signature was proposed which showed potent prognostic capability in stratifying patients into low- and high-risk subgroups (P < .0001). Furthermore, time-dependent receiver operating characteristics curves confirmed that the final AS signature was effective and robust in predicting overall survival for pRCC patients with the area under the curve above 0.9 from 1 to 5 years. In addition, splicing correlation network was built to uncover the potential regulatory pattern among prognostic splicing factors and candidate AS events. Besides, gene set enrichment analysis revealed the involvement of these candidates AS events in tumor-related pathways including extracellular matrix organization, oxidative phosphorylation, and P53 signaling pathways. Taken together, our results could contribute to elucidating the underlying mechanism of AS in the oncogenesis process and broaden the novel field of prognostic and clinical application of molecule-targeted approaches in pRCC.     

Physical maturity in common bottlenose dolphins (Tursiops truncatus) from Sarasota Bay,Florida

Cetacean physical maturity is defined by growth cessation and complete fusion of epiphyses to vertebral bodies indicated by invisible sutures. Many studies have shown epiphyseal fusion is highly variable among individuals. In-depth examinations into fusion variability are lacking. We analyzed vertebrae of 37 (n = 21 female, n = 16 male) stranded common bottlenose dolphins (Tursiops truncatus) from the well-studied Gulf of Mexico, Sarasota Bay community. For each specimen, vertebrae were examined by vertebral region for degree of fusion anteriorly and posteriorly of each centrum and categorized from unfused to fused in five degrees. An ordinal logistic regression was used to estimate degree of fusion probability for each epiphysis. The model had fixed effects for age, number of offspring, sex, sexual maturity, and a random effect for epiphysis. Results show that age/reproductive status significantly explains an individual's degree of fusion. Adult females with fewer calves had more fusion than those with more reproductive experience across multiple ages. Access to long-term observational and sample data on the dolphins residing in the area served by Mote Marine Laboratory's Stranding Investigations Program offers a unique opportunity to examine the relationship between energetic demands of reproduction (calcium production/reproductive output) versus preconceived definitions of physical maturity (skeletal fusion) more closely.     

Comparison of complex modeling strategies for prediction of a binary outcome based on a few,highly correlated predictors

Motivated by a clinical prediction problem, a simulation study was performed to compare different approaches for building risk prediction models. Robust prediction models for hospital survival in patients with acute heart failure were to be derived from three highly correlated blood parameters measured up to four times, with predictive ability having explicit priority over interpretability. Methods that relied only on the original predictors were compared with methods using an expanded predictor space including transformations and interactions. Predictors were simulated as transformations and combinations of multivariate normal variables which were fitted to the partly skewed and bimodally distributed original data in such a way that the simulated data mimicked the original covariate structure. Different penalized versions of logistic regression as well as random forests and generalized additive models were investigated using classical logistic regression as a benchmark. Their performance was assessed based on measures of predictive accuracy, model discrimination, and model calibration. Three different scenarios using different subsets of the original data with different numbers of observations and events per variable were investigated. In the investigated setting, where a risk prediction model should be based on a small set of highly correlated and interconnected predictors, Elastic Net and also Ridge logistic regression showed good performance compared to their competitors, while other methods did not lead to substantial improvements or even performed worse than standard logistic regression. Our work demonstrates how simulation studies that mimic relevant features of a specific data set can support the choice of a good modeling strategy.     

A prognostic eight-gene expression signature for patients with breast cancer receiving adjuvant chemotherapy

Breast cancer is a popularly diagnosed malignant tumor. Genomic profiling studies suggest that breast cancer is a disease with heterogeneity. Chemotherapy is one of the chief means to treat breast cancer, while its responses and clinical outcomes vary largely due to the conventional clinicopathological factors and inherent chemosensitivity of breast cancer. Using the least absolute shrinkage and selection operator (LASSO) Cox regression model, our study established a multi-mRNA-based signature model and constructed a relative nomogram in predicting distant-recurrence-free survival for patients receiving surgery and following chemotherapy. We constructed a signature of eight mRNAs (IPCEF1, SYNDIG1, TIGIT, SPESP1, C2CD4A, CLCA2, RLN2, and CCL19) with the LASSO model, which was employed to separate subjects into groups with high- and low-risk scores. Obvious differences of distant-recurrence-free survival were found between these two groups. This eight-mRNA-based signature was independently associated with the prognosis and had better prognostic value than classical clinicopathologic factors according to multivariate Cox regression results. Receiver operating characteristic results demonstrated excellent performance in diagnosing 3-year distant-recurrence by the eight-mRNA signature. A nomogram that combined both the eight-mRNA-based signature and clinicopathological risk factors was constructed. Comparing with an ideal model, the nomograms worked well both in the training and validation sets. Through the results that the eight-mRNA signature effectively classified patients into low- and high-risk of distant recurrence, we concluded that this eight-mRNA-based signature played a promising predictive role in prognosis and could be clinically applied in breast cancer patients receiving adjuvant chemotherapy.     

Structures composing protein domains

This review summarizes available data concerning intradomain structures (IS) such as functionally important amino acid residues, short linear motifs, conserved or disordered regions, peptide repeats, broadly occurring secondary structures or folds, etc. IS form structural features (units or elements) necessary for interactions with proteins or non-peptidic ligands, enzyme reactions and some structural properties of proteins. These features have often been related to a single structural level (e.g. primary structure) mostly requiring certain structural context of other levels (e.g. secondary structures or supersecondary folds) as follows also from some examples reported or demonstrated here. In addition, we deal with some functionally important dynamic properties of IS (e.g. flexibility and different forms of accessibility), and more special dynamic changes of IS during enzyme reactions and allosteric regulation. Selected notes concern also some experimental methods, still more necessary tools of bioinformatic processing and clinically interesting relationships.     

The telocentric tandem repeat at the p-arm is not conserved in Mus musculus subspecies

Mouse chromosomes, with the exception of the Y chromosome, are telocentric. The telomere at the p-arm is separated from the centromere by the tL1 sequence and TLC tandem repeats. A previous report showed that the TLC array was also conserved in other strains of the subgenus Mus. These results suggest that the TLC arrays promote the stable evolutionary maintenance of a telocentric karyotype in the subgenus Mus. In this study, we investigated the degree of conservation of TLC arrays among a variety of wild-derived inbred strains, all of which are descendants of wild mice captured in several areas of the world. Genomic PCR analysis indicates that the sequential order of telomere-tL1 is highly conserved in all strains, whereas tL1-TLC is not. Next, Southern blot analysis of DNAs isolated from a panel of mouse subspecies showed both Mus musculus domesticus and Mus musculus castaneus subspecies possess TLC arrays. Unexpectedly, this repeat appears to be lost in almost all Mus musculus musculus and Mus musculus molossinus subspecies, which show a clear geographic divide. These results indicate that either other unknown sequences were replaced by the TLC repeat or almost all M. m. musculus and M. m. molossinus subspecies do not have any sequence between the telomere and minor satellites. Our observation suggests that the TLC array might be evolutionarily unstable and not essential for murine chromosomal conformation. This is the first example of the subspecies-specific large genome alterations in mice.