Combined immune checkpoint inhibition with durvalumab and tremelimumab with and without radiofrequency ablation in patients with advanced biliary tract carcinoma.

BACKGROUND: Current standard of care for advanced biliary tract cancer (BTC) is gemcitabine, cisplatin plus anti-PD1/PD-L1, but response rates are modest. The purpose of this study was to explore the efficacy and safety of durvalumab (anti-PD-L1) and tremelimumab (anti-CTLA-4), with and without an interventional radiology (IR) procedure in advanced BTC. METHODS: Eligible patients with advanced BTC who had received or refused at least one prior line of systemic therapy were treated with tremelimumab and durvalumab for four combined doses followed by monthly durvalumab alone with and without an IR procedure until the progression of disease or unacceptable toxicity. Objective response was assessed through CT or MRI by Response Evaluation Criteria in Solid Tumors (RECIST, version 1.1) every 8 weeks. Adverse events (AEs) were recorded and managed. The primary endpoint was 6-month progression-free survival (PFS). RESULTS: Twenty-three patients with advanced BTC were enrolled; 17 patients were assigned to treatment with durvalumab and tremelimumab (Durva/Treme); and 6 patients were treated with the combination of durvalumab, tremelimumab plus IR procedure (Durva/Treme + IR). The best clinical responses in the Durva/Treme arm were partial response (n = 1), stable disease (n = 5), progressive disease (n = 5), and in the Durva/Treme + IR arm: partial response (n = 0), stable disease (n = 3), progressive disease (n = 3). The median PFS was 2.2 months (95% CI: 1.3-3.1 months) in the Durva/Treme arm and 2.9 months (95% CI: 1.9-4.7 months) in the Durva/Treme + IR arm (p = 0.27). The median OS was 5.1 months (95% CI: 2.5-6.9 months) in the Durva/Treme arm and 5.8 months (95% CI: 2.9-40.1 months) in the Durva/Treme + IR arm (p = 0.31). The majority of AEs were grades 1-2. CONCLUSION: Durva/Treme and Durva/Treme + IR showed similar efficacy. With a manageable safety profile. Larger studies are needed to fully characterize the efficacy of Durva/Treme ± IR in advanced BTC.

Development of a bicistronic anti-CD19/CD20 CAR construct including abrogation of unexpected nucleic acid sequence deletions.

To address CD19 loss from lymphoma after anti-CD19 chimeric antigen receptor (CAR) T cell therapy, we designed a bicistronic construct encoding an anti-CD19 CAR and an anti-CD20 CAR. We detected deletions from the expected bicistronic construct sequence in a minority of transcripts by mRNA sequencing. Loss of bicistronic construct transgene DNA was also detected. Deletions of sequence were present at much higher frequencies in transduced T cell mRNA versus gamma-retroviral vector RNA. We concluded that these deletions were caused by intramolecular template switching of the reverse transcriptase enzyme during reverse transcription of gamma-retroviral vector RNA into transgene DNA of transduced T cells. Intramolecular template switching was driven by repeated regions of highly similar nucleic acid sequence within CAR sequences. We optimized the sequence of the bicistronic CAR construct to reduce repeated regions of highly similar sequences. This optimization nearly eliminated sequence deletions. This work shows that repeated regions of highly similar nucleic acid sequence must be avoided in complex CAR constructs. We further optimized the bicistronic construct by lengthening the linker of the anti-CD20 single-chain variable fragment. This modification increased CD20-specific interleukin-2 release and reduced CD20-specific activation-induced cell death. We selected an optimized anti-CD19/CD20 bicistronic construct for clinical development.

CD47-Dependent Regulation of Immune Checkpoint Gene Expression and MYCN mRNA Splicing in Murine CD8 and Jurkat T Cells.

