Search

A limited number of cell lines have fueled the majority of preclinical Prostate cancer (PCa) research. Despite tremendous effort in characterizing their molecular profiles, comprehensive whole genome sequencing with allelic phasing of somatic genome alterations has not been undertaken to date. Here, we utilized whole genome Linked-read sequencing to obtain haplotype information from the seven most commonly used PCa cell lines (PC3, LNCaP, DU145, CWR22Rv1, VCaP, LAPC4, MDA-PCa-2b), four castrate resistant (CR) subclones (LNCaP_Abl, LNCaP_C42b, VCaP-CR, LAPC4-CR), and an immortalized prostate epithelial line RWPE-1. Phasing of mutations allowed derivation of “Gene-level Haplotype” to assess whether a gene harbored heterozygous mutations in one or both alleles, providing a comprehensive catalogue of mono or bi-allelically inactivated genes. Phased structural variant analysis allowed identification of complex rearrangement chains consistent with chromothripsis and chromoplexy, with breakpoints occurred across a single allele, providing further evidence that complex SVs occurred in a concerted event, rather than through accumulation of multiple independent rearrangements. Additionally, comparison of parental and CR subclones revealed previously known and novel genomic alterations associated with the CR clones. This study therefore comprehensively characterized phased genomic alterations in the commonly used PCa cell lines and provided a useful resource for future cancer research.

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

Motivation DNA sequencing archives have grown to enormous scales in recent years, and thousands of human genomes have already been sequenced. The size of these data sets has made searching the raw read data infeasible without high-performance data-query technology. Additionally, it is challenging to search a repository of short-read data using relational logic and to apply that logic across samples from multiple whole-genome sequencing samples. Results We have built a compact, efficiently-indexed database that contains the raw read data for over 250 human genomes, encompassing trillions of bases of DNA, and that allows users to search these data in real-time. The Terabase Search Engine enables retrieval from this database of all the reads for any genomic location in a matter of seconds. Users can search using a range of positions or a specific sequence that is aligned to the genome on the fly. Availability and implementation Public access to the Terabase Search Engine database is available at http://tse.idies.jhu.edu. Supplementary information Supplementary data are available at Bioinformatics online.

Chromatin alterations mediate mutations and gene expression changes in cancer. Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) has been utilized to study genome-wide chromatin structure in human cancer cell lines, yet numerous technical challenges limit comparable analyses in primary tumors. Here we have developed a new whole-genome analytic pipeline to optimize ChIP-Seq protocols on patient-derived xenografts from human papillomavirus–related (HPV+) head and neck squamous cell carcinoma (HNSCC) samples. We further associated chromatin aberrations with gene expression changes from a larger cohort of the tumor and normal samples with RNA-Seq data. We detect differential histone enrichment associated with tumor-specific gene expression variation, sites of HPV integration in the human genome, and HPV-associated histone enrichment sites upstream of cancer driver genes, which play central roles in cancer-associated pathways. These comprehensive analyses enable unprecedented characterization of the complex network of molecular changes resulting from chromatin alterations that drive HPV-related tumorigenesis. Cancer Res; 77(23); 6538–50. ©2017 AACR.

Candida glabratais an opportunistic pathogen in humans, responsible for approximately 20% of disseminated candidiasis. Candida glabrata's ability to adhere to host tissue is mediated by GPI-anchored cell wall proteins (GPI-CWPs); the corresponding genes contain long tandem repeat regions. These repeat regions resulted in assembly errors in the reference genome. Here, we performed a de novo assembly of the C. glabrata type strain CBS138 using long single-molecule real-time reads, with short read sequences (Illumina) for refinement, and constructed telomere-to-telomere assemblies of all 13 chromosomes. Our assembly has excellent agreement overall with the current reference genome, but we made substantial corrections within tandem repeat regions. Specifically, we removed 62 genes of which 45 were scrambled due to misassembly in the reference. We annotated 31 novel ORFs of which 24 ORFs are GPI-CWPs. In addition, we corrected the tandem repeat structure of an additional 21 genes. Our corrections to the genome were substantial, with the length of new genes and tandem repeat corrections amounting to approximately 3.8% of the ORFeome length. As most corrections were within the coding regions of GPI-CWP genes, our genome assembly establishes a high-quality reference set of genes and repeat structures for the functional analysis of these cell surface proteins.

