Addressing age measurement errors in fish growth estimation from length-stratified samples.

Fish growth models are crucial for fisheries stock assessments and are commonly estimated using fish length-at-age data. This data is widely collected using length-stratified age sampling (LSAS), a cost-effective two-phase response-selective sampling method. The data may contain age measurement errors (MEs). We propose a methodology that accounts for both LSAS and age MEs to accurately estimate fish growth. The proposed methods use empirical proportion likelihood methodology for LSAS and the structural errors in variables methodology for age MEs. We provide a measure of uncertainty for parameter estimates and standardized residuals for model validation. To model the age distribution, we employ a continuation ratio-logit model that is consistent with the random nature of the true age distribution. We also apply a discretization approach for age and length distributions, which significantly improves computational efficiency and is consistent with the discrete age and length data typically encountered in practice. Our simulation study shows that neglecting age MEs can lead to significant bias in growth estimation, even with small but non-negligible age MEs. However, our new approach performs well regardless of the magnitude of age MEs and accurately estimates SEs of parameter estimators. Real data analysis demonstrates the effectiveness of the proposed model validation device. Computer codes to implement the methodology are provided.

Multifactor dimensionality reduction method identifies novel SNP interactions in the WNT protein interaction networks that are associated with recurrence risk in colorectal cancer.

Background: Interactions among genetic variants are rarely studied but may explain a part of the variability in patient outcomes. Objectives: In this study, we aimed to identify 1 to 3 way interactions among SNPs from five Wnt protein interaction networks that predict the 5-year recurrence risk in a cohort of stage I-III colorectal cancer patients. Methods: 423 patients recruited to the Newfoundland Familial Colorectal Cancer Registry were included. Five Wnt family member proteins (Wnt1, Wnt2, Wnt5a, Wnt5b, and Wnt11) were selected. The BioGRID database was used to identify the proteins interacting with each of these proteins. Genotypes of the SNPs located in the interaction network genes were retrieved from a genome-wide SNP genotype data previously obtained in the patient cohort. The GMDR 0.9 program was utilized to examine 1-, 2-, and 3-SNP interactions using a 5-fold cross validation step. Top GMDR 0.9 models were assessed by permutation testing and, if significant, prognostic associations were verified by multivariable logistic regression models. Results: GMDR 0.9 has identified novel 1, 2, and 3-way SNP interactions associated with 5-year recurrence risk in colorectal cancer. Nine of these interactions were multi loci interactions (2-way or 3-way). Identified interaction models were able to distinguish patients based on their 5-year recurrence-free status in multivariable regression models. The significance of interactions was the highest in the 3-SNP models. Several of the identified SNPs were eQTLs, indicating potential biological roles of the genes they were associated with in colorectal cancer recurrence. Conclusions: We identified novel interacting genetic variants that associate with 5-year recurrence risk in colorectal cancer. A significant portion of the genes identified were previously linked to colorectal cancer pathogenesis or progression. These variants and genes are of interest for future functional and prognostic studies. Our results provide further evidence for the utility of GMDR models in identifying novel prognostic biomarkers and the biological importance of the Wnt pathways in colorectal cancer.

Examining SNP-SNP interactions and risk of clinical outcomes in colorectal cancer using multifactor dimensionality reduction based methods.