Elevated expression of CD47 in some cancers is associated with poor survival related to its function as an innate immune checkpoint when expressed on tumor cells. In contrast, elevated CD47 expression in cutaneous melanomas is associated with improved survival. Previous studies implicated protective functions of CD47 expressed by immune cells in the melanoma tumor microenvironment. RNA sequencing analysis of responses induced by CD3 and CD28 engagement on wild type and CD47-deficient Jurkat T lymphoblast cells identified additional regulators of T cell function that were also CD47-dependent in mouse CD8 T cells. MYCN mRNA expression was upregulated in CD47-deficient cells but downregulated in CD47-deficient cells following activation. CD47 also regulated alternative splicing that produces two N-MYC isoforms. The CD47 ligand thrombospondin-1 inhibited expression of these MYCN mRNA isoforms, as well as induction of the oncogenic decoy MYCN opposite strand (MYCNOS) RNA during T cell activation. Analysis of mRNA expression data for melanomas in The Cancer Genome Atlas identified a significant coexpression of MYCN with CD47 and known regulators of CD8 T cell function. Thrombospondin-1 inhibited the induction of TIGIT, CD40LG, and MCL1 mRNAs following T cell activation in vitro. Increased mRNA expression of these T cell transcripts and MYCN in melanomas was associated with improved overall survival.

The epigenetic regulator LSH maintains fork protection and genomic stability via MacroH2A deposition and RAD51 filament formation.

The Immunodeficiency Centromeric Instability Facial Anomalies (ICF) 4 syndrome is caused by mutations in LSH/HELLS, a chromatin remodeler promoting incorporation of histone variant macroH2A. Here, we demonstrate that LSH depletion results in degradation of nascent DNA at stalled replication forks and the generation of genomic instability. The protection of stalled forks is mediated by macroH2A, whose knockdown mimics LSH depletion and whose overexpression rescues nascent DNA degradation. LSH or macroH2A deficiency leads to an impairment of RAD51 loading, a factor that prevents MRE11 and EXO1 mediated nascent DNA degradation. The defect in RAD51 loading is linked to a disbalance of BRCA1 and 53BP1 accumulation at stalled forks. This is associated with perturbed histone modifications, including abnormal H4K20 methylation that is critical for BRCA1 enrichment and 53BP1 exclusion. Altogether, our results illuminate the mechanism underlying a human syndrome and reveal a critical role of LSH mediated chromatin remodeling in genomic stability.

LSH mediates gene repression through macroH2A deposition.

The human Immunodeficiency Centromeric Instability Facial Anomalies (ICF) 4 syndrome is a severe disease with increased mortality caused by mutation in the LSH gene. Although LSH belongs to a family of chromatin remodeling proteins, it remains unknown how LSH mediates its function on chromatin in vivo. Here, we use chemical-induced proximity to rapidly recruit LSH to an engineered locus and find that LSH specifically induces macroH2A1.2 and macroH2A2 deposition in an ATP-dependent manner. Tethering of LSH induces transcriptional repression and silencing is dependent on macroH2A deposition. Loss of LSH decreases macroH2A enrichment at repeat sequences and results in transcriptional reactivation. Likewise, reduction of macroH2A by siRNA interference mimicks transcriptional reactivation. ChIP-seq analysis confirmed that LSH is a major regulator of genome-wide macroH2A distribution. Tethering of ICF4 mutations fails to induce macroH2A deposition and ICF4 patient cells display reduced macroH2A deposition and transcriptional reactivation supporting a pathogenic role for altered marcoH2A deposition. We propose that LSH is a major chromatin modulator of the histone variant macroH2A and that its ability to insert marcoH2A into chromatin and transcriptionally silence is disturbed in the ICF4 syndrome.

Lsh/HELLS is required for B lymphocyte development and immunoglobulin class switch recombination.

