CBCRisk-Black: a personalized contralateral breast cancer risk prediction model for black women.

PURPOSE: Black breast cancer (BC) survivors have a higher risk of developing contralateral breast cancer (CBC) than Whites. Existing CBC risk prediction tools are developed based on mostly White women. To address this racial disparity, it is crucial to develop tools tailored for Black women to help them inform about their actual risk of CBC. METHODS: We propose an absolute risk prediction model, CBCRisk-Black, specifically for Black BC patients. It uses data on Black women from two sources: Breast Cancer Surveillance Consortium (BCSC) and Surveillance, Epidemiology, and End Results (SEER). First, a matched lasso logistic regression model for estimating relative risks (RR) is developed. Then, it is combined with relevant hazard rates and attributable risks to obtain absolute risks. Six-fold cross-validation is used to internally validate CBCRisk-Black. We also compare CBCRisk-Black with CBCRisk, an existing CBC risk prediction model. RESULTS: The RR model uses data from BCSC on 744 Black women (186 cases). CBCRisk-Black has four risk factors (RR compared to baseline): breast density (2.13 for heterogeneous/extremely dense), family history of BC (2.28 for yes), first BC tumor size (2.14 for T3/T4, 1.56 for TIS), and age at first diagnosis of BC (1.41 for < 40). The area under the receiver operating characteristic curve (AUC) for 3- and 5-year predictions are 0.72 and 0.65 for CBCRisk-Black while those are 0.65 and 0.60 for CBCRisk. CONCLUSION: CBCRisk-Black may serve as a useful tool to clinicians in counseling Black BC patients by providing a more accurate and personalized CBC risk estimate.

A segmented measurement error model for modeling and analysis of method comparison data.

Method comparison studies are concerned with estimating relationship between two clinical measurement methods. The methods often exhibit a structural change in the relationship over the measurement range. Ignoring this change would lead to an inaccurate estimate of the relationship. Motivated by a study of two digoxin assays where such a change occurs, this article develops a statistical methodology for appropriately analyzing such studies. Specifically, it proposes a segmented extension of the classical measurement error model to allow a piecewise linear relationship between the methods. The changepoint at which the transition takes place is treated as an unknown parameter in the model. An expectation-maximization-type algorithm is developed to fit the model and appropriate extensions of the existing measures are proposed for segment-specific evaluation of similarity and agreement. Bootstrapping and large-sample theory of maximum likelihood estimators are employed to perform the relevant inferences. The proposed methodology is evaluated by simulation and is illustrated by analyzing the digoxin data.

A Bayesian Hierarchical Framework for Pathway Analysis in Genome-Wide Association Studies.

BACKGROUND: Pathway analysis allows joint consideration of multiple SNPs belonging to multiple genes, which in turn belong to a biologically defined pathway. This type of analysis is usually more powerful than single-SNP analyses for detecting joint effects of variants in a pathway. METHODS: We develop a Bayesian hierarchical model by fully modeling the 3-level hierarchy, namely, SNP-gene-pathway that is naturally inherent in the structure of the pathways, unlike the currently used ad hoc ways of combining such information. We model the effects at each level conditional on the effects of the levels preceding them within the generalized linear model framework. To deal with the high dimensionality, we regularize the regression coefficients through an appropriate choice of priors. The model is fit using a combination of iteratively weighted least squares and expectation-maximization algorithms to estimate the posterior modes and their standard errors. A normal approximation is used for inference. RESULTS: We conduct simulations to study the proposed method and find that our method has higher power than some standard approaches in several settings for identifying pathways with multiple modest-sized variants. We illustrate the method by analyzing data from two genome-wide association studies on breast and renal cancers. CONCLUSION: Our method can be helpful in detecting pathway association.

Dose-dependent effect of mammographic breast density on the risk of contralateral breast cancer.

