RESUMEN
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.
Asunto(s)
Simulación por Computador , Marcadores Genéticos/genética , Estudio de Asociación del Genoma Completo/métodos , Modelos Genéticos , Enfermedades Raras/genética , Estudios de Asociación Genética , Variación Genética/genética , Humanos , FenotipoRESUMEN
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.
Asunto(s)
Estudios de Asociación Genética/métodos , Presión Sanguínea/genética , Factores de Confusión Epidemiológicos , Conjuntos de Datos como Asunto , Variación Genética , Humanos , Modelos Genéticos , Fenotipo , Polimorfismo de Nucleótido Simple/genética , Análisis de Regresión , Proyectos de Investigación , Programas InformáticosRESUMEN
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.
Asunto(s)
Supervivientes de Cáncer/estadística & datos numéricos , Neoplasias Colorrectales/mortalidad , Adulto , Anciano , Estudios de Cohortes , Neoplasias Colorrectales/epidemiología , Neoplasias Colorrectales/genética , Supervivencia sin Enfermedad , Femenino , Estudios de Seguimiento , Humanos , Masculino , Inestabilidad de Microsatélites , Persona de Mediana Edad , Mutación Missense , Recurrencia Local de Neoplasia/epidemiología , Recurrencia Local de Neoplasia/genética , Fenotipo , Pronóstico , Proteínas Proto-Oncogénicas B-raf/genética , Análisis de Supervivencia , Factores de Tiempo , Adulto JovenRESUMEN
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.
RESUMEN
Lymphatic invasion (LVI) is associated with disease recurrence in axillary node-negative (ANN) breast cancer. Using gene expression profiling of 105 ANN tumors, we found that podocalyxin (PODXL) was more highly expressed in tumors with LVI (LVI+) than in those without LVI (LVI-). Differences in PODXL expression were validated using real-time polymerase chain reaction as well as by immunohistochemistry in an independent set of 652 tumors on tissue microarrays. Disease-free survival (DFS) analyses were conducted for association of high PODXL protein expression with risk of distant recurrence overall and within breast cancer subtypes using both Cox and cure-rate models. High PODXL expression was associated with poor prognosis features including large tumor size, high histological grade, estrogen and progesterone receptor negativity, and with clinical alterations characteristic of the basal-like breast cancer phenotype. Surprisingly, despite having other poor prognosis characteristics, women with high PODXL expressing tumors had better long-term DFS in multivariate analysis with traditional clinicopathologic factors including LVI and HER2 status (P = 0.001). PODXL has the potential to be a useful biomarker for identifying good prognosis patients in characteristically poor prognosis breast cancer groups and may impact treatment of women with this disease.
Asunto(s)
Neoplasias de la Mama/metabolismo , Neoplasias de la Mama/patología , Ganglios Linfáticos/patología , Fenotipo , Sialoglicoproteínas/metabolismo , Axila , Biomarcadores de Tumor/metabolismo , Neoplasias de la Mama/genética , Neoplasias de la Mama/mortalidad , Femenino , Expresión Génica , Humanos , Metástasis Linfática , Pronóstico , ARN Mensajero/genética , ARN Mensajero/metabolismo , Sialoglicoproteínas/genética , Carga TumoralRESUMEN
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.
RESUMEN
Group 14 of Genetic Analysis Workshop 17 examined several issues related to analysis of complex traits using DNA sequence data. These issues included novel methods for analyzing rare genetic variants in an aggregated manner (often termed collapsing rare variants), evaluation of various study designs to increase power to detect effects of rare variants, and the use of machine learning approaches to model highly complex heterogeneous traits. Various published and novel methods for analyzing traits with extreme locus and allelic heterogeneity were applied to the simulated quantitative and disease phenotypes. Overall, we conclude that power is (as expected) dependent on locus-specific heritability or contribution to disease risk, large samples will be required to detect rare causal variants with small effect sizes, extreme phenotype sampling designs may increase power for smaller laboratory costs, methods that allow joint analysis of multiple variants per gene or pathway are more powerful in general than analyses of individual rare variants, population-specific analyses can be optimal when different subpopulations harbor private causal mutations, and machine learning methods may be useful for selecting subsets of predictors for follow-up in the presence of extreme locus heterogeneity and large numbers of potential predictors.
Asunto(s)
Predisposición Genética a la Enfermedad/genética , Epidemiología Molecular/métodos , Polimorfismo de Nucleótido Simple/genética , Análisis de Regresión , Inteligencia Artificial , Interpretación Estadística de Datos , Minería de Datos , Exoma , Variación Genética , Proyecto Genoma Humano , Humanos , Metaanálisis como Asunto , Análisis de SecuenciaRESUMEN
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.
RESUMEN
Sequentially observed survival times are of interest in many studies but there are difficulties in analyzing such data using nonparametric or semiparametric methods. First, when the duration of followup is limited and the times for a given individual are not independent, induced dependent censoring arises for the second and subsequent survival times. Non-identifiability of the marginal survival distributions for second and later times is another issue, since they are observable only if preceding survival times for an individual are uncensored. In addition, in some studies a significant proportion of individuals may never have the first event. Fully parametric models can deal with these features, but robustness is a concern. We introduce a new approach to address these issues. We model the joint distribution of the successive survival times by using copula functions, and provide semiparametric estimation procedures in which copula parameters are estimated without parametric assumptions on the marginal distributions. This provides more robust estimates and checks on the fit of parametric models. The methodology is applied to a motivating example involving relapse and survival following colon cancer treatment.
Asunto(s)
Modelos Estadísticos , Análisis de Varianza , Neoplasias del Colon/mortalidad , Neoplasias del Colon/terapia , Humanos , Funciones de Verosimilitud , Recurrencia , Tasa de SupervivenciaRESUMEN
Copula models for multivariate lifetimes have become widely used in areas such as biomedicine, finance and insurance. This paper fills some gaps in existing methodology for copula parameters and model assessment. We consider procedures based on likelihood and pseudolikelihood ratio statistics and introduce semiparametric maximum likelihood estimation leading to semiparametric versions. For cases where standard asymptotic approximations do not hold, we propose an efficient simulation technique for obtaining p-values. We apply these methods to tests for a copula model, based on embedding it in a larger copula family. It is shown that the likelihood and pseudolikelihood ratio tests are consistent even when the expanded copula model is misspecified. Power comparisons with two other tests of fit indicate that model expansion provides a convenient, powerful and robust approach. The methods are illustrated on an application concerning the time to loss of vision in the two eyes of an individual.
Asunto(s)
Biometría/métodos , Funciones de Verosimilitud , Simulación por Computador , Retinopatía Diabética/terapia , Humanos , Fotocoagulación/normas , Modelos Biológicos , Modelos Estadísticos , Agudeza VisualRESUMEN
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.
Asunto(s)
Neoplasias Colorrectales , Polimorfismo de Nucleótido Simple , Neoplasias Colorrectales/patología , Variaciones en el Número de Copia de ADN/genética , Proteínas de Unión al ADN/genética , Predisposición Genética a la Enfermedad , Estudio de Asociación del Genoma Completo , Humanos , Polimorfismo de Nucleótido Simple/genética , Sitios de Carácter Cuantitativo/genética , Factores de Empalme de ARN/genéticaRESUMEN
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.
RESUMEN
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.
Asunto(s)
Neoplasias Colorrectales/genética , Proteínas de Unión al ADN/genética , Metionina/genética , Metástasis de la Neoplasia , Treonina/genética , Timidilato Sintasa/genética , Neoplasias Colorrectales/patología , Humanos , Polimorfismo de Nucleótido SimpleRESUMEN
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.
Asunto(s)
Marcadores Genéticos , Sitios de Carácter Cuantitativo , Humanos , Polimorfismo de Nucleótido SimpleRESUMEN
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.
Asunto(s)
Neoplasias Colorrectales/genética , Variaciones en el Número de Copia de ADN , Mutación INDEL , Anciano , Supervivencia sin Enfermedad , Femenino , Predisposición Genética a la Enfermedad , Humanos , Masculino , Pronóstico , RecurrenciaRESUMEN
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.
RESUMEN
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.
RESUMEN
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.
Asunto(s)
Biomarcadores de Tumor/análisis , Neoplasias de la Mama/genética , Marcadores Genéticos/genética , Modelos Estadísticos , Adulto , Anciano , Axila , Neoplasias de la Mama/patología , Estudios de Cohortes , Progresión de la Enfermedad , Supervivencia sin Enfermedad , Femenino , Humanos , Estimación de Kaplan-Meier , Queratina-5/análisis , Ganglios Linfáticos/patología , Metástasis Linfática/diagnóstico , Persona de Mediana Edad , Recurrencia Local de Neoplasia/diagnóstico , Valor Predictivo de las Pruebas , Pronóstico , Receptor ErbB-2/análisis , Receptores de Estrógenos/análisis , Receptores de Progesterona/análisis , Recurrencia , Factores de TiempoRESUMEN
Use of trait-dependent sampling designs in whole-genome association studies of sequence data can reduce total sequencing costs with modest losses of statistical efficiency. In a quantitative trait (QT) analysis of data from the Genetic Analysis Workshop 17 mini-exome for unrelated individuals in the Asian subpopulation, we investigate alternative designs that sequence only 50% of the entire cohort. In addition to a simple random sampling design, we consider extreme-phenotype designs that are of increasing interest in genetic association analysis of QTs, especially in studies concerned with the detection of rare genetic variants. We also evaluate a novel sampling design in which all individuals have a nonzero probability of being selected into the sample but in which individuals with extreme phenotypes have a proportionately larger probability. We take differential sampling of individuals with informative trait values into account by inverse probability weighting using standard survey methods which thus generalizes to the source population. In replicate 1 data, we applied the designs in association analysis of Q1 with both rare and common variants in the FLT1 gene, based on knowledge of the generating model. Using all 200 replicate data sets, we similarly analyzed Q1 and Q4 (which is known to be free of association with FLT1) to evaluate relative efficiency, type I error, and power. Simulation study results suggest that the QT-dependent selection designs generally yield greater than 50% relative efficiency compared to using the entire cohort, implying cost-effectiveness of 50% sample selection and worthwhile reduction of sequencing costs.