Single-Cell Transcriptomics of Bone Marrow Stromal Cells in Diversity Outbred Mice: A Model for Population-Level scRNA-Seq Studies

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics data sets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five diversity outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell types. SCENIC was used to reconstruct gene regulatory networks (GRNs), and we observed that cell types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically informative model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells in large populations. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease

A major fraction of loci identified by genome-wide association studies (GWASs) lead to alterations in alternative splicing, but interpretation of how such alterations impact proteins is hindered by the technical limitations of short-read RNA-seq, which cannot directly link splicing events to full-length transcript or protein isoforms. Long-read RNA-seq represents a powerful tool to define and quantify transcript isoforms, and recently, infer protein isoform existence. Here we present a novel approach that integrates information from GWAS, splicing QTL (sQTL), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes which colocalized with BMD associations (H 4 PP ≥ 0.75). We generated deep coverage PacBio long-read RNA-seq data (N=∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were novel. By casting the colocalized sQTLs directly onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Using these data, we created one of the first proteome-scale resources defining full-length isoforms impacted by colocalized sQTLs. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense mediated decay (NMD) and 190 that potentially resulted in the expression of new protein isoforms. Finally, we identified colocalizing sQTLs in TPM2 for splice junctions between two mutually exclusive exons, and two different transcript termination sites, making it impossible to interpret without long-read RNA-seq data. siRNA mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization. We expect our approach to be widely generalizable across diverse clinical traits and accelerate system-scale analyses of protein isoform activities modulated by GWAS loci.

Intramembranous bone regeneration in diversity outbred mice is heritable

There are over one million cases of failed bone repair in the U.S. annually, resulting in substantial patient morbidity and societal costs. Multiple candidate genes affecting bone traits such as bone mineral density have been identified in human subjects and animal models using genome-wide association studies (GWAS). This approach for understanding the genetic factors affecting bone repair is impractical in human subjects but could be performed in a model organism if there is sufficient variability and heritability in the bone regeneration response. Diversity Outbred (DO) mice, which have significant genetic diversity and have been used to examine multiple intact bone traits, would be an excellent possibility. Thus, we sought to evaluate the phenotypic distribution of bone regeneration, sex effects and heritability of intramembranous bone regeneration on day 7 following femoral marrow ablation in 47 12-week old DO mice (23 males, 24 females). Compared to a previous study using 4 inbred mouse strains, we found similar levels of variability in the amount of regenerated bone (coefficient of variation of 86 % v. 88 %) with approximately the same degree of heritability (0.42 v. 0.49). There was a trend toward more bone regeneration in males than females. The amount of regenerated bone was either weakly or not correlated with bone mass at intact sites, suggesting that the genetic factors responsible for bone regeneration and intact bone phenotypes are at least partially independent. In conclusion, we demonstrate that DO mice exhibit variation and heritability of intramembranous bone regeneration that will be suitable for future GWAS.

Fentanyl-induced acute and conditioned behaviors in two inbred mouse lines: Potential role for Glyoxalase

An increase in opioid-overdose deaths was evident before the COVID-19 pandemic, and has escalated since its onset. Fentanyl, a highly potent synthetic opioid, is the primary driver of these recent trends. The current study used two inbred mouse strains, C57BL/6 J and A/J, to investigate the genetics of behavioral responses to fentanyl. Mice were tested for conditioned place preference and fentanyl-induced locomotor activity. C57BL/6J mice formed a conditioned place preference to fentanyl injections and fentanyl increased their activity. Neither effect was noted in A/J mice. We conducted RNA-sequencing on the nucleus accumbens of mice used for fentanyl-induced locomotor activity. Surprisingly, we noted few differentially expressed genes using treatment as the main factor. However many genes differed between strains. We validated differences in two genes: suppressor APC domain containing 1 (Sapcd1) and Glyoxalase 1 (Glo1), with quantitative PCR on RNA from the nucleus accumbens and prefrontal cortex (). In both regions A/J mice had significantly higher expression of both genes than did C57BL/6 J. In prefrontal cortex, fentanyl treatment decreased Glo1 mRNA. Glyoxalase 1 catalyzes the detoxification of reactive alpha-oxoaldehydes such as glyoxal and methylglyoxal, is associated with anxiety and activity levels, and its inhibition reduces alcohol intake. We suggest that future studies assess the ability of Glo1 and related metabolites to modify opioid intake.

Genomic variants within chromosome 14q32.32 regulate bone mass through MARK3 signaling in osteoblasts

Bone mineral density (BMD) is a highly heritable predictor of osteoporotic fracture. GWAS have identified hundreds of loci influencing BMD, but few have been functionally analyzed. In this study, we show that SNPs within a BMD locus on chromosome 14q32.32 alter splicing and expression of PAR-1a/microtubule affinity regulating kinase 3 (MARK3), a conserved serine/threonine kinase known to regulate bioenergetics, cell division, and polarity. Mice lacking Mark3 either globally or selectively in osteoblasts have increased bone mass at maturity. RNA profiling from Mark3-deficient osteoblasts suggested changes in the expression of components of the Notch signaling pathway. Mark3-deficient osteoblasts exhibited greater matrix mineralization compared with controls that was accompanied by reduced Jag1/Hes1 expression and diminished downstream JNK signaling. Overexpression of Jag1 in Mark3-deficient osteoblasts both in vitro and in vivo normalized mineralization capacity and bone mass, respectively. Together, these findings reveal a mechanism whereby genetically regulated alterations in Mark3 expression perturb cell signaling in osteoblasts to influence bone mass.

Genetic analysis of osteoblast activity identifies Zbtb40 as a regulator of osteoblast activity and bone mass

Osteoporosis is a genetic disease characterized by progressive reductions in bone mineral density (BMD) leading to an increased risk of fracture. Over the last decade, genome-wide association studies (GWASs) have identified over 1000 associations for BMD. However, as a phenotype BMD is challenging as bone is a multicellular tissue affected by both local and systemic physiology. Here, we focused on a single component of BMD, osteoblast-mediated bone formation in mice, and identified associations influencing osteoblast activity on mouse Chromosomes (Chrs) 1, 4, and 17. The locus on Chr. 4 was in an intergenic region between Wnt4 and Zbtb40, homologous to a locus for BMD in humans. We tested both Wnt4 and Zbtb40 for a role in osteoblast activity and BMD. Knockdown of Zbtb40, but not Wnt4, in osteoblasts drastically reduced mineralization. Additionally, loss-of-function mouse models for both genes exhibited reduced BMD. Our results highlight that investigating the genetic basis of in vitro osteoblast mineralization can be used to identify genes impacting bone formation and BMD.

Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10⁻¹²) BMD locus on Chromosome 3@52.5 Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp^-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp^-/- mice displayed increased osteogenic differentiation. Lhfp^-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10⁻⁵) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass.

