Genetic architecture of microRNA expression: implications for the transcriptome and complex traits.

We sought to comprehensively and systematically characterize the relationship between genetic variation, miRNA expression, and mRNA expression. Genome-wide expression profiling of samples of European and African ancestry identified in each population hundreds of miRNAs whose increased expression is correlated with correspondingly reduced expression of target mRNAs. We scanned 3' UTR SNPs with a potential functional effect on miRNA binding for cis-acting expression quantitative trait loci (eQTLs) for the corresponding proximal target genes. To extend sequence-based, localized analyses of SNP effect on miRNA binding, we proceeded to dissect the genetic basis of miRNA expression variation; we mapped miRNA expression levels-as quantitative traits-to loci in the genome as miRNA eQTLs, demonstrating that miRNA expression is under significant genetic control. We found that SNPs associated with miRNA expression are significantly enriched with those SNPs already shown to be associated with mRNA. Moreover, we discovered that many of the miRNA-associated genetic variations identified in our study are associated with a broad spectrum of human complex traits from the National Human Genome Research Institute catalog of published genome-wide association studies. Experimentally, we replicated miRNA-induced mRNA expression inhibition and the cis-eQTL relationship to the target gene for several identified relationships among SNPs, miRNAs, and mRNAs in an independent set of samples; furthermore, we conducted miRNA overexpression and inhibition experiments to functionally validate the miRNA-mRNA relationships. This study extends our understanding of the genetic regulation of the transcriptome and suggests that genetic variation might underlie observed relationships between miRNAs and mRNAs more commonly than has previously been appreciated.

Integrative analyses of genetic variation, epigenetic regulation, and the transcriptome to elucidate the biology of platinum sensitivity.

BACKGROUND: Using genome-wide genetic, gene expression, and microRNA expression (miRNA) data, we developed an integrative approach to investigate the genetic and epigenetic basis of chemotherapeutic sensitivity. RESULTS: Through a sequential multi-stage framework, we identified genes and miRNAs whose expression correlated with platinum sensitivity, mapped these to genomic loci as quantitative trait loci (QTLs), and evaluated the associations between these QTLs and platinum sensitivity. A permutation analysis showed that top findings from our approach have a much lower false discovery rate compared to those from a traditional GWAS of drug sensitivity. Our approach identified five SNPs associated with 10 miRNAs and the expression level of 15 genes, all of which were associated with carboplatin sensitivity. Of particular interest was one SNP (rs11138019), which was associated with the expression of both miR-30d and the gene ABCD2, which were themselves correlated with both carboplatin and cisplatin drug-specific phenotype in the HapMap samples. Functional study found that knocking down ABCD2 in vitro led to increased apoptosis in ovarian cancer cell line SKOV3 after cisplatin treatment. Over-expression of miR-30d in vitro caused a decrease in ABCD2 expression, suggesting a functional relationship between the two. CONCLUSIONS: We developed an integrative approach to the investigation of the genetic and epigenetic basis of human complex traits. Our approach outperformed standard GWAS and provided hints at potential biological function. The relationships between ABCD2 and miR-30d, and ABCD2 and platin sensitivity were experimentally validated, suggesting a functional role of ABCD2 and miR-30d in sensitivity to platinating agents.

Chemotherapeutic-induced apoptosis: a phenotype for pharmacogenomics studies.

