DNA supercoiling is an inherent consequence of twisting DNA and is critical for regulating gene expression and DNA replication. However, DNA supercoiling at a genomic scale in human cells is uncharacterized. To map supercoiling we used biotinylated-trimethylpsoralen as a DNA structure probe to show the genome is organized into supercoiling domains. Domains are formed and remodeled by RNA polymerase and topoisomerase activities and are flanked by GC-AT boundaries and CTCF binding sites. Under-wound domains are transcriptionally active, enriched in topoisomerase I, “open” chromatin fibers and DNaseI sites, but are depleted of topoisomerase II. Furthermore DNA supercoiling impacts on additional levels of chromatin compaction as under-wound domains are cytologically decondensed, topologically constrained, and decompacted by transcription of short RNAs. We suggest that supercoiling domains create a topological environment that facilitates gene activation providing an evolutionary purpose for clustering genes along chromosomes.

Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges, including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures. Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice, very useful profiles can be obtained with~0.13 genome coverage. We improve on previous methods by first implementing a combined correction for sequence mappability and GC content, and second, by applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions was previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1000 samples, most of which were obtained from the fixed tissue archives of more than 25 institutions. We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches and better copy number data than high-resolution microarrays at a substantially lower cost.

Despite developments in targeted gene sequencing and whole-genome analysis techniques, the robust detection of all genetic variation, including structural variants, in and around genes of interest and in an allele-specific manner remains a challenge. Here we present targeted locus amplification (TLA), a strategy to selectively amplify and sequence entire genes on the basis of the crosslinking of physically proximal sequences. We show that, unlike other targeted re-sequencing methods, TLA works without detailed prior locus information, as one or a few primer pairs are sufficient for sequencing tens to hundreds of kilobases of surrounding DNA. This enables robust detection of single nucleotide variants, structural variants and gene fusions in clinically relevant genes, including BRCA1 and BRCA2, and enables haplotyping. We show that TLA can also be used to uncover insertion sites and sequences of integrated transgenes and viruses. TLA therefore promises to be a useful method in genetic research and diagnostics when comprehensive or allele-specific genetic information is needed.

Screening of therapeutics relies on representative cancer models. The representation of human glioblastoma by in vitro cell culture models is questionable. We obtained genomic profiles by array comparative genomic hybridization of both short-and long-term primary cell and spheroid cultures, derived from seven glioblastomas and one anaplastic oligodendroglioma. Chromosomal copy numbers were compared between cell cultures and spheroids and related to the parental gliomas using unsupervised hierarchical clustering and correlation coefficient. In seven out of eight short-term cell cultures, the genomic profiles clustered further apart from their parental tumors than spheroid cultures. In four out of eight samples, the genetic changes in cell culture were substantial. The average correlation coefficient between parental tumors and spheroid profiles was 0.89 (range: 0.79-0.97), whereas that between parental tumors and cell cultures was 0.62 (range: 0.10-0.96). In two out of three long-term cell cultures progressive genetic changes had developed, whereas the spheroid cultures were genetically stable. It is concluded that genomic profiles of primary cell cultures from glioblastoma are frequently deviant from parental tumor profiles, whereas spheroids are genetically more representative of the glioblastoma. This implies that glioma cell culture data have to be handled with the highest caution.

Gemcitabine is a deoxycytidine (dCyd) analogue with activity against several solid cancers. Gemcitabine is activated by dCyd kinase (dCK) and interferes, as its triphosphate dFdCTP, with tumor growth through incorporation into DNA. Alternatively, the metabolite gemcitabine diphosphate (dFdCDP) can interfere with DNA synthesis and thus tumor growth through inhibition of ribonucleotide reductase. Gemcitabine can be inactivated by the enzyme dCyd deaminase (dCDA). In most in vitro models, resistance to gemcitabine was associated with a decreased dCK activity. In all these models, resistance was established using continuous exposure to gemcitabine with increasing concentrations; however, these in vitro models have limited clinical relevance. To develop in vivo resistance to gemcitabine, we treated mice bearing a moderately sensitive tumor Colon 26-A (T/C = 0.25) with a clinically relevant schedule (120 mg/kg every 3 days). By repeated transplant of the most resistant tumor and continuation of gemcitabine treatment for >1 year, the completely resistant tumor Colon 26-G (T/C = 0.96) was created. Initial studies focused on resistance mechanisms known from in vitro studies. In Colon 26-G, dCK activity was 1.7-fold decreased; dCDA and DNA polymerase were not changed; and Colon 26-G accumulated 1.5-fold less dFdCTP, 6 hours after a gemcitabine injection, than the parental tumor. Based on in vitro studies, these relative minor changes were considered insufficient to explain the completely resistant phenotype. Therefore, an expression microarray was done with Colon 26-A versus Colon 26-G. Using independently grown nonresistant and resistant tumors, a striking increase in expression of the RRM1 subunit gene was found in Colon 26-G. The expression of RRM1 mRNA was 25-fold increased in the resistant tumor, as measured by real-time PCR, which was confirmed by Western blotting. In contrast, RRM2 mRNA was 2-fold decreased. However, ribonucleotide reductase enzyme activity was only moderately increased in Colon 26-G. In conclusion, this is the first model with in vivo induced resistance to gemcitabine. In contrast to most in vitro studies, dCK activity was not the most important determinant of gemcitabine resistance. Expression microarray identified RRM1 as the gene with the highest increase in expression in the Colon 26-G, which might clarify its complete gemcitabine-resistant phenotype. This study is the first in vivo evidence for a key role for RRM1 in acquired gemcitabine resistance. (Cancer Res 2005; 65(20): 9510-6)

