Mobile genes in the human microbiome are structured from global to individual scales.

Recent work has underscored the importance of the microbiome in human health, and has largely attributed differences in phenotype to differences in the species present among individuals. However, mobile genes can confer profoundly different phenotypes on different strains of the same species. Little is known about the function and distribution of mobile genes in the human microbiome, and in particular whether the gene pool is globally homogenous or constrained by human population structure. Here, we investigate this question by comparing the mobile genes found in the microbiomes of 81 metropolitan North Americans with those of 172 agrarian Fiji islanders using a combination of single-cell genomics and metagenomics. We find large differences in mobile gene content between the Fijian and North American microbiomes, with functional variation that mirrors known dietary differences such as the excess of plant-based starch degradation genes found in Fijian individuals. Notably, we also observed differences between the mobile gene pools of neighbouring Fijian villages, even though microbiome composition across villages is similar. Finally, we observe high rates of recombination leading to individual-specific mobile elements, suggesting that the abundance of some genes may reflect environmental selection rather than dispersal limitation. Together, these data support the hypothesis that human activities and behaviours provide selective pressures that shape mobile gene pools, and that acquisition of mobile genes is important for colonizing specific human populations.

A phylogenetic analysis using full-length viral genomes of South American dengue serotype 3 in consecutive Venezuelan outbreaks reveals a novel NS5 mutation.

Dengue virus currently causes 50-100 million infections annually. Comprehensive knowledge about the evolution of Dengue in response to selection pressure is currently unavailable, but would greatly enhance vaccine design efforts. In the current study, we sequenced 187 new dengue virus serotype 3 (DENV-3) genotype III whole genomes isolated from Asia and the Americas. We analyzed them together with previously-sequenced isolates to gain a more detailed understanding of the evolutionary adaptations existing in this prevalent American serotype. In order to analyze the phylogenetic dynamics of DENV-3 during outbreak periods; we incorporated datasets of 48 and 11 sequences spanning two major outbreaks in Venezuela during 2001 and 2007-2008, respectively. Our phylogenetic analysis of newly sequenced viruses shows that subsets of genomes cluster primarily by geographic location, and secondarily by time of virus isolation. DENV-3 genotype III sequences from Asia are significantly divergent from those from the Americas due to their geographical separation and subsequent speciation. We measured amino acid variation for the E protein by calculating the Shannon entropy at each position between Asian and American genomes. We found a cluster of seven amino acid substitutions having high variability within E protein domain III, which has previously been implicated in serotype-specific neutralization escape mutants. No novel mutations were found in the E protein of sequences isolated during either Venezuelan outbreak. Shannon entropy analysis of the NS5 polymerase mature protein revealed that a G374E mutation, in a region that contributes to interferon resistance in other flaviviruses by interfering with JAK-STAT signaling was present in both the Asian and American sequences from the 2007-2008 Venezuelan outbreak, but was absent in the sequences from the 2001 Venezuelan outbreak. In addition to E, several NS5 amino acid changes were unique to the 2007-2008 epidemic in Venezuela and may give additional insight into the adaptive response of DENV-3 at the population level.

Genome variation in Cryptococcus gattii, an emerging pathogen of immunocompetent hosts.

Cryptococcus gattii recently emerged as the causative agent of cryptococcosis in healthy individuals in western North America, despite previous characterization of the fungus as a pathogen in tropical or subtropical regions. As a foundation to study the genetics of virulence in this pathogen, we sequenced the genomes of a strain (WM276) representing the predominant global molecular type (VGI) and a clinical strain (R265) of the major genotype (VGIIa) causing disease in North America. We compared these C. gattii genomes with each other and with the genomes of representative strains of the two varieties of Cryptococcus neoformans that generally cause disease in immunocompromised people. Our comparisons included chromosome alignments, analysis of gene content and gene family evolution, and comparative genome hybridization (CGH). These studies revealed that the genomes of the two representative C. gattii strains (genotypes VGI and VGIIa) are colinear for the majority of chromosomes, with some minor rearrangements. However, multiortholog phylogenetic analysis and an evaluation of gene/sequence conservation support the existence of speciation within the C. gattii complex. More extensive chromosome rearrangements were observed upon comparison of the C. gattii and the C. neoformans genomes. Finally, CGH revealed considerable variation in clinical and environmental isolates as well as changes in chromosome copy numbers in C. gattii isolates displaying fluconazole heteroresistance.

Human genome anatomy: BACs integrating the genetic and cytogenetic maps for bridging genome and biomedicine.