AIM: To determine whether cellular apoptosis is a suitable phenotypic trait for pharmacogenomics studies by evaluating caspase 3/7-mediated activity in lymphoblastoid cell lines after treatment with six chemotherapeutic agents: 5'-deoxyfluorouridine, pemetrexed, cytarabine, paclitaxel, carboplatin, and cisplatin. MATERIALS AND METHODS: Using monozygotic twin pair and sibling pair lymphoblastoid cell lines, we identified conditions for measurement of caspase 3/7 activity in lymphoblastoid cell lines. Genome-wide association studies were performed with over 2 million single nucleotide polymorphisms (SNPs) and cisplatin-induced apoptosis in HapMap CEU cell lines (n=77). RESULTS: Although treatment with 5'-deoxyfluorouridine and pemetrexed for up to 24 h resulted in low levels of apoptosis or interindividual variation in caspase-dependent cell death; paclitaxel, cisplatin, carboplatin, and cytarabine treatment for 24 h resulted in 9.4-fold, 9.1-fold, 7.0-fold, and 6.0-fold increases in apoptosis relative to control, respectively. There was a weak correlation between caspase activity and cytotoxicity (r(2)=0.03-0.29) demonstrating that cytotoxicity and apoptosis are two distinct phenotypes that may produce independent genetic associations. Estimated heritability (h(2)) for apoptosis was 0.57 and 0.29 for cytarabine (5 and 40 µmol/l, respectively), 0.22 for paclitaxel (12.5 nmol/l), and 0.34 for cisplatin (5 µmol/l). In the genome-wide association study using the HapMap CEU panel, we identified a significant enrichment of cisplatin-induced cytotoxicity SNPs within the significant cisplatin-induced apoptosis SNPs and an enrichment of expression quantitative trait loci (eQTL). Among these eQTLs, we identified several eQTLs with known function related to apoptosis and/or cytotoxicity. CONCLUSION: Our study identifies apoptosis as a phenotype for pharmacogenomic studies in lymphoblastoid cell lines after treatment with paclitaxel, cisplatin, carboplatin, and cytarabine that may have utility for discovering biomarkers to predict response to certain chemotherapeutics.

Population differences in microRNA expression and biological implications.

Population differences observed for complex traits may be attributed to the combined effect of socioeconomic, environmental, genetic and epigenetic factors. To better understand population differences in complex traits, genome-wide genetic and gene expression differences among ethnic populations have been studied. Here we set out to evaluate population differences in small non-coding RNAs through an evaluation of microRNA (miRNA) baseline expression in HapMap lymphoblastoid cell lines (LCLs) derived from 53 CEU (Utah residents with northern and western European ancestry) and 54 YRI (African from Ibadan, Nigeria). Using the Exiqon miRCURYTM LNA arrays, we found that 16% of all miRNAs evaluated in our study differ significantly between these 2 ethnic groups (pBonferroni corrected< 0.05). Furthermore, we explored the potential biological function of these observed differentially expressed miRNAs by comprehensively examining their effect on the transcriptome and their relationship with cellular sensitivity drug phenotypes. After multiple testing adjustment (false discovery rate (FDR)< 0.1), we found that 55% and 88% of the differentially expressed miRNAs were significantly and inversely correlated with an mRNA expression phenotype in the CEU and YRI samples, respectively. Interestingly, a substantial proportion (64%) of these miRNAs correlated with cellular sensitivity to chemotherapeutic agents (FDR< 0.05). Lastly, upon performing a genome-wide association study between SNPs and miRNA expression, we identified a large number of SNPs exhibiting different allele frequencies that affect the expression of these differentially expressed miRNAs, suggesting the role of genetic variants in mediating the observed population differences.

Utility of patient-derived lymphoblastoid cell lines as an ex vivo capecitabine sensitivity prediction model for breast cancer patients.

Capecitabine is commonly used in treating breast cancer; however, therapeutic response varies among patients and there is no clinically validated model to predict individual outcomes. Here, we investigated whether drug sensitivity quantified in ex vivo patients' blood-derived cell lines can predict response to capecitabine in vivo. Lymphoblastoid cell lines (LCLs) were established from a cohort of metastatic breast cancer patients (n = 53) who were prospectively monitored during treatment with single agent capecitabine at 2000 mg/m2/day. LCLs were treated with increasing concentrations of 5'-DFUR, a major capecitabine metabolite, to assess patients' ex vivo sensitivity to this drug. Subsequently, ex vivo phenotype was compared to observed patient disease response and drug induced-toxicities. We acquired an independent cohort of breast cancer cell lines and LCLs derived from the same donors from ATCC, compared their sensitivity to 5'-DFUR. As seen in the patient population, we observed large inter-individual variability in response to 5'-DFUR treatment in patient-derived LCLs. Patients whose LCLs were more sensitive to 5'-DFUR had a significantly longer median progression free survival (9-month vs 6-month, log rank p-value = 0.017). In addition, this significant positive correlation for 5'-DFUR sensitivity was replicated in an independent cohort of 8 breast cancer cell lines and LCLs derived from the same donor. Our data suggests that at least a portion of the individual sensitivity to capecitabine is shared between germline tissue and tumor tissue. It also supports the utility of patient-derived LCLs as a predictive model for capecitabine treatment efficacy in breast cancer patients.

