Integrated genomic and functional analyses of human skin-associated Staphylococcus reveal extensive inter- and intra-species diversity

Human skin is stably colonized by a distinct microbiota that functions together with epidermal cells to maintain a protective physical barrier. Staphylococcus, a prominent genus of the skin microbiota, participates in colonization resistance, tissue repair, and host immune regulation in strain-specific manners. To unlock the potential of engineering skin microbial communities, we aim to characterize the diversity of this genus within the context of the skin environment. We reanalyzed an extant 16S rRNA amplicon dataset obtained from distinct body sites of healthy volunteers, providing a detailed biogeographic depiction of staphylococcal species that colonize our skin. S. epidermidis, S. capitis, and S. hominis were the most abundant staphylococcal species present in all volunteers and were detected at all body sites. Pan-genome analysis of isolates from these three species revealed that the genus-core was dominated by central metabolism genes. Species-restricted-core genes encoded known host colonization functions. The majority (~68%) of genes were detected only in a fraction of isolate genomes, underscoring the immense strain-specific gene diversity. Conspecific genomes grouped into phylogenetic clades, exhibiting body site preference. Each clade was enriched for distinct gene sets that are potentially involved in site tropism. Finally, we conducted gene expression studies of select isolates showing variable growth phenotypes in skin-like medium. In vitro expression revealed extensive intra- and inter-species gene expression variation, substantially expanding the functional diversification within each species. Our study provides an important resource for future ecological and translational studies to examine the role of shared and strain-specific staphylococcal genes within the skin environment.

Expanded microbiome niches of RAG-deficient patients

The complex interplay between microbiota and immunity is important to human health. To explore how altered adaptive immunity influences the microbiome, we characterize skin, nares, and gut microbiota of patients with recombination-activating gene (RAG) deficiency-a rare genetically defined inborn error of immunity (IEI) that results in a broad spectrum of clinical phenotypes. Integrating de novo assembly of metagenomes from RAG-deficient patients with reference genome catalogs provides an expansive multi-kingdom view of microbial diversity. RAG-deficient patient microbiomes exhibit inter-individual variation, including expansion of opportunistic pathogens (e.g., Corynebacterium bovis, Haemophilus influenzae), and a relative loss of body site specificity. We identify 35 and 27 bacterial species derived from skin/nares and gut microbiomes, respectively, which are distinct to RAG-deficient patients compared to healthy individuals. Underscoring IEI patients as potential reservoirs for viral persistence and evolution, we further characterize the colonization of eukaryotic RNA viruses (e.g., Coronavirus 229E, Norovirus GII) in this patient population.

Integrated genomic and functional analyses of human skin-associated Staphylococcus reveals extensive inter- and intra-species diversity

Human skin is stably colonized by a distinct microbiota that functions together with epidermal cells to maintain a protective physical barrier. Staphylococcus, a prominent genus of the skin microbiota, participates in colonization resistance, tissue repair, and host immune regulation in strain specific manners. To unlock the potential of engineering skin microbial communities, we aim to fully characterize the functional diversity of this genus within the context of the skin environment. We conducted metagenome and pan-genome analyses of isolates obtained from distinct body sites of healthy volunteers, providing a detailed biogeographic depiction of staphylococcal species that colonize our skin. S. epidermidis, S. capitis, and S. hominis were the most abundant species present in all volunteers and were detected at all body sites. Pan-genome analysis of these three species revealed that the genus-core was dominated by central metabolism genes. Species-specific core genes were enriched in host colonization functions. The majority (~68%) of genes were detected only in a fraction of isolate genomes, underscoring the immense strain-specific gene diversity. Conspecific genomes grouped into phylogenetic clades, exhibiting body site preference. Each clade was enriched for distinct gene-sets that are potentially involved in site tropism. Finally, we conducted gene expression studies of select isolates showing variable growth phenotypes in skin-like medium. In vitro expression revealed extensive intra- and inter-species gene expression variation, substantially expanding the functional diversification within each species. Our study provides an important resource for future ecological and translational studies to examine the role of shared and strain-specific staphylococcal genes within the skin environment.