Objectives: To validate next-generation sequencing (NGS) technology for clinical diagnosis and to determine appropriate read depth. Methods: We validated the KRAS , BRAF , and EGFR genes within the Ion AmpliSeq Cancer Hotspot Panel using the Ion Torrent Personal Genome Machine (Life Technologies, Carlsbad, CA). Results: We developed a statistical model to determine the read depth needed for a given percent tumor cellularity and number of functional genomes. Bottlenecking can result from too few input genomes. By using 16 formalin-fixed, paraffin-embedded (FFPE) cancer-free specimens and 118 cancer specimens with known mutation status, we validated the six traditional analytic performance characteristics recommended by the Next-Generation Sequencing: Standardization of Clinical Testing Working Group. Baseline noise is consistent with spontaneous and FFPE-induced C:G→T:A deamination mutations. Conclusions: Redundant bioinformatic pipelines are essential, since a single analysis pipeline gave false-negative and false-positive results. NGS is sufficiently robust for the clinical detection of gene mutations, with attention to potential artifacts.

MotivationGenomics features with similar genomewide distributions are generally hypothesized to be functionally related, for example, co-localization of histones and transcription start sites indicate chromatin regulation of transcription factor activity. Therefore, statistical algorithms to perform spatial, genomewide correlation among genomic features are required.ResultsHere, we propose a method, StereoGene, that rapidly estimates genomewide correlation among pairs of genomic features. These features may represent high throughput data mapped to reference genome or sets of genomic annotations in that reference genome. StereoGene enables correlation of continuous data directly, avoiding the data binarization and subsequent data loss. Correlations are computed among neighboring genomic positions using kernel correlation. Representing the correlation as a function of the genome position, StereoGene outputs the local correlation track as part of the analysis. StereoGene also accounts for confounders such as input DNA by partial correlation. We apply our method to numerous comparisons of ChIP-Seq datasets from the Human Epigenome Atlas and FANTOM CAGE to demonstrate its wide applicability. We observe the changes in the correlation between epigenomic features across developmental trajectories of several tissue types consistent with known biology, and find a novel spatial correlation of CAGE clusters with donor splice sites and with poly(A) sites. These analyses provide examples for the broad applicability of StereoGene for regulatory genomics.AvailabilityThe StereoGene C++ source code, program documentation, Galaxy integration scripts and examples are available from the project homepage http://stereogene.bioinf.fbb.msu.ru/Contactfavorov@sensi.orgSupplementary informationSupplementary data are available online.

Characterizing structural variants in the human genome is of great importance, but a genome wide analysis to detect interspersed repeats has not been done. Thus, the degree to which mobile DNAs contribute to genetic diversity, heritable disease, and oncogenesis remains speculative. We perform transposon insertion profiling by microarray (TIP-chip) to map human L1(Ta) retrotransposons (LINE-1 s) genome-wide. This identified numerous novel human L1(Ta) insertional polymorphisms with highly variant allelic frequencies. We also explored TIP-chip's usefulness to identify candidate alleles associated with different phenotypes in clinical cohorts. Our data suggest that the occurrence of new insertions is twice as high as previously estimated, and that these repeats are under-recognized as sources of human genomic and phenotypic diversity. We have just begun to probe the universe of human L1(Ta) polymorphisms, and as TIP-chip is applied to other insertions such as Alu SINEs, it will expand the catalog of genomic variants even further.

Long interspersed element type 1 (L1) retrotransposons are ubiquitous mammalian mobile elements and potential tools for in vivo mutagenesis; however, native L1 elements are relatively inactive in mice when introduced as transgenes. We have previously described a synthetic L1 element, ORFeus , containing two synonymously recoded ORFs relative to mouse L1. It is significantly more active for retrotransposition in cell culture than all native L1 elements tested. To study its activity in vivo , we developed a transgenic mouse model in which ORFeus expression was controlled by a constitutive heterologous promoter, and we established definitive evidence for ORFeus retrotransposition activity both in germ line and somatic tissues. Germ line retrotransposition frequencies resulting in 0.33 insertions per animal are seen among progeny of ORFeus donor element heterozygotes derived from a single founder, representing a >20-fold increase over native L1 elements. We observe somatic transposition events in 100% of the ORFeus donor-containing animals, and an average of 17 different insertions are easily recovered from each animal; modeling suggests that the number of somatic insertions per animal exceeds this number by perhaps several orders of magnitude. Nearly 200 insertions were precisely mapped, and their distribution in the mouse genome appears random relative to transcription units and guanine-cytosine content. The results suggest that ORFeus may be developed into useful tools for in vivo mutagenesis.

