In addition to protein coding sequence, the human genome contains a significant amount of regulatory DNA, the identification of which is proving somewhat recalcitrant to both in silico and functional methods. An approach that has been used with some success is comparative sequence analysis, whereby equivalent genomic regions from different organisms are compared in order to identify both similarities and differences. In general, similarities in sequence between highly divergent organisms imply functional constraint. We have used a whole-genome comparison between humans and the pufferfish, Fugu rubripes, to identify nearly 1,400 highly conserved non-coding sequences. Given the evolutionary divergence between these species, it is likely that these sequences are found in, and furthermore are essential to, all vertebrates. Most, and possibly all, of these sequences are located in and around genes that act as developmental regulators. Some of these sequences are over 90% identical across more than 500 bases, being more highly conserved than coding sequence between these two species. Despite this, we cannot find any similar sequences in invertebrate genomes. In order to begin to functionally test this set of sequences, we have used a rapid in vivo assay system using zebrafish embryos that allows tissue-specific enhancer activity to be identified. Functional data is presented for highly conserved non-coding sequences associated with four unrelated developmental regulators (SOX21, PAX6, HLXB9, and SHH), in order to demonstrate the suitability of this screen to a wide range of genes and expression patterns. Of 25 sequence elements tested around these four genes, 23 show significant enhancer activity in one or more tissues. We have identified a set of non-coding sequences that are highly conserved throughout vertebrates. They are found in clusters across the human genome, principally around genes that are implicated in the regulation of development, including many transcription factors. These highly conserved non-coding sequences are likely to form part of the genomic circuitry that uniquely defines vertebrate development.

Successful host defense against numerous pulmonary infections depends on bacterial clearance by polymorphonuclear leukocytes (PMNs); however, excessive PMN accumulation can result in life-threatening lung injury. Local expression of CXC chemokines is critical for PMN recruitment. The impact of chemokinedependent PMN recruitment during pulmonary Mycobacterium tuberculosis infection is not fully understood. Here, we analyzed expression of genes encoding CXC chemokines in M. tuberculosis-infected murine lung tissue and found that M. tuberculosis infection promotes upregulation of Cxcr2 and its ligand Cxcl5. To determine the contribution of CXCL5 in pulmonary PMN recruitment, we generated Cxcl5 -/-mice and analyzed their immune response against M. tuberculosis. Both Cxcr2 -/-mice and Cxcl5 -/-mice, which are deficient for only one of numerous CXCR2 ligands, exhibited enhanced survival compared with that of WT mice following high-dose M. tuberculosis infection. The resistance of Cxcl5 -/-mice to M. tuberculosis infection was not due to heightened M. tuberculosis clearance but was the result of impaired PMN recruitment, which reduced pulmonary inflammation. Lung epithelial cells were the main source of CXCL5 upon M. tuberculosis infection, and secretion of CXCL5 was reduced by blocking TLR2 signaling. Together, our data indicate that TLR2-induced epithelial-derived CXCL5 is critical for PMN-driven destructive inflammation in pulmonary tuberculosis.

Determining the underlying haplotypes of individual human genomes is an essential, but currently difficult, step toward a complete understanding of genome function. Fosmid pool-based next-generation sequencing allows genome-wide generation of 40-kb haploid DNA segments, which can be phased into contiguous molecular haplotypes computationally by Single Individual Haplotyping (SIH). Many SIH algorithms have been proposed, but the accuracy of such methods has been difficult to assess due to the lack of real benchmark data. To address this problem, we generated whole genome fosmid sequence data from a HapMap trio child, NA12878, for which reliable haplotypes have already been produced. We assembled haplotypes using eight algorithms for SIH and carried out direct comparisons of their accuracy, completeness and efficiency. Our comparisons indicate that fosmid-based haplotyping can deliver highly accurate results even at low coverage and that our SIH algorithm, ReFHap, is able to efficiently produce high-quality haplotypes. We expanded the haplotypes for NA12878 by combining the current haplotypes with our fosmid-based haplotypes, producing near-to-complete new gold-standard haplotypes containing almost 98% of heterozygous SNPs. This improvement includes notable fractions of disease-related and GWA SNPs. Integrated with other molecular biological data sets, this phase information will advance the emerging field of diploid genomics.

The nutrients consumed by intracellular pathogens are mostly unknown. This is mainly due to the challenge of disentangling host and pathogen metabolism sharing the majority of metabolic pathways and hence metabolites. Here, we investigated the metabolic changes of Mycobacterium tuberculosis, the causative agent of tuberculosis, and its human host cell during early infection. To this aim, we combined gene expression data of both organisms and metabolite changes during the course of infection through integration into a genome-wide metabolic network. This led to the identification of infection-specific metabolic alterations, which we further exploited to model host-pathogen interactions quantitatively by flux balance analysis. These in silico data suggested that tubercle bacilli consume up to 33 different nutrients during early macrophage infection, which the bacteria utilize to generate energy and biomass to establish intracellular growth. Such multisubstrate fueling strategy renders the pathogen’s metabolism robust toward perturbations, such as innate immune responses or antibiotic treatments.