Plasmid Dissemination and Selection of a Multidrug-Resistant Klebsiella pneumoniae Strain during Transplant-Associated Antibiotic Therapy

Antibiotics, which are used both to prevent and to treat infections, are a mainstay therapy for lifesaving procedures such as transplantation. For this reason, and many others, increased antibiotic resistance among human-associated pathogens, such as the carbapenem-resistant Enterobacteriaceae species, is of grave concern. In this study, we report on a hematopoietic stem cell transplant recipient in whom cultures detected the emergence of carbapenem resistance and spread across five strains of bacteria that persisted for over a year. Carbapenem resistance in Citrobacter freundii, Enterobacter cloacae, Klebsiella aerogenes, and Klebsiella pneumoniae was linked to a pair of plasmids, each carrying the Klebsiella pneumoniae carbapenemase gene (blaKPC). Surveillance cultures identified a carbapenem-susceptible strain of Citrobacter freundii that may have become resistant through horizontal gene transfer of these plasmids. Selection of a multidrug-resistant Klebsiella pneumoniae strain was also detected following combination antibiotic therapy. Here we report a plasmid carrying the blaKPC gene with broad host range that poses the additional threat of spreading to endogenous members of the human gut microbiome.IMPORTANCE Antibiotic-resistant bacteria are a serious threat to medically fragile patient populations. The spread of antibiotic resistance through plasmid-mediated mechanisms is of grave concern as it can lead to the conversion of endogenous patient-associated strains to difficult-to-treat pathogens.

Investigation of a Cluster of Sphingomonas koreensis Infections

BACKGROUND

Plumbing systems are an infrequent but known reservoir for opportunistic microbial pathogens that can infect hospitalized patients. In 2016, a cluster of clinical sphingomonas infections prompted an investigation.

METHODS

We performed whole-genome DNA sequencing on clinical isolates of multidrug-resistant Sphingomonas koreensis identified from 2006 through 2016 at the National Institutes of Health (NIH) Clinical Center. We cultured S. koreensis from the sinks in patient rooms and performed both whole-genome and shotgun metagenomic sequencing to identify a reservoir within the infrastructure of the hospital. These isolates were compared with clinical and environmental S. koreensis isolates obtained from other institutions.

RESULTS

The investigation showed that two isolates of S. koreensis obtained from the six patients identified in the 2016 cluster were unrelated, but four isolates shared more than 99.92% genetic similarity and were resistant to multiple antibiotic agents. Retrospective analysis of banked clinical isolates of sphingomonas from the NIH Clinical Center revealed the intermittent recovery of a clonal strain over the past decade. Unique single-nucleotide variants identified in strains of S. koreensis elucidated the existence of a reservoir in the hospital plumbing. Clinical S. koreensis isolates from other facilities were genetically distinct from the NIH isolates. Hospital remediation strategies were guided by results of microbiologic culturing and fine-scale genomic analyses.

CONCLUSIONS

This genomic and epidemiologic investigation suggests that S. koreensis is an opportunistic human pathogen that both persisted in the NIH Clinical Center infrastructure across time and space and caused health care-associated infections. (Funded by the NIH Intramural Research Programs.).

Expanded skin virome in DOCK8-deficient patients

Human microbiome studies have revealed the intricate interplay of host immunity and bacterial communities to achieve homeostatic balance. Healthy skin microbial communities are dominated by bacteria with low viral representation1-3, mainly bacteriophage. Specific eukaryotic viruses have been implicated in both common and rare skin diseases, but cataloging skin viral communities has been limited. Alterations in host immunity provide an opportunity to expand our understanding of microbial-host interactions. Primary immunodeficient patients manifest with various viral, bacterial, fungal, and parasitic infections, including skin infections4. Dedicator of cytokinesis 8 (DOCK8) deficiency is a rare primary human immunodeficiency characterized by recurrent cutaneous and systemic infections, as well as atopy and cancer susceptibility5. DOCK8, encoding a guanine nucleotide exchange factor highly expressed in lymphocytes, regulates actin cytoskeleton, which is critical for migration through collagen-dense tissues such as skin6. Analyzing deep metagenomic sequencing data from DOCK8-deficient skin samples demonstrated a notable increase in eukaryotic viral representation and diversity compared with healthy volunteers. De novo assembly approaches identified hundreds of novel human papillomavirus genomes, illuminating microbial dark matter. Expansion of the skin virome in DOCK8-deficient patients underscores the importance of immune surveillance in controlling eukaryotic viral colonization and infection.