PURPOSE: Increased mammographic breast density is a significant risk factor for breast cancer. It is not clear if it is also a risk factor for the development of contralateral breast cancer. METHODS: The data were obtained from Breast Cancer Surveillance Consortium and included women diagnosed with invasive breast cancer or ductal carcinoma in situ between ages 18 and 88 and years 1995 and 2009. Each case of contralateral breast cancer was matched with three controls based on year of first breast cancer diagnosis, race, and length of follow-up. A total of 847 cases and 2541 controls were included. The risk factors included in the study were mammographic breast density, age of first breast cancer diagnosis, family history of breast cancer, anti-estrogen treatment, hormone replacement therapy, menopausal status, and estrogen receptor status, all from the time of first breast cancer diagnosis. Both univariate analysis and multivariate conditional logistic regression analysis were performed. RESULTS: In the final multivariate model, breast density, family history of breast cancer, and anti-estrogen treatment remained significant with p values less than 0.01. Increasing breast density had a dose-dependent effect on the risk of contralateral breast cancer. Relative to 'almost entirely fat' category of breast density, the adjusted odds ratios (and p values) in the multivariate analysis for 'scattered density,' 'heterogeneously dense,' and 'extremely dense' categories were 1.65 (0.036), 2.10 (0.002), and 2.32 (0.001), respectively. CONCLUSION: Breast density is an independent and significant risk factor for development of contralateral breast cancer. This risk factor should contribute to clinical decision making.

A model for individualized risk prediction of contralateral breast cancer.

PURPOSE: Patients diagnosed with invasive breast cancer (BC) or ductal carcinoma in situ are increasingly choosing to undergo contralateral prophylactic mastectomy (CPM) to reduce their risk of contralateral BC (CBC). This is a particularly disturbing trend as a large proportion of these CPMs are believed to be medically unnecessary. Many BC patients tend to substantially overestimate their CBC risk. Thus, there is a pressing need to educate patients effectively on their CBC risk. We develop a CBC risk prediction model to aid physicians in this task. METHODS: We used data from two sources: Breast Cancer Surveillance Consortium and Surveillance, Epidemiology, and End Results to build the model. The model building steps are similar to those used in developing the BC risk assessment tool (popularly known as Gail model) for counseling women on their BC risk. Our model, named CBCRisk, is exclusively designed for counseling women diagnosed with unilateral BC on the risk of developing CBC. RESULTS: We identified eight factors to be significantly associated with CBC-age at first BC diagnosis, anti-estrogen therapy, family history of BC, high-risk pre-neoplasia status, estrogen receptor status, breast density, type of first BC, and age at first birth. Combining the relative risk estimates with the relevant hazard rates, CBCRisk projects absolute risk of developing CBC over a given period. CONCLUSIONS: By providing individualized CBC risk estimates, CBCRisk may help in counseling of BC patients. In turn, this may potentially help alleviate the rate of medically unnecessary CPMs.

Semiparametric modeling and analysis of longitudinal method comparison data.

Studies comparing two or more methods of measuring a continuous variable are routinely conducted in biomedical disciplines with the primary goal of measuring agreement between the methods. Often, the data are collected by following a cohort of subjects over a period of time. This gives rise to longitudinal method comparison data where there is one observation trajectory for each method on every subject. It is not required that observations from all methods be available at each observation time. The multiple trajectories on the same subjects are dependent. We propose modeling the trajectories nonparametrically through penalized regression splines within the framework of mixed-effects models. The model also uses random effects of subjects and their interactions to capture dependence in observations from the same subjects. It additionally allows the within-subject errors of each method to be correlated. It is fit using the method of maximum likelihood. Agreement between the methods is evaluated by performing inference on measures of agreement, such as concordance correlation coefficient and total deviation index, which are functions of parameters of the assumed model. Simulations indicate that the proposed methodology performs reasonably well for 30 or more subjects. Its application is illustrated by analyzing a dataset of percentage body fat measurements. Copyright © 2017 John Wiley & Sons, Ltd.

Tolerance bands for functional data.

Often the object of inference in biomedical applications is a range that brackets a given fraction of individual observations in a population. A classical estimate of this range for univariate measurements is a "tolerance interval." This article develops its natural extension for functional measurements, a "tolerance band," and proposes a methodology for constructing its pointwise and simultaneous versions that incorporates both sparse and dense functional data. Assuming that the measurements are observed with noise, the methodology uses functional principal component analysis in a mixed model framework to represent the measurements and employs bootstrapping to approximate the tolerance factors needed for the bands. The proposed bands also account for uncertainty in the principal components decomposition. Simulations show that the methodology has, generally, acceptable performance unless the data are quite sparse and unbalanced, in which case the bands may be somewhat liberal. The methodology is illustrated using two real datasets, a sparse dataset involving CD4 cell counts and a dense dataset involving core body temperatures.