Integrative "omic" analysis for tamoxifen sensitivity through cell based models.

It has long been observed that tamoxifen sensitivity varies among breast cancer patients. Further, ethnic differences of tamoxifen therapy between Caucasian and African American have also been reported. Since most studies have been focused on Caucasian people, we sought to comprehensively evaluate genetic variants related to tamoxifen therapy in African-derived samples. An integrative "omic" approach developed by our group was used to investigate relationships among endoxifen (an active metabolite of tamoxifen) sensitivity, SNP genotype, mRNA and microRNA expressions in 58 HapMap YRI lymphoblastoid cell lines. We identified 50 SNPs that associate with cellular sensitivity to endoxifen through their effects on 34 genes and 30 microRNA expression. Some of these findings are shared in both Caucasian and African samples, while others are unique in the African samples. Among gene/microRNA that were identified in both ethnic groups, the expression of TRAF1 is also correlated with tamoxifen sensitivity in a collection of 44 breast cancer cell lines. Further, knock-down TRAF1 and over-expression of hsa-let-7i confirmed the roles of hsa-let-7i and TRAF1 in increasing tamoxifen sensitivity in the ZR-75-1 breast cancer cell line. Our integrative omic analysis facilitated the discovery of pharmacogenomic biomarkers that potentially affect tamoxifen sensitivity.

Genetic variation that predicts platinum sensitivity reveals the role of miR-193b* in chemotherapeutic susceptibility.

Platinum agents are the backbone of cancer chemotherapy. Recently, we identified and replicated the role of a single nucleotide polymorphism (SNP, rs1649942) in predicting platinum sensitivity both in vitro and in vivo. Using the CEU samples from the International HapMap Project, we found the same SNP to be a master regulator of multiple gene expression phenotypes, prompting us to investigate whether rs1649942-mediated regulation of miRNAs may in part contribute to variation in platinum sensitivity. To these ends, 60 unrelated HapMap CEU I/II samples were used for our discovery-phase study using high-throughput genome-wide miRNA and gene expression profiling. Examining the relationships among rs1649942, its gene expression targets, genome-wide miRNA expression, and cellular sensitivity to carboplatin and cisplatin, we identified 2 platinum-associated miRNAs (miR-193b* and miR-320) that inhibit the expression of 5 platinum-associated genes (CRIM1, IFIT2, OAS1, KCNMA1, and GRAMD1B). We further replicated the relationship between the expression of miR-193b*, CRIM1, IFIT2, KCNMA1, and GRAMD1B, and platinum sensitivity in a separate HapMap CEU III dataset. We then showed that overexpression of miR-193b* in a randomly selected HapMap cell line results in resistance to both carboplatin and cisplatin. This relationship was also found in 7 ovarian cancer cell lines from NCI60 dataset and confirmed in an OVCAR-3 that overexpression of miR-193b* leads to increased resistance to carboplatin. Our findings highlight a potential mechanism of action for a previously observed genotype-survival outcome association. Further examination of miR-193b* in platinum sensitivity in ovarian cancer is warranted.

Comprehensive evaluation of the contribution of X chromosome genes to platinum sensitivity.

Using a genome-wide gene expression data set generated from Affymetrix GeneChip Human Exon 1.0ST array, we comprehensively surveyed the role of 322 X chromosome gene expression traits on cellular sensitivity to cisplatin and carboplatin. We identified 31 and 17 X chromosome genes whose expression levels are significantly correlated (after multiple testing correction) with sensitivity to carboplatin and cisplatin, respectively, in the combined HapMap CEU (Utah residents with ancestry from northern and western Europe) and YRI (Yoruba in Ibahan, Nigeria) populations (false discovery rate, FDR < 0.05). Of those, 14 overlap for both cisplatin and carboplatin. Using an independent gene expression quantification method, the Illumina Sentrix Human-6 Expression BeadChip, measured on the same HapMap cell lines, we found that 4 and 2 of these genes are significantly associated with carboplatin and cisplatin sensitivity, respectively, in both analyses. Two genes, CTPS2 and DLG3, were identified by both genome-wide gene expression analyses as correlated with cellular sensitivity to both platinating agents. The expression of DLG3 gene was also found to correlate with cellular sensitivity to platinating agents in NCI-60 cancer cell lines. In addition, we evaluated whether the expression of X chromosome genes contributed to the observed differences in sensitivity to the platinums between CEU and YRI-derived cell lines. Of the 34 distinct genes significantly correlated with either carboplatin or cisplatin sensitivity, 14 are differentially expressed (defined as P < 0.05) between CEU and YRI. Thus, sex chromosome genes play a role in cellular sensitivity to platinating agents and differences in the expression level of these genes are an important source of variation that should be included in comprehensive pharmacogenomic studies.

