High-density Oligonucleotide Arrays Research Articles

Comparisons of Gene Expression Indexes for Oligonucleotide Arrays

High density oligonucleotide arrays have become a standard research tool to monitor the expression of thousands of genes simultaneously. Affymetrix GeneChip arrays are the most popular. They use short oligonucleotides to probe for genes in an RNA sample. However, important challenges remain in estimating expression level from raw hybridization in tensities on the array. In this paper, we deal with the problem of estimating gene expression based on a statistical model. The present method is like Li and Wong model (2001a), but assumes more generality. More precisely, we show how the model introduced by Li and Wong can be generalized to provide new measure of gene expression. Moreover, we provide a comparison between these two models.

Expression of autophagy-modulating genes in peripheral blood mononuclear cells from familial clustering patients with chronic hepatitis B virus infection.

We previously found that genetic factors are associated with a familial predisposition for developing liver cirrhosis and hepatocellular carcinoma during chronic hepatitis B virus (HBV) infection. Autophagy has been shown to play a role in HBV replication and the course of disease. More than 190 host genes have been identified that modify the process of autophagy, but which of these genes are involved in chronicity of HBV infection and how this occurs remains unclear. Chronic hepatitis B (CHB) patients were recruited to investigate the expression of autophagy-modulating genes in peripheral blood mononuclear cells (PBMCs). mRNA prepared from PBMCs from members of two families with clustering HBV infection, including 11 CHB patients and nine healthy spouses, was hybridized to high-density oligonucleotide arrays. Immunoblot analysis was used to determine the level of autophagy. Of the 192 autophagy-modulating genes, 18 were found to be differently expressed. Of these, 11 displayed decreased expression in CHB patients, while seven displayed increased expression compared to those in healthy controls. Functional analysis showed that these genes are closely involved in initiation, nucleation, elongation of phagophores, formation of autophagosomes, transportation to lysosomes, and the process of degradation. Western blot analysis revealed inhibited autophagy in PBMCs based on decreased lipidation of LC3II. A differential expression profile of autophagy-modulating genes was observed, and decreased autophagy in PBMCs could be closely associated with chronicity of HBV infection, suggesting a novel strategy for the treatment of patients with chronic HBV infection.

The autoimmune targets in IPEX are dominated by gut epithelial proteins

The autoimmune targets in IPEX are dominated by gut epithelial proteins

Transfer of embryos with segmental mosaicism is associated with a significant reduction in live-birth rate

To evaluate the impact of segmental mosaicism on pregnancy outcomes from the transfer of embryos previously designated as euploid. Retrospective cohort analysis. Single, private, high-volume fertility center. Three hundred and twenty-seven women who underwent 377 frozen single euploid embryo transfers. Trophectoderm biopsy of embryos cultured to the blastocyst stage, where all transferred embryos were designated euploid by high-density oligonucleotide array comparative genomic hybridization (aCGH); after ascertaining all outcomes, revaluation of aCGH results for evidence of segmental mosaicism (defined as mosaicism on a portion of a chromosome). Live-birth rate and spontaneous abortion rate. Of the 377 embryos transferred, 357 were euploid with no mosaicism, and 20 embryos had segmental mosaicism. Segmental mosaics had a statistically significantly lower live-birth rate compared with euploid controls (30.0% vs. 53.8%). When controlling for age and day of Trophectoderm biopsy, the odds for live birth after transfer of segmental mosaics were reduced by 66% compared with euploid controls (0.34; 95% confidence interval, 0.13-0.92). The spontaneous abortion rate was statistically significantly higher after transfer of segmental mosaics compared with euploid controls (40.0% vs. 18.2%). Blastocysts with segmental mosaicism have reduced reproductive potential but retain the ability to result in live birth. These results support reporting segmental mosaicism to optimize selection of a single embryo for transfer that will maximize the chance of life birth.

Whole-genome sequencing analysis of CNV using low-coverage and paired-end strategies is efficient and outperforms array-based CNV analysis