Background: SNP interactions may explain the variable outcome risk among colorectal cancer patients. Examining SNP interactions is challenging, especially with large datasets. Multifactor Dimensionality Reduction (MDR)-based programs may address this problem. Objectives: 1) To compare two MDR-based programs for their utility; and 2) to apply these programs to sets of MMP and VEGF-family gene SNPs in order to examine their interactions in relation to colorectal cancer survival outcomes. Methods: This study applied two data reduction methods, Cox-MDR and GMDR 0.9, to study one to three way SNP interactions. Both programs were run using a 5-fold cross validation step and the top models were verified by permutation testing. Prognostic associations of the SNP interactions were verified using multivariable regression methods. Eight datasets, including SNPs from MMP family genes (n = 201) and seven sets of VEGF-family interaction networks (n = 1,517 SNPs) were examined. Results: ∼90 million potential interactions were examined. Analyses in the MMP and VEGF gene family datasets found several novel 1- to 3-way SNP interactions. These interactions were able to distinguish between the patients with different outcome risks (regression p-values 0.03-2.2E-09). The strongest association was detected for a 3-way interaction including CHRM3.rs665159_EPN1.rs6509955_PTGER3.rs1327460 variants. Conclusion: Our work demonstrates the utility of data reduction methods while identifying potential prognostic markers in colorectal cancer.

A comprehensive analysis of SNPs and CNVs identifies novel markers associated with disease outcomes in colorectal cancer.

We aimed to examine the associations of a genome-wide set of single nucleotide polymorphisms (SNPs) and 254 copy number variations (CNVs) and/or insertion/deletions (INDELs) with clinical outcomes in colorectal cancer patients (n = 505). We also aimed to investigate whether their associations changed (e.g., appeared, diminished) over time. Multivariable Cox proportional hazards and piece-wise Cox regression models were used to examine the associations. The Cancer Genome Atlas (TCGA) datasets were used for replication purposes and to examine the gene expression differences between tumor and nontumor tissue samples. A common SNP (WBP11-rs7314075) was associated with disease-specific survival with P-value of 3.2 × 10-8 . Association of this region with disease-specific survival was also detected in the TCGA patient cohort. Two expression quantitative trait loci (eQTLs) were identified in this locus that were implicated in the regulation of ERP27 expression. Interestingly, expression levels of ERP27 and WBP11 were significantly different between colorectal tumors and nontumor tissues. Three SNPs predicted the risk of recurrent disease only after 5 years postdiagnosis. Overall, our study identified novel variants, one of which also showed an association in the TCGA dataset, but no CNVs/INDELs, that associated with outcomes in colorectal cancer. Three SNPs were candidate predictors of long-term recurrence/metastasis risk.

Powerful rare variant association testing in a copula-based joint analysis of multiple phenotypes.

In genetic association studies of rare variants, the low power of association tests is one of the main challenges. In this study, we propose a new single-marker association test called C-JAMP (Copula-based Joint Analysis of Multiple Phenotypes), which is based on a joint model of multiple phenotypes given genetic markers and other covariates. We evaluated its performance and compared its empirical type I error and power with existing univariate and multivariate single-marker and multi-marker rare-variant tests in extensive simulation studies. C-JAMP yielded unbiased genetic effect estimates and valid type I errors with an adjusted test statistic. When strongly dependent traits were jointly analyzed, C-JAMP had the highest power in all scenarios except when a high percentage of variants were causal with moderate/small effect sizes. When traits with weak or moderate dependence were analyzed, whether C-JAMP or competing approaches had higher power depended on the effect size. When C-JAMP was applied with a misspecified copula function, it still achieved high power in some of the scenarios considered. In a real-data application, we analyzed sequencing data using C-JAMP and performed the first genome-wide association studies of high-molecular-weight and medium-molecular-weight adiponectin plasma concentrations. C-JAMP identified 20 rare variants with p-values smaller than 10-5 , while all other tests resulted in the identification of fewer variants with higher p-values. In summary, the results indicate that C-JAMP is a powerful, flexible, and robust method for association studies, and we identified novel candidate markers for adiponectin. C-JAMP is implemented as an R package and freely available from https://cran.r-project.org/package=CJAMP.

The long-term survival characteristics of a cohort of colorectal cancer patients and baseline variables associated with survival outcomes with or without time-varying effects.