Osteoporosis is a common, chronic disease characterized by low bone mineral density (BMD) that puts millions of Americans at high risk of fracture. Variation in BMD in the general population is, in large part, determined by genetic factors. To identify novel genes influencing BMD, we performed a genome-wide association study in a panel of inbred mouse strains. We identified a locus on Chromosome 3 strongly associated with BMD. Using a combination of systems genetics approaches, we connected the expression of the Lhfp gene with BMD-associated genetic variants and predicted it influenced BMD by altering the activity of bone-forming osteoblasts. Using mice deficient in Lhfp, we demonstrated that Lhfp negatively regulates bone formation and BMD. These data suggest that inhibiting Lhfp may represent a novel therapeutic strategy to increase BMD and decrease the risk of fracture.

Integrating GWAS and Co-expression Network Data Identifies Bone Mineral Density Genes SPTBN1 and MARK3 and an Osteoblast Functional Module

Bone mineral density (BMD) is a highly heritable predictor of osteoporotic fracture. Genome-wide association studies (GWAS) for BMD have identified dozens of associations; yet, the genes responsible for most associations remain elusive. Here, we used a bone co-expression network to predict causal genes at BMD GWAS loci based on the premise that genes underlying a disease are often functionally related and functionally related genes are often co-expressed. By mapping genes implicated by BMD GWAS onto a bone co-expression network, we predicted and inferred the function of causal genes for 30 of 64 GWAS loci. We experimentally confirmed that two of the genes predicted to be causal, SPTBN1 and MARK3, are potentially responsible for the effects of GWAS loci on chromosomes 2p16.2 and 14q32.32, respectively. This approach provides a roadmap for the dissection of additional BMD GWAS associations. Furthermore, it should be applicable to GWAS data for a wide range of diseases.

Purification of restriction fragments containing replication intermediates from complex genomes for 2-D gel analysis

In order to perform 2-D gel analyses on restriction fragments from higher eukaryotic genomes, it is necessary to remove most of the linear, nonreplicating, fragments from the starting DNA preparation. This is so because the replication intermediates in a single-copy locus constitute such a minute fraction of all of the restriction fragments in a standard DNA preparation-whether isolated from synchronized or asynchronous cultures. Furthermore, the very long DNA strands that characterize higher eukaryotic genomes are inordinately subject to branch migration and shear. We have developed a method that results in significant enrichment of replicating fragments that largely maintain their branched intermediates. The method depends upon two important factors: (1) replicating fragments in higher eukaryotic nuclei appear to be attached to the nuclear matrix in a supercoiled fashion, and (2) partially single-stranded fragments (e.g., those containing replication forks) are selectively adsorbed to benzoylated naphthoylated DEAE (BND)-cellulose in high salt concentrations. By combining matrix-enrichment and BND-cellulose chromatography, it is possible to obtain preparations that are enriched as much as 200-fold over the starting genomic DNA, and are thus suitable for analysis on 2-D gels.