Mutation of HELLS (Helicase, Lymphoid-Specific)/Lsh in human DNA causes a severe immunodeficiency syndrome, but the nature of the defect remains unknown. We assessed here the role of Lsh in hematopoiesis using conditional Lsh knockout mice with expression of Mx1 or Vav Cre-recombinase. Bone marrow transplantation studies revealed that Lsh depletion in hematopoietic stem cells severely reduced B cell numbers and impaired B cell development in a hematopoietic cell-autonomous manner. Lsh-deficient mice without bone marrow transplantation exhibited lower Ig levels in vivo compared to controls despite normal peripheral B cell numbers. Purified B lymphocytes proliferated normally but produced less immunoglobulins in response to in vitro stimulation, indicating a reduced capacity to undergo class switch recombination (CSR). Analysis of germline transcripts, examination of double-stranded breaks using biotin-labeling DNA break assay, and End-seq analysis indicated that the initiation of the recombination process was unscathed. In contrast, digestion-circularization PCR analysis and high-throughput sequencing analyses of CSR junctions and a chromosomal break repair assay indicated an impaired ability of the canonical end-joining pathway in Lsh-deficient B cells. Our data suggest a hematopoietic cell-intrinsic role of Lsh in B cell development and in CSR providing a potential target for immunodeficiency therapy.

The chromatin remodeling protein Lsh alters nucleosome occupancy at putative enhancers and modulates binding of lineage specific transcription factors.

Dynamic regulation of chromatin accessibility is a key feature of cellular differentiation during embryogenesis, but the precise factors that control access to chromatin remain largely unknown. Lsh/HELLS is critical for normal development and mutations of Lsh in human cause the ICF (Immune deficiency, Centromeric instability, Facial anomalies) syndrome, a severe immune disorder with multiple organ deficiencies. We report here that Lsh, previously known to regulate DNA methylation level, has a genome wide chromatin remodeling function. Using micrococcal nuclease (MNase)-seq analysis, we demonstrate that Lsh protects MNase accessibility at transcriptional regulatory regions characterized by DNase I hypersensitivity and certain histone 3 (H3) tail modifications associated with enhancers. Using an auxin-inducible degron system, allowing proteolytical degradation of Lsh, we show that Lsh mediated changes in nucleosome occupancy are independent of DNA methylation level and are characterized by reduced H3 occupancy. While Lsh mediated nucleosome occupancy prevents binding sites for transcription factors in wild type cells, depletion of Lsh leads to an increase in binding of ectopically expressed tissue specific transcription factors to their respective binding sites. Our data suggests that Lsh mediated chromatin remodeling can modulate nucleosome positioning at a subset of putative enhancers contributing to the preservation of cellular identity through regulation of accessibility.

Genomic Heterogeneity Within Individual Prostate Cancer Foci Impacts Predictive Biomarkers of Targeted Therapy.

BACKGROUND: Most lethal prostate cancers progress from relapse of aggressive primary disease. Recently, the most significant advances in survival benefit from systemic therapy have come from moving the administration of therapy to an earlier disease state. There is movement toward using biomarkers from the intraprostatic index lesion to guide early systemic therapy. OBJECTIVE: To determine the genomic heterogeneity, including the heterogeneity of predictive biomarkers, within the index focus of treatment-naïve prostate cancer. DESIGN, SETTING, AND PARTICIPANTS: Ten patients with treatment-naïve prostate cancer underwent prostatectomy. DNA was extracted from 70 spatially distinct regions of the 10 index foci. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Single nucleotide mutations, small indels, and copy number changes were identified. Intrafocal genomic heterogeneity and heterogeneity of alterations that predict response to therapy was determined. RESULTS AND LIMITATIONS: Exome sequencing and copy number estimates demonstrate branched evolution with >75% of point mutations being subclonal, including numerous pathways associated with castrate-resistant prostate cancer. Seven of 10 patients harbor alterations in one of five genes that predict response to targeted therapies with survival benefit in prostate cancer. Within biomarker-positive cases, 25% of intraprostatic regions are biomarker negative, with discordance between intraprostatic regions and lymph node metastases. CONCLUSIONS: Treatment-naïve, nonmetastatic prostate cancer has marked intrafocal heterogeneity. Numerous alterations in pathways associated with castration-resistant prostate cancer are present in subclonal populations, including biomarkers predictive of response to targeted therapy. PATIENT SUMMARY: Untreated patients' tumors have alterations that predict response to targeted therapies, but the presence of a biomarker is dependent on what region of the tumor was evaluated.