Evidence that synthetic lethality underlies the mutual exclusivity of oncogenic KRAS and EGFR mutations in lung adenocarcinoma

Human lung adenocarcinomas (LUAD) contain mutations in EGFR in ∼15% of cases and in KRAS in ∼30%, yet no individual adenocarcinoma appears to carry activating mutations in both genes, a finding we have confirmed by re-analysis of data from over 600 LUAD. Here we provide evidence that co-occurrence of mutations in these two genes is deleterious. In transgenic mice programmed to express both mutant oncogenes in the lung epithelium, the resulting tumors express only one oncogene. We also show that forced expression of a second oncogene in human cancer cell lines with an endogenous mutated oncogene is deleterious. The most prominent features accompanying loss of cell viability were vacuolization, other changes in cell morphology, and increased macropinocytosis. Activation of ERK, p38 and JNK in the dying cells suggests that an overly active MAPK signaling pathway may mediate the phenotype. Together, our findings indicate that mutual exclusivity of oncogenic mutations may reveal unexpected vulnerabilities and therapeutic possibilities.

DOI:http://dx.doi.org/10.7554/eLife.06907.001

A person develops cancer when changes in a cell's DNA (called mutations) allow the cell to grow rapidly and spread around the body. The mutated genes are often involved in controlling the growth of cells, such as two genes called EGFR and KRAS, which are associated with forms of lung cancer.

In a type of lung cancer called adenocarcinoma, the KRAS gene is mutated in about one-third of tumors and the EGFR gene is mutated in about 15%. However, the two mutations rarely or never occur in the same tumor. This could be because the effects of the mutations overlap, so that cells with both mutations have no advantages over cells with just one. Alternatively, it is possible that having both mutations may be harmful to tumor cells.

Here, Unni, Lockwood et al. analyzed genetic data from over 600 lung tumors and confirmed that none of them have cancer-causing mutations in both KRAS and EGFR. Then, Unni, Lockwood et al. carried out experiments using genetically engineered mice with mutated forms of both KRAS and EGFR that are activated by a drug called doxycycline. As expected, the mice developed lung tumors when exposed to the drug, but these tumors didn't grow any faster than mouse tumors that had mutations in only one of the genes. In the mice with both mutant genes, only one of the two genes was actually active in most of the tumor cells.

Unni, Lockwood et al. manipulated human lung tumor cells in the laboratory so that the cells had mutated versions of both genes. These cells developed serious abnormalities and died, which may be due to the over-activation of a communication pathway within the cells called MAPK signaling. The next challenges are to understand why the combination of these two mutant genes kills these cancer cells and to look for other combinations of mutations that can be toxic to cancer cells. In the future, it might be possible to develop drugs that can mimic the effects of these gene mutations to treat cancers.

DOI:http://dx.doi.org/10.7554/eLife.06907.002

Identification of candidate genes involved in coronary artery calcification by transcriptome sequencing of cell lines

Background

Massively-parallel cDNA sequencing (RNA-Seq) is a new technique that holds great promise for cardiovascular genomics. Here, we used RNA-Seq to study the transcriptomes of matched coronary artery disease cases and controls in the ClinSeq® study, using cell lines as tissue surrogates.

Results