A heteroscedastic measurement error model for method comparison data with replicate measurements.

Measurement error models offer a flexible framework for modeling data collected in studies comparing methods of quantitative measurement. These models generally make two simplifying assumptions: (i) the measurements are homoscedastic, and (ii) the unobservable true values of the methods are linearly related. One or both of these assumptions may be violated in practice. In particular, error variabilities of the methods may depend on the magnitude of measurement, or the true values may be nonlinearly related. Data with these features call for a heteroscedastic measurement error model that allows nonlinear relationships in the true values. We present such a model for the case when the measurements are replicated, discuss its fitting, and explain how to evaluate similarity of measurement methods and agreement between them, which are two common goals of data analysis, under this model. Model fitting involves dealing with lack of a closed form for the likelihood function. We consider estimation methods that approximate either the likelihood or the model to yield approximate maximum likelihood estimates. The fitting methods are evaluated in a simulation study. The proposed methodology is used to analyze a cholesterol dataset.

Measuring agreement in method comparison studies with heteroscedastic measurements.

We propose a methodology for evaluation of agreement between two methods of measuring a continuous variable whose variability changes with magnitude. This problem routinely arises in method comparison studies that are common in health-related disciplines. Assuming replicated measurements, we first model the data using a heteroscedastic mixed-effects model, wherein a suitably defined true measurement serves as the variance covariate. Fitting this model poses some computational difficulties as the likelihood function is not available in a closed form. We deal with this issue by suggesting four estimation methods to obtain approximate maximum likelihood estimates. Two of these methods are based on numerical approximation of the likelihood, and the other two are based on approximation of the model. Next, we extend the existing agreement evaluation methodology designed for homoscedastic data to work under the proposed heteroscedastic model. This methodology can be used with any scalar measure of agreement. Simulations show that the suggested inference procedures generally work well for moderately large samples. They are illustrated by analyzing a data set of cholesterol measurements.

L-statistics for Repeated Measurements Data With Application to Trimmed Means, Quantiles and Tolerance Intervals.

The L-statistics form an important class of estimators in nonparametric statistics. Its members include trimmed means and sample quantiles and functions thereof. This article is devoted to theory and applications of L-statistics for repeated measurements data, wherein the measurements on the same subject are dependent and the measurements from different subjects are independent. This article has three main goals: (a) Show that the L-statistics are asymptotically normal for repeated measurements data. (b) Present three statistical applications of this result, namely, location estimation using trimmed means, quantile estimation and construction of tolerance intervals. (c) Obtain a Bahadur representation for sample quantiles. These results are generalizations of similar results for independently and identically distributed data. The practical usefulness of these results is illustrated by analyzing a real data set involving measurement of systolic blood pressure. The properties of the proposed point and interval estimators are examined via simulation.

Breastfed infants metabolize perchlorate.

Bifidobacteria are the dominant intestinal bacteria in breastfed infants. It is known that they can reduce nitrate. Although no direct experiments have been conducted until now, inferred pathways for Bifidobacterium bifidum include perchlorate reduction via perchlorate reductase. We show that when commercially available strains of bifidobacteria are cultured in milk, spiked with perchlorate, perchlorate is consumed. We studied 13 breastfed infant-mother pairs who provided 43 milk samples and 39 infant urine samples, and 5 formula-fed infant-mother pairs who provided 21 formula samples and 21 infant urine samples. Using iodine as a conservative tracer, we determined the average urinary iodine (UI) to milk iodine (MI) concentration ratio to be 2.87 for the breastfed infants. For the same samples, the corresponding perchlorate concentration ratio was 1.37 (difference significant, p < 0.001), indicating that perchlorate is lost. For the formula fed infant group the same ratios were 1.20 and 1.58; the difference was not significant (p = 0.68). However, the small number of subjects in the latter group makes it more difficult to conclude definitively whether perchlorate reduction does or does not occur.