BACKGROUND: Colorectal cancer is the third most common cancer in the world. In this study, we assessed the long-term survival characteristics and prognostic associations and potential time-varying effects of clinico-demographic variables and two molecular markers (microsatellite instability (MSI) and BRAF Val600Glu mutation) in a population-based patient cohort followed up to ~ 19 years. METHODS: The patient cohort included 738 incident cases diagnosed between 1999 and 2003. Cox models were used to analyze the association between the variables and a set of survival outcome measures (overall survival (OS), disease-specific survival (DSS), recurrence-free survival (RFS), metastasis-free survival (MFS), recurrence/metastasis-free survival (RMFS), and event-free survival (EFS)). Cox proportional hazard (PH) assumption was tested for all variables, and Cox models with time-varying effects were used if any departure from the PH assumption was detected. RESULTS: During the follow-up, ~ 61% patients died from any cause, ~ 26% died from colorectal cancer, and ~ 10% and ~ 20% experienced recurrences and distant metastases, respectively. Stage IV disease and post-diagnostic recurrence or metastasis were strongly linked to risk of death from colorectal cancer. If a patient had survived the first 6 years without any disease-related event (i.e., recurrence, metastasis, or death from colorectal cancer), their risks became very minimal after this time period. Distinct sets of markers were associated with different outcome measures. In some cases, the effects by variables were constant throughout the follow-up. For example, MSI-high tumor phenotype and older age at diagnosis predicted longer MFS times consistently over the follow-up. However, in some other cases, the effects of the variables varied with time. For example, adjuvant radiotherapy treatment was associated with increased risk of metastasis in patients who received this treatment after 5.5 years post-diagnosis, but not before that. CONCLUSIONS: This study describes the long-term survival characteristics of a prospective cohort of colorectal cancer patients, relationships between baseline variables and a detailed set of patient outcomes over a long time, and time-varying effects of a group of variables. The results presented advance our understanding of the long-term prognostic characteristics in colorectal cancer and are expected to inspire future studies and clinical care strategies.

Correction to: A genome-wide association study identifies single nucleotide polymorphisms associated with time-to-metastasis in colorectal cancer.

.

A genome-wide association study identifies single nucleotide polymorphisms associated with time-to-metastasis in colorectal cancer.

BACKGROUND: Differentiating between cancer patients who will experience metastasis within a short time and who will be long-term survivors without metastasis is a critical aim in healthcare. The microsatellite instability (MSI)-high tumor phenotype is such a differentiator in colorectal cancer, as patients with these tumors are unlikely to experience metastasis. Our aim in this study was to determine if germline genetic variations could further differentiate colorectal cancer patients based on the long-term risk and timing of metastasis. METHODS: The patient cohort consisted of 379 stage I-III Caucasian colorectal cancer patients with microsatellite stable or MSI-low tumors. We performed univariable analysis on 810,622 common single nucleotide polymorphisms (SNPs) under different genetic models. Depending on the long-term metastasis-free survival probability estimates, we applied a mixture cure model, Cox proportional hazards regression model, or log-rank test. For SNPs reaching Bonferroni-corrected significance (p < 6.2 × 10- 8) having valid genetic models, multivariable analysis adjusting for significant baseline characteristics was conducted. RESULTS: After adjusting for significant baseline characteristics, specific genotypes of ten polymorphisms were significantly associated with time-to-metastasis. These polymorphisms are three intergenic SNPs, rs5749032 (p = 1.28 × 10- 10), rs2327990 (p = 9.59 × 10- 10), rs1145724 (p = 3 × 10- 8), and seven SNPs within the non-coding sequences of three genes: FHIT (p = 2.59 × 10- 9), EPHB1 (p = 8.23 × 10- 9), and MIR7515 (p = 4.87 × 10- 8). CONCLUSIONS: Our results suggest novel associations of specific genotypes of SNPs with early metastasis in Caucasian colorectal cancer patients. These associations, once replicated in other patient cohorts, could assist in the development of personalized treatment strategies for colorectal cancer patients.

Functional assessment of rat pulmonary flavin-containing monooxygenase activity.