Myocardial right ventricular (RV) hypertrophy due to pulmonary hypertension is aimed at normalizing ventricular wall stress. Depending on the degree of pressure overload, RV hypertrophy may progress to a state of impaired contractile function and heart failure, but this cannot be discerned during the early stages of ventricular remodeling. We tested whether critical differences in gene expression profiles exist between ventricles before the ultimate development of either a compensated or decompensated hypertrophic phenotype. Both phenotypes were selectively induced in Wistar rats by a single subcutaneous injection of either a low or a high dose of the pyrrolizidine alkaloid monocrotaline (MCT). Spotted oligonucleotide microarrays were used to investigate pressure-dependent cardiac gene expression profiles at 2 wk after the MCT injections, between control rats and rats that would ultimately develop either compensated or decompensated hypertrophy. Clustering of significantly regulated genes revealed specific expression profiles for each group, although the degree of hypertrophy was still similar in both. The ventricles destined to progress to failure showed activation of pro-apoptotic pathways, particularly related to mitochondria, whereas the group developing compensated hypertrophy showed blocked pro-death effector signaling via p38-MAPK, through upregulation of MAPK phosphatase-1. In summary, we show that, already at an early time point, pivotal differences in gene expression exist between ventricles that will ultimately develop either a compensated or a decompensated phenotype, depending on the degree of pressure overload. These data reveal genes that may provide markers for the early prediction of clinical outcome as well as potential targets for early intervention.

Five individual colonies of Haemophilus influenzae were isolated from each of one to three cultures of sputum collected from 18 patients with chronic obstructive pulmonary disease (COPD). The isolates were studied to investigate whether the major outer membrane proteins (MOMPs) changed during persistence. The relationship between isolates was analyzed by fingerprinting their chromosomal DNA. The fingerprints of eight strains (isolated from eight patients) with various MOMP compositions were different, whereas fingerprints of isolates with identical MOMP compositions were indistinguishable. In 12 patients, two or more strains with different MOMP compositions were found; one strain was isolated from the sputum samples of each of the six remaining patients. In seven of the 12 patients, strains with different MOMPs but with indistinguishable fingerprints were found. The differences were found in proteins b,c (five patients) and d (five patients). In patients with COPD, the MOMPs of H. influenzae are subject to changes that may enable this bacterium to escape immunological defense mechanisms.

MicroRNAs (miRs) have been recognized as promising biomarkers. It is unknown to what extent tumor-derived miRs are differentially expressed between primary colorectal cancers (pCRCs) and metastatic lesions, and to what extent the expression profiles of tumor tissue differ from the surrounding normal tissue. Next-generation sequencing (NGS) of 220 fresh-frozen samples, including paired primary and metastatic tumor tissue and non-tumorous tissue from 38 patients, revealed expression of 2245 known unique mature miRs and 515 novel candidate miRs. Unsupervised clustering of miR expression profiles of pCRC tissue with paired metastases did not separate the two entities, whereas unsupervised clustering of miR expression profiles of pCRC with normal colorectal mucosa demonstrated complete separation of the tumor samples from their paired normal mucosa. Two hundred and twenty-two miRs differentiated both pCRC and metastases from normal tissue samples (false discovery rate (FDR) <0.05). The highest expressed tumor-specific miRs were miR-21 and miR-92a, both previously described to be involved in CRC with potential as circulating biomarker for early detection. Only eight miRs, 0.5% of the analysed miR transcriptome, were differentially expressed between pCRC and the corresponding metastases (FDR <0.1), consisting of five known miRs (miR-320b, miR-320d, miR-3117, miR-1246 and miR-663b) and three novel candidate miRs (chr 1-2552-5p, chr 8-20656-5p and chr 10-25333-3p). These results indicate that previously unrecognized candidate miRs expressed in advanced CRC were identified using NGS. In addition, miR expression profiles of pCRC and metastatic lesions are highly comparable and may be of similar predictive value for prognosis or response to treatment in patients with advanced CRC.