BackgroundCopy number variation (CNV) analysis is an integral component of the study of human genomes in both research and clinical settings. Array-based CNV analysis is the current first-tier approach in...

HDR-del: A tool based on Hamming distance for prioritizing pathogenic chromosomal deletions in exome sequencing.

High-density oligonucleotide arrays have widely been used to detect pathogenic chromosomal deletions. In addition to high-density oligonucleotide arrays, programs using whole-exome sequencing have become available for estimating copy-number variations using depth of coverage. Here, we propose a new statistical method, HDR-del, to prioritize pathogenic chromosomal deletions based on Hamming distance in exome sequencing. In vcf (variant call format) files generated from exome sequencing, hemizygous chromosomal deletion regions lack heterozygous variants and lead to apparent long runs of homozygosity (ROH). In our Hamming distance ratio (HDR)-del approach, we calculate the "difference" in heterozygous status between an affected individual and control individuals using the HDR over all candidate chromosomal deletion regions defined as ROH longer than 1Mbp. Using a suitable test statistic, which is expected to be large for a true pathogenic deletion region, we prioritize candidate chromosomal deletion regions based on this statistic. In our approach, we were able to considerably narrow down true pathogenic chromosomal deletion regions, which were confirmed by high-density oligonucleotide arrays in four mitochondrial disease patients. Our HDR-del approach represents an easy method for detecting chromosomal deletions.

Transfer of embryos with segmental mosaicism is associated with a significant reduction in live birth rate compared to euploid controls

Segmental mosaicism is a complex issue with minimal to no data on established thresholds for labeling an embryo euploid versus aneuploid, or on pregnancy outcomes on the transfer of embryos with segmental mosaicism. We aim to evaluate the impact of segmental mosaicism on pregnancy and live birth rates (LBR) from the transfer of embryos previously designated as euploid. Retrospective cohort analysis. All PGS-screened frozen single embryo transfers from January 2015 - December 2015 from a high-volume fertility center were included. Transferred embryos were all designated euploid by high-density oligonucleotide array comparative genomic hybridization (aCGH) PGS by the reporting laboratory. Embryos with segmental mosaicism, designated as euploid, were classified by the reporting laboratory based on the size, location, and percentage of mosaicism of the duplication or deletion. The detailed aCGH results from all transferred euploid embryos were subsequently unmasked to determine the impact of segmental mosaicism on pregnancy outcomes. Chi-square analyses were used to detect differences between euploid versus segmental mosaic embryos. Logistic regression was performed to control for age and day of blastocyst biopsy. A total of 327 women underwent 377 frozen single euploid embryo transfers. After unmasking PGS results, 357 embryos were euploid with no mosaicism and 20 embryos had segmental mosaicism. Average patient age was similar between groups. LBR was significantly lower in the segmental mosaic group (30.0% vs 53.8%, p= 0.04). The odds of a live birth after transfer of embryos with segmental mosaicism was reduced by 66% compared to euploid embryos, controlled for age and day of biopsy (0.34 [95% CI, 0.13- 0.92]) (Table 1). Spontaneous miscarriage, defined as a failed pregnancy following a positive pregnancy test, was significantly higher in the segmental mosaic group compared to euploid controls (40.0% vs. 18.2%, p= 0.04). Blastocysts with segmental mosaicism result in significantly decreased LBR and higher SAB rates compared to euploid controls. These results support reporting all segmental mosaicism, regardless of degree, to optimize embryo selection and LBR.Tabled 1Transfer outcomes of embryos with segmental mosaicism compared to euploid controlsEuploid Embryos N=357Segmental Mosaicism Detected N=20Adjusted Odds Ratio (95% CI)Live Birth N (%)192 (53.8%)6 (30.0%)0.34 (0.13-0.92)Spontaneous Miscarriage N (%)65 (18.2%)8 (40.0%)3.02 (1.18-7.76) Open table in a new tab

Extracting genotype information of Arabidopsis thaliana recombinant inbred lines from transcript profiles established with high-density oligonucleotide arrays.