Joint risk prediction for hazardous use of alcohol, cannabis, and tobacco among adolescents: A preliminary study using statistical and machine learning.

For some, substance use during adolescence may be a stepping stone on the way to substance use disorders in adulthood. Risk prediction models may help identify adolescent users at elevated risk for hazardous substance use. This preliminary analysis used cross-sectional data (n = 270, ages 13-18) from the baseline dataset of a randomized controlled trial intervening with adolescent alcohol and/or cannabis use. Models were developed for jointly predicting quantitative scores on three measures of hazardous substance use (Rutgers Alcohol Problems Index, Adolescent Cannabis Problem Questionnaire, and Hooked on Nicotine Checklist) based on personal risk factors using two statistical and machine learning methods: multivariate covariance generalized linear models (MCGLM) and penalized multivariate regression with a lasso penalty. The predictive accuracy of a model was evaluated using root mean squared error computed via leave-one-out cross-validation. The final proposed model was an MCGLM model. It has eleven risk factors: age, early life stress, age of first tobacco use, age of first cannabis use, lifetime use of other substances, age of first use of other substances, maternal education, parental attachment, family cigarette use, family history of hazardous alcohol use, and family history of hazardous cannabis use. Different subsets of these risk factors feature in the three outcome-specific components of this joint model. The quantitative risk estimate provided by the proposed model may help identify adolescent substance users of cannabis, alcohol, and tobacco who may be at an elevated risk of developing hazardous substance use.

Krüppel-like Factor 4 activates HBG gene expression in primary erythroid cells.

The SP1/Krüppel-like Factor (SP1/KLF) family of transcription factors plays a role in diverse cellular processes, including proliferation, differentiation and control of gene transcription. The discovery of KLF1 (EKLF), a key regulator of HBB (ß-globin) gene expression, expanded our understanding of the role of KLFs in erythropoiesis. In this study, we investigated a mechanism of HBG (γ-globin) regulation by KLF4. siRNA-mediated gene silencing and enforced expression of KLF4 in K562 cells substantiated the ability of KLF4 to positively regulate endogenous HBG gene transcription. The physiological significance of this finding was confirmed in primary erythroid cells, where KLF4 knockdown at day 11 significantly attenuated HBG mRNA levels and enforced expression at day 28 stimulated the silenced HBG genes. In vitro binding characterization using the γ-CACCC and ß-CACCC probes demonstrated KLF4 preferentially binds the endogenous γ-CACCC, while CREB binding protein (CREBBP) binding was not selective. Co-immunoprecipitation studies confirmed protein-protein interaction between KLF4 and CREBBP. Furthermore, sequential chromatin immunoprecipitation assays showed co-localization of both factors in the γ-CACCC region. Subsequent luciferase reporter studies demonstrated that KLF4 trans-activated HBG promoter activity and that CREBBP enforced expression resulted in gene repression. Our data supports a model of antagonistic interaction of KLF4/CREBBP trans-factors in HBG regulation.

A preliminary risk prediction model for cannabis use disorder.

The ongoing trend toward legalization of cannabis for medicinal/recreational purposes is expected to increase the prevalence of cannabis use disorder (CUD). Thus, it is imperative to be able to predict the quantitative risk of developing CUD for a cannabis user based on their personal risk factors. Yet no such model currently exists. In this study, we perform preliminary analysis toward building such a model. The data come from n = 94 regular cannabis users recruited from Albuquerque, New Mexico during 2007-2010. As the data are cross-sectional, we only consider risk factors that remain relatively stable over time. We apply statistical and machine learning classification techniques that allow n to be small relative to the number of predictors. We use predictive accuracy estimated using leave-one-out-cross-validation to evaluate model performance. The final model is a LASSO logistic regression model consisting of the following seven risk factors: age; level of enjoyment from initial cigarette smoking; total score on Impulsive Sensation-Seeking Scale questionnaire; score on cognitive instability factor of Barratt Impulsivity Scale questionnaire; and scores on neuroticism, openness, and conscientiousness personality traits of Neuroticism, Extraversion, and Openness inventory. This model has an overall accuracy of 0.66 and the area under its receiver operating characteristic curve is 0.65. In summary, a preliminary relative risk model for predicting the quantitative risk of CUD is developed. It can be employed to identify users at high risk of CUD who may be provided with early intervention.

