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 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.

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.

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.

A novel instrument for separating large DNA molecules with pulsed homogeneous electric fields.

A new instrument has been developed for the electrophoretic separation of large DNA molecules that can independently regulate the voltage of each of 24 electrodes and allow the magnitude, orientation, homogeneity, and duration of the electric field to be precisely controlled. Each parameter can be varied at any time during the electrophoretic process. Thus distinct sets of conditions can be combined to optimize the separation of various fragment sizes in a single run. Independent control of electrode voltage allows all of the fields to be generated with electrodes arranged in a closed contour, independent of a particular geometry. This device increases both the resolution in any size range and the speed of separation, especially for DNA molecules larger than 3 megabases.

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.

Butyrate-induced changes in nuclease sensitivity of chromatin cannot be correlated with transcriptional activation.

We examined in the H4IIE rat hepatoma cell line the relationship between butyrate-induced changes in the nuclease sensitivity of chromatin and changes in transcriptional activity of specific genes. The butyrate-inducible metallothionein I (MT-I) gene underwent a dramatic increase in DNase I sensitivity after 3 h of butyrate treatment. However, genes not transcribed in H4IIE cells underwent the same changes in DNase I sensitivity. Thus, butyrate-induced increases in DNase I sensitivity are not sufficient for the transcriptional activation of a gene. Butyrate treatment has also been reported to alter the sensitivity of sequences to micrococcal nuclease (MNase) in a manner reflecting their tissue-specific expression. Butyrate exposure caused increased digestion of the MT-I gene by MNase. However, butyrate-induced MNase sensitivity also occurred for genes which are neither transcribed in untreated cells nor butyrate inducible. Moreover, cadmium, a potent transcriptional activator of the MT-I gene, does not alter the sensitivity of the MT-I gene to MNase. Thus, the butyrate-induced alterations in MNase sensitivity are neither sufficient for, necessary for, nor indicative of transcriptional activation.

Rat metallothionein-1 structural gene and three pseudogenes, one of which contains 5'-regulatory sequences.

As shown by Southern blot analysis, the metallothionein-1 (MT-1) genes in rats comprise a multigene family. We present the sequence of the MT-1 structural gene and compare its features with other metallothionein genes. Three MT-1 pseudogenes which we sequenced apparently arose by reverse transcription of processed mRNA transcripts. Two of these, MT-1 psi a and MT-1 psi c, are retrogenes which derive from the MT-1 mRNA, having diverged from the MT-1 gene 6.9 and 2.6 million years ago, respectively. The third, MT-1 psi b, differs from the MT-1 cDNA by only three nucleotide alterations. Surprisingly, MT-1 psi b also preserves sequence homology for 142 base pairs 5' to the transcription initiation site of the parent gene; it contains a promoter sequence sufficient for specifying metal ion induction. We identified, by S1 nuclease mapping, an RNA polymerase II initiation site 432 base pairs 5' of the MT-1 transcription initiation site of the MT-1 structural gene which could explain the formation of the mRNA precursor to this pseudogene. We were unable to detect MT-1 psi b transcripts, either in liver tissue or after transfection. We conclude that the absence of detectable transcripts from this pseudogene is due to either a reduced level of transcription or the formation of unstable transcripts as a consequence of the lack of a consensus sequence normally found 3' of transcription termination in the MT-1 structural gene.

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.

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.

Regulation of the rat metallothionein-I gene by sodium butyrate.

Sodium butyrate selectively induces accumulation of metallothionein-I (MT-I) RNA in H4IIE rat hepatoma cells. The induction is rapid; significant elevation in cytoplasmic MT-I RNA can be observed within three hours after exposure to 5 mM butyrate. Maximal levels of MT-I RNA are obtained after eight hours. Butyrate stimulates MT RNA accumulation in the absence of de novo protein synthesis, indicating that MT induction by butyrate is not a distal step in a cascade of gene activation events. Butyrate blocks the induction of tyrosine amino transferase by dexamethasone. In contrast, butyrate and dexamethasone induced MT RNA elevations are additive. Butyrate induced MT-I RNA transcripts initiate at the correct start site. Measurements of the transcriptional activity of the MT-I gene indicate that butyrate stimulates MT-I transcription. The rapid, direct nature of the induction of MT-I by butyrate, combined with the extensive characterization of the metallothionein gene, provide an excellent system in which to study the effects of butyrate on a small, well-defined, responsive region of chromatin.

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.

