Differential spatial and structural organization of the X chromosome underlies dosage compensation in C. elegans.

The adjustment of X-linked gene expression to the X chromosome copy number (dosage compensation [DC]) has been widely studied as a model of chromosome-wide gene regulation. In Caenorhabditis elegans, DC is achieved by twofold down-regulation of gene expression from both Xs in hermaphrodites. We show that in males, the single X chromosome interacts with nuclear pore proteins, while in hermaphrodites, the DC complex (DCC) impairs this interaction and alters X localization. Our results put forward a structural model of DC in which X-specific sequences locate the X chromosome in transcriptionally active domains in males, while the DCC prevents this in hermaphrodites.

Melanotic tumors of the nervous system are characterized by distinct mutational, chromosomal and epigenomic profiles.

Melanotic tumors of the nervous system show overlapping histological characteristics but differ substantially in their biological behavior. In order to achieve a better delineation of such tumors, we performed an in-depth molecular characterization. Eighteen melanocytomas, 12 melanomas, and 14 melanotic and 14 conventional schwannomas (control group) were investigated for methylome patterns (450k array), gene mutations associated with melanotic tumors and copy number variants (CNVs). The methylome fingerprints assigned tumors to entity-specific groups. Methylation groups also showed a substantial overlap with histology-based diagnosis suggesting that they represent true biological entities. On the molecular level, melanotic schwannomas were characterized by a complex karyotype with recurrent monosomy of chromosome 22q and variable whole chromosomal gains and recurrent losses commonly involving chromosomes 1, 17p and 21. Melanocytomas carried GNAQ/11 mutations and presented with CNV involving chromosomes 3 and 6. Melanomas were frequently mutated in the TERT promoter, harbored additional oncogene mutations and showed recurrent chromosomal losses involving chromosomes 9, 10 and 6q, as well as gains of 22q. Together, melanotic nervous system tumors have several distinct mutational and chromosomal alterations and can reliably be distinguished by methylome profiling.

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

Segregation of viral plasmids depends on tethering to chromosomes and is regulated by phosphorylation

Eukaryotic viruses can maintain latency in dividing cells as extrachromosomal nuclear plasmids. Segregation and nuclear retention of DNA is, therefore, a key issue in retaining copy number. The E2 enhancer protein of the papillomaviruses is required for viral DNA replication and transcription. Viral mutants that prevent phosphorylation of the bovine papillomavirus type 1 (BPV) E2 protein are transformation-defective, despite normal viral gene expression and replication function. Cell colonies harboring such mutants show sectoring of viral DNA and are unable to maintain the episome. We find that transforming viral DNA attaches to mitotic chromosomes, in contrast to the mutant genome encoding the E2 phosphorylation mutant. Second-site suppressor mutations were uncovered in both E1 and E2 genes that allow for transformation, maintenance, and chromosomal attachment. E2 protein was also found to colocalize to mitotic chromosomes, whereas the mutant did not, suggesting a direct role for E2 in viral attachment to chromosomes. Such viral hitch-hiking onto cellular chromosomes is likely to provide a general mechanism for maintaining nuclear plasmids.

Dominant transformation by mutated human ras genes in vitro requires more than 100 times higher expression than is observed in cancers