Mobile elements are important components of our genomes, with diverse and significant effects on phenotype. Not only can transposons inactivate genes by direct disruption and shuffle the genome through recombination, they can also alter gene expression subtly or powerfully. Currently active transposons are highly polymorphic in host populations, including, among hundreds of others, L1 and Alu elements in humans and Ty1 elements in yeast. For this reason, we wished to develop a simple genome-wide method for identifying all transposons in any given sample. We have designed a transposon insertion site profiling chip (TIP-chip), a microarray intended for use as a high-throughput technique for mapping transposon insertions. By selectively amplifying transposon flanking regions and hybridizing them to the array, we can locate all transposons present in a sample. We have tested the TIP-chip extensively to map Ty1 retrotransposon insertions in yeast and have achieved excellent results in two laboratory strains as well as in evolved Ty1 high-copy strains. We are able to identify all of the theoretically detectable transposons in the FY2 lab strain, with essentially no false positives. In addition, we mapped many new transposon copies in the high-copy Ty1 strain and determined its Ty1 insertion pattern.

The L1 retrotransposon is the most highly successful autonomous retrotransposon in mammals. This prolific genome parasite may on occasion benefit its host through genome rearrangements or adjustments of host gene expression. In examining possible effects of L1 elements on host gene expression, we investigated whether a full-length L1 element inserted in the antisense orientation into an intron of a cellular gene may actually split the gene's transcript into two smaller transcripts: (1) a transcript containing the upstream exons and terminating in the major antisense polyadenylation site (MAPS) of the L1, and (2) a transcript derived from the L1 antisense promoter (ASP) that includes the downstream exons of the gene. Bioinformatic analysis and experimental follow-up provide evidence for this L1 “gene-breaking” hypothesis. We identified three human genes apparently “broken” by L1 elements, as well as 12 more candidate genes. Most of the inserted L1 elements in our 15 candidate genes predate the human/chimp divergence. If indeed split, the transcripts of these genes may in at least one case encode potentially interacting proteins, and in another case may encode novel proteins. Gene-breaking represents a new mechanism through which L1 elements remodel mammalian genomes.

Single-nucleotide polymorphisms (SNPs) constitute the great majority of variations in the human genome, and as heritable variable landmarks they are useful markers for disease mapping and resolving population structure. Redundant coverage in overlaps of large-insert genomic clones, sequenced as part of the Human Genome Project, comprises a quarter of the genome, and it is representative in terms of base compositional and functional sequence features. We mined these regions to produce 500,000 high-confidence SNP candidates as a uniform resource for describing nucleotide diversity and its regional variation within the genome. Distributions of marker density observed at different overlap length scales under a model of recombination and population size change show that the history of the population represented by the public genome sequence is one of collapse followed by a recent phase of mild size recovery. The inferred times of collapse and recovery are Upper Paleolithic, in agreement with archaeological evidence of the initial modern human colonization of Europe.