Bicc1 is a genetic determinant of osteoblastogenesis and bone mineral density

Patient bone mineral density (BMD) predicts the likelihood of osteoporotic fracture. While substantial progress has been made toward elucidating the genetic determinants of BMD, our understanding of the factors involved remains incomplete. Here, using a systems genetics approach in the mouse, we predicted that bicaudal C homolog 1 (Bicc1), which encodes an RNA-binding protein, is responsible for a BMD quantitative trait locus (QTL) located on murine chromosome 10. Consistent with this prediction, mice heterozygous for a null allele of Bicc1 had low BMD. We used a coexpression network-based approach to determine how Bicc1 influences BMD. Based on this analysis, we inferred that Bicc1 was involved in osteoblast differentiation and that polycystic kidney disease 2 (Pkd2) was a downstream target of Bicc1. Knock down of Bicc1 and Pkd2 impaired osteoblastogenesis, and Bicc1 deficiency-dependent osteoblast defects were rescued by Pkd2 overexpression. Last, in 2 human BMD genome-wide association (GWAS) meta-analyses, we identified SNPs in BICC1 and PKD2 that were associated with BMD. These results, in both mice and humans, identify Bicc1 as a genetic determinant of osteoblastogenesis and BMD and suggest that it does so by regulating Pkd2 transcript levels.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Cryptococcus neoformans var. grubii is a major human pathogen responsible for deadly meningoencephalitis in immunocompromised patients. Here, we report the sequencing and annotation of its genome. Evidence for extensive intron splicing, antisense transcription, non-coding RNAs, and alternative polyadenylation indicates the potential for highly intricate regulation of gene expression in this opportunistic pathogen. In addition, detailed molecular, genetic, and genomic studies were performed to characterize structural features of the genome, including centromeres and origins of replication. Finally, the phenotypic and genome re-sequencing analysis of a collection of isolates of the reference H99 strain resulting from laboratory passage revealed that microevolutionary processes during in vitro culturing of pathogenic fungi can impact virulence.

Bubble-seq analysis of the human genome reveals distinct chromatin-mediated mechanisms for regulating early- and late-firing origins

We have devised a method for isolating virtually pure and comprehensive libraries of restriction fragments that contained replication initiation sites (bubbles) in vivo. We have now sequenced and mapped the bubble-containing fragments from GM06990, a near-normal EBV-transformed lymphoblastoid cell line, and have compared origin distributions with a comprehensive replication timing study recently published for this cell line. We find that early-firing origins, which represent ∼32% of all origins, overwhelmingly represent zones, associate only marginally with active transcription units, are localized within large domains of open chromatin, and are significantly associated with DNase I hypersensitivity. Origin "density" falls from early- to mid-S-phase, but rises again in late S-phase to levels only 17% lower than in early S-phase. Unexpectedly, late origin density calculated on the 1-Mb scale increases as a function of increasing chromatin compaction. Furthermore, the median efficiency of origins in late-replicating, heterochromatic domains is only 25% lower than in early-replicating euchromatic loci. Thus, the activation of early- and late-firing origins must be regulated by quintessentially different mechanisms. The aggregate data can be unified into a model in which initiation site selection is driven almost entirely by epigenetic factors that fashion both the long-range and local chromatin environments, with underlying DNA sequence and local transcriptional activity playing only minor roles. Importantly, the comprehensive origin map we have prepared for GM06990 overlaps moderately well with origin maps recently reported for the genomes of four different human cell lines based on the distributions of small nascent strands.

Integrity of chromatin and replicating DNA in nuclei released from fission yeast by semi-automated grinding in liquid nitrogen

Background

Studies of nuclear function in many organisms, especially those with tough cell walls, are limited by lack of availability of simple, economical methods for large-scale preparation of clean, undamaged nuclei.

Findings

Here we present a useful method for nuclear isolation from the important model organism, the fission yeast, Schizosaccharomyces pombe. To preserve in vivo molecular configurations, we flash-froze the yeast cells in liquid nitrogen. Then we broke their tough cell walls, without damaging their nuclei, by grinding in a precision-controlled motorized mortar-and-pestle apparatus. The cryo-ground cells were resuspended and thawed in a buffer designed to preserve nuclear morphology, and the nuclei were enriched by differential centrifugation. The washed nuclei were free from contaminating nucleases and have proven well-suited as starting material for genome-wide chromatin analysis and for preparation of fragile DNA replication intermediates.