The expression of flavin-containing monooxygenase (FMO) varies extensively between human and commonly used preclinical species such as rat and mouse. The aim of this study was to investigate the pulmonary FMO activity in rat using benzydamine. Furthermore, the contribution of rat lung to the clearance of benzydamine was investigated using an in vivo pulmonary extraction model. Benzydamine N-oxygenation was observed in lung microsomes and lung slices. Thermal inactivation of FMO and CYP inhibition suggested that rat pulmonary N-oxygenation is predominantly FMO mediated while any contribution from CYPs is negligible. The predicted lung clearance (CLlung) estimated from microsomes and slices was 16 ± 0.6 and 2.1 ± 0.3 mL/min/kg, respectively. The results from in vivo pulmonary extraction indicated no pulmonary extraction following intravenous and intra-arterial dosing to rats. Interestingly, the predicted CLlung using rat lung microsomes corresponded to approximately 35% of rat CLliver suggesting that the lung makes a smaller contribution to the whole body clearance of benzydamine. Although benzydamine clearance in rat appears to be predominantly mediated by hepatic metabolism, the data suggest that the lung may also make a smaller contribution to its whole body clearance.

Comparison of Rat and Human Pulmonary Metabolism Using Precision-cut Lung Slices (PCLS).

BACKGROUND: Although the liver is the primary organ of drug metabolism, the lungs also contain drug-metabolizing enzymes and may, therefore, contribute to the elimination of drugs. In this investigation, the Precision-cut Lung Slice (PCLS) technique was standardized with the aims of characterizing and comparing rat and human pulmonary drug metabolizing activity. METHOD: Due to the limited availability of human lung tissue, standardization of the PCLS method was performed with rat lung tissue. Pulmonary enzymatic activity was found to vary significantly with rat age and rat strain. The Dynamic Organ Culture (DOC) system was superior to well-plates for tissue incubations, while oxygen supply appeared to have a limited impact within the 4h incubation period used here. RESULTS: The metabolism of a range of phase I and phase II probe substrates was assessed in rat and human lung preparations. Cytochrome P450 (CYP) activity was relatively low in both species, whereas phase II activity appeared to be more significant. CONCLUSION: PCLS is a promising tool for the investigation of pulmonary drug metabolism. The data indicates that pulmonary CYP activity is relatively low and that there are significant differences in enzyme activity between rat and human lung.

Associations of single nucleotide polymorphisms with mucinous colorectal cancer: genome-wide common variant and gene-based rare variant analyses.

BACKGROUND: Colorectal cancer has significant impact on individuals and healthcare systems. Many genes have been identified to influence its pathogenesis. However, the genetic basis of mucinous tumor histology, an aggressive subtype of colorectal cancer, is currently not well-known. This study aimed to identify common and rare genetic variations that are associated with the mucinous tumor phenotype. METHODS: Genome-wide single nucleotide polymorphism (SNP) data was investigated in a colorectal cancer patient cohort (n = 505). Association analyses were performed for 729,373 common SNPs and 275,645 rare SNPs. Common SNP association analysis was performed using univariable and multivariable logistic regression under different genetic models. Rare-variant association analysis was performed using a multi-marker test. RESULTS: No associations reached the traditional genome-wide significance. However, promising genetic associations were identified. The identified common SNPs significantly improved the discriminatory accuracy of the model for mucinous tumor phenotype. Specifically, the area under the receiver operating characteristic curve increased from 0.703 (95% CI: 0.634-0.773) to 0.916 (95% CI: 0.873-0.960) when considering the most significant SNPs. Additionally, the rare variant analysis identified a number of genetic regions that potentially contain causal rare variants associated with the mucinous tumor phenotype. CONCLUSIONS: This is the first study applying both common and rare variant analyses to identify genetic associations with mucinous tumor phenotype using a genome-wide genotype data. Our results suggested novel associations with mucinous tumors. Once confirmed, these results will not only help us understand the biological basis of mucinous histology, but may also help develop targeted treatment options for mucinous tumors.