This project develops a novel experimental technique to perform ChIP-Seq (chromatin immunoprecipitation with massively parallel DNA sequencing) analysis of chromatin structure in primary tumor tissues from high risk HPV-related head and neck squamous cell carcinomas (HPV+ HNSCC). Recent data suggest that chromatin structure is the central regulator and predictor of cancer-specific expression and mutagenesis landscape of diseased cells. Genome-wide gene expression dysregulation in many tumors, including HPV+ HNSCC, are incompletely described by current knowledge. Methods for study of chromatin structure in primary tumor tissue are needed to better understand the role global epigenetic changes may play in these tumors. However, ChIP-Seq, which is the state-of-the-art method of elucidating chromatin structure, until now, has not been reliably performed on any HNSCC samples. Because chromatin structure is disrupted at room temperature, ChIP-Seq is especially complicated for primary patient tissues, which are primarily obtained as surgical waste after pathology review. Snap freezing of leftover waste surgical tissues and further tissue thawing for the analysis decreases chromatin structure integrity necessary for highly sensitive ChIP-Seq methodology, especially for tumor samples with chromatin structure deformed during carcinogenesis. To improve the chromatin structure integrity in tumor sample we added a xenografting step and minimized the exposure of cancer tissue to room temperature conditions after mouse surgery. We also minimized patient non-cancer tissue preservation at ambient temperature after patient surgery. We successfully performed ChIP-Seq for H3K4me3, H3K9me3, and H3K9ac on frozen uvulopalatopharyngoplasty (UPPP) primary tissues, frozen patient derived xenograft tissues, and freshly-cultured head and neck squamous cell carcinoma cell lines, revealing comparable success rates between tissue type and sample preservation techniques. ChIP-Seq techniques were performed and cross validated using tried and true qRT-PCR methods to demonstrate data reproducibility. The biological relevance of the ChIP-Seq data was confirmed through massive RNA-Seq analysis of 47 HPV+ HNSCC samples and 25 non-cancer controls. Analysis revealed that most H3K9ac and H3K9me3 enrichment is similar in primary tissues, regardless of disease status. Only small portion of them showed differential histone enrichment, which correlated with differential expression of corresponding genes. On the other hand, H3K4me3 showed strong tissue specificity and were found differentially enriched especially in tumor samples. The proposed experimental pipeline demonstrates high reproducibility between biological replicates, diversity of tissue models, and low dependence of ChIP-Seq analysis on tissue preservation techniques. Citation Format: Dylan Z. Kelley, Emily L. Flam, Hildegard A. Wulf, Theresa Guo, Evgeny Izumchenko, Dzov A. Singman, Ludmila V. Danilova, Elena D. Stavrovskaya, Michael Considine, Justin A. Bishop, William H. Westra, Zubair Khan, Wayne M. Koch, David Sidransky, Sarah Wheelan, Joseph A. Califano, Alexander V. Favorov, Elana J. Fertig, Daria A. Gaykalova. The in-parallel whole-genome ChIP-Seq analysis of primary tissues, patient derived xenografts, and cancer cell lines from HPV-relative HNSCC samples [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2424. doi:10.1158/1538-7445.AM2017-2424

Carbone, Lucia et al.-- This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported licence., Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

We have developed a computer program that aligns spliced sequences to genomic sequences, using local alignment algorithms and heuristics to put together a global spliced alignment. Spidey can produce reliable alignments quickly, even when confronted with noise from alternative splicing, polymorphisms, sequencing errors, or evolutionary divergence. We show how Spidey was used to align reference sequences to known genomic sequences and then to the draft human genome, to align mRNAs to gene clusters, and to align mouse mRNAs to human genomic sequence. We compared Spidey to two other spliced alignment programs; Spidey generally performed quite well in a very reasonable amount of time.

Human identity testing is critical to the fields of forensics, paternity, and hematopoietic stem cell transplantation. Most bone marrow (BM) engraftment testing currently uses microsatellites or short tandem repeats that are resolved by capillary electrophoresis. Single-nucleotide polymorphisms (SNPs) are theoretically a better choice among polymorphic DNA; however, ultrasensitive detection of SNPs using next-generation sequencing is currently not possible because of its inherently high error rate. We circumvent this problem by analyzing blocks of closely spaced SNPs, or haplotypes. As proof-of-principle, we chose the HLA-A locus because it is highly polymorphic and is already genotyped to select proper donors for BM transplant recipients. We aligned common HLA-A alleles and identified a region containing 18 closely spaced SNPs, flanked by nonpolymorphic DNA for primer placement. Analysis of cell line mixtures shows that the assay is accurate and precise, and has a lower limit of detection of approximately 0.01%. The BM from a series of hematopoietic stem cell transplantation patients who tested as all donor by short tandem repeat analysis demonstrated 0% to 1.5% patient DNA. Comprehensive analysis of the human genome using the 1000 Genomes database identified many additional loci that could be used for this purpose. This assay may prove useful to identify hematopoietic stem cell transplantation patients destined to relapse, microchimerism associated with solid organ transplantation, forensic applications, and possibly patient identification.

The transposon piggyBac is being used increasingly for genetic studies. Here, we describe modified versions of piggyBac transposase that have potentially wide-ranging applications, such as reversible transgenesis and modified targeting of insertions. piggyBac is distinguished by its ability to excise precisely, restoring the donor site to its pretransposon state. This characteristic makes piggyBac useful for reversible transgenesis, a potentially valuable feature when generating induced pluripotent stem cells without permanent alterations to genomic sequence. To avoid further genome modification following piggyBac excision by reintegration, we generated an excision competent/integration defective (Exc(+)Int(-)) transposase. Our findings also suggest the position of a target DNA-transposase interaction. Another goal of genome engineering is to develop reagents that can guide transgenes to preferred genomic regions. Others have shown that piggyBac transposase can be active when fused to a heterologous DNA-binding domain. An Exc(+)Int(-) transposase, the intrinsic targeting of which is defective, might also be a useful intermediate in generating a transposase whose integration activity could be rescued and redirected by fusion to a site-specific DNA-binding domain. We show that fusion to two designed zinc finger proteins rescued the Int(-) phenotype. Successful guided transgene integration into genomic DNA would have broad applications to gene therapy and molecular genetics. Thus, an Exc(+)Int(-) transposase is a potentially useful reagent for genome engineering and provides insight into the mechanism of transposase-target DNA interaction.

We have created a statistically grounded tool for determining the correlation of genomewide data with other datasets or known biological features, intended to guide biological exploration of high-dimensional datasets, rather than providing immediate answers. The software enables several biologically motivated approaches to these data and here we describe the rationale and implementation for each approach. Our models and statistics are implemented in an R package that efficiently calculates the spatial correlation between two sets of genomic intervals (data and/or annotated features), for use as a metric of functional interaction. The software handles any type of pointwise or interval data and instead of running analyses with predefined metrics, it computes the significance and direction of several types of spatial association; this is intended to suggest potentially relevant relationships between the datasets. Availability and implementation: The package, GenometriCorr, can be freely downloaded at http://genometricorr.sourceforge.net/. Installation guidelines and examples are available from the sourceforge repository. The package is pending submission to Bioconductor.

Background Although nearly half of the human genome is comprised of repetitive sequences, the expression profile of these elements remains largely uncharacterized. Recently developed high throughput sequencing technologies provide us with a powerful new set of tools to study repeat elements. Hence, we performed whole transcriptome sequencing to investigate the expression of repetitive elements in human frontal cortex using postmortem tissue obtained from the Stanley Medical Research Institute. Results We found a significant amount of reads from the human frontal cortex originate from repeat elements. We also noticed that Alu elements were expressed at levels higher than expected by random or background transcription. In contrast, L1 elements were expressed at lower than expected amounts. Conclusions Repetitive elements are expressed abundantly in the human brain. This expression pattern appears to be element specific and can not be explained by random or background transcription. These results demonstrate that our knowledge about repetitive elements is far from complete. Further characterization is required to determine the mechanism, the control, and the effects of repeat element expression.

Human cancer cell lines and xenografts are valuable samples for whole-genome sequencing of human cancer. Tumors can be maintained by serial xenografting in athymic (nude) or severe combined immunodeficient (SCID) mice. In the current study, we developed a molecular assay to quantify the relative contributions of human and mouse in mixed DNA samples. The assay was designed based on deletion/insertion variation between human and mouse genomes. The percentage of mouse DNA was calculated according to the relative peak heights of PCR products analyzed by capillary electrophoresis. Three markers from chromosomes 9 and 10 accurately predicted the mouse genome ratio and were combined into a multiplex PCR reaction. We used the assay to quantify the relative DNA amounts of 93 mouse xenografts used for a recently reported integrated genomic analysis of human pancreatic cancer. Of the 93 xenografts, the mean percentage of contaminating mouse DNA was 47%, ranging from 17% to 73%, with 43% of samples having >50% mouse DNA. We then comprehensively compared the human and mouse genomes to identify 370 additional candidate gene loci demonstrating human-mouse length variation. With increasing whole-genome sequencing of human cancers, this assay should be useful to monitor strategies to enrich human cancer cells from mixed human-mouse cell xenografts. Finally, we discuss how contaminating mouse DNA affects next-generation DNA sequencing.

Transposable elements are ubiquitous among sequenced genomes. The host genomes roughly subdivide into two types: 1) streamlined, that is, small,with little space between genes and lacking large introns, or 2) bulky, with lots of space between genes and many large introns. Most microorganisms, along with selected vertebrates like the pufferfish, fall into the first class, whereas mammals and most plants fall into the second class. As can be seen from Fig. 1, transposable element abundance mirrors the genome type of the host, with mobile elements comprising half or more of many of these bulky genomes! Mobile elements are of two basic types: DNA transposons, which predominantly mobilize via a cut and paste mechanism, and retrotransposons, which move by a copy and paste mechanism involving reverse transcription of an RNA intermediate (Fig. 1 right panel; Curcio and Derbyshire 2003). Retrotransposons are found in virtually all eukaryotes, from yeast (Kim et al. 1998) to human (Lander et al. 2001). Remarkably, in a yeast cell, the number of retrotransposon copies can be changed rather dramatically without a major impact on the phenotype of the host. The change in copy number can be seen using a new tiling array technique by which it is possible to comprehensively map the unique genomic regions adjacent to all transposable element copies probed (Fig. 2; Wheelan et al. 2006). The ability of yeast strains to tolerate very high copy numbers of transposons is due in part to the fact that, in yeast, most insertions are targeted to non-essential genomic regions, even though most of the genome is protein-coding (Chalker and Sandmeyer 1990; Devine and Boeke 1996; Ji et al. 1993; Zou et al. 1996). This property and many others suggest that retrotransposons are highly coevolved with their hosts.