The gene-mutation-cancer hypothesis holds that mutated cellular protooncogenes, such as point-mutated proto-ras, “play a dominant part in cancer,” because they are sufficient to transform transfected mouse cell lines in vitro [Alberts, B., Bray, D., Lewis, J., Raff, M., Roberts, K. & Watson, J. D. (1994) Molecular Biology of the Cell (Garland, New York)]. However, in cells transformed in vitro mutated human ras genes are expressed more than 100-fold than in the cancers from which they are isolated. In view of the discrepancy between the very low levels of ras transcription in cancers and the very high levels in cells transformed in vitro, we have investigated the minimal level of human ras expression for transformation in vitro. Using point-mutated human ras genes recombined with different promoters from either human metallothionein-IIA or human fibronectin or from retroviruses we found dominant in vitro transformation of the mouse C3H cell line only with ras genes linked to viral promoters. These ras genes were expressed more than 120-fold higher than are native ras genes of C3H cells. The copy number of transfected ras genes ranged from 2–6 in our system. In addition, nondominant transformation was observed in a small percentage (2–7%) of C3H cells transfected with ras genes that are expressed less than 20 times higher than native C3H ras genes. Because over 90% of cells expressing ras at this moderately enhanced level were untransformed, transformation must follow either a nondominant ras mechanism or a non-ras mechanism. We conclude that the mutated, but normally expressed, ras genes found in human and animal cancers are not likely to “play a dominant part in cancer.” The conclusion that mutated ras genes are not sufficient or dominant for cancer is directly supported by recent discoveries of mutated ras in normal animals, and in benign human tissue, “which has little potential to progress” [Jen, J., Powell, S. M., Papadopoulos, N., Smith, K. J., Hamilton, S. R., Vogelstein, B. & Kinzler, K. W. (1994) Cancer Res. 54, 5523–5526]. Even the view that mutated ras is necessary for cancer is hard to reconcile with (i) otherwise indistinguishable cancers with and without ras mutations, (ii) metastases of the same human cancers with and without ras mutations, (iii) retroviral ras genes that are oncogenic without point mutations, and (iv) human tumor cells having spontaneously lost ras mutation but not tumorigencity.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

Molecular characterization of a mutable pigmentation phenotype and isolation of the first active transposable element from Sorghum bicolor

Accumulation of red phlobaphene pigments in sorghum grain pericarp is under the control of the Y gene. A mutable allele of Y, designated as y-cs (y-candystripe), produces a variegated pericarp phenotype. Using probes from the maize p1 gene that cross-hybridize with the sorghum Y gene, we isolated the y-cs allele containing a large insertion element. Our results show that the Y gene is a member of the MYB-transcription factor family. The insertion element, named Candystripe1 (Cs1), is present in the second intron of the Y gene and shares features of the CACTA superfamily of transposons. Cs1 is 23,018 bp in size and is bordered by 20-bp terminal inverted repeat sequences. It generated a 3-bp target site duplication upon insertion within the Y gene and excised from y-cs, leaving a 2-bp footprint in two cases analyzed. Reinsertion of the excised copy of Cs1 was identified by Southern hybridization in the genome of each of seven red pericarp revertant lines tested. Cs1 is the first active transposable element isolated from sorghum. Our analysis suggests that Cs1-homologous sequences are present in low copy number in sorghum and other grasses, including sudangrass, maize, rice, teosinte, and sugarcane. The low copy number and high transposition frequency of Cs1 imply that this transposon could prove to be an efficient gene isolation tool in sorghum.

Mechanisms of Human Mitochondrial DNA Maintenance: The Determining Role of Primary Sequence and Length over Function

Although the regulation of mitochondrial DNA (mtDNA) copy number is performed by nuclear-coded factors, very little is known about the mechanisms controlling this process. We attempted to introduce nonhuman ape mtDNA into human cells harboring either no mtDNA or mutated mtDNAs (partial deletion and tRNA gene point mutation). Unexpectedly, only cells containing no mtDNA could be repopulated with nonhuman ape mtDNA. Cells containing a defective human mtDNA did not incorporate or maintain ape mtDNA and therefore died under selection for oxidative phosphorylation function. On the other hand, foreign human mtDNA was readily incorporated and maintained in these cells. The suicidal preference for self-mtDNA showed that functional parameters associated with oxidative phosphorylation are less relevant to mtDNA maintenance and copy number control than recognition of mtDNA self-determinants. Non–self-mtDNA could not be maintained into cells with mtDNA even if no selection for oxidative phosphorylation was applied. The repopulation kinetics of several mtDNA forms after severe depletion by ethidium bromide treatment showed that replication and maintenance of mtDNA in human cells are highly dependent on molecular features, because partially deleted mtDNA molecules repopulated cells significantly faster than full-length mtDNA. Taken together, our results suggest that mtDNA copy number may be controlled by competition for limiting levels of trans-acting factors that recognize primarily mtDNA molecular features. In agreement with this hypothesis, marked variations in mtDNA levels did not affect the transcription of nuclear-coded factors involved in mtDNA replication.