Chromatic: WebAssembly-Based Cancer Genome Viewer.

Chromatic is a novel web-browser tool that enables researchers to visually inspect genomic variations identified through next-generation sequencing of cancer data sets to determine whether such calls are, in fact, valid. It is the first cancer bioinformatics tool developed using WebAssembly technology, which comprises a portable, low-level byte code format that provides for the rapid execution of programs within supported web browsers. It has been designed expressly for ease of use by scientists without extensive expertise in bioinformatics.

Transcriptional alterations in hereditary and sporadic nonfunctioning pancreatic neuroendocrine tumors according to genotype.

BACKGROUND: Nonfunctioning pancreatic neuroendocrine tumors (NFPanNETs) may be sporadic or inherited because of germline mutations associated with von Hippel-Lindau disease (VHL) or multiple endocrine neoplasia type 1 (MEN1). The clinical behavior of NFPanNETs is difficult to predict, even in tumors of the same stage and grade. The authors analyzed genotype-specific patterns of transcriptional messenger RNA (mRNA) levels of NFPanNETs to understand the molecular features that determine PanNET phenotype. METHODS: Thirty-two samples were included for genome-wide mRNA gene expression analysis (9 VHL-associated, 10 MEN1-associated, and 9 sporadic NFPanNETs and 4 purified normal islet cell [NIC] samples). Validation of genes was performed by real-time polymerase chain reaction analysis and immunohistochemistry. Gene expression profiles were analyzed by tumor genotype, and pathway analysis was curated. RESULTS: Consensus clustering of mRNA expression revealed separate clustering of NICs, VHL-associated NFPanNETs, and MEN1-associated NFPanNETs; whereas some sporadic tumors clustered with MEN1. Four of 5 MEN1-like sporadic PanNET subtypes had loss of heterozygosity at the MEN1 gene locus. Pathway analysis demonstrated subtype-specific pathway activation, comprising angiogenesis and immune response in VHL; neuronal development in MEN1; protein ubiquitination in the new MEN1/sporadic subtype; and cytokinesis and cilium/microtubule development in sporadic NFPanNETs. Among many genes, platelet-derived growth factor receptor ß (PDGFRB), lymphoid enhancer-binding factor-1 (Lef-1), cyclin-dependent kinase 4 (CDK4), and CDK6 were upregulated in VHL or MEN1 NFPanNETs, providing potential subtype-specific treatment targets. CONCLUSIONS: Distinct mRNA expression patterns were identified in sporadic-associated, VHL-associated, and MEN1-associated NFPanNETs. The current results uncover new pathways involved in NFPanNETs that are subtype-specific and provide potential new diagnostic or therapeutic targets based on tumor subtype. Cancer 2018;124:636-47. © 2017 American Cancer Society.

Alview: Portable Software for Viewing Sequence Reads in BAM Formatted Files.

The name Alview is a contraction of the term Alignment Viewer. Alview is a compiled to native architecture software tool for visualizing the alignment of sequencing data. Inputs are files of short-read sequences aligned to a reference genome in the SAM/BAM format and files containing reference genome data. Outputs are visualizations of these aligned short reads. Alview is written in portable C with optional graphical user interface (GUI) code written in C, C++, and Objective-C. The application can run in three different ways: as a web server, as a command line tool, or as a native, GUI program. Alview is compatible with Microsoft Windows, Linux, and Apple OS X. It is available as a web demo at https://cgwb.nci.nih.gov/cgi-bin/alview. The source code and Windows/Mac/Linux executables are available via https://github.com/NCIP/alview.

Generating a focused view of disease ontology cancer terms for pan-cancer data integration and analysis.