Conclusions

We have developed a simple, reproducible, economical procedure for large-scale preparation of endogenous-nuclease-free, morphologically intact nuclei from fission yeast. With appropriate modifications, this procedure may well prove useful for isolation of nuclei from other organisms with, or without, tough cell walls.

Bubble-chip analysis of human origin distributions demonstrates on a genomic scale significant clustering into zones and significant association with transcription

We have used a novel bubble-trapping procedure to construct nearly pure and comprehensive human origin libraries from early S- and log-phase HeLa cells, and from log-phase GM06990, a karyotypically normal lymphoblastoid cell line. When hybridized to ENCODE tiling arrays, these libraries illuminated 15.3%, 16.4%, and 21.8% of the genome in the ENCODE regions, respectively. Approximately half of the origin fragments cluster into zones, and their signals are generally higher than those of isolated fragments. Interestingly, initiation events are distributed about equally between genic and intergenic template sequences. While only 13.2% and 14.0% of genes within the ENCODE regions are actually transcribed in HeLa and GM06990 cells, 54.5% and 25.6% of zonal origin fragments overlap transcribed genes, most with activating chromatin marks in their promoters. Our data suggest that cell synchronization activates a significant number of inchoate origins. In addition, HeLa and GM06990 cells activate remarkably different origin populations. Finally, there is only moderate concordance between the log-phase HeLa bubble map and published maps of small nascent strands for this cell line.

Preferential association of irreversibly silenced E2F-target genes with pericentromeric heterochromatin in differentiated muscle cells

The heterochromatin-associated H3K9 tri-methylase Suv39h1 is involved in the permanent silencing of E2F target genes in differentiating but not in quiescent cells. Here, we tested the hypothesis that permanent silencing of E2F target genes is associated with their subnuclear positioning close to the pericentromeric heterochromatin compartment, enriched in Suv39h1. Using fluorescence in situ hybridization, we analyzed the subnuclear localization of three E2F target genes relative to the pericentromeric heterochromatin, in cycling fibroblasts or differentiating myoblasts. We observed that all three E2F-target genes have a tendency to relocate closer to the pericentromeric heterochromatin, only when cells differentiate and undergo an irreversible cell cycle withdrawal. These data suggest that repression of E2F target genes in cycling or in differentiating cells is achieved through distinct mechanisms. In differentiating cells, permanent silencing is driven by a Suv39h1-dependent H3K9 tri-methylation and positioning close to the heterochromatin compartment, whereas repression in cycling cells seems independent from subnuclear positioning and requires distinct H3K9 methylation levels.

Purification of restriction fragments containing replication intermediates from complex genomes for 2-D gel analysis

In order to perform 2-D gel analyses on restriction fragments from higher eukaryotic genomes, it is necessary to remove most of the linear, nonreplicating, fragments from the starting DNA preparation. This is so because the replication intermediates in a single-copy locus constitute such a minute fraction of all of the restriction fragments in a standard DNA preparation - whether isolated from synchronized or asynchronous cultures. Furthermore, the very long DNA strands that characterize higher eukaryotic genomes are inordinately subject to branch migration and shear. We have developed a method that results in significant enrichment of replicating fragments that largely maintain their branched intermediates. The method depends upon two important factors: (1) replicating fragments in higher eukaryotic nuclei appear to be attached to the nuclear matrix in a supercoiled fashion, and (2) partially single-stranded fragments (e.g., those containing replication forks) are selectively adsorbed to benzoylated napthoylated DEAE (BND)-cellulose in high salt conCentrations. By combining matrix-enrichment and BND-cellulose chromatography, it is possible to obtain preparations that are enriched as much as 200-fold over the starting genomic DNA and are thus suitable for analysis on 2-D gels.

Isolation of restriction fragments containing origins of replication from complex genomes

The identification and isolation of origins of replication from mammalian genomes has been a demanding task owing to the great complexity of these genomes. However, two methods have been refined in recent years each of which allows significant enrichment of recently activated origins of replication from asynchronous cell cultures. In one of these, nascent strands are melted from the long template DNA, and the small, origin-centered strands are isolated on sucrose gradients. The second method involves the selective entrapment of bubble-containing fragments in gelling agarose and their subsequent recovery and isolation by molecular cloning. Libraries prepared by this method from Chinese hamster and human cells have been shown to be extremely pure, and provide a renewable resource of origins that can be used as probes on microarrays or sequenced by high-throughput techniques to localize them within the genomic source. The bubble-trapping method is described here for asynchronous mammalian cells that grow with reasonable doubling times and from which nuclear matrices can be reliably prepared. The method for nuclear matrix preparation and enrichment of replication intermediates is described in an accompanying chapter entitled, "Purification of Restriction Fragments Containing Replication Intermediates from Mammalian Cells for 2-D Gel Analysis").