Polymorphic probes identified via a sequence-based approach are suitable to infer the genotypes of recombinant inbred lines from hybridisation intensities of GeneChip ® transcript profiling experiments. The sequences of the probes of the ATH1 GeneChip® exactly match transcript sequences of the Arabidopsis thaliana reference genome Col-0, whereas nucleotide differences and/or insertions/deletions may be observed for transcripts of other A.thaliana accessions. Individual probes of the GeneChip® that show sequence polymorphisms between different A.thaliana accessions may serve as single-feature polymorphism (SFP) markers, provided that the sequence changes cause differences in hybridisation intensity for the accessions of interest. A sequence-based approach identified features on the high-density oligonucleotide array that showed sequence polymorphisms between A.thaliana accessions Col-0 and C24. Hybridisation intensities of polymorphic probes were extracted from genome-wide transcript profiles of Col-0/C24 and C24/Col-0 recombinant inbred lines and assessed after standardisation via sliding window analyses to identify SFP markers. The genotypes of the recombinant inbred lines were determined with the SFP markers and the resulting data were integrated with information, which had been established previously with single nucleotide polymorphism and insertion/deletion markers, to enrich the linkage map of the Col-0/C24 and C24/Col-0 recombinant inbred populations. Congruence between the molecular marker map and the sequence maps of the A.thaliana Col-0 chromosomes proved the reliability of the genotype information which was deduced from the transcript profiles of the Col-0/C24 and C24/Col-0 recombinant inbred lines.

Tailoring high-density oligonucleotide arrays for transcript profiling of different Arabidopsis thaliana accessions using a sequence-based approach.

Excluding polymorphic probes from GeneChip ® transcript profiling experiments via a sequence-based approach results in improved detection of differentially expressed genes in developing seeds of Arabidopsis thaliana accessions Col-0 and C24. GeneChip® arrays represent a powerful tool for transcript profiling experiments. The ATH1 GeneChip® has been designed based on the sequence of the Arabidopsis thaliana reference genome Col-0, hence the features on the array exactly match the sequences of Col-0 transcripts. In contrast, transcripts of other A. thaliana accessions or related species may show nucleotide differences and/or insertions/deletions when compared to the corresponding Col-0 transcripts, therefore, comparisons of transcript abundance involving different A. thaliana accessions or related species may be compromised for a certain number of transcripts. To tackle this limitation, a sequence-based strategy was developed. Only features on the array that were identical in sequence for the specimen to be compared were considered for transcript profiling. The impact of the proposed strategy was evaluated for transcript profiles that were established for developing seeds of A. thaliana accessions Col-0 and C24.

Comparisons of normalization methods for relative quantitative in real-time polymerase chain reaction

ABSTRACTThe real-time polymerase chain reaction (rtPCR) provides sensitive and accurate quantitative results and becomes a widespread technique in analyzing gene expressions. House-keeping genes are required as references to normalize data of target genes, which may be unstable. This normalization process is similar to the normalization in analyzing high-density oligonucleotide arrays. This article evaluates the feasibility of normalizations for high-density oligonucleotide arrays to normalize data collected in rtPCR experiments. Since data features are different, simulations are used to evaluate the performance of these normalizations to rtPCR data based on five indices. Their feasibilities are illustrated by a rtPCR data.

MFAP5 and TNNC1: Potential markers for predicting occult cervical lymphatic metastasis and prognosis in early stage tongue cancer.

The purpose of this study is to identify candidate genes that could predict prognosis of early-stage tongue squamous cell carcinoma (TSCC) and its occult cervical lymphatic metastasis by large-scale gene expression profiling. Tumor tissue and matched normal mucosa samples were collected from patients with TSCC and analyzed with Affymetrix HTA2.0 high-density oligonucleotide array. Differentially expressed genes in TSCC with cervical lymph node metastasis (CLNM) were further analyzed with Gene Ontology and Kyoto Encyclopedia of Genes and Genomes for their functions and related pathways. A total of 107 differentially expressed genes (p < 0.05) were identified by microarray in TSCC samples with CLNM (n = 6) compared to those without CLNM (n = 6). Genes involved in the cell-matrix adherens junction and migration function including MFAP5, TNNC1, MGP, FBFBP1 and FBXO32 were selected and validated by RT-PCR in TSCC samples (n = 32). Of the five genes, MFAP5 and TNCC1 expressions were further validated by immohistochemistry (n = 61). The significant positive correlation between MFAP5 and TNNC1 expression (p<0.001) was observed. Notably, over-expression of MFAP5 and TNNC1 were correlated with CLNM, metastasis relapse-free survival and overall survival. Our findings indicated that MFAP5 and TNNC1 may be potential markers for predicting occult cervical lymphatic metastasis and prognosis of oral tongue carcinoma.