Bio-ontologies provide terminologies for the scientific community to describe biomedical entities in a standardized manner. There are multiple initiatives that are developing biomedical terminologies for the purpose of providing better annotation, data integration and mining capabilities. Terminology resources devised for multiple purposes inherently diverge in content and structure. A major issue of biomedical data integration is the development of overlapping terms, ambiguous classifications and inconsistencies represented across databases and publications. The disease ontology (DO) was developed over the past decade to address data integration, standardization and annotation issues for human disease data. We have established a DO cancer project to be a focused view of cancer terms within the DO. The DO cancer project mapped 386 cancer terms from the Catalogue of Somatic Mutations in Cancer (COSMIC), The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium, Therapeutically Applicable Research to Generate Effective Treatments, Integrative Oncogenomics and the Early Detection Research Network into a cohesive set of 187 DO terms represented by 63 top-level DO cancer terms. For example, the COSMIC term 'kidney, NS, carcinoma, clear_cell_renal_cell_carcinoma' and TCGA term 'Kidney renal clear cell carcinoma' were both grouped to the term 'Disease Ontology Identification (DOID):4467 / renal clear cell carcinoma' which was mapped to the TopNodes_DOcancerslim term 'DOID:263 / kidney cancer'. Mapping of diverse cancer terms to DO and the use of top level terms (DO slims) will enable pan-cancer analysis across datasets generated from any of the cancer term sources where pan-cancer means including or relating to all or multiple types of cancer. The terms can be browsed from the DO web site (http://www.disease-ontology.org) and downloaded from the DO's Apache Subversion or GitHub repositories. Database URL: http://www.disease-ontology.org

Genome-wide transcriptional sequencing identifies novel mutations in metabolic genes in human hepatocellular carcinoma.

We report on next-generation transcriptome sequencing results of three human hepatocellular carcinoma tumor/tumor-adjacent pairs. This analysis robustly examined ∼12,000 genes for both expression differences and molecular alterations. We observed 4,513 and 1,182 genes demonstrating 2-fold or greater increase or decrease in expression relative to their normal, respectively. Network analysis of expression data identified the Aurora B signaling, FOXM1 transcription factor network and Wnt signaling pathways pairs being altered in HCC. We validated as differential gene expression findings in a large data set containing of 434 liver normal/tumor sample pairs. In addition to known driver mutations in TP53 and CTNNB1, our mutation analysis identified non-synonymous mutations in genes implicated in metabolic diseases, i.e. diabetes and obesity: IRS1, HMGCS1, ATP8B1, PRMT6 and CLU, suggesting a common molecular etiology for HCC of alternative pathogenic origin.

Discovery analysis of TCGA data reveals association between germline genotype and survival in ovarian cancer patients.

BACKGROUND: Ovarian cancer remains a significant public health burden, with the highest mortality rate of all the gynecological cancers. This is attributable to the late stage at which the majority of ovarian cancers are diagnosed, coupled with the low and variable response of advanced tumors to standard chemotherapies. To date, clinically useful predictors of treatment response remain lacking. Identifying the genetic determinants of ovarian cancer survival and treatment response is crucial to the development of prognostic biomarkers and personalized therapies that may improve outcomes for the late-stage patients who comprise the majority of cases. METHODS: To identify constitutional genetic variations contributing to ovarian cancer mortality, we systematically investigated associations between germline polymorphisms and ovarian cancer survival using data from The Cancer Genome Atlas Project (TCGA). Using stage-stratified Cox proportional hazards regression, we examined >650,000 SNP loci for association with survival. We additionally examined whether the association of significant SNPs with survival was modified by somatic alterations. RESULTS: Germline polymorphisms at rs4934282 (AGAP11/C10orf116) and rs1857623 (DNAH14) were associated with stage-adjusted survival (p= 1.12e-07 and 1.80e-07, FDR q= 1.2e-04 and 2.4e-04, respectively). A third SNP, rs4869 (C10orf116), was additionally identified as significant in the exome sequencing data; it is in near-perfect LD with rs4934282. The associations with survival remained significant when somatic alterations. CONCLUSIONS: Discovery analysis of TCGA data reveals germline genetic variations that may play a role in ovarian cancer survival even among late-stage cases. The significant loci are located near genes previously reported as having a possible relationship to platinum and taxol response. Because the variant alleles at the significant loci are common (frequencies for rs4934282 A/C alleles = 0.54/0.46, respectively; rs1857623 A/G alleles = 0.55/0.45, respectively) and germline variants can be assayed noninvasively, our findings provide potential targets for further exploration as prognostic biomarkers and individualized therapies.

Genetic alterations activating kinase and cytokine receptor signaling in high-risk acute lymphoblastic leukemia.

Genomic profiling has identified a subtype of high-risk B-progenitor acute lymphoblastic leukemia (B-ALL) with alteration of IKZF1, a gene expression profile similar to BCR-ABL1-positive ALL and poor outcome (Ph-like ALL). The genetic alterations that activate kinase signaling in Ph-like ALL are poorly understood. We performed transcriptome and whole genome sequencing on 15 cases of Ph-like ALL and identified rearrangements involving ABL1, JAK2, PDGFRB, CRLF2, and EPOR, activating mutations of IL7R and FLT3, and deletion of SH2B3, which encodes the JAK2-negative regulator LNK. Importantly, several of these alterations induce transformation that is attenuated with tyrosine kinase inhibitors, suggesting the treatment outcome of these patients may be improved with targeted therapy.

Bambino: a variant detector and alignment viewer for next-generation sequencing data in the SAM/BAM format.

SUMMARY: Bambino is a variant detector and graphical alignment viewer for next-generation sequencing data in the SAM/BAM format, which is capable of pooling data from multiple source files. The variant detector takes advantage of SAM-specific annotations, and produces detailed output suitable for genotyping and identification of somatic mutations. The assembly viewer can display reads in the context of either a user-provided or automatically generated reference sequence, retrieve genome annotation features from a UCSC genome annotation database, display histograms of non-reference allele frequencies, and predict protein-coding changes caused by SNPs. AVAILABILITY: Bambino is written in platform-independent Java and available from https://cgwb.nci.nih.gov/goldenPath/bamview/documentation/index.html, along with documentation and example data. Bambino may be launched online via Java Web Start or downloaded and run locally.

Genetic variations at loci involved in the immune response are risk factors for hepatocellular carcinoma.

UNLABELLED: Primary liver cancer is the third most common cause of cancer-related death worldwide, with a rising incidence in Western countries. Little is known about the genetic etiology of this disease. To identify genetic factors associated with hepatocellular carcinoma (HCC) and liver cirrhosis (LC), we conducted a comprehensive, genome-wide variation analysis in a population of unrelated Asian individuals. Copy number variation (CNV) and single nucleotide polymorphisms (SNPs) were assayed in peripheral blood with the high-density Affymetrix SNP6.0 microarray platform. We used a two-stage discovery and replication design to control for overfitting and to validate observed results. We identified a strong association with CNV at the T-cell receptor gamma and alpha loci (P < 1 × 10(-15)) in HCC cases when contrasted with controls. This variation appears to be somatic in origin, reflecting differences between T-cell receptor processing in lymphocytes from individuals with liver disease and healthy individuals that is not attributable to chronic hepatitis virus infection. Analysis of constitutional variation identified three susceptibility loci including the class II MHC complex, whose protein products present antigen to T-cell receptors and mediate immune surveillance. Statistical analysis of biologic networks identified variation in the "antigen presentation and processing" pathway as being highly significantly associated with HCC (P = 1 × 10(-11)). SNP analysis identified two variants whose allele frequencies differ significantly between HCC and LC. One of these (P = 1.74 × 10(-12)) lies in the PTEN homolog TPTE2. CONCLUSION: Combined analysis of CNV, individual SNPs, and pathways suggest that HCC susceptibility is mediated by germline factors affecting the immune response and differences in T-cell receptor processing.

Systematic analysis of genetic alterations in tumors using Cancer Genome WorkBench (CGWB).

Systematic investigations of genetic changes in tumors are expected to lead to greatly improved understanding of cancer etiology. To meet the analytical challenges presented by such studies, we developed the Cancer Genome WorkBench (http://cgwb.nci.nih.gov), the first computational platform to integrate clinical tumor mutation profiles with the reference human genome. A novel heuristic algorithm, IndelDetector, was developed to automatically identify insertion/deletion (indel) polymorphisms as well as indel somatic mutations with high sensitivity and accuracy. It was incorporated into an automated pipeline that detects genetic alterations and annotates their effects on protein coding and 3D structure. The ability of the system to facilitate identifying genetic alterations is illustrated in three projects with publicly accessible data. Mutagenesis in tumor DNA replication leading to complex genetic changes in the EGFR kinase domain is suggested by a novel deletion-insertion combination observed in paired tumor-normal lung cancer resequencing data. Automated analysis of 152 genes resequenced by the SeattleSNPs group was able to identify 91% of the 1251 indel polymorphisms discovered by SeattleSNPs. In addition, our system discovered 518 novel indels in this data set, 451 of which were found to be valid by manual inspection of sequence traces. Our experience demonstrates that CGWB not only greatly improves the productivity and the accuracy of mutation identification, but also, through its data integration and visualization capabilities, facilitates identification of underlying genetic etiology.

A high-resolution multistrain haplotype analysis of laboratory mouse genome reveals three distinctive genetic variation patterns.

Understanding of the structure and the origin of genetic variation patterns in the laboratory inbred mouse provides insight into the utility of the mouse model for studying human complex diseases and strategies for disease gene mapping. In order to address this issue, we have constructed a multistrain, high-resolution haplotype map for the 99-Mb mouse Chromosome 16 using approximately 70,000 single nucleotide polymorphism (SNP) markers derived from whole-genome shotgun sequencing of five laboratory inbred strains. We discovered that large polymorphic blocks (i.e., regions where only two haplotypes, thus one SNP conformation, are found in the five strains), large monomorphic blocks (i.e., regions where the five strains share the same haplotype), and fragmented blocks (i.e., regions of greater complexity not resembling at all the first two categories) span 50%, 18%, and 32% of the chromosome, respectively. The haplotype map has 98% accuracy in predicting mouse genotypes in two other studies. Its predictions are also confirmed by experimental results obtained from resequencing of 40-kb genomic sequences at 21 distinct genomic loci in 13 laboratory inbred strains and 12 wild-derived strains. We demonstrate that historic recombination, intra-subspecies variations and inter-subspecies variations have all contributed to the formation of the three distinctive genetic signatures. The results suggest that the controlled complexity of the laboratory inbred strains may provide a means for uncovering the biological factors that have shaped genetic variation patterns.

Detecting false expression signals in high-density oligonucleotide arrays by an in silico approach.

High-density oligonucleotide arrays have become a popular assay for concurrent measurement of mRNA expression at the genome scale. Much effort has been devoted to the development of statistical analysis tools aimed at reducing experimental noise and normalizing experimental variation in gene expression analysis. However, these investigations do not detect or catalog systematic problems associated with specific oligonucleotide probes. Here, we present an investigation of problematic probes that yield consistent but inaccurate signals across multiple experiments. By evaluating data integrity among gene, probe sequence, and genomic structure we identified a total of 20,696 (10.5%) nonspecific probes that could cross-hybridize to multiple genes and a total of 18,363 (9.3%) probes that miss the target transcript sequences on the Affymetrix GeneChip U95A/Av2 array. The numbers of nonspecific and mistargeted probes on the U133A array are 29,405 (12.1%) and 19,717 (8.0%), respectively. The poor performance of the mistargeted probes was confirmed in two GeneChip experiments, in which these probes showed a 20-30% decrease in detecting present signals compared with normal probes. Comparison of qualitative expression signals obtained from SAGE and EST data with those from GeneChip arrays showed that the consistency of the two platforms is 30% lower in problematic probes than in normal probes. A Web application was developed to apply our results for improving the accuracy of expression analysis.