Lymphoblastoid cell lines (LCLs) from 16 cases and controls representing phenotypic extremes for coronary calcification were cultured and analyzed using RNA-Seq. All cell lines were then independently re-cultured and along with another set of 16 independent cases and controls, were profiled with Affymetrix microarrays to perform a technical validation of the RNA-Seq results. Statistically significant changes (p < 0.05) were detected in 186 transcripts, many of which are expressed at extremely low levels (5–10 copies/cell), which we confirmed through a separate spike-in control RNA-Seq experiment. Next, by fitting a linear model to exon-level RNA-Seq read counts, we detected signals of alternative splicing in 18 transcripts. Finally, we used the RNA-Seq data to identify differential expression (p < 0.0001) in eight previously unannotated regions that may represent novel transcripts. Overall, differentially expressed genes showed strong enrichment (p = 0.0002) for prior association with cardiovascular disease. At the network level, we found evidence for perturbation in pathways involving both cardiovascular system development and function as well as lipid metabolism.

Conclusions

We present a pilot study for transcriptome involvement in coronary artery calcification and demonstrate how RNA-Seq analyses using LCLs as a tissue surrogate may yield fruitful results in a clinical sequencing project. In addition to canonical gene expression, we present candidate variants from alternative splicing and novel transcript detection, which have been unexplored in the context of this disease.

Balancing selection maintains a form of ERAP2 that undergoes nonsense-mediated decay and affects antigen presentation

A remarkable characteristic of the human major histocompatibility complex (MHC) is its extreme genetic diversity, which is maintained by balancing selection. In fact, the MHC complex remains one of the best-known examples of natural selection in humans, with well-established genetic signatures and biological mechanisms for the action of selection. Here, we present genetic and functional evidence that another gene with a fundamental role in MHC class I presentation, endoplasmic reticulum aminopeptidase 2 (ERAP2), has also evolved under balancing selection and contains a variant that affects antigen presentation. Specifically, genetic analyses of six human populations revealed strong and consistent signatures of balancing selection affecting ERAP2. This selection maintains two highly differentiated haplotypes (Haplotype A and Haplotype B), with frequencies 0.44 and 0.56, respectively. We found that ERAP2 expressed from Haplotype B undergoes differential splicing and encodes a truncated protein, leading to nonsense-mediated decay of the mRNA. To investigate the consequences of ERAP2 deficiency on MHC presentation, we correlated surface MHC class I expression with ERAP2 genotypes in primary lymphocytes. Haplotype B homozygotes had lower levels of MHC class I expressed on the surface of B cells, suggesting that naturally occurring ERAP2 deficiency affects MHC presentation and immune response. Interestingly, an ERAP2 paralog, endoplasmic reticulum aminopeptidase 1 (ERAP1), also shows genetic signatures of balancing selection. Together, our findings link the genetic signatures of selection with an effect on splicing and a cellular phenotype. Although the precise selective pressure that maintains polymorphism is unknown, the demonstrated differences between the ERAP2 splice forms provide important insights into the potential mechanism for the action of selection.

It has long been known that the extremely high levels of genetic diversity present in the major histocompatibility locus (MHC) are due to balancing selection, a type of natural selection that maintains advantageous genetic diversity in populations. The MHC encodes for molecules required for a type of antigen presentation that mediates detection of infected and cancerous cells by the immune system; the genetic diversity of the MHC thus ensures an adequate response to the wide variety of pathogens that humans encounter. Here, we show that other genes involved in the same antigen-presentation pathway are also subject to balancing selection in humans. Specifically, we show that balancing selection acts to maintain two forms of the endoplasmic reticulum aminopeptidase 2 gene (ERAP2), which encodes a protein also involved in antigen presentation. Although the two ERAP2 forms are present in a similar frequency (close to 0.5), they are associated with differences with respect to the levels of MHC molecules on the cell surface of immune cells. In summary, our findings show that natural selection maintains variants of ERAP2 that affect immune surveillance; they also establish ERAP2 as one of the few examples of balancing selection in humans where the selected variant, its functional consequences, and its influence in interpersonal diversity are known.