The MITE family Heartbreaker (Hbr): Molecular markers in maize

Transposable elements are ubiquitous in plant genomes, where they frequently comprise the majority of genomic DNA. The maize genome, which is believed to be structurally representative of large plant genomes, contains single genes or small gene islands interspersed with much longer blocks of retrotransposons. Given this organization, it would be desirable to identify molecular markers preferentially located in genic regions. In this report, the features of a newly described family of miniature inverted repeat transposable elements (MITEs) (called Heartbreaker), including high copy number and polymorphism, stability, and preference for genic regions, have been exploited in the development of a class of molecular markers for maize. To this end, a modification of the AFLP procedure called transposon display was used to generate and display hundreds of genomic fragments anchored in Hbr elements. An average of 52 markers were amplified for each primer combination tested. In all, 213 polymorphic fragments were reliably scored and mapped in 100 recombinant inbred lines derived from a cross between the maize inbreds B73 × Mo17. In this mapping population, Hbr markers are distributed evenly across the 10 maize chromosomes. This procedure should be of general use in the development of markers for other MITE families in maize and in other plant and animal species where MITEs have been identified.

Mendel-GFDb and Mendel-ESTS: databases of plant gene families and ESTs annotated with gene family numbers and gene family names

There is no control over the information provided with sequences when they are deposited in the sequence databases. Consequently mistakes can seed the incorrect annotation of other sequences. Grouping genes into families and applying controlled annotation overcomes the problems of incorrect annotation associated with individual sequences. Two databases (http://www.mendel.ac.uk) were created to apply controlled annotation to plant genes and plant ESTs: Mendel-GFDb is a database of plant protein (gene) families based on gapped-BLAST analysis of all sequences in the SWISS-PROT family of databases. Sequences are aligned (ClustalW) and identical and similar residues shaded. The families are visually curated to ensure that one or more criteria, for example overall relatedness and/or domain similarity relate all sequences within a family. Sequence families are assigned a ‘Gene Family Number’ and a unified description is developed which best describes the family and its members. If authority exists the gene family is assigned a ‘Gene Family Name’. This information is placed in Mendel-GFDb. Mendel-ESTS is primarily a database of plant ESTs, which have been compared to Mendel-GFDb, completely sequenced genomes and domain databases. This approach associated ESTs with individual sequences and the controlled annotation of gene families and protein domains; the information being placed in Mendel-ESTS. The controlled annotation applied to genes and ESTs provides a basis from which a plant transcription database can be developed.

Antisense-mediated decrease in DNA ligase III expression results in reduced mitochondrial DNA integrity

The human DNA ligase III gene encodes both nuclear and mitochondrial proteins. Abundant evidence supports the conclusion that the nuclear DNA ligase III protein plays an essential role in both base excision repair and homologous recombination. However, the role of DNA ligase III protein in mitochondrial genome dynamics has been obscure. Human tumor-derived HT1080 cells were transfected with an antisense DNA ligase III expression vector and clones with diminished levels of DNA ligase III activity identified. Mitochondrial protein extracts prepared from these clones had decreased levels of DNA ligase III relative to extracts from cells transfected with a control vector. Analysis of these clones revealed that the DNA ligase III antisense mRNA-expressing cells had reduced mtDNA content compared to control cells. In addition, the residual mtDNA present in these cells had numerous single-strand nicks that were not detected in mtDNA from control cells. Cells expressing antisense ligase III also had diminished capacity to restore their mtDNA to pre-irradiation levels following exposure to γ-irradiation. An antisense-mediated reduction in cellular DNA ligase IV had no effect on the copy number or integrity of mtDNA. This observaion, coupled with other evidence, suggests that DNA ligase IV is not present in the mitochondria and does not play a role in maintaining mtDNA integrity. We conclude that DNA ligase III is essential for the proper maintenance of mtDNA in cultured mammalian somatic cells.

Directed evolution of polymerase function by compartmentalized self-replication

We describe compartmentalized self-replication (CSR), a strategy for the directed evolution of enzymes, especially polymerases. CSR is based on a simple feedback loop consisting of a polymerase that replicates only its own encoding gene. Compartmentalization serves to isolate individual self-replication reactions from each other. In such a system, adaptive gains directly (and proportionally) translate into genetic amplification of the encoding gene. CSR has applications in the evolution of polymerases with novel and useful properties. By using three cycles of CSR, we obtained variants of Taq DNA polymerase with 11-fold higher thermostability than the wild-type enzyme or with a >130-fold increased resistance to the potent inhibitor heparin. Insertion of an extra stage into the CSR cycle before the polymerase reaction allows its application to enzymes other than polymerases. We show that nucleoside diphosphate kinase and Taq polymerase can form such a cooperative CSR cycle based on reciprocal catalysis, whereby nucleoside diphosphate kinase produces the substrates required for the replication of its own gene. We also find that in CSR the polymerase genes themselves evolve toward more efficient replication. Thus, polymerase genes and their encoded polypeptides cooperate to maximize postselection copy number. CSR should prove useful for the directed evolution of enzymes, particularly DNA or RNA polymerases, as well as for the design and study of in vitro self-replicating systems mimicking prebiotic evolution and viral replication.

Mechanism of DNA polymerase II-mediated frameshift mutagenesis

Escherichia coli possesses three SOS-inducible DNA polymerases (Pol II, IV, and V) that were recently found to participate in translesion synthesis and mutagenesis. Involvement of these polymerases appears to depend on the nature of the lesion and its local sequence context, as illustrated by the bypass of a single N-2-acetylaminofluorene adduct within the NarI mutation hot spot. Indeed, error-free bypass requires Pol V (umuDC), whereas mutagenic (−2 frameshift) bypass depends on Pol II (polB). In this paper, we show that purified DNA Pol II is able in vitro to generate the −2 frameshift bypass product observed in vivo at the NarI sites. Although the ΔpolB strain is completely defective in this mutation pathway, introduction of the polB gene on a low copy number plasmid restores the −2 frameshift pathway. In fact, modification of the relative copy number of polB versus umuDC genes results in a corresponding modification in the use of the frameshift versus error-free translesion pathways, suggesting a direct competition between Pol II and V for the bypass of the same lesion. Whether such a polymerase competition model for translesion synthesis will prove to be generally applicable remains to be confirmed.

Single-copy transgenic mice with chosen-site integration.

We describe a general way of introducing transgenes into the mouse germ line for comparing different sequences without the complications of variation in copy number and insertion site. The method uses homologous recombination in embryonic stem (ES) cells to generate mice having a single copy of a transgene integrated into a chosen location in the genome. To test the method, a single copy murine bcl-2 cDNA driven by either a chicken beta-actin promoter or a human beta-actin promoter has been inserted immediately 5' to the X-linked hypoxanthine phosphoribosyltransferase locus by a directly selectable homologous recombination event. The level of expression of the targeted bcl-2 transgene in ES cells is identical in independently isolated homologous recombinants having the same promoter yet varies between the different promoters. In contrast, the expression of bcl-2 transgenes having the same (chicken beta-actin) promoter varies drastically when they are independently integrated at random insertion sites. Both promoters direct broad expression of the single-copy transgene in mice derived from the respective targeted ES cells. In vitro and in vivo, the human beta-actin promoter consistently directed a higher level of transgene expression than the chicken beta-actin promoter.