Modern molecular biology has brought many new tools to the geneticist as well as an exponentially expanding database of genomes and new genes for study. Of particular use in the analysis of these genes is the synthetic gene, a nucleotide sequence designed to the specifications of the investigator. Typically, synthetic genes encode the same product as the gene of interest, but the synthetic nucleotide sequence for that protein may contain modifications affecting expression or base composition. Other desirable changes typically involve the revision of restriction sites. Designing synthetic genes by hand is a time-consuming and error-prone process that may involve several computer programs. We have developed a tools environment that combines many modules to provide a platform for rapid synthetic gene design for multikilobase sequences. We have used GeneDesign to successfully design a synthetic Ty1 element and a large variety of other synthetic sequences. GeneDesign has been implemented as a publicly accessible Web-based resource and can be found at http://slam.bs.jhmi.edu/gd.

A major barrier to conceptual advances in understanding the mechanisms and regulation of imprinting of a genomic region is our relatively poor understanding of the overall organization of genes and of the potentially important cis-acting regulatory sequences that lie in the nonexonic segments that make up 97% of the genome. Interspecies sequence comparison offers an effective approach to identify sequence from conserved functional elements. In this article we describe the successful use of this approach in comparing a approximately 1-Mb imprinted genomic domain on mouse chromosome 7 to its orthologous region on human 11p15.5. Within the region, we identified 112 exons of known genes as well as a novel gene identified uniquely in the mouse region, termed Msuit, that was found to be imprinted. In addition to these coding elements, we identified 33 CpG islands and 49 orthologous nonexonic, nonisland sequences that met our criteria as being conserved, and making up 4.1% of the total sequence. These conserved noncoding sequence elements were generally clustered near imprinted genes and the majority were between Igf2 and H19 or within Kvlqt1. Finally, the location of CpG islands provided evidence that suggested a two-island rule for imprinted genes. This study provides the first global view of the architecture of an entire imprinted domain and provides candidate sequence elements for subsequent functional analyses.