Identification of neural crest and glial enhancers at the mouse Sox10 locus through transgenesis in zebrafish

Sox10 is a dynamically regulated transcription factor gene that is essential for the development of neural crest-derived and oligodendroglial populations. Developmental genes often require multiple regulatory sequences that integrate discrete and overlapping functions to coordinate their expression. To identify Sox10 cis-regulatory elements, we integrated multiple model systems, including cell-based screens and transposon-mediated transgensis in zebrafish, to scrutinize mammalian conserved, noncoding genomic segments at the mouse Sox10 locus. We demonstrate that eight of 11 Sox10 genomic elements direct reporter gene expression in transgenic zebrafish similar to patterns observed in transgenic mice, despite an absence of observable sequence conservation between mice and zebrafish. Multiple segments direct expression in overlapping populations of neural crest derivatives and glial cells, ranging from pan-Sox10 and pan-neural crest regulatory control to the modulation of expression in subpopulations of Sox10-expressing cells, including developing melanocytes and Schwann cells. Several sequences demonstrate overlapping spatial control, yet direct expression in incompletely overlapping developmental intervals. We were able to partially explain neural crest expression patterns by the presence of head to head SoxE family binding sites within two of the elements. Moreover, we were able to use this transcription factor binding site signature to identify the corresponding zebrafish enhancers in the absence of overall sequence homology. We demonstrate the utility of zebrafish transgenesis as a high-fidelity surrogate in the dissection of mammalian gene regulation, especially those with dynamically controlled developmental expression.

Functional analyses of glycyl-tRNA synthetase mutations suggest a key role for tRNA-charging enzymes in peripheral axons

Charcot-Marie-Tooth disease type 2D (CMT2D) and distal spinal muscular atrophy type V (dSMA-V) are axonal neuropathies characterized by a phenotype that is more severe in the upper extremities. We previously implicated mutations in the gene encoding glycyl-tRNA synthetase (GARS) as the cause of CMT2D and dSMA-V. GARS is a member of the family of aminoacyl-tRNA synthetases responsible for charging tRNA with cognate amino acids; GARS ligates glycine to tRNA(Gly). Here, we present functional analyses of disease-associated GARS mutations and show that there are not any significant mutation-associated changes in GARS expression levels; that the majority of identified GARS mutations modeled in yeast severely impair viability; and that, in most cases, mutant GARS protein mislocalizes in neuronal cells. Indeed, four of the five mutations studied show loss-of-function features in at least one assay, suggesting that tRNA-charging deficits play a role in disease pathogenesis. Finally, we detected endogenous GARS-associated granules in the neurite projections of cultured neurons and in the peripheral nerve axons of normal human tissue. These data are particularly important in light of the recent identification of CMT-associated mutations in another tRNA synthetase gene [YARS (tyrosyl-tRNA synthetase gene)]. Together, these findings suggest that tRNA-charging enzymes play a key role in maintaining peripheral axons.

Deletion of long-range sequences at Sox10 compromises developmental expression in a mouse model of Waardenburg-Shah (WS4) syndrome

The transcription factor SOX10 is mutated in the human neurocristopathy Waardenburg-Shah syndrome (WS4), which is characterized by enteric aganglionosis and pigmentation defects. SOX10 directly regulates genes expressed in neural crest lineages, including the enteric ganglia and melanocytes. Although some SOX10 target genes have been reported, the mechanisms by which SOX10 expression is regulated remain elusive. Here, we describe a transgene-insertion mutant mouse line (Hry) that displays partial enteric aganglionosis, a loss of melanocytes, and decreased Sox10 expression in homozygous embryos. Mutation analysis of Sox10 coding sequences was negative, suggesting that non-coding regulatory sequences are disrupted. To isolate the Hry molecular defect, Sox10 genomic sequences were collected from multiple species, comparative sequence analysis was performed and software was designed (ExactPlus) to identify identical sequences shared among species. Mutation analysis of conserved sequences revealed a 15.9 kb deletion located 47.3 kb upstream of Sox10 in Hry mice. ExactPlus revealed three clusters of highly conserved sequences within the deletion, one of which shows strong enhancer potential in cultured melanocytes. These studies: (i) present a novel hypomorphic Sox10 mutation that results in a WS4-like phenotype in mice; (ii) demonstrate that a 15.9 kb deletion underlies the observed phenotype and likely removes sequences essential for Sox10 expression; (iii) combine a novel in silico method for comparative sequence analysis with in vitro functional assays to identify candidate regulatory sequences deleted in this strain. These studies will direct further analyses of Sox10 regulation and provide candidate sequences for mutation detection in WS4 patients lacking a SOX10-coding mutation.