Pulsed field gel electrophoresis techniques for separating 1- to 50-kilobase DNA fragments.

Conventional agarose gel electrophoresis separates DNA using a static electric field. The maximum size limit for separation of DNA by this method is about 20 kilobase pairs (kb). A number of new electrophoretic techniques which employ periodic reorientation of electric fields permit separation of DNA well beyond this size limit. We sought to determine whether the use of very fast (millisecond) field switching could improve separation of DNA in the size range of 1 to 50 kb. Additionally, we have compared the resolution obtained with each of the different field switching regimens for DNA in this size range. Switching intervals of from 0.2 to 900 ms were used with unidirectional pulsing of a single electric field, with pulsed field gels, and with field inversion gel electrophoresis. Plotting the mobility of DNA as a function of size demonstrates that under the conditions used, each of these techniques offers comparable resolution. We also have examined the separation obtained when field inversion gels are run with forward and reverse fields of equal voltage and different durations, versus using fields of equal duration and different voltages. Field inversion which uses forward and reverse fields of different voltages yields resolution which is superior to the other methods examined.

Distinct forms of the beta subunit of GTP-binding regulatory proteins identified by molecular cloning.

Two distinct beta subunits of guanine nucleotide-binding regulatory proteins have been identified by cDNA cloning and are referred to as beta 1 and beta 2 subunits. The bovine transducin beta subunit (beta 1) has been cloned previously. We have now isolated and analyzed cDNA clones that encode the beta 2 subunit from bovine adrenal, bovine brain, and a human myeloid leukemia cell line, HL-60. The 340-residue Mr 37,329 beta 2 protein is 90% identical with beta 1 in predicted amino acid sequence, and it is also organized as a series of repetitive homologous segments. The major mRNA that encodes the bovine beta 2 subunit is 1.7 kilobases in length. It is expressed at lower levels than beta 1 subunit mRNA in all tissues examined. The beta 1 and beta 2 messages are expressed in cloned human cell lines. Hybridization of cDNA probes to bovine DNA showed that beta 1 and beta 2 are encoded by separate genes. The amino acid sequences for the bovine and human beta 2 subunit are identical, as are the amino acid sequences for the bovine and human beta 1 subunit. This evolutionary conservation suggests that the two beta subunits have different roles in the signal transduction process.

Metallothionein synthesis in foetal, neonatal and maternal rat liver.

The synthesis of hepatic metallothionein relative to other cytosol proteins was measured by [35S]cysteine incorporation in foetal, neonatal and pregnant rats. The relative rate of hepatic metallothionein synthesis reached a maximum in foetal liver on days 18-21 of gestation. Metallothionein synthesis then declined until weaning, when adult levels were established. The rate of metallothionein synthesis was greater in pregnant rats at term than in nulliparous rats. To determine if circulating inducing agents could play a role in the regulation of metallothionein synthesis in foetal liver we treated pregnant rats with inducers at a time prior to the normal rise in foetal liver metallothionein synthesis. Injections of copper, cadmium or hydrocortisone to 17-day-pregnant dams failed to induce foetal metallothionein synthesis. In contrast, zinc injection to the dam was an effective inducer in the foetuses. Maternal laparotomy (performed to expose the foetus for direct injection of inducers) induced foetal metallothionein synthesis. Metallothionein synthesis in the livers of 17-day-gestation dams was induced by all metal injections and laparotomy but, surprisingly, not by hydrocortisone injection. Maternal adrenalectomy did not influence the subsequent normal elevation in foetal or maternal metallothionein synthesis. These results, in conjunction with previous reports, suggest that mobilization of zinc in serum during late gestation may regulate foetal and maternal changes in metallothionein synthesis.