Human cancer cell lines and xenografts are valuable samples for whole-genome sequencing of human cancer. Tumors can be maintained by serial xenografting in athymic (nude) or severe combined immunodeficient (SCID) mice. Critical to next generation based sequencing is the relative proportions of human and mouse DNA. In the current study, we developed a molecular assay to quantify the relative contributions of human and mouse in mixed DNA samples. The assay was designed based on deletion/insertion variation between human and mouse genomes. The percentage of mouse DNA was calculated according to the relative peak heights of PCR products analyzed by capillary electrophoresis. Three markers from chromosomes 9 and 10 accurately predicted the mouse genome ratio and were combined into a multiplex PCR reaction. We used the assay to quantify the relative DNA amounts of 93 mouse xenografts used for a recently reported integrated genomic analysis of human pancreatic cancer. Of the 93 xenografts, the mean % of contaminating mouse DNA was 47%, ranging from 17% to 73%, with 43% of samples having greater than 50% mouse DNA. We then comprehensively compared the human and mouse genomes to identify 370 additional candidate gene loci demonstrating human-mouse length variation. With increasing use of whole genome sequencing of human cancers, this assay should be useful to monitor strategies to enrich human cancer cells from mixed human-mouse cell xenografts. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 2210.

Summary Female eutherian mammals use X chromosome inactivation (XCI) to epigenetically regulate gene expression from ∼4% of the genome. To quantitatively map the topography of XCI for defined cell types at single cell resolution, we have generated female mice that carry X-linked, Cre-activated, and nuclear-localized fluorescent reporters—GFP on one X chromosome and tdTomato on the other. Using these reporters in combination with different Cre drivers, we have defined the topographies of XCI mosaicism for multiple CNS cell types and of retinal vascular dysfunction in a model of Norrie disease. Depending on cell type, fluctuations in the XCI mosaic are observed over a wide range of spatial scales, from neighboring cells to left versus right sides of the body. These data imply a major role for XCI in generating female-specific, genetically directed, stochastic diversity in eutherian mammals on spatial scales that would be predicted to affect CNS function within and between individuals.