Glycyl tRNA synthetase mutations in Charcot-Marie-Tooth disease type 2D and distal spinal muscular atrophy type V

Charcot-Marie-Tooth disease type 2D (CMT2D) and distal spinal muscular atrophy type V (dSMA-V) are axonal peripheral neuropathies inherited in an autosomal dominant fashion. Our previous genetic and physical mapping efforts localized the responsible gene(s) to a well-defined region on human chromosome 7p. Here, we report the identification of four disease-associated missense mutations in the glycyl tRNA synthetase gene in families with CMT2D and dSMA-V. This is the first example of an aminoacyl tRNA synthetase being implicated in a human genetic disease, which makes genes that encode these enzymes relevant candidates for other inherited neuropathies and motor neuron diseases.

Glycyl tRNA synthetase mutations in Charcot-Marie-Tooth disease type 2D and distal spinal muscular atrophy type V

Charcot-Marie-Tooth disease type 2D (CMT2D) and distal spinal muscular atrophy type V (dSMA-V) are axonal peripheral neuropathies inherited in an autosomal dominant fashion. Our previous genetic and physical mapping efforts localized the responsible gene(s) to a well-defined region on human chromosome 7p. Here, we report the identification of four disease-associated missense mutations in the glycyl tRNA synthetase gene in families with CMT2D and dSMA-V. This is the first example of an aminoacyl tRNA synthetase being implicated in a human genetic disease, which makes genes that encode these enzymes relevant candidates for other inherited neuropathies and motor neuron diseases.

Pericentromeric duplications in the laboratory mouse

Duplications have long been postulated to be an important mechanism by which genomes evolve. Interspecies genomic comparisons are one method by which the origin and molecular mechanism of duplications can be inferred. By comparative mapping in human, mouse, and rat, we previously found evidence for a recent chromosome-fission event that occurred in the mouse lineage. Cytogenetic mapping revealed that the genomic segments flanking the fission site appeared to be duplicated, with copies residing near the centromere of multiple mouse chromosomes. Here we report the mapping and sequencing of the regions of mouse chromosomes 5 and 6 involved in this chromosome-fission event as well as the results of comparative sequence analysis with the orthologous human and rat genomic regions. Our data indicate that the duplications associated with mouse chromosomes 5 and 6 are recent and that the resulting duplicated segments share significant sequence similarity with a series of regions near the centromeres of the mouse chromosomes previously identified by cytogenetic mapping. We also identified pericentromeric duplicated segments shared between mouse chromosomes 5 and 1. Finally, novel mouse satellite sequences as well as putative chimeric transcripts were found to be associated with the duplicated segments. Together, these findings demonstrate that pericentromeric duplications are not restricted to primates and may be a common mechanism for genome evolution in mammals.

Parallel construction of orthologous sequence-ready clone contig maps in multiple species