XRCC3 Thr241Met and TYMS variable number tandem repeat polymorphisms are associated with time-to-metastasis in colorectal cancer.

BACKGROUND: Metastasis is a major cause of mortality in cancer. Identifying prognostic factors that distinguish patients who will experience metastasis in the short-term and those that will be free of metastasis in the long-term is of particular interest in current medical research. The objective of this study was to examine if select genetic polymorphisms can differentiate colorectal cancer patients based on timing and long-term risk of metastasis. METHODS: The patient cohort consisted of 402 stage I-III colorectal cancer patients with microsatellite instability (MSI)-low (MSI-L) or microsatellite stable (MSS) tumors. We applied multivariable mixture cure model, which is the proper model when there is a substantial group of patients who remain free of metastasis in the long-term, to 26 polymorphisms. Time-dependent receiver operator characteristic (ROC) curve analysis was performed to determine the change in discriminatory accuracy of the models when the significant SNPs were included. RESULTS: After adjusting for significant baseline characteristics, two polymorphisms were significantly associated with time-to-metastasis: TT and TC genotypes of the XRCC3 Thr241Met (p = 0.042) and the 3R/3R genotype of TYMS 5'-UTR variable number tandem repeat (VNTR) (p = 0.009) were associated with decreased time-to-metastasis. ROC curves showed that the discriminatory accuracy of the model is increased slightly when these polymorphisms were added to the significant baseline characteristics. CONCLUSIONS: Our results indicate XRCC3 Thr241Met and TYMS 5'-UTR VNTR polymorphisms are associated with time-to-metastasis, and may have potential biological roles in expediting the metastatic process. Once replicated, these associations could contribute to the development of precision medicine for colorectal cancer patients.

Assessment of the pulmonary CYP1A1 metabolism of mavoglurant (AFQ056) in rat.

1. AFQ056 phenotyping results indicate that CYP1A1 is responsible for the formation of the oxidative metabolite, M3. In line with the predominant assumption that CYP1A1 is mainly expressed in extrahepatic tissues, only traces of M3 were detected in hepatic systems. The aim of this study was to investigate the pulmonary CYP1A1 mediated metabolism of AFQ056 in rat. 2. Western blot analysis confirmed that CYP1A1 is expressed in rat lung albeit at low levels. M3 formation was clearly observed in recombinant rat CYP1A1, lung microsomes and lung tissue slices and was strongly inhibited by ketoconazole in the incubations. As CYP3A4 and CYP2C9 metabolites were only observed at trace levels, we concluded that the reduced M3 formation was due to CYP1A1 inhibition. 3. AFQ056 lung clearance (CLlung) as estimated from in vitro data was predicted to be negligible (<1% pulmonary blood flow). This was confirmed by in vivo experiments where intravenous and intra-arterial dosing to rats failed to show significant pulmonary extraction. 4. While rat lung may make a contribution to the formation of M3, it is unlikely to be the only organ involved in this process and further experiments are required to investigate the potential metabolic elimination routes for AFQ056.

Estimating and testing direct genetic effects in directed acyclic graphs using estimating equations.

In genetic association studies, it is important to distinguish direct and indirect genetic effects in order to build truly functional models. For this purpose, we consider a directed acyclic graph setting with genetic variants, primary and intermediate phenotypes, and confounding factors. In order to make valid statistical inference on direct genetic effects on the primary phenotype, it is necessary to consider all potential effects in the graph, and we propose to use the estimating equations method with robust Huber-White sandwich standard errors. We evaluate the proposed causal inference based on estimating equations (CIEE) method and compare it with traditional multiple regression methods, the structural equation modeling method, and sequential G-estimation methods through a simulation study for the analysis of (completely observed) quantitative traits and time-to-event traits subject to censoring as primary phenotypes. The results show that CIEE provides valid estimators and inference by successfully removing the effect of intermediate phenotypes from the primary phenotype and is robust against measured and unmeasured confounding of the indirect effect through observed factors. All other methods except the sequential G-estimation method for quantitative traits fail in some scenarios where their test statistics yield inflated type I errors. In the analysis of the Genetic Analysis Workshop 19 dataset, we estimate and test genetic effects on blood pressure accounting for intermediate gene expression phenotypes. The results show that CIEE can identify genetic variants that would be missed by traditional regression analyses. CIEE is computationally fast, widely applicable to different fields, and available as an R package.