Monocyte protein signatures of disease severity in sickle cell anemia.

Using two-dimensional difference gel electrophoresis (2D DIGE) we have analyzed monocytes derived from 10 sickle cell disease patients (5 males and 5 females ages 12-18) to generate hypotheses regarding signature proteins that appear most positively and negatively correlated with vasoocclusive event rate. Signature proteins have been identified by tandem mass spectrometry. Based on the limited number of samples analyzed, the most negatively correlated proteins related to crises rate were transketolase and coronin in the membrane fraction and heat shock 70 kDa protein cognate 4, and adenylate kinase isoenzyme 2, mitochondrial found in the cytosolic fraction. The protein spots that were most positively correlated with crisis rate in the cytoplasmic fraction were far upstream element-binding protein and Alpha actinin 1 or Alpha actinin 4. Utilizing StepSIM analysis, vinculin was able to classify all samples from the combined set and the membrane-only set, and cytosolic leukotriene A-4 hydrolase and phosphoglycerate kinase were also identified as important indicators for differentiating between low and high vasoocclusive event rates.

Pharmaco-proteomic study of hydroxyurea-induced modifications in the sickle red blood cell membrane proteome.

Hydroxyurea (HU) is an effective oral drug for the management of homozygous sickle cell anemia (SS) in part because it increases fetal hemoglobin (HbF) levels within sickle red blood cells (RBCs) and thus reduces sickling. However, results from the Multicenter Study of HU suggested that clinical symptoms often improved before a significant increase in HbF levels occurred. This indicated that HU may be acting through the modification of additional cellular mechanisms that are yet to be identified. Hence, in this study, we focused on the analysis of the sickle RBC membrane proteome +/- HU treatment. 2D-DIGE (Two Dimensional Difference In-Gel Electrophoresis) technology and tandem mass spectrometry has been used to determine quantitative differences between sickle cell membrane proteins in the presence and absence of a clinically relevant concentration of HU. In vitro protein profiling of 13 sickle RBC membrane samples +/- 50 muM HU identified 10 statistically significant protein spots. Of these, the most remarkable class of proteins to show a statistically significant increase was the anti-oxidant enzymes-catalase, thioredoxin peroxidase and biliver-din reductase and the chaperonin containing TCP1 complex assisting in the folding of RBC cytoskeletal proteins. Interestingly, catalase immunoblots showed an increase in the acidic forms of the enzyme within sickle RBC membranes on incubation with 50 muM HU. We further identified this modification in catalase to be phosphorylation and demonstrated that HU exposed SS RBC membranes showed a 2-fold increase in tyrosine phosphorylation of catalase as compared to counterparts not exposed to HU. These results present an attractive model for HU-induced post-translational modification and potential activation of catalase in mature sickle RBCs. These findings also identify protein targets of HU other than fetal hemoglobin and enhance the understanding of the drug mechanism.

Variation in the monocyte proteome.

We have determined the variability of the monocyte proteome and identified those proteins that demonstrate the greatest variation in the general control population. Monocytes were isolated from 18 healthy (9 male and 9 female) donors ages 18-50 and with no known genetic or blood disorder. A combination of Ficoll-Paque PLUS density centrifugation of cells found in the buffy coat and positive selection with monoclonal antibodies against CD14, coupled to magnetic beads, led to >98% purity of monocytes. A 100,000 g microsomal membrane fraction or 100,000 g supernatant fraction from a control subject was compared to the equivalent fractions from a distinct control subject by two-dimensional differential gel electrophoresis (2D DIGE). Those protein spots that demonstrated Cy3-/Cy5- ratios greater than 2.5-fold in at least one experiment were selected for further statistical analysis. We determined variability for 31 cytosolic and 12 membrane protein spots. Proteins have been identified for 27 of the cytosolic protein spots and 9 of the microsomal membrane protein spots by in-gel digestion with trypsin followed by reverse-phase high-performance liquid chromatography in line with tandem mass spectrometry. We identified 24 distinct monocyte proteins that demonstrated the greatest variability in this general control population. The proteins demonstrating the greatest variance in the cytosolic fraction were enolase-1 and WD (tryptophan-aspartate) repeat-containing protein 1, and in the membrane fraction they were lamin B1 and L-plastin. This study demonstrates the importance of considering variance in the control population when performing future protein profiling comparisons of monocytes derived from disease versus control populations.