A revisionist replicon model for higher eukaryotic genomes

The replicon model devised to explain replication control in bacteria has served as the guiding paradigm in the search for origins of replication in the more complex genomes of eukaryotes. In Saccharomyces cerevisiae, this model has proved to be extremely useful, leading to the identification of specific genetic elements (replicators) and the interacting initiator proteins that activate them. However, replication control in organisms ranging from Schizosaccharomyces pombe to mammals is far more fluid: only a small number of origins seem to represent classic replicators, while the majority correspond to zones of inefficient, closely spaced start sites none of which are indispensable for origin activity. In addition, it is apparent that the epigenetic state of a given sequence largely determines its ability to be used as a replication initiation site. These conclusions were arrived at over a period of three decades, and required the development of several novel replicon mapping techniques, as well as new ways of examining the chromatin architecture of any sequence of interest. Recently, methods have been elaborated for isolating all of the active origins in the genomes of higher eukaryotes en masse. Microarray analyses and more recent high-throughput sequencing technology will allow all the origins to be mapped onto the chromosomes of any organism whose genome has been sequenced. With the advent of whole-genome studies on gene expression and chromatin composition, the field is now positioned to define both the genetic and epigenetic rules that govern origin activity.

Isolating apparently pure libraries of replication origins from complex genomes

Because of the complexity of higher eukaryotic genomes and the lack of a reliable autonomously replicating sequence (ARS) assay for isolating potential replicators, the identification of origins has proven to be extremely challenging and time consuming. We have developed a new origin-trapping method based on the partially circular nature of restriction fragments containing replication bubbles and have prepared a library of approximately 1,000 clones from early S phase CHO cells. When 15 randomly selected clones were analyzed by a stringent two-dimensional (2D) gel replicon mapping method, all were shown to correspond to active, early firing origins. Furthermore, most of these appear to derive from broad zones of potential sites, and the five that were analyzed in a time-course study are all inefficient. This bubble-trapping scheme will allow the construction of comprehensive origin libraries from any complex genome so that their natures and distributions vis-a-vis other chromosomal markers can be established.

Reproducible doxycycline-inducible transgene expression at specific loci generated by Cre-recombinase mediated cassette exchange

Comparative analysis of mutants using transfection is complicated by clones exhibiting variable levels of gene expression due to copy number differences and genomic position effects. Recombinase-mediated cassette exchange (RMCE) can overcome these problems by introducing the target gene into pre-determined chromosomal loci, but recombination between the available recombinase targeting sites can reduce the efficiency of targeted integration. We developed a new LoxP site (designated L3), which when used with the original LoxP site (designated L2), allows highly efficient and directional replacement of chromosomal DNA with incoming DNA. A total of six independent LoxP integration sites introduced either by homologous recombination or retroviral delivery were analyzed; 70-80% of the clones analyzed in hamster and human cells were correct recombinants. We combined the RMCE strategy with a new, tightly regulated tetracycline induction system to produce a robust, highly reliable system for inducible transgene expression. We observed stable inducible expression for over 1 month, with uniform expression in the cell population and between clones derived from the same integration site. This system described should find significant applications for studies requiring high level and regulated transgene expression and for determining the effects of various stresses or oncogenic conditions in vivo and in vitro.

Specific signals at the 3' end of the DHFR gene define one boundary of the downstream origin of replication