Functional Genomics Using the Saccharomyces cerevisiae Yeast Deletion Collections

Constructed by a consortium of 16 laboratories, the Saccharomyces genome-wide deletion collections have, for the past decade, provided a powerful, rapid, and inexpensive approach for functional profiling of the yeast genome. Loss-of-function deletion mutants were systematically created using a polymerase chain reaction (PCR)-based gene deletion strategy to generate a start-to-stop codon replacement of each open reading frame by homologous recombination. Each strain carries two molecular barcodes that serve as unique strain identifiers, enabling their growth to be analyzed in parallel and the fitness contribution of each gene to be quantitatively assessed by hybridization to high-density oligonucleotide arrays or through the use of next-generation sequencing technologies. Functional profiling of the deletion collections, using either strain-by-strain or parallel assays, provides an unbiased approach to systematically survey the yeast genome. The Saccharomyces yeast deletion collections have proved immensely powerful in contributing to the understanding of gene function, including functional relationships between genes and genetic pathways in response to diverse genetic and environmental perturbations.

A Comprehensive Genomic Analysis Reveals the Genetic Landscape of Mitochondrial Respiratory Chain Complex Deficiencies.

Mitochondrial disorders have the highest incidence among congenital metabolic disorders characterized by biochemical respiratory chain complex deficiencies. It occurs at a rate of 1 in 5,000 births, and has phenotypic and genetic heterogeneity. Mutations in about 1,500 nuclear encoded mitochondrial proteins may cause mitochondrial dysfunction of energy production and mitochondrial disorders. More than 250 genes that cause mitochondrial disorders have been reported to date. However exact genetic diagnosis for patients still remained largely unknown. To reveal this heterogeneity, we performed comprehensive genomic analyses for 142 patients with childhood-onset mitochondrial respiratory chain complex deficiencies. The approach includes whole mtDNA and exome analyses using high-throughput sequencing, and chromosomal aberration analyses using high-density oligonucleotide arrays. We identified 37 novel mutations in known mitochondrial disease genes and 3 mitochondria-related genes (MRPS23, QRSL1, and PNPLA4) as novel causative genes. We also identified 2 genes known to cause monogenic diseases (MECP2 and TNNI3) and 3 chromosomal aberrations (6q24.3-q25.1, 17p12, and 22q11.21) as causes in this cohort. Our approaches enhance the ability to identify pathogenic gene mutations in patients with biochemically defined mitochondrial respiratory chain complex deficiencies in clinical settings. They also underscore clinical and genetic heterogeneity and will improve patient care of this complex disorder.

The Screening of Genes Sensitive to Long-Term, Low-Level Microwave Exposure and Bioinformatic Analysis of Potential Correlations to Learning and Memory