Comparison is a fundamental tool for analyzing DNA sequence. Interspecies sequence comparison is particularly powerful for inferring genome function and is based on the simple premise that conserved sequences are likely to be important. Thus, the comparison of a genomic sequence with its orthologous counterpart from another species is increasingly becoming an integral component of genome analysis. In ideal situations, such comparisons are performed with orthologous sequences from multiple species. To facilitate multispecies comparative sequence analysis, a robust and scalable strategy for simultaneously constructing sequence-ready bacterial artificial chromosome (BAC) contig maps from targeted genomic regions has been developed. Central to this approach is the generation and utilization of "universal" oligonucleotide-based hybridization probes ("overgo" probes), which are designed from sequences that are highly conserved between distantly related species. Large collections of these probes are used en masse to screen BAC libraries from multiple species in parallel, with the isolated clones assembled into physical contig maps. To validate the effectiveness of this strategy, efforts were focused on the construction of BAC-based physical maps from multiple mammalian species (chimpanzee, baboon, cat, dog, cow, and pig). Using available human and mouse genomic sequence and a newly developed computer program to design the requisite probes, sequence-ready maps were constructed in all species for a series of targeted regions totaling approximately 16 Mb in the human genome. The described approach can be used to facilitate the multispecies comparative sequencing of targeted genomic regions and can be adapted for constructing BAC contig maps in other vertebrates.

Comparative physical mapping of targeted regions of the rat genome

The comparative mapping and sequencing of vertebrate genomes is now a key priority for the Human Genome Project. In addition to finishing the human genome sequence and generating a 'working draft' of the mouse genome sequence, significant attention is rapidly turning to the analysis of other model organisms, such as the laboratory rat (Rattus norvegicus). As a complement to genome-wide mapping and sequencing efforts, it is often important to generate detailed maps and sequence data for specific regions of interest. Using an adaptation of our previously described approach for constructing mouse comparative and physical maps, we have established a general strategy for targeted mapping of the rat genome. Specifically, we constructed a framework comparative map of human Chromosome (Chr) 7 and the orthologous regions of the rat genome, as well as two large (>1-Mb) P1-derived artificial chromosome (PAC)-based physical maps. Generation of these physical maps involved the use of mouse-derived probes that cross-hybridized with rat PAC clones. The first PAC map encompasses the cystic fibrosis transmembrane conductance regulator gene (Cftr), while the second map allows a three-species comparison of a genomic region containing intra- and inter-chromosomal evolutionary rearrangements. The studies reported here further demonstrate that cross-species hybridization between related animals, such as rat and mouse, can be readily used for the targeted construction of clone-based physical maps, thereby accelerating the analysis of biologically interesting regions of vertebrate genomes.

Comparative genome mapping in the sequence-based era: early experience with human chromosome 7

The success of the ongoing Human Genome Project has resulted in accelerated plans for completing the human genome sequence and the earlier-than-anticipated initiation of efforts to sequence the mouse genome. As a complement to these efforts, we are utilizing the available human sequence to refine human-mouse comparative maps and to assemble sequence-ready mouse physical maps. Here we describe how the first glimpses of genomic sequence from human chromosome 7 are directly facilitating these activities. Specifically, we are actively enhancing the available human-mouse comparative map by analyzing human chromosome 7 sequence for the presence of orthologs of mapped mouse genes. Such orthologs can then be precisely positioned relative to mapped human STSs and other genes. The chromosome 7 sequence generated to date has allowed us to more than double the number of genes that can be placed on the comparative map. The latter effort reveals that human chromosome 7 is represented by at least 20 orthologous segments of DNA in the mouse genome. A second component of our program involves systematically analyzing the evolving human chromosome 7 sequence for the presence of matching mouse genes and expressed-sequence tags (ESTs). Mouse-specific hybridization probes are designed from such sequences and used to screen a mouse bacterial artificial chromosome (BAC) library, with the resulting data used to assemble BAC contigs based on probe-content data. Nascent contigs are then expanded using probes derived from newly generated BAC-end sequences. This approach produces BAC-based sequence-ready maps that are known to contain a gene(s) and are homologous to segments of the human genome for which sequence is already available. Our ongoing efforts have thus far resulted in the isolation and mapping of >3,800 mouse BACs, which have been assembled into >100 contigs. These contigs include >250 genes and represent approximately 40% of the mouse genome that is homologous to human chromosome 7. Together, these approaches illustrate how the availability of genomic sequence directly facilitates studies in comparative genomics and genome evolution.