Comparison of single-marker and multi-marker tests in rare variant association studies of quantitative traits.

In genetic association studies of rare variants, low statistical power and potential violations of established estimator properties are among the main challenges of association tests. Multi-marker tests (MMTs) have been proposed to target these challenges, but any comparison with single-marker tests (SMTs) has to consider that their aim is to identify causal genomic regions instead of variants. Valid power comparisons have been performed for the analysis of binary traits indicating that MMTs have higher power, but there is a lack of conclusive studies for quantitative traits. The aim of our study was therefore to fairly compare SMTs and MMTs in their empirical power to identify the same causal loci associated with a quantitative trait. The results of extensive simulation studies indicate that previous results for binary traits cannot be generalized. First, we show that for the analysis of quantitative traits, conventional estimation methods and test statistics of single-marker approaches have valid properties yielding association tests with valid type I error, even when investigating singletons or doubletons. Furthermore, SMTs lead to more powerful association tests for identifying causal genes than MMTs when the effect sizes of causal variants are large, and less powerful tests when causal variants have small effect sizes. For moderate effect sizes, whether SMTs or MMTs have higher power depends on the sample size and percentage of causal SNVs. For a more complete picture, we also compare the power in studies of quantitative and binary traits, and the power to identify causal genes with the power to identify causal rare variants. In a genetic association analysis of systolic blood pressure in the Genetic Analysis Workshop 19 data, SMTs yielded smaller p-values compared to MMTs for most of the investigated blood pressure genes, and were least influenced by the definition of gene regions.

Germline INDELs and CNVs in a cohort of colorectal cancer patients: their characteristics, associations with relapse-free survival time, and potential time-varying effects on the risk of relapse.

INDELs and CNVs are structural variations that may play roles in cancer susceptibility and patient outcomes. Our objectives were a) to computationally detect and examine the genome-wide INDEL/CNV profiles in a cohort of colorectal cancer patients, and b) to examine the associations of frequent INDELs/CNVs with relapse-free survival time. We also identified unique variants in 13 Familial Colorectal Cancer Type X (FCCX) cases. The study cohort consisted of 495 colorectal cancer patients. QuantiSNP and PennCNV algorithms were utilized to predict the INDELs/CNVs using genome-wide signal intensity data. Duplex PCR was used to validate predictions for 10 variants. Multivariable Cox regression models were used to test the associations of 106 common variants with relapse-free survival time. Score test and the multivariable Cox proportional hazards models with time-varying coefficients were applied to identify the variants with time-varying effects on the relapse-free survival time. A total of 3486 distinct INDELs/CNVs were identified in the patient cohort. The majority of these variants were rare (83%) and deletion variants (81%). The results of the computational predictions and duplex PCR results were highly concordant (93-100%). We identified four promising variants significantly associated with relapse-free survival time (P < 0.05) in the multivariable Cox proportional hazards regression models after adjustment for clinical factors. More importantly, two additional variants were identified to have time-varying effects on the risk of relapse. Finally, 58 rare variants were identified unique to the FCCX cases; none of them were detected in more than one patient. This is one of the first genome-wide analyses that identified the germline INDEL/CNV profiles in colorectal cancer patients. Our analyses identified novel variants and genes that can biologically affect the risk of relapse in colorectal cancer patients. Additionally, for the first time, we identified germline variants that can potentially be early-relapse markers in colorectal cancer.