ObjectiveTo gain a better understanding of gene expression changes in the brain following microwave exposure in mice. This study hopes to reveal mechanisms contributing to microwave-induced learning and memory dysfunction. MethodsMice were exposed to whole body 2100 MHz microwaves with specific absorption rates (SARs) of 0.45 W/kg, 1.8 W/kg, and 3.6 W/kg for 1 hour daily for 8 weeks. Differentially expressing genes in the brains were screened using high-density oligonucleotide arrays, with genes showing more significant differences further confirmed by RT-PCR. ResultsThe gene chip results demonstrated that 41 genes (0.45 W/kg group), 29 genes (1.8 W/kg group), and 219 genes (3.6 W/kg group) were differentially expressed. GO analysis revealed that these differentially expressed genes were primarily involved in metabolic processes, cellular metabolic processes, regulation of biological processes, macromolecular metabolic processes, biosynthetic processes, cellular protein metabolic processes, transport, developmental processes, cellular component organization, etc. KEGG pathway analysis showed that these genes are mainly involved in pathways related to ribosome, Alzheimer's disease, Parkinson's disease, long-term potentiation, Huntington's disease, and Neurotrophin signaling. Construction of a protein interaction network identified several important regulatory genes including synbindin (sbdn), Crystallin (CryaB), PPP1CA, Ywhaq, Psap, Psmb1, Pcbp2, etc., which play important roles in the processes of learning and memory. ConclusionLong-term, low-level microwave exposure may inhibit learning and memory by affecting protein and energy metabolic processes and signaling pathways relating to neurological functions or diseases.

Mouse model of proximal colon-specific tumorigenesis driven by microsatellite instability-induced Cre-mediated inactivation of Apc and activation of Kras.

KRAS gene mutations are found in 40-50% of colorectal cancer cases, but their functional contribution is not fully understood. To address this issue, we generated genetically engineered mice with colon tumors expressing an oncogenic Kras(G12D) allele in the context of the Adenomatous polyposis coli (Apc) deficiency to compare them to tumors harboring Apc deficiency alone. CDX2P9.5-G22Cre (referred to as G22Cre) mice showing inducible Cre recombinase transgene expression in the proximal colon controlled under the CDX2 gene promoter were intercrossed with Apc (flox/flox) mice and LSL-Kras (G12D) mice carrying loxP-flanked Apc and Lox-Stop-Lox oncogenic Kras(G12D) alleles, respectively, to generate G22Cre; Apc(flox/flox); Kras(G12D) and G22Cre; Apc(flox/flox); KrasWT mice. Gene expression profiles of the tumors were analyzed using high-density oligonucleotide arrays. Morphologically, minimal difference in proximal colon tumor was observed between the two mouse models. Consistent with previous findings in vitro, Glut1 transcript and protein expression was up-regulated in the tumors of G22Cre;Apc (flox/flox) ; Kras(G12D) mice. Immunohistochemical staining analysis revealed that GLUT1 protein expression correlated with KRAS mutations in human colorectal cancer. Microarray analysis identified 11 candidate genes upregulated more than fivefold and quantitative PCR analysis confirmed that Aqp8, Ttr, Qpct, and Slc26a3 genes were upregulated 3.7- to 30.2-fold in tumors with mutant Kras. These results demonstrated the validity of the G22Cre; Apc(flox/flox) ;Kras (G12D) mice as a new mouse model with oncogenic Kras activation. We believe that this model can facilitate efforts to define novel factors that contribute to the pathogenesis of human colorectal cancer with KRAS mutations.

Abstract 2853: Investigating DNAPK as a biomarker and a novel therapeutic target in aggressive prostate cancer