Human genome sequencing is accelerating rapidly. Multiple genome maps link this sequence to problems in biology and clinical medicine. Because each map represents a different aspect of the structure, content, and behavior of human chromosomes, these fundamental properties must be integrated with the genome to understand disease genes, cancer instability, and human evolution. Cytogenetic maps use 400-850 visible band landmarks and are the primary means for defining prenatal defects and novel cancer breakpoints, thereby providing simultaneous examination of the entire genome. Recent genetic, physical, and transcript maps use PCR-based landmarks called sequence-tagged sites (STSs). We have integrated these genome maps by anchoring the human cytogenetic to the STS-based genetic and physical maps with 1021 STS-BAC pairs at an average spacing of approximately 1 per 3 Mb. These integration points are represented by 872 unique STSs, including 642 polymorphic markers and 957 bacterial artificial chromosomes (BACs), each of which was localized on high resolution fluorescent banded chromosomes. These BACs constitute a resource that bridges map levels and provides the tools to seamlessly translate questions raised by genomic change seen at the chromosomal level into answers based at the molecular level. We show how the BACs provide molecular links for understanding human genomic duplications, meiosis, and evolution, as well as reagents for conducting genome-wide prenatal diagnosis at the molecular level and for detecting gene candidates associated with novel cancer breakpoints.

A novel member of the F-box/WD40 gene family, encoding dactylin, is disrupted in the mouse dactylaplasia mutant.

Early outgrowth of the vertebrate embryonic limb requires signalling by the apical ectodermal ridge (AER) to the progress zone (PZ), which in response proliferates and lays down the pattern of the presumptive limb in a proximal to distal progression. Signals from the PZ maintain the AER until the anlagen for the distal phalanges have been formed. The semidominant mouse mutant dactylaplasia (Dac) disrupts the maintenance of the AER, leading to truncation of distal structures of the developing footplate, or autopod. Adult Dac homozygotes thus lack hands and feet except for malformed single digits, whereas heterozygotes lack phalanges of the three middle digits. Dac resembles the human autosomal dominant split hand/foot malformation (SHFM) diseases. One of these, SHFM3, maps to chromosome 10q24 (Refs 6,7), which is syntenic to the Dac region on chromosome 19, and may disrupt the orthologue of Dac. We report here the positional cloning of Dac and show that it belongs to the F-box/WD40 gene family, which encodes adapters that target specific proteins for destruction by presenting them to the ubiquitination machinery. In conjuction with recent biochemical studies, this report demonstrates the importance of this gene family in vertebrate embryonic development.

A YAC-based physical map of the mouse genome.

A physical map of the mouse genome is an essential tool for both positional cloning and genomic sequencing in this key model system for biomedical research. Indeed, the construction of a mouse physical map with markers spaced at an average interval of 300 kb is one of the stated goals of the Human Genome Project. Here we report the results of a project at the Whitehead Institute/MIT Center for Genome Research to construct such a physical map of the mouse. We built the map by screening sequenced-tagged sites (STSs) against a large-insert yeast artificial chromosome (YAC) library and then integrating the STS-content information with a dense genetic map. The integrated map shows the location of 9,787 loci, providing landmarks with an average spacing of approximately 300 kb and affording YAC coverage of approximately 92% of the mouse genome. We also report the results of a project at the MRC UK Mouse Genome Centre targeted at chromosome X. The project produced a YAC-based map containing 619 loci (with 121 loci in common with the Whitehead map and 498 additional loci), providing especially dense coverage of this sex chromosome. The YAC-based physical map directly facilitates positional cloning of mouse mutations by providing ready access to most of the genome. More generally, use of this map in addition to a newly constructed radiation hybrid (RH) map provides a comprehensive framework for mouse genomic studies.

Mouse molecular cytogenetic resource: 157 BACs link the chromosomal and genetic maps.