Protein profiling of sickle cell versus control RBC core membrane skeletons by ICAT technology and tandem mass spectrometry.

A proteomic approach using a cleavable ICAT reagent and nano-LC ESI tandem mass spectrometry was used to perform protein profiling of core RBC membrane skeleton proteins between sickle cell patients (SS) and controls (AA), and determine the efficacy of this technology. The data was validated through Peptide/Protein Prophet and protein ratios were calculated through ASAPratio. Through an ANOVA test, it was determined that there is no significant difference in the mean ratios from control populations (AA1/AA2) and sickle cell versus control populations (AA/SS). The mean ratios were not significantly different from 1.0 in either comparison for the core skeleton proteins (alpha spectrin, beta spectrin, band 4.1 and actin). On the natural-log scale, the variation (standard deviation) of the method was determined to be 14.1% and the variation contributed by the samples was 13.8% which together give a total variation of 19.7% in the ratios.

Characterization of Cement Particles Found in Peri-implantitis-Affected Human Biopsy Specimens.

PURPOSE: Peri-implantitis is a disease characterized by soft tissue inflammation and continued loss of supporting bone, which can result in implant failure. Peri-implantitis is a multifactorial disease, and one of its triggering factors may be the presence of excess cement in the soft tissues surrounding an implant. This descriptive study evaluated the composition of foreign particles from 36 human biopsy specimens with 19 specimens selected for analysis. The biopsy specimens were obtained from soft tissues affected by peri-implantitis around cement-retained implant crowns and compared with the elemental composition of commercial luting cement. MATERIALS AND METHODS: Nineteen biopsy specimens were chosen for the comparison, and five test cements (TempBond, Telio, Premier Implant Cement, Intermediate Restorative Material, and Relyx) were analyzed using scanning electron microscopy equipped with energy dispersive x-ray spectroscopy. This enabled the identification of the chemical composition of foreign particles embedded in the tissue specimens and the composition of the five cements. Statistical analysis was conducted using classification trees to pair the particles present in each specimen with the known cements. RESULTS: The particles in each biopsy specimen could be associated with one of the commercial cements with a level of probability ranging between .79 and 1. TempBond particles were found in one biopsy specimen, Telio particles in seven, Premier Implant Cement particles in four, Relyx particles in four, and Intermediate Restorative Material particles in three. CONCLUSION: Particles found in human soft tissue biopsy specimens around implants affected by peri-implant disease were associated with five commercially available dental cements.

Early dysregulation of the mitochondrial proteome in a mouse model of Alzheimer's disease.

Mitochondrial structural and functional alterations appear to play to an important role in the pathogenesis of Alzheimer's disease (AD). In the present study, we used a quantitative comparative proteomic profiling approach to analyze changes in the mitochondrial proteome in AD. A triple transgenic mouse model of AD (3xTg-AD) which harbors mutations in three human transgenes, APP(Swe), PS1(M146V) and Tau(P301L), was used in these experiments. Quantitative differences in the mitochondrial proteome between the cerebral cortices of 6-month-old male 3xTg-AD and non-transgenic mice were determined by using two-dimensional difference gel electrophoresis (2D-DIGE) and tandem mass spectrometry. We identified 23 different proteins whose expression levels differed significantly between triple transgenic and non-transgenic mitochondria. Both down-regulated and up-regulated mitochondrial proteins were observed in transgenic AD cortices. Proteins which were dysregulated in 3xTg-AD cortices functioned in a wide variety of metabolic pathways, including the citric acid cycle, oxidative phosphorylation, pyruvate metabolism, glycolysis, oxidative stress, fatty acid oxidation, ketone body metabolism, ion transport, apoptosis, and mitochondrial protein synthesis. These alterations in the mitochondrial proteome of the cerebral cortices of triple transgenic AD mice occurred before the development of significant amyloid plaque and neurofibrillary tangles, indicating that mitochondrial dysregulation is an early event in AD.