Abstract Background/purpose: Despite recent advances in the development of androgen deprivation therapies (ADT), prostate cancer (PCa) remains the second leading cause of cancer death in American men. Thus, there is a critical need to improve treatment outcomes for men with aggressive PCa. The purpose of this study was to identify and functionally characterize kinases that can serve as both prognostic biomarkers and therapeutic targets in PCa. Methods: To identify kinases associated with metastatic progression of PCa, we utilized high-density oligonucleotide arrays to interrogate the expression of all kinases in 545 prostatectomy samples, obtained from high-risk patients with long-term clinical follow-up (&amp;gt;10 years). We ranked all kinases by their fold change in expression between PCas that subsequently metastasized versus those that did not. We then characterized the top ranked kinase in preclinical PCa models with in vitro mechanistic experiments and in vivo therapeutic studies. Results: We nominated DNA-dependent protein kinase (DNAPK) as the top kinase associated with metastatic PCa progression. Interrogation of pathways associated with DNAPK knockdown in PC3, DU145, VCaP and C4-2B, identified canonical Wnt signaling as the top DNAPK-activated pathway. In PCa cells (PC3, LNCaP C.S. {LNCaP grown under charcoal-stripped serum condition}, C4-2B and LNCaP-AR) DNAPK inhibition significantly abrogates Wnt signaling and abolishes cell proliferation, migration and invasion. Interestingly, we found that DNAPK directly modulates the Wnt signaling independent of the androgen receptor (AR), highlighting the therapeutic potential of DNAPK inhibition in cancers resistant to ADT. Indeed, in VCaP xenograft models, pharmacologic inhibition of DNAPK with NU7441 demonstrated 3.4 fold (p&amp;lt;0.01) reduction of tumor growth at non-toxic doses. Using an independent prostate cancer tissue cohort, we validate DNAPK expression as a biomarker for metastatic disease progression (p = 6.4 e−06, hazard ratio = 2). Conclusion: In summary, we nominate DNAPK as a potential biomarker of disease progression and as a novel therapeutic target in aggressive prostate cancer. Our data suggests that DNAPK mediates PCa progression by upregulating Wnt signaling, and demonstrates that DNAPK inhibition results in significant tumor responses in PCa xenografts. Citation Format: Vishal Kothari, Jonathan F. Goodwin, Shuang Zhao, Elai Davicioni, Jeffrey R. Karnes, Robert B. Den, Rohit Mehra, Karen E. Knudsen, Felix Y. Feng. Investigating DNAPK as a biomarker and a novel therapeutic target in aggressive prostate cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2853. doi:10.1158/1538-7445.AM2015-2853

Alu-mediated diverse and complex pathogenic copy-number variants within human chromosome 17 at p13.3.

Alu repetitive elements are known to be major contributors to genome instability by generating Alu-mediated copy-number variants (CNVs). Most of the reported Alu-mediated CNVs are simple deletions and duplications, and the mechanism underlying Alu-Alu-mediated rearrangement has been attributed to non-allelic homologous recombination (NAHR). Chromosome 17 at the p13.3 genomic region lacks extensive low-copy repeat architecture; however, it is highly enriched for Alu repetitive elements, with a fraction of 30% of total sequence annotated in the human reference genome, compared with the 10% genome-wide and 18% on chromosome 17. We conducted mechanistic studies of the 17p13.3 CNVs by performing high-density oligonucleotide array comparative genomic hybridization, specifically interrogating the 17p13.3 region with ∼150 bp per probe density; CNV breakpoint junctions were mapped to nucleotide resolution by polymerase chain reaction and Sanger sequencing. Studied rearrangements include 5 interstitial deletions, 14 tandem duplications, 7 terminal deletions and 13 complex genomic rearrangements (CGRs). Within the 17p13.3 region, Alu-Alu-mediated rearrangements were identified in 80% of the interstitial deletions, 46% of the tandem duplications and 50% of the CGRs, indicating that this mechanism was a major contributor for formation of breakpoint junctions. Our studies suggest that Alu repetitive elements facilitate formation of non-recurrent CNVs, CGRs and other structural aberrations of chromosome 17 at p13.3. The common observation of Alu-mediated rearrangement in CGRs and breakpoint junction sequences analysis further demonstrates that this type of mechanism is unlikely attributed to NAHR, but rather may be due to a recombination-coupled DNA replicative repair process.

Involvement of TG-interacting factor in microglial activation during experimental traumatic brain injury.

Traumatic brain injury (TBI) is a complex injury involving several physiological alterations, potentially leading to neurological impairment. Previous mouse studies using high-density oligonucleotide array analysis have confirmed the upregulation of transforming growth-interacting factor (TGIF) mRNA in TBI. TGIF is a transcriptional corepressor of transforming growth factor beta (TGF-β) signaling which plays a protective role in TBI. However, the functional roles of TGIF in TBI are not well understood. In this study, we used confocal microscopy after immunofluorescence staining to demonstrate the increase of TGIF levels in the activated microglia of the pericontusional cortex of rats with TBI. Intracerebral knockdown of TGIF in the pericontusional cortex significantly downregulated TGIF expression, attenuated microglial activation, reduced the volume of damaged brain tissue, and facilitated recovery of limb motor function. Collectively, our results indicate that TGIF is involved in TBI-induced microglial activation, resulting in secondary brain injury and motor dysfunction. This study investigated the roles of transforming growth-interacting factor (TGIF) in a traumatic brain injury (TBI)-rat model. We demonstrated the increase of TGIF levels in the activated microglia of the pericontusional cortex of rats with TBI. Intracerebral knockdown of TGIF in the pericontusional cortex of the TBI rats significantly attenuated micoglial activation, reduced the volume of damaged brain tissue, and facilitated recovery of limb motor function. We suggest that inhibition of TGIF might provide a promising therapeutic strategy for TBI.

Efficacy of ribavirin against malignant glioma cell lines.

Ribavirin (1-β-D-ribofuranosy-1,2,4-triazole-3-carboxamide) has been widely administered as an antiviral agent against RNA and DNA viruses. Ribavirin, in combination with interferon, has predominantly been applied in the treatment of the hepatitis C virus infection and its potential antitumor efficacy has recently become a point of interest. The aim of the present study was to evaluate the effect of ribavirin on the growth of malignant glioma cells, to identify novel predictive genes in malignant glioma cells (by analyzing gene expression profiles) and to assess the influence of ribavirin on the cell cycle of malignant glioma cells. The present study evaluated the antitumor efficacy of ribavirin against various malignant glioma cell lines (A-172, AM-38, T98G, U-87MG, U-138MG, U-251MG and YH-13). After culturing the cells in ribavirin-containing culture medium (final concentration, 0–1,000 μM) for 72 h, the viable proliferated cells were harvested and counted. The half maximal inhibitory concentration of ribavirin, with regard to the growth of the malignant glioma cell lines, was determined from the concentration of ribavirin required for 50% growth inhibition in comparison to the untreated control cells. Furthermore, the current study identified the genes in which the gene expression levels correlated with the ribavirin sensitivity of the malignant glioma cells lines, using a high-density oligonucleotide array. Finally, cell cycle analysis was performed on the U-87MG cell line. It was identified that ribavirin inhibited the growth of all of the malignant glioma cell lines in a dose-dependent manner, although the ribavirin sensitivity varied between each cell line. Of the extracted genes, PDGFRA demonstrated the strongest positive correlation between gene expression level and ribavirin sensitivity. Cell cycle analysis of the U-87MG cell line demonstrated that ribavirin treatment induces G0/G1 arrest and thus may be an effective agent for inhibiting malignant glioma cell growth. Therefore, the results of the current study indicate that ribavirin may have potential as a therapeutic agent in the treatment of malignant gliomas.

Identification and Validation of Genetic Variants that Influence Transcription Factor and Cell Signaling Protein Levels

Many genetic variants associated with human disease have been found to be associated with alterations in mRNA expression. Although it is commonly assumed that mRNA expression changes will lead to consequent changes in protein levels, methodological challenges have limited our ability to test the degree to which this assumption holds true. Here, we further developed the micro-western array approach and globally examined relationships between human genetic variation and cellular protein levels. We collected more than 250,000 protein level measurements comprising 441 transcription factor and signaling protein isoforms across 68 Yoruba (YRI) HapMap lymphoblastoid cell lines (LCLs) and identified 12 cis and 160 trans protein level QTLs (pQTLs) at a false discovery rate (FDR) of 20%. Whereas up to two thirds of cis mRNA expression QTLs (eQTLs) were also pQTLs, many pQTLs were not associated with mRNA expression. Notably, we replicated and functionally validated a trans pQTL relationship between the KARS lysyl-tRNA synthetase locus and levels of the DIDO1 protein. This study demonstrates proof of concept in applying an antibody-based microarray approach to iteratively measure the levels of human proteins and relate these levels to human genome variation and other genomic data sets. Our results suggest that protein-based mechanisms might functionally buffer genetic alterations that influence mRNA expression levels and that pQTLs might contribute phenotypic diversity to a human population independently of influences on mRNA expression.