We have established a collection of strong molecular cytogenetic markers that span the mouse autosomes and X chromosome at an average spacing of one per 19 Mb and identify 127 distinct band landmarks. In addition, this Mouse Molecular Cytogenetic Resource relates the ends of the genetic maps to their chromosomal locations. The resource consists of 157 bacterial artificial chromosome (BAC) clones, each of which identifies specific mouse chromosome bands or band borders, and 42 of which are linked to genetic markers that define the centromeric and telomeric ends of the Whitehead/MIT recombinational maps. In addition, 108 randomly selected and 6 STS-linked BACs have been assigned to single chromosome bands. We have also developed a high-resolution fluorescent reverse-banding technique for mouse chromosomes that allows simultaneous localization of probes by fluorescence in situ hybridization (FISH) with respect to the cytogenetic landmarks. This approach integrates studies of the entire mouse genome. Moreover, these reagents will simplify gene mapping and analyses of genomic fragments in fetal and adult mouse models. As shown with the MMU16 telomeric marker for the trisomy 16 mouse model of Down syndrome, these clones can obviate the need for metaphase analyses. The potential contribution of this resource and associated methods extends well beyond mapping and includes clues to understanding mouse chromosomes and their rearrangements in cancers and evolution. Finally it will facilitate the development of an integrated view of the mouse genome by providing anchor points from the genetic to the cytogenetic and functional maps of the mouse as we attempt to understand mutations, their biological consequences, and gene function.

The mouse pudgy mutation disrupts Delta homologue Dll3 and initiation of early somite boundaries.

Pudgy (pu) homozygous mice exhibit clear patterning defects at the earliest stages of somitogenesis, resulting in adult mice with severe vertebral and rib deformities. By positional cloning and complementation, we have determined that the pu phenotype is caused by a mutation in the delta-like 3 gene (Dll3), which is homologous to the Notch-ligand Delta in Drosophila. Histological and molecular marker analyses show that the pu mutation disrupts the proper formation of morphological borders in early somite formation and of rostral-caudal compartment boundaries within somites. Viability analysis also indicates an important role in early development. The results point to a key role for a Notch-signalling pathway in the initiation of patterning of vertebrate paraxial mesoderm.

Serrate2 is disrupted in the mouse limb-development mutant syndactylism.

The mouse syndactylism (sm) mutation impairs some of the earliest aspects of limb development and leads to subsequent abnormalities in digit formation. In sm homozygotes, the apical ectodermal ridge (AER) is hyperplastic by embryonic day 10.5, leading to abnormal dorsoventral thickening of the limb bud, subsequent merging of the skeletal condensations that give rise to cartilage and bone in the digits, and eventual fusion of digits. The AER hyperplasia and its effect on early digital patterning distinguish sm from many other syndactylies that result from later failure of cell death in the interdigital areas. Here we use positional cloning to show that the gene mutated in sm mice encodes the putative Notch ligand Serrate. The results provide direct evidence that a Notch signalling pathway is involved in the earliest stages of limb-bud patterning and support the idea that an ancient genetic mechanism underlies both AER formation in vertebrates and wing-margin formation in flies. In addition to cloning the sm gene, we have mapped three modifiers of sm, for which we suggest possible candidate genes.

A bacterial artificial chromosome-based framework contig map of human chromosome 22q.

We have constructed a physical map of human chromosome 22q using bacterial artificial chromosome (BAC) clones. The map consists of 613 chromosome 22-specific BAC clones that have been localized and assembled into contigs using 452 landmarks, 346 of which were previously ordered and mapped to specific regions of the q arm of the chromosome by means of chromosome 22-specific yeast artificial chromosome clones. The BAC-based map provides immediate access to clones that are stable and convenient for direct genome analysis. The approach to rapidly developing marker-specific BAC contigs is relatively straightforward and can be extended to generate scaffold BAC contig maps of the rest of the chromosomes. These contigs will provide substrates for sequencing the entire human genome. We discuss how to efficiently close contig gaps using the end sequences of BAC clone inserts.

Construction and characterization of a human bacterial artificial chromosome library.

We have constructed an arrayed human genomic BAC library with approximately 4x coverage that is represented by 96,000 BAC clones with average insert size of nearly 140 kb. A new BAC vector that allows color-based positive screening to identify transformants with inserts has increased BAC cloning efficiency. The library was gridded onto hybridization filters at high density for efficient identification of BAC clones by colony hybridization. The library was also formulated into characteristic DNA pools to allow for PCR screening of the library for STS content. We have characterized the library mainly by screening with more than 300 different landmarks that include cDNA, STSs, and cosmid clones. We describe methods for using BAC clones and discuss the implications for genome characterization, mapping, and sequencing.

A human chromosome 22 fosmid resource: mapping and analysis of 96 clones.