Mutations in the gene encoding KRIT1, a Krev-1/rap1a binding protein, cause cerebral cavernous malformations (CCM1)

Cerebral cavernous malformations (CCM) are congenital vascular anomalies of the brain that can cause significant neurological disabilities, including intractable seizures and hemorrhagic stroke. One locus for autosomal dominant CCM ( CCM1 ) maps to chromosome 7q21-q22. Recombination events in linked family members define a critical region of approximately 2 Mb and a shared disease haplotype associated with a presumed founder effect in families of Mexican-American descent points to a potentially smaller region of interest. Using a genomic sequence-based positional cloning strategy, we have identified KRIT1, encoding a protein that interacts with the Krev-1/rap1a tumor suppressor, as the CCM1 gene. Seven different KRIT1 mutations have been identified in 23 distinct CCM1 families. The identical mutation is present in 16 of 21 Mexican-American families analyzed, substantiating a founder effect in this population. Other Mexican-American and non-Hispanic Caucasian CCM1 kindreds harbor other KRIT1 mutations. Identification of a common Mexican-American mutation has potential clinical significance for presymptomatic diagnosis of CCM in this population. In addition, these data point to a key role for the Krev-1/rap1a signaling pathway in angiogenesis and cerebrovascular disease.

Transplantation of transduced nonhuman primate CD34+ cells using a gibbon ape leukemia virus vector: restricted expression of the gibbon ape leukemia virus receptor to a subset of CD34+ cells

The transduction efficiencies of immunoselected rhesus macaque (Macaca mulatta) CD34+ cells and colony-forming progenitor cells based on polymerase chain reaction (PCR) analysis were comparable for an amphotropic Moloney murine leukemia virus (MLV) retroviral vector and a retroviral vector derived from the gibbon ape leukemia virus (GaLV) packaging cell line, PG13. On performing autologous transplantation studies using immunoselected CD34+ cells transduced with the GaLV envelope (env) retroviral vector, less than 1% of peripheral blood (PB) contained provirus. This was true whether bone marrow (BM) or cytokine-mobilized PB immunoselected CD34+ cells were reinfused. This level of marking was evident in two animals whose platelet counts never fell below 50,000/microliter and whose leukocyte counts had recovered by days 8 and 10 after having received 1.7 x 10(7) or greater of cytokine-mobilized CD34+ PB cells/kg. Reverse transcriptase(RT)-PCR analysis of CD34+ subsets for both the GaLV and amphotropic receptor were performed. The expression of the GaLV receptor was determined to be restricted to CD34+ Thy-1+ cells, and both CD34+ CD38+ and CD34+ CD38dim cells, while the amphotropic receptor was present on all CD34+ cell subsets examined. Our findings suggest that, in rhesus macaques, PG13-derived retroviral vectors may only be able to transduce a subset of CD34+ cells as only CD34+ Thy-1+ cells express the GaLV receptor.

Cloning of the BamHI methyl transferase gene from Bacillus amyloliquefaciens

We wish to report the initial characterization of a recombinant clone containing the BamHI methylase gene. Genomic chromosomal DNA purified from Bacillus amyloliquefaciens was partially cleaved with HindIII, fractionated by size, and cloned into pSP64. Plasmid DNA from this library was challenged with BamHI endonuclease and transformed into Escherichia coli HB101. A recombinant plasmid pBamM6.5 and a subclone pBamM2.5 were shown to contain the BamHI methylase gene based on three independent observations. Both plasmids were found to be resistant to BamHI endonuclease cleavage, and chromosomal DNA isolated from E.coli HB101 cells harboring either of the plasmids pBamM6.5 or pBamM2.5 was resistant to cleavage by BamHI endonuclease. In addition, DNA isolated from lambda phage passaged through E.coli HB101 containing either plasmid was also resistant to BamHI cleavage. Expression of the BamHI methylase gene is dependent on orientation in pSP64. In these clones preliminary evidence indicates that methylase gene expression may be under the direction of the plasmid encoded LacZ promoter.