Germline polymorphisms discovered via a cell-based, genome-wide approach predict platinum response in head and neck cancers.

Identifying patients prior to treatment who are more likely to benefit from chemotherapeutic agents or more likely to experience adverse events is an aim of personalized medicine. Pharmacogenomics offers a potential means of achieving this goal through the discovery of predictive germline genetic biomarkers. When applied particularly to the treatment of head and neck cancers, such information could offer significant benefit to patients as a means of potentially reducing morbidity associated with platinum-based chemotherapy. We developed a genome-wide, cell-based approach to identify single nucleotide polymorphisms (SNPs) associated with platinum susceptibility and then evaluated these SNPs as predictors for response and toxicity in head and neck cancer patients treated with platinum-based therapy as part of a phase II clinical trial. Sixty head and neck cancer patients were evaluated. Of 45 genome-wide SNPs examined, we found that 2 SNPs, rs6870861 (P=0.004; false discovery rate [FDR] <0.05) and rs2551038 (P=0.005; FDR <0.05), were associated significantly with overall response to carboplatin-based induction chemotherapy when incorporated into a model along with total carboplatin exposure. Interestingly, these 2 SNPs are associated strongly with the baseline expression of >20 genes (all P ≤10(-4)), and that 2 genes (SLC22A5 and SLCO4C1) are important organic cation/anion transporters known to affect platinum uptake and clearance. Several other SNPs were associated nominally with carboplatin-related hematologic toxicities. These findings demonstrate importantly that a genome-wide, cell-based model can identify novel germline genetic biomarkers of platinum susceptibility, which are replicable in a clinical setting with treated cancer patients and seem clinically meaningful for potentially enabling future personalization of care in such patients.

Platinum sensitivity-related germline polymorphism discovered via a cell-based approach and analysis of its association with outcome in ovarian cancer patients.

PURPOSE: Cell-based approaches were used to identify genetic markers predictive of patients' risk for poor response prior to chemotherapy. EXPERIMENTAL DESIGN: We conducted genome-wide association studies (GWAS) to identify single-nucleotide polymorphisms (SNP) associated with cellular sensitivity to carboplatin through their effects on mRNA expression using International HapMap lymphoblastoid cell lines (LCL) and replicated them in additional LCLs. SNPs passing both stages of the cell-based study were tested for association with progression-free survival (PFS) in patients. Phase 1 validation was based on 377 ovarian cancer patients receiving at least four cycles of carboplatin and paclitaxel from the Australian Ovarian Cancer Study (AOCS). Positive associations were then assessed in phase 2 validation analysis of 1,326 patients from the Ovarian Cancer Association Consortium and The Cancer Genome Atlas. RESULTS: In the initial GWAS, 342 SNPs were associated with carboplatin-induced cytotoxicity, of which 18 unique SNPs were retained after assessing their association with gene expression. One SNP (rs1649942) was replicated in an independent LCL set (Bonferroni adjusted P < 0.05). It was found to be significantly associated with decreased PFS in phase 1 AOCS patients (P(per-allele) = 2 × 10(-2)), with a stronger effect in the subset of women with optimally debulked tumors (P(per-allele) = 4 × 10(-3)). rs1649942 was also associated with poorer overall survival in women with optimally debulked tumors (P(per-allele) = 9 × 10(-3)). However, this SNP was not significant in phase 2 validation analysis with patients from numerous cohorts. CONCLUSION: This study shows the potential of cell-based, genome-wide approaches to identify germline predictors of treatment outcome and highlights the need for extensive validation in patients to assess their clinical effect.