We have created a resource for chromosome 22 consisting of 96 unique, well-characterized Fosmids. The Fosmid vector permits efficient cloning of DNA fragments averaging 40 kb in a single-copy vector based on the F factor of Escherichia coli. We have found that Fosmid clones from human chromosome 22 show remarkable stability and are useful for a wide variety of applications in genome analysis. These 96 clones have been localized by FISH, using high-resolution fluorescent banding and multicolor mapping techniques, and their position on the chromosome was correlated with their content of a number of common repeated sequence elements. We identified a subset of clones likely to contain genes by restriction analysis using the enzymes NotI, MluI, SacII, and BssHII. This collection of cytogenetically anchored clones, representing nearly 7% of the chromosome, is of immediate value for detecting chromosomal rearrangements, for use in gene isolation, and as a framework for physical mapping.

Disruption of the nuclear hormone receptor RORalpha in staggerer mice.

Homozygous staggerer (sg) mice show a characteristic severe cerebellar ataxia due to a cell-autonomous defect in the development of Purkinje cells. These cells show immature morphology, synaptic arrangement, biochemical properties and gene expression, and are reduced in numbers. In addition, sg heterozygotes show accelerated dendritic atrophy and cell loss, suggesting that sg has a role in mature Purkinje cells. Effects of this mutation on cerebellar development have been studied for 25 years, but its molecular basis has remained unknown. We have genetically mapped staggerer to an interval of 160 kilobases on mouse chromosome 9 which was found to contain the gene encoding RORalpha, a member of the nuclear hormone-receptor superfamily. Staggerer mice were found to carry a deletion within the RORalpha gene that prevents translation of the ligand-binding homology domain. We propose a model based on these results, in which RORalpha interacts with the thyroid hormone signalling pathway to induce Purkinje-cell maturation.

An STS-based map of the human genome.

A physical map has been constructed of the human genome containing 15,086 sequence-tagged sites (STSs), with an average spacing of 199 kilobases. The project involved assembly of a radiation hybrid map of the human genome containing 6193 loci and incorporated a genetic linkage map of the human genome containing 5264 loci. This information was combined with the results of STS-content screening of 10,850 loci against a yeast artificial chromosome library to produce an integrated map, anchored by the radiation hybrid and genetic maps. The map provides radiation hybrid coverage of 99 percent and physical coverage of 94 percent of the human genome. The map also represents an early step in an international project to generate a transcript map of the human genome, with more than 3235 expressed sequences localized. The STSs in the map provide a scaffold for initiating large-scale sequencing of the human genome.

Use of secondary pulsed field gel electrophoresis in separation of large DNA.

It has been reported that secondary pulsed field gel (SPFG) electrophoresis can dramatically increase the speed of separation of large DNA molecules without a decrease in resolution (Zhang, T. Y., Smith, C. L., and Cantor, C. R. (1991) Nucleic Acids Res. 19, 1291-1296). However, our attempts to duplicate previous SPFG conditions were unsuccessful. We therefore sought to more precisely define the effects of secondary pulsing on the separation of large DNA and to determine the value of the technique in separating molecules up to 1100 kb. Here we report on two of the key SPFG parameters, namely the frequency and duration of the secondary pulse on the migration and resolution of DNA in SPFG. We found that the size range of separation is determined by the sum of the duration of the primary and secondary pulses. Under optimal conditions, a 25-70% increase in velocity can be achieved without loss in resolution.

The basis of high resolution separation of small DNAs by asymmetric-voltage field inversion electrophoresis and its application to DNA sequencing gels.

We have previously shown that asymmetric-voltage field inversion electrophoresis produces more uniform separation for fragments between 1 and 50 kilobases (kb) than other modes of pulsed field gel electrophoresis. We now report on the basis of this phenomenon. As in conventional electrophoresis, the pulsed field mobility of DNAs between 1 and 50 kb varies with voltage in a size dependent manner. The complex migration pattern obtained with asymmetric-voltage field inversion electrophoresis reflects the difference between the mobilities of each sized fragment under the conditions used for the forward and reverse fields. We have applied this technique to DNA sequencing gels and find improvement in resolution for single-stranded fragments in polyacrylamide gels.

Pulsed field gel electrophoresis: studies of DNA migration made with the programmable, autonomously-controlled electrode electrophoresis system.

We have studied the migration of DNA in pulsed field agarose gels under a variety of electrophoresis conditions. We have made use of an instrument which can generate electric fields of any orientation, magnitude, or duration to compare different separation techniques for DNA molecules of from 1 to several thousand kilobase pairs. We discuss the capabilities of the system and present results of gel runs in which electrophoresis conditions were changed individually or in combination. The mobility of DNA in pulsed field gels is shown to reflect a number of interdependent physical parameters.