The Chinese hamster dihydrofolate reductase (DHFR) origin of replication consists of a 55-kb zone of potential initiation sites lying between the convergently transcribed DHFR and 2BE2121 genes. Two subregions within this zone (ori-beta/ori-beta' and ori-gamma) are preferred. In the DHFR-deficient variant, DR8, which has deleted a 14-kb sequence straddling the 3' end of the DHFR gene, early-firing origin activity in the downstream ori-beta/ori-beta' and ori-gamma regions is completely suppressed. We show that the critical deleted sequences reside within a 168-bp segment encompassing the intron 5/exon 6 boundary, exon 6, 54 bp of the 3' untranslated region (UTR), but not the three natural polyA sites. In wild-type cells, this sequence efficiently arrests transcription in a region a few kilobases downstream, which coincides with the 5' boundary of the replication initiation zone. In DR8, DHFR-specific transcripts efficiently use an alternative sixth exon (6c) and polyA signals near the middle of the former intergenic region to process primary transcripts. However, transcription proceeds to a position almost 35 kb downstream from these signals, and replication initiation can only be detected beyond this point. When the wild-type 168-bp 3' element is inserted into DR8 at the same position as alternative exon 6c, transcription is arrested efficiently and initiations occur almost immediately downstream. Thus, the normal 3' end of the DHFR gene constitutes a boundary element not only for the gene but also for the local origin of replication.

The promoter of the Chinese hamster ovary dihydrofolate reductase gene regulates the activity of the local origin and helps define its boundaries

The dihydrofolate reductase (DHFR) and 2BE2121 genes in the Chinese hamster are convergently transcribed in late G1 and ea ly S phase, and bracket an early-firing origin of replication that consists of a 55-kb zone of potential initiation sites. To test whether transcription through the DHFR gene is required to activate this origin in early S phase, we examined the two-dimension (2D) gel patterns of replication intermediates from several variants in which parts or all of the DHFR promote had been deleted. In those variants in which transcription was undetectable, initiation in the intergenic space was markedly suppressed (but not eliminated) in early S phase. Further more, replication of the locus required virtually the entire S period, as opposed to the usual 3-4 h. However, restoration of transcription with either the wild-type Chinese hamster promote or a Drosophila-based construct restored origin activity to the wild-type pattern. Surprisingly, 2D gel analysis of promote less variants revealed that initiation occurs at a low level in ea ly S phase not only in the intergenic region, but also in the body of the DHFR gene. The latter phenomenon has never been observed in the wild-type locus. These studies suggest that transcription through the gene normally increases the efficiency of origin firing in early S phase, but also suppresses initiation in the body of the gene, thus helping to define the boundaries of the downstream origin.

The matrix attachment region in the Chinese hamster dihydrofolate reductase origin of replication may be required for local chromatid separation

Centered in the Chinese hamster dihydrofolate reductase origin of replication is a prominent nuclear matrix attachment region (MAR). Indirect lines of evidence suggested that this MAR might be required for origin activation in early S phase. To test this possibility, we have deleted the MAR from a Chinese hamster ovary variant harboring a single copy of the dihydrofolate reductase locus. However, 2D gel replicon mapping shows that removal of the MAR has no significant effect either on the frequency or timing of initiation in this locus. Rather, fluorescence in situ hybridization studies on cells swollen under either neutral or alkaline conditions show that deletion of the MAR interferes with local separation of daughter chromatids. This surprising result provides direct genetic evidence that at least a subset of MARs performs an important biological function, possibly related to chromatid cohesion and separation.

The dihydrofolate reductase origin of replication does not contain any nonredundant genetic elements required for origin activity

The Chinese hamster dihydrofolate reductase (DHFR) origin of replication consists of a broad zone of potential initiation sites scattered throughout a 55-kb intergenic spacer, with at least three sites being preferred (ori-beta, ori-beta', and ori-gamma). We previously showed that deletion of the most active site or region (ori-beta) has no demonstrable effect on initiation in the remainder of the intergenic spacer nor on the time of replication of the DHFR locus as a whole. In the present study, we have now deleted ori-beta', both ori-beta and ori-beta', an 11-kb region just downstream from the DHFR gene, or the central approximately 40-kb core of the spacer. The latter two deletions together encompass >95% of the initiation sites that are normally used in this locus. Two-dimensional gel analysis shows that initiation still occurs in the early S phase in the remainder of the intergenic spacer in each of these deletion variants. Even removal of the 40-kb core fails to elicit a significant effect on the time of replication of the DHFR locus in the S period; indeed, in the truncated spacer that remains, the efficiency of initiation actually appears to increase relative to the corresponding region in the wild-type locus. Thus, if replicators control the positions of nascent strand start sites in this complex origin, either (i) there must be a very large number of redundant elements in the spacer, each of which regulates initiation only in its immediate environment, or (ii) they must lie outside the central core in which the vast majority of nascent strand starts occur.

Amplification of the human dihydrofolate reductase gene via double minutes is initiated by chromosome breaks

DNA sequence amplification is one of the most frequent manifestations of genomic instability in human tumors. We have shown previously that amplification of the dihydrofolate reductase (DHFR) gene in Chinese hamster cells is initiated by chromosome breaks, followed by bridge-breakage-fusion cycles that generate large intrachromosomal repeats; these are ultimately trimmed by an unknown process to smaller, more homogenous units manifested as homogenously staining chromosome regions (HSRs). However, in most human tumor cells, amplified DNA sequences are borne on unstable, extrachromosomal double minutes (DMs), which suggests the operation of a different amplification mechanism. In this study, we have isolated a large number of independent methotrexate-resistant human cell lines, all of which contained DHFR-bearing DMs. Surprisingly, all but one of these also had suffered partial or complete loss of one of the parental DHFR-bearing chromosomes. Cells in a few populations displayed what could be transient intermediates in the amplification process, including an initial HSR, its subsequent breakage, the appearance of DHFR-containing fragments, and, finally, DMs. Our studies suggest that HSRs and DMs both are initiated by chromosome breaks, but that cell types differ in how the extra sequences ultimately are processed and/or maintained.

Phosphoaminoglycosides inhibit SWI2/SNF2 family DNA-dependent molecular motor domains

Members of the SWI2/SNF2 family of proteins participate in an array of nucleic acid metabolic functions, including chromatin remodeling and transcription. The present studies identify a novel strategy to specifically inhibit the functional DNA-dependent adenosinetriphosphatase (ATPase) motor domain common to SWI2/SNF2 family members. We have identified preparations of phosphoaminoglycosides, which are natural products of aminoglycoside-resistant bacteria, as inhibitors of the in vitro activities of three SWI2/SNF2 family members. These compounds inhibit the ATPase activity of the active DNA-dependent ATPase A domain (ADAAD) by competing with respect to DNA and thus have no effect on DNA-independent ATPases or on RNA-dependent ATPases. Within the superfamily of DNA-dependent ATPases, these compounds are most potent toward SWI2/SNF2 family members and less potent toward other DNA-dependent ATPases. We demonstrate that it is feasible to target DNA-dependent ATPases of a particular type without affecting the function of other ATPases. As the SWI2/SNF2 proteins have been proposed to function in all aspects of DNA metabolism, this paper provides an archetype for development of DNA metabolic inhibitors.

Dispersive initiation of replication in the Chinese hamster rhodopsin locus

Several higher eukaryotic replication origins appear to be composed of broad zones of potential nascent strand start sites, while others are more circumscribed, resembling those of yeast, bacteria, and viruses. The most delocalized origin identified so far is approximately 55 kb in length and lies between the convergently transcribed dihydrofolate reductase (DHFR) and the 2BE2121 genes on chromosome 2 in the Chinese hamster genome. In some of our studies, we have utilized the rhodopsin origin as an early replicating internal standard for assessing the effects of deleting various parts of the DHFR locus on DHFR origin activity. However, it had not been previously established that the rhodopsin locus was located at a site far enough away to be immune to such deletions, nor had the mechanism of initiation at this origin been characterized. In the present study, we have localized the rhodopsin domain to a pair of small metacentric chromosomes and have used neutral/neutral 2-D gel replicon mapping to show that initiation in this origin is also highly delocalized, encompassing a region more than 50 kb in length that includes the nontranscribed rhodopsin gene itself. The initiation zone is flanked at least on one end by an actively transcribed gene that does not support initiation. Thus, the DHFR and rhodopsin origins belong to a class of complex, polydisperse origins that appears to be unique to higher eukaryotic cells.

A eukaryotic SWI2/SNF2 domain, an exquisite detector of double-stranded to single-stranded DNA transition elements

Many members of the SWI2/SNF2 family of adenosine triphosphatases participate in the assembly/disassembly of multiprotein complexes involved in the DNA metabolic processes of transcription, recombination, repair, and chromatin remodeling. The DNA molecule serves as an essential effector or catalyst for most of the members of this particular class of proteins, and the structure of the DNA may be more important than the nucleotide sequence. Inspection of the DNA structure at sites where multiprotein complexes are assembled/disassembled for these various DNA metabolic processes reveals the presence of a common element: a double-stranded to single-stranded transition region. We now show that this DNA element is crucial for the ATP hydrolytic function of an SWI2/SNF2 family member: DNA-dependent ATPase A. We further demonstrate that a domain containing the seven helicase-related motifs that are common to the SWI2/SNF2 family of proteins mediates the interaction with the DNA, yielding specific DNA structural recognition. This study forms a primary step toward understanding the physico-biochemical nature of the interaction between a particular class of DNA-dependent ATPase and their DNA effectors. Furthermore, this study provides a foundation for development of mechanisms to specifically target this class of DNA-dependent ATPases.

Distal sequences, but not ori-beta/OBR-1, are essential for initiation of DNA replication in the Chinese hamster DHFR origin

In the Chinese hamster dihydrofolate reductase replication initiation zone, the ori-beta locus is preferred over other start sites. To test the hypothesis that ori-beta contains a genetic replicator, we restored a deletion in the 3' end of the DHFR gene with a cosmid that provides the missing sequence and simultaneously knocks out the downstream ori-beta locus. Replication initiates normally in ori-beta knockout cell lines, and the DHFR domain is still synthesized in early S phase. However, initiation is completely suppressed in the starting deletion variant lacking the 3' end of the gene. We conclude that ori-beta does not contain an essential replicator, but that distant sequence elements have profound effects on origin activity in this locus.

DNA-dependent adenosinetriphosphatase A: immunoaffinity purification and characterization of immunological reagents

We describe an immunoaffinity purification of DNA-dependent ATPase A from fetal calf thymus. The rapid purification increases the yield of enzymatically active enzyme approximately 4-fold, with up to a 7-fold increase in specific activity, and significantly improves the yield of a higher molecular weight species of ATPase A. In the presence of a denatured calf thymus DNA effector, the immunoaffinity-purified enzyme has a specific activity that is more than 10-fold higher than reported for any other eukaryotic DNA-dependent ATPase and 100-fold higher than most others. The improvement in yield has allowed several polypeptides to be identified using monoclonal antibodies, and these polypeptides are demonstrated to be structurally related by partial peptide mapping with N-chlorosuccinimide. The preferred DNA effector for ATP hydrolysis continues to be a DNA primer-template junction with an adjacent stretch of single-stranded DNA. We have used the immunoaffinity-purified enzyme to develop additional stable murine hybridoma monoclones, resulting in a bank of antibodies that recognize a number of different epitopes. All of the monoclonal antibodies react with both calf thymus DNA-dependent ATPase A and bacteriophage T4 gene 44 protein, a DNA-dependent ATPase essential for DNA replication in the bacteriophage T4 system. These monoclonal antibodies should facilitate the development of our understanding with respect to the role and regulation of DNA-dependent ATPases in eukaryotic DNA replication.

Laser cross-linking of proteins to nucleic acids. II. Interactions of the bacteriophage T4 DNA replication polymerase accessory proteins complex with DNA

In this paper we examine the interactions of the polymerase accessory proteins subassembly of the bacteriophage T4 DNA replication complex, using single-pulse ultraviolet laser excitation to induce protein-nucleic acid cross-links. The laser-induced cross-linking permits effective "freezing" of the instantaneous equilibrium state of the complex and thus provides a mechanism to dissect the individual protein-nucleic acid interactions involved in complex assembly. We find that the binding of the gene 44, 62, and 45 proteins is dependent not only on the presence of each of the other proteins, but also on the presence of adenine nucleotide cofactors. We find that the nonhydrolyzable analogs of ATP often behave more like ADP than ATP in these experiments. Gene 45 protein is able to induce an increase in cross-linking of the gp44/62 complex to nucleic acids, and this increased cross-linking correlates with changes in the apparent Km of the gp44/62 complex for polynucleotides and with changes in Vmax during ATP hydrolysis. Our results suggest that the enhanced DNA binding is predominately through the gene 62 protein and not the ATPase catalytic subunit (gene 44 protein). Thus the gene 62 protein seems to play an integral role in gp45-mediated enhancement of the ATP hydrolytic activity of gp44. These results are summarized and integrated in the form of a model for the multiple interactions of the accessory proteins with DNA and one another in the presence of mononucleotide cofactors and substrates.

Probing the energetics of oligo(dT).poly(dA) by laser cross-linking

Experimentally determined changes in free energy (delta G(o)) for thymine-thymine interactions occurring in oligo(dT).poly(dA) are dependent on the method used for preparation of the double-stranded template. A rapid laser cross-linking technique was used to examine the equilibrium between oligomers of (dT) bound to either poly(dA) or poly(rA). The single-pulse (4-6 nsec) ultraviolet laser excitation of these polynucleotides causes pyrimidine dimer formation between contiguous oligo(dT) molecules, resulting in a "ligation" of the oligomers. Analysis of the resulting data using standard binding isotherms allowed determination of the degree of cooperativity existing between oligomers. Using the cooperativity, delta G(o), delta H(o), and delta S(o) are calculated, thereby providing thermodynamic parameters for this interaction. The measured cooperativity of oligo(dT) molecule interactions allows direct calculation of the number of 3' ends available as nicked structures or the number of 3' ends associated with gaps for oligo(dT).poly(dA) when used as a substrate for DNA synthesis.

DNA-dependent adenosinetriphosphatase A is the eukaryotic analogue of the bacteriophage T4 gene 44 protein: immunological identity of DNA replication-associated ATPases

We report the construction of three stable murine hybridomas that secrete monoclonal antibodies which recognize calf thymus DNA-dependent adenosinetriphosphatase A. All three of the antibodies react specifically with calf thymus ATPase A and the gene 44 protein from the bacteriophage T4 DNA-dependent ATPase. Each of the three anti-ATPase A antibodies appears to recognize a different epitope and none of the antibodies inhibit DNA-dependent ATP hydrolysis by ATPase A. Furthermore, one of the antibodies has been shown to react with two different preparations of HeLa cell DNA-dependent ATPases and a yeast DNA-dependent ATPase, all of which have been implicated in the enzymology of DNA replication. These findings provide strong evidence for the role of ATPase A in DNA replication. These observations lead us to conclude that, apart from the nucleotide binding sites, there are at least three epitopes common to both the bacteriophage and eukaryotic DNA-dependent ATPases that we have examined and that the different preparations of the eukaryotic ATPases contain the same DNA-dependent ATPase.