A genome-wide association study in human lymphoblastoid cells supports safety of mitochondrial complex I inhibitor.
Ontology highlight
ABSTRACT:
SUBMITTER: Gao H
PROVIDER: S-EPMC8743030 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Radiation therapy is used to treat half of all cancer patients. Response to radiation therapy varies widely among patients. Therefore, we performed a genome-wide association study (GWAS) to identify biomarkers to help predict radiation response using 277 ethnically defined human lymphoblastoid cell lines (LCLs). Basal gene expression levels and 1.3 million genome-wide single nucleotide polymorphism (SNP) markers from both Affymetrix and Illumina platforms were assayed for all 277 human LCLs. MTS [3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium] assays for radiation cytotoxicity were also performed to obtain area under the curve (AUC) as a radiation response phenotype for use in the association studies. Functional validation of candidate genes, selected from an integrated analysis that used SNP, expression, and AUC data, was performed with multiple cancer cell lines using specific siRNA knockdown, followed by MTS and colony-forming assays. A total of 27 loci, each containing at least two SNPs within 50 kb with P-values less than 10(-4) were associated with radiation AUC. A total of 270 expression probe sets were associated with radiation AUC with P < 10(-3). The integrated analysis identified 50 SNPs in 14 of the 27 loci that were associated with both AUC and the expression of 39 genes, which were also associated with radiation AUC (P < 10(-3)). Functional validation using siRNA knockdown in multiple tumor cell lines showed that C13orf34, MAD2L1, PLK4, TPD52, and DEPDC1B each significantly altered radiation sensitivity in at least two cancer cell lines. Studies performed with LCLs can help to identify novel biomarkers that might contribute to variation in response to radiation therapy and enhance our understanding of mechanisms underlying that variation.
| S-EPMC2963812 | biostudies-literature
Project description:Association mapping with lymphoblastoid cell lines (LCLs) is a promising approach in pharmacogenomics research, and in the current study we utilized LCLs to perform association mapping for 29 chemotherapy drugs.Currently, we use LCLs to perform genome-wide association mapping of the cytotoxic response of 520 European-Americans to 29 different anticancer drugs; the largest LCL study to date. A novel association approach using a multivariate analysis of covariance design was employed with the software program MAGWAS, testing for differences in the dose-response profiles between genotypes without making assumptions about the response curve or the biologic mode of association. Additionally, by classifying 25 of the 29 drugs into eight families according to structural and mechanistic relationships, MAGWAS was used to test for associations that were shared across each drug family. Finally, a unique algorithm using multivariate responses and multiple linear regressions across pairs of response curves was used for unsupervised clustering of drugs.Among the single-drug studies, suggestive associations were obtained for 18 loci, 12 within/near genes. Three of these, MED12L, CHN2 and MGMT, have been previously implicated in cancer pharmacogenomics. The drug family associations resulted in four additional suggestive loci (three contained within/near genes). One of these genes, HDAC4, associated with the DNA alkylating agents, shows possible clinical interactions with temozolomide. For the drug clustering analysis, 18 of 25 drugs clustered into the appropriate family.This study demonstrates the utility of LCLs in identifying genes that have clinical importance in drug response and for assigning unclassified agents to specific drug families, and proposes new candidate genes for follow-up in a large number of chemotherapy drugs.
| S-EPMC4055088 | biostudies-literature
Project description:The extent to which RNA stability differs between individuals and its contribution to the interindividual expression variation remain unknown. We conducted a genome-wide analysis of RNA stability in seven human HapMap lymphoblastoid cell lines (LCLs) and analyzed the effect of DNA sequence variation on RNA half-life differences. Twenty-six percent of the expressed genes exhibited RNA half-life differences between LCLs at a false discovery rate (FDR) < 0.05, which accounted for ~ 37% of the gene expression differences between individuals. Nonsense polymorphisms were associated with reduced RNA half-lives. In genes presenting interindividual RNA half-life differences, higher coding GC3 contents (G and C percentages at the third-codon positions) were correlated with increased RNA half-life. Consistently, G and C alleles of single nucleotide polymorphisms (SNPs) in protein coding sequences were associated with enhanced RNA stability. These results suggest widespread interindividual differences in RNA stability related to DNA sequence and composition variation.
| S-EPMC3576867 | biostudies-literature
Project description:BackgroundTaxane is one of the first line treatments of lung cancer. In order to identify novel single nucleotide polymorphisms (SNPs) that might contribute to taxane response, we performed a genome-wide association study (GWAS) for two taxanes, paclitaxel and docetaxel, using 276 lymphoblastoid cell lines (LCLs), followed by genotyping of top candidate SNPs in 874 lung cancer patient samples treated with paclitaxel.MethodsGWAS was performed using 1.3 million SNPs and taxane cytotoxicity IC50 values for 276 LCLs. The association of selected SNPs with overall survival in 76 small or 798 non-small cell lung cancer (SCLC, NSCLC) patients were analyzed by Cox regression model, followed by integrated SNP-microRNA-expression association analysis in LCLs and siRNA screening of candidate genes in SCLC (H196) and NSCLC (A549) cell lines.Results147 and 180 SNPs were associated with paclitaxel or docetaxel IC50s with p-values <10-4 in the LCLs, respectively. Genotyping of 153 candidate SNPs in 874 lung cancer patient samples identified 8 SNPs (p-value < 0.05) associated with either SCLC or NSCLC patient overall survival. Knockdown of PIP4K2A, CCT5, CMBL, EXO1, KMO and OPN3, genes within 200 kb up-/downstream of the 3 SNPs that were associated with SCLC overall survival (rs1778335, rs2662411 and rs7519667), significantly desensitized H196 to paclitaxel. SNPs rs2662411 and rs1778335 were associated with mRNA expression of CMBL or PIP4K2A through microRNA (miRNA) hsa-miR-584 or hsa-miR-1468.ConclusionsGWAS in an LCL model system, joined with clinical translational and functional studies, might help us identify genetic variations associated with overall survival of lung cancer patients treated paclitaxel.
| S-EPMC3573965 | biostudies-literature
Project description:Recently, lymphoblastoid cell lines (LCLs) have emerged as an innovative model system for mapping gene variants that predict the dose response to chemotherapy drugs.In the current study, this strategy was expanded to the in-vitro genome-wide association approach, using 516 LCLs derived from a White cohort to assess the cytotoxic response to temozolomide.Genome-wide association analysis using ?2.1 million quality-controlled single-nucleotide polymorphisms (SNPs) identified a statistically significant association (P<10(-8)) with SNPs in the O(6)-methylguanine-DNA methyltransferase (MGMT) gene. We also show that the primary SNP in this region is significantly associated with the differential gene expression of MGMT (P<10(-26)) in LCLs and differential methylation in glioblastoma samples from The Cancer Genome Atlas.The previously documented clinical and functional relationships between MGMT and temozolomide response highlight the potential of well-powered genome-wide association studies of the LCL model system to identify meaningful genetic associations.
| S-EPMC3691078 | biostudies-literature
Project description:BACKGROUND: Previous studies of network properties of human disease genes have mainly focused on monogenic diseases or cancers and have suffered from discovery bias. Here we investigated the network properties of complex disease genes identified by genome-wide association studies (GWAs), thereby eliminating discovery bias. PRINCIPAL FINDINGS: We derived a network of complex diseases (n = 54) and complex disease genes (n = 349) to explore the shared genetic architecture of complex diseases. We evaluated the centrality measures of complex disease genes in comparison with essential and monogenic disease genes in the human interactome. The complex disease network showed that diseases belonging to the same disease class do not always share common disease genes. A possible explanation could be that the variants with higher minor allele frequency and larger effect size identified using GWAs constitute disjoint parts of the allelic spectra of similar complex diseases. The complex disease gene network showed high modularity with the size of the largest component being smaller than expected from a randomized null-model. This is consistent with limited sharing of genes between diseases. Complex disease genes are less central than the essential and monogenic disease genes in the human interactome. Genes associated with the same disease, compared to genes associated with different diseases, more often tend to share a protein-protein interaction and a Gene Ontology Biological Process. CONCLUSIONS: This indicates that network neighbors of known disease genes form an important class of candidates for identifying novel genes for the same disease.
| S-EPMC2779513 | biostudies-literature
Project description:The mammalian target of rapamycin (mTOR) inhibitors, a set of promising potential anti-cancer agents, has shown response variability among individuals. This study aimed to identify novel biomarkers and mechanisms that might influence the response to Rapamycin and Everolimus. Genome-wide association (GWA) analyses involving single nucleotide polymorphisms (SNPs), mRNA, and microRNAs microarray data were assessed for association with area under the cytotoxicity dose response curve (AUC) of two mTOR inhibitors in 272 human lymphoblastoid cell lines (LCLs). Integrated analysis among SNPs, expression data, microRNA data and AUC values were also performed to help select candidate genes for further functional characterization. Functional validation of candidate genes using siRNA screening in multiple cell lines followed by MTS assays for the two mTOR inhibitors were performed. We found that 16 expression probe sets (genes) that overlapped between the two drugs were associated with AUC values of two mTOR inhibitors. One hundred and twenty seven and one hundred SNPs had P < 10(-4), while 8 and 10 SNPs had P < 10(-5) with Rapamycin and Everolimus AUC, respectively. Functional studies indicated that 13 genes significantly altered cell sensitivity to either one or both drugs in at least one cell line. Additionally, one microRNA, miR-10a, was significantly associated with AUC values for both drugs and was shown to repress expression of genes that were associated with AUC and desensitize cells to both drugs. In summary, this study identified genes and a microRNA that might contribute to response to mTOR inhibitors.
| S-EPMC3757297 | biostudies-literature
Project description:Metformin is currently considered as a promising anticancer agent in addition to its anti-diabetic effect. To better individualize metformin therapy and explore novel molecular mechanisms in cancer treatment, we conducted a pharmacogenomic study using 266 lymphoblastoid cell lines (LCLs). Metformin cytotoxicity assay was performed using the MTS assay. Genome-wide association (GWA) analyses were performed in LCLs using 1.3 million SNPs, 485k DNA methylation probes, 54k mRNA expression probe sets, and metformin cytotoxicity (IC50s). Top candidate genes were functionally validated using siRNA screening, followed by MTS assay in breast cancer cell lines. Further study of one top candidate, STUB1, was performed to elucidate the mechanisms by which STUB1 might contribute to metformin action. GWA analyses in LCLs identified 198 mRNA expression probe sets, 12 SNP loci, and 5 DNA methylation loci associated with metformin IC50 with P-values?<10?4 or?<10?5. Integrated SNP/methylation loci-expression-IC50 analyses found 3 SNP loci or 5 DNA methylation loci associated with metformin IC50 through trans-regulation of expression of 11 or 26 genes with P-value?<10?4. Functional validation of top 61 candidate genes in 4 IPA networks indicated down regulation of 14 genes significantly altered metformin sensitivity in two breast cancer cell lines. Mechanistic studies revealed that the E3 ubiquitin ligase, STUB1, could influence metformin response by facilitating proteasome-mediated degradation of cyclin A. GWAS using a genomic data-enriched LCL model system, together with functional and mechanistic studies using cancer cell lines, help us to identify novel genetic and epigenetic biomarkers involved in metformin anticancer response.
| S-EPMC6078562 | biostudies-literature
Project description:BackgroundExpression level of many genes shows abundant natural variation in human populations. The variations in gene expression are believed to contribute to phenotypic differences. Emerging evidence has shown that microRNAs (miRNAs) are one of the key regulators of gene expression. However, past studies have focused on the miRNA target genes and used loss- or gain-of-function approach that may not reflect natural association between miRNA and mRNAs.Methodology/principal findingsTo examine miRNA regulatory effect on global gene expression under endogenous condition, we performed pair-wise correlation coefficient analysis on expression levels of 366 miRNAs and 14,174 messenger RNAs (mRNAs) in 90 immortalized lymphoblastoid cell lines, and observed significant correlations between the two species of RNA transcripts. We identified a total of 7,207 significantly correlated miRNA-mRNA pairs (false discovery rate q<0.01). Of those, 4,085 pairs showed positive correlations while 3,122 pairs showed negative correlations. Gene ontology analyses on the miRNA-correlated genes revealed significant enrichments in several biological processes related to cell cycle, cell communication and signal transduction. Individually, each of three miRNAs (miR-331, -98 and -33b) demonstrated significant correlation with the genes in cell cycle-related biological processes, which is consistent with important role of miRNAs in cell cycle regulation.Conclusions/significanceThis study demonstrates feasibility of using naturally expressed transcript profiles to identify endogenous correlation between miRNA and miRNA. By applying this genome-wide approach, we have identified thousands of miRNA-correlated genes and revealed potential role of miRNAs in several important cellular functions. The study results along with accompanying data sets will provide a wealth of high-throughput data to further evaluate the miRNA-regulated genes and eventually in phenotypic variations of human populations.
| S-EPMC2691578 | biostudies-literature
Project description:We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.
| S-EPMC8742999 | biostudies-literature