The application of crystallization in the presence of additives to enable drug-in-capsule technology.

Crystallization in the presence of additives such as surfactants, polymers or impurities has been widely investigated in the pharmaceutical and chemical industries in order to change the crystal habit or to obtain a more desirable polymorph, affect crystal growth and influence dissolution. However, in this study, we investigated the concept of crystallization in the presence of surfactants in order to incorporate into the crystal lattice, a small amount (less than 1% w/w) of surfactant, sodium lauryl sulfate (SLS). The resulting crystals were further compared to crystals coated with SLS using a washing procedure; in order to assess whether either procedure generates improvements in the apparent solubility and dissolution of a poorly soluble drug so it can be filled directly into a capsule without the need of a complex formulation process.

Genetic association analysis based on a joint model of gene expression and blood pressure.

Recent work on genetic association studies suggests that much of the heritable variation in complex traits is unexplained, which indicates a need for using more biologically meaningful modeling approaches and appropriate statistical methods. In this study, we propose a biological framework and a corresponding statistical model incorporating multilevel biological measures, and illustrate it in the analysis of the real data provided by the Genetic Analysis Workshop (GAW) 19, which contains whole genome sequence (WGS), gene expression (GE), and blood pressure (BP) data. We investigate the direct effect of single-nucleotide variants (SNVs) on BP and GE, while considering the non-directional dependence between BP and GE, by using copula functions to jointly model BP and GE conditional on SNVs. We implement the method for analysis on a genome-wide scale, and illustrate it within an association analysis of 68,727 SNVs on chromosome 19 that lie in or around genes with available GE measures. Although there is no indication for inflated type I errors under the proposed method, our results show that the association tests have smaller p values than tests under univariate models for common and rare variants using single-variant tests and gene-based multimarker tests. Hence, considering multilevel biological measures and modeling the dependence structure between these measures by using a plausible graphical approach may lead to more informative findings than standard univariate tests of common variants and well-recognized gene-based rare variant tests.

Bivariate genetic association analysis of systolic and diastolic blood pressure by copula models.

We conduct genetic association analysis in the subset of unrelated individuals from the San Antonio Family Studies pedigrees, applying a two-stage approach to take account of the dependence between systolic and diastolic blood pressure (SBP and DBP). In the first stage, we adjust blood pressure for the effects of age, sex, smoking, and use of antihypertensive medication based on a novel modification of censored regression. In the second stage, we model the bivariate distribution of the adjusted SBP and DBP phenotypes by a copula function with interpretable SBP-DBP correlation parameters. This allows us to identify genetic variants associated with each of the adjusted blood pressures, as well as variants that explain the association between the two phenotypes. Within this framework, we define a pleiotropic variant as one that reduces the SBP-DBP correlation. Our results for whole genome sequence variants in the gene ULK4 on chromosome 3 suggest that inference obtained from a copula model can be more informative than findings from the SBP-specific and DBP-specific univariate models alone.

Insights from mixture cure modeling of molecular markers for prognosis in breast cancer.

With the ultimate aim of improving clinical management of breast cancer, investigators have sought to identify molecular genetic markers that stratify newly diagnosed patients into subtypes differing in short- or long-term prognosis. Conventional survival models can fail to describe adequately the relationship between subtype and disease recurrence, particularly when there is a substantial proportion of long-term disease-free survivors. The observed patterns of disease-free survival in an undifferentiated patient cohort may be explained by an underlying mixture of two subgroups: patients who will remain free of disease in the long term (ie, cured), and those who will experience disease recurrence within their lifetime (ie, susceptible.) In this article, we review the concepts and methods of the mixture cure models and apply them in the analysis of molecular genetic prognostic factors for disease-free survival and time to disease recurrence in a cohort of patients with axillary lymph node-negative breast cancer.