Cell type-specific network analysis in Diversity Outbred mice identifies genes potentially responsible for human bone mineral density GWAS associations.

Genome-wide association studies (GWASs) have identified many sources of genetic variation associated with bone mineral density (BMD), a clinical predictor of fracture risk and osteoporosis. Aside from the identification of causal genes, other difficult challenges to informing GWAS include characterizing the roles of predicted causal genes in disease and providing additional functional context, such as the cell type predictions or biological pathways in which causal genes operate. Leveraging single-cell transcriptomics (scRNA-seq) can assist in informing BMD GWAS by linking disease-associated variants to genes and providing a cell type context for which these causal genes drive disease. Here, we use large-scale scRNA-seq data from bone marrow-derived stromal cells cultured under osteogenic conditions (BMSC-OBs) from Diversity Outbred (DO) mice to generate cell type-specific networks and contextualize BMD GWAS-implicated genes. Using trajectories inferred from the scRNA-seq data, we identify networks enriched with genes that exhibit the most dynamic changes in expression across trajectories. We discover 21 network driver genes, which are likely to be causal for human BMD GWAS associations that colocalize with expression/splicing quantitative trait loci (eQTL/sQTL). These driver genes, including Fgfrl1 and Tpx2, along with their associated networks, are predicted to be novel regulators of BMD via their roles in the differentiation of mesenchymal lineage cells. In this work, we showcase the use of single-cell transcriptomics from mouse bone-relevant cells to inform human BMD GWAS and prioritize genetic targets with potential causal roles in the development of osteoporosis.

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.

Systems genetics in diversity outbred mice inform BMD GWAS and identify determinants of bone strength.

Genome-wide association studies (GWASs) for osteoporotic traits have identified over 1000 associations; however, their impact has been limited by the difficulties of causal gene identification and a strict focus on bone mineral density (BMD). Here, we use Diversity Outbred (DO) mice to directly address these limitations by performing a systems genetics analysis of 55 complex skeletal phenotypes. We apply a network approach to cortical bone RNA-seq data to discover 66 genes likely to be causal for human BMD GWAS associations, including the genes SERTAD4 and GLT8D2. We also perform GWAS in the DO for a wide-range of bone traits and identify Qsox1 as a gene influencing cortical bone accrual and bone strength. In this work, we advance our understanding of the genetics of osteoporosis and highlight the ability of the mouse to inform human genetics.

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-12) 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-5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass.

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.

Network Analysis Implicates Alpha-Synuclein (Snca) in the Regulation of Ovariectomy-Induced Bone Loss.

The postmenopausal period in women is associated with decreased circulating estrogen levels, which accelerate bone loss and increase the risk of fracture. Here, we gained novel insight into the molecular mechanisms mediating bone loss in ovariectomized (OVX) mice, a model of human menopause, using co-expression network analysis. Specifically, we generated a co-expression network consisting of 53 gene modules using expression profiles from intact and OVX mice from a panel of inbred strains. The expression of four modules was altered by OVX, including module 23 whose expression was decreased by OVX across all strains. Module 23 was enriched for genes involved in the response to oxidative stress, a process known to be involved in OVX-induced bone loss. Additionally, module 23 homologs were co-expressed in human bone marrow. Alpha synuclein (Snca) was one of the most highly connected "hub" genes in module 23. We characterized mice deficient in Snca and observed a 40% reduction in OVX-induced bone loss. Furthermore, protection was associated with the altered expression of specific network modules, including module 23. In summary, the results of this study suggest that Snca regulates bone network homeostasis and ovariectomy-induced bone loss.

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.

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" (Chapter 16 ).

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.

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.

A winding road to origin discovery.

Studies in our laboratory over the last three decades have shown that the Chinese hamster dihydrofolate reductase (DHFR) origin of replication corresponds to a broad zone of inefficient initiation sites distributed throughout the spacer between the convergently transcribed DHFR and 2BE2121 genes. It is clear from mutational analysis that none of these sites is genetically required for controlling origin activity. However, the integrity of the promoter of the DHFR gene is needed to activate the downstream origin, while the 3' processing signals prevent invasion and inactivation of the downstream origin by transcription forks. Several other origins in metazoans have been shown to correspond to zones of inefficient sites, while a different subset appears to be similar to the fixed replicators that characterize origins in S. cerevisiae and lower organisms. These observations have led us to suggest a model in which the mammalian genome is dotted with a hierarchy of degenerate, redundant, and inefficient replicators at intervals of a kilobase or less, some of which may have evolved to be highly circumscribed and efficient. The activities of initiation sites are proposed to be largely regulated by local transcription and chromatin architecture. Recently, we and others have devised strategies for identifying active origins on a genome-wide scale in order to define their distributions between fixed and dispersive origin types and to detect relationships among origins, genes, and epigenetic markers. The global pictures emerging are suggestive but far from complete and appear to be plagued by some of the same uncertainties that have led to conflicting views of individual origins in the past (particularly DHFR). In this paper, we will trace the history of origin discovery in mammalian genomes, primarily using the well-studied DHFR origin as a model, because it has been analyzed by nearly every available origin mapping technique in several different laboratories, while many origins have been identified by only one. We will address the strengths and shortcomings of the various methods utilized to identify and characterize origins in complex genomes and will point out how we and others were sometimes led astray by false assumptions and biases, as well as insufficient information. The goal is to help guide future experiments that will provide a truly comprehensive and accurate portrait of origins and their regulation. After all, in the words of George Santayana, "Those who do not learn from history are doomed to repeat it."