Genomic characterization of Salmonella Cerro ST367, an emerging Salmonella subtype in cattle in the United States.

BACKGROUND: Within the last decade, Salmonella enterica subsp. enterica serovar Cerro (S. Cerro) has become one of the most common serovars isolated from cattle and dairy farm environments in the northeastern US. The fact that this serovar is commonly isolated from subclinically infected cattle and is rarely associated with human disease, despite its frequent isolation from cattle, has led to the hypothesis that this emerging serovar may be characterized by reduced virulence. We applied comparative and population genomic approaches to (i) characterize the evolution of this recently emerged serovar and to (ii) gain a better understanding of genomic features that could explain some of the unique epidemiological features associated with this serovar. RESULTS: In addition to generating a de novo draft genome for one Salmonella Cerro strain, we also generated whole genome sequence data for 26 additional S. Cerro isolates, including 16 from cattle operations in New York (NY) state, 2 from human clinical cases from NY in 2008, and 8 from diverse animal sources (7 from Washington state and 1 from Florida). All isolates sequenced in this study represent sequence type ST367. Population genomic analysis showed that isolates from the NY cattle operations form a well-supported clade within S. Cerro ST367 (designated here "NY bovine clade"), distinct from isolates from Washington state, Florida and the human clinical cases. A molecular clock analysis indicates that the most recent common ancestor of the NY bovine clade dates back to 1998, supporting the recent emergence of this clone.Comparative genomic analyses revealed several relevant genomic features of S. Cerro ST367, that may be responsible for reduced virulence of S. Cerro, including an insertion creating a premature stop codon in sopA. In addition, patterns of gene deletion in S. Cerro ST367 further support adaptation of this clone to a unique ecological or host related niche. CONCLUSIONS: Our results indicate that the increase in prevalence of S. Cerro ST367 is caused by a highly clonal subpopulation and that S. Cerro ST367 is characterized by unique genomic deletions that may indicate adaptation to specific ecological niches and possibly reduced virulence in some hosts.

Mutability of Y-chromosomal microsatellites: rates, characteristics, molecular bases, and forensic implications.

Nonrecombining Y-chromosomal microsatellites (Y-STRs) are widely used to infer population histories, discover genealogical relationships, and identify males for criminal justice purposes. Although a key requirement for their application is reliable mutability knowledge, empirical data are only available for a small number of Y-STRs thus far. To rectify this, we analyzed a large number of 186 Y-STR markers in nearly 2000 DNA-confirmed father-son pairs, covering an overall number of 352,999 meiotic transfers. Following confirmation by DNA sequence analysis, the retrieved mutation data were modeled via a Bayesian approach, resulting in mutation rates from 3.78 × 10(-4) (95% credible interval [CI], 1.38 × 10(-5) - 2.02 × 10(-3)) to 7.44 × 10(-2) (95% CI, 6.51 × 10(-2) - 9.09 × 10(-2)) per marker per generation. With the 924 mutations at 120 Y-STR markers, a nonsignificant excess of repeat losses versus gains (1.16:1), as well as a strong and significant excess of single-repeat versus multirepeat changes (25.23:1), was observed. Although the total repeat number influenced Y-STR locus mutability most strongly, repeat complexity, the length in base pairs of the repeated motif, and the father's age also contributed to Y-STR mutability. To exemplify how to practically utilize this knowledge, we analyzed the 13 most mutable Y-STRs in an independent sample set and empirically proved their suitability for distinguishing close and distantly related males. This finding is expected to revolutionize Y-chromosomal applications in forensic biology, from previous male lineage differentiation toward future male individual identification.

Escherichia coli serotype O55:H7 diversity supports parallel acquisition of bacteriophage at Shiga toxin phage insertion sites during evolution of the O157:H7 lineage.

Enteropathogenic Escherichia coli (EPEC) continues to be a leading cause of mortality and morbidity in children around the world. Two EPEC genomes have been fully sequenced: those of EPEC O127:H6 strain E2348/69 (United Kingdom, 1969) and EPEC O55:H7 strain CB9615 (Germany, 2003). The O55:H7 serotype is a recent precursor to the virulent enterohemorrhagic E. coli O157:H7. To explore the diversity of O55:H7 and better understand the clonal evolution of O157:H7, we fully sequenced EPEC O55:H7 strain RM12579 (California, 1974), which was collected 1 year before the first U.S. isolate of O157:H7 was identified in California. Phage-related sequences accounted for nearly all differences between the two O55:H7 strains. Additionally, O55:H7 and O157:H7 strains were tested for the presence and insertion sites of Shiga toxin gene (stx)-containing bacteriophages. Analysis of non-phage-associated genes supported core elements of previous O157:H7 stepwise evolutionary models, whereas phage composition and insertion analyses suggested a key refinement. Specifically, the placement and presence of lambda-like bacteriophages (including those containing stx) should not be considered stable evolutionary markers or be required in placing O55:H7 and O157:H7 strains within the stepwise evolutionary models. Additionally, we suggest that a 10.9-kb region (block 172) previously believed unique to O55:H7 strains can be used to identify early O157:H7 strains. Finally, we defined two subsets of O55:H7 strains that share an as-yet-unobserved or extinct common ancestor with O157:H7 strains. Exploration of O55:H7 diversity improved our understanding of the evolution of E. coli O157:H7 and suggested a key revision to accommodate existing and future configurations of stx-containing bacteriophages into current models.

PDX1 regulation of FABP1 and novel target genes in human intestinal epithelial Caco-2 cells.

The transcription factor pancreatic and duodenal homeobox 1 (PDX1) plays an essential role in pancreatic development and in maintaining proper islet function via target gene regulation. Few intestinal PDX1 targets, however, have been described. We sought to define novel PDX1-regulated intestinal genes. Caco-2 human intestinal epithelial cells were engineered to overexpress PDX1 and gene expression profiles relative to control cells were assessed. Expression of 80 genes significantly increased while that of 49 genes significantly decreased more than 4-fold following PDX1 overexpression in differentiated Caco-2 cells. Analysis of the differentially regulated genes with known functional annotations revealed genes encoding transcription factors, growth factors, kinases, digestive glycosidases, nutrient transporters, nutrient binding proteins, and structural components. The gene for fatty acid binding protein 1, liver, FABP1, is repressed by PDX1 in Caco-2 cells. PDX1 overexpression in Caco-2 cells also results in repression of promoter activity driven by the 0.6kb FABP1 promoter. PDX1 regulation of promoter activity is consistent with the decrease in FABP1 RNA abundance resulting from PDX1 overexpression and identifies FABP1 as a candidate PDX1 target. PDX1 repression of FABP1, LCT, and SI suggests a role for PDX1 in patterning anterior intestinal development.

A whole-genome single nucleotide polymorphism-based approach to trace and identify outbreaks linked to a common Salmonella enterica subsp. enterica serovar Montevideo pulsed-field gel electrophoresis type.

In this study, we report a whole-genome single nucleotide polymorphism (SNP)-based evolutionary approach to study the epidemiology of a multistate outbreak of Salmonella enterica subsp. enterica serovar Montevideo. This outbreak included 272 cases that occurred in 44 states between July 2009 and April 2010. A case-control study linked the consumption of salami made with contaminated black and red pepper to the outbreak. We sequenced, on the SOLiD System, 47 isolates with XbaI PFGE pattern JIXX01.0011, a common pulsed-field gel electrophoresis (PFGE) pattern associated with isolates from the outbreak. These isolates represented 20 isolates collected from human sources during the period of the outbreak and 27 control isolates collected from human, food, animal, and environmental sources before the outbreak. Based on 253 high-confidence SNPs, we were able to reconstruct a tip-dated molecular clock phylogeny of the isolates and to assign four human isolates to the actual outbreak. We developed an SNP typing assay to rapidly discriminate between outbreak-related cases and non-outbreak-related cases and tested this assay on an extended panel of 112 isolates. These results suggest that only a very small percentage of the human isolates with the outbreak PFGE pattern and obtained during the outbreak period could be attributed to the actual pepper-related outbreak (20%), while the majority (80%) of the putative cases represented background cases. This study demonstrates that next-generation-based SNP typing provides the resolution and accuracy needed for outbreak investigations of food-borne pathogens that cannot be distinguished by currently used subtyping methods.

Single nucleotide polymorphisms and haplotypes in Native American populations.

Autosomal DNA polymorphisms can provide new information and understanding of both the origins of and relationships among modern Native American populations. At the same time that autosomal markers can be highly informative, they are also susceptible to ascertainment biases in the selection of the markers to use. Identifying markers that can be used for ancestry inference among Native American populations can be considered separate from identifying markers to further the quest for history. In the current study, we are using data on nine Native American populations to compare the results based on a large haplotype-based dataset with relatively small independent sets of single nucleotide polymorphisms. We are interested in what types of limited datasets an individual laboratory might be able to collect are best for addressing two different questions of interest. First, how well can we differentiate the Native American populations and/or infer ancestry by assigning an individual to her population(s) of origin? Second, how well can we infer the historical/evolutionary relationships among Native American populations and their Eurasian origins? We conclude that only a large comprehensive dataset involving multiple autosomal markers on multiple populations will be able to answer both questions; different small sets of markers are able to answer only one or the other of these questions. Using our largest dataset, we see a general increasing distance from Old World populations from North to South in the New World except for an unexplained close relationship between our Maya and Quechua samples.

SNPs for a universal individual identification panel.

An efficient method to uniquely identify every individual would have value in quality control and sample tracking of large collections of cell lines or DNA as is now often the case with whole genome association studies. Such a method would also be useful in forensics. SNPs represent the best markers for such purposes. We have developed a globally applicable resource of 92 SNPs for individual identification (IISNPs) with extremely low probabilities of any two unrelated individuals from anywhere in the world having identical genotypes. The SNPs were identified by screening over 500 likely/candidate SNPs on samples of 44 populations representing the major regions of the world. All 92 IISNPs have an average heterozygosity [0.4 and the F(st) values are all\0.06 on our 44 populations making these a universally applicable panel irrespective of ethnicity or ancestry. No significant linkage disequilibrium (LD) occurs for all unique pairings of 86 of the 92 IISNPs (median LD = 0.011) in all of the 44 populations. The remaining 6 IISNPs show strong LD in most of the 44 populations for a small subset (7) of the unique pairings in which they occur due to close linkage. 45 of the 86 SNPs are spread across the 22 human autosomes and show very loose or no genetic linkage with each other. These 45 IISNPs constitute an excellent panel for individual identification including paternity testing with associated probabilities of individual genotypes less than 10(-15), smaller than achieved with the current panels of forensic markers. This panel also improves on an interim panel of 40 IISNPs previously identified using 40 population samples. The unlinked status of the subset of 45 SNPs we have identified also makes them useful for situations involving close biological relationships. Comparisons with random sets of SNPs illustrate the greater discriminating power, efficiency, and more universal applicability of this IISNP panel to populations around the world. The full set of 86 IISNPs that do not show LD can be used to provide even smaller genotype match probabilities in the range of 10(-31)-10(-35) based on the 44 population samples studied.

Pdx1 inactivation restricted to the intestinal epithelium in mice alters duodenal gene expression in enterocytes and enteroendocrine cells.

Null mutant mice lacking the transcription factor pancreatic and duodenal homeobox 1 (Pdx1) are apancreatic and survive only a few days after birth. The role of Pdx1 in regulating intestinal gene expression has therefore yet to be determined in viable mice with normal pancreatic development. We hypothesized that conditional inactivation of Pdx1 restricted to the intestinal epithelium would alter intestinal gene expression and cell differentiation. Pdx1(flox/flox);VilCre mice with intestine-specific Pdx1 inactivation were generated by crossing a transgenic mouse strain expressing Cre recombinase, driven by a mouse villin 1 gene promoter fragment, with a mutant mouse strain homozygous for loxP site-flanked Pdx1. Pdx1 protein is undetectable in all epithelial cells in the intestinal epithelium of Pdx1(flox/flox);VilCre mice. Goblet cell number and mRNA abundance for mucin 3 and mucin 13 genes in the proximal small intestine are comparable between Pdx1(flox/flox);VilCre and control mice. Similarly, Paneth cell number and expression of Paneth cell-related genes Defa1, Defcr-rs1, and Mmp7 in the proximal small intestine remain statistically unchanged by Pdx1 inactivation. Although the number of enteroendocrine cells expressing chromogranin A/B, gastric inhibitory polypeptide (Gip), or somatostatin (Sst) is unaffected in the Pdx1(flox/flox);VilCre mice, mRNA abundance for Gip and Sst is significantly reduced in the proximal small intestine. Conditional Pdx1 inactivation attenuates intestinal alkaline phosphatase (IAP) activity in the duodenal epithelium, consistent with an average 91% decrease in expression of the mouse enterocyte IAP gene, alkaline phosphatase 3 (a novel Pdx1 target candidate), in the proximal small intestine following Pdx1 inactivation. We conclude that Pdx1 is necessary for patterning appropriate gene expression in enterocytes and enteroendocrine cells of the proximal small intestine.

Comprehensive mutation analysis of 17 Y-chromosomal short tandem repeat polymorphisms included in the AmpFlSTR Yfiler PCR amplification kit.

The Y-chromosomal short tandem repeat (Y-STR) polymorphisms included in the AmpFlSTR Yfiler polymerase chain reaction amplification kit have become widely used for forensic and evolutionary applications where a reliable knowledge on mutation properties is necessary for correct data interpretation. Therefore, we investigated the 17 Yfiler Y-STRs in 1,730-1,764 DNA-confirmed father-son pairs per locus and found 84 sequence-confirmed mutations among the 29,792 meiotic transfers covered. Of the 84 mutations, 83 (98.8%) were single-repeat changes and one (1.2%) was a double-repeat change (ratio, 1:0.01), as well as 43 (51.2%) were repeat gains and 41 (48.8%) repeat losses (ratio, 1:0.95). Medians from Bayesian estimation of locus-specific mutation rates ranged from 0.0003 for DYS448 to 0.0074 for DYS458, with a median rate across all 17 Y-STRs of 0.0025. The mean age (at the time of son's birth) of fathers with mutations was with 34.40 (+/-11.63) years higher than that of fathers without ones at 30.32 (+/-10.22) years, a difference that is highly statistically significant (p < 0.001). A Poisson-based modeling revealed that the Y-STR mutation rate increased with increasing father's age on a statistically significant level (alpha = 0.0294, 2.5% quantile = 0.0001). From combining our data with those previously published, considering all together 135,212 meiotic events and 331 mutations, we conclude for the Yfiler Y-STRs that (1) none had a mutation rate of >1%, 12 had mutation rates of >0.1% and four of <0.1%, (2) single-repeat changes were strongly favored over multiple-repeat ones for all loci but 1 and (3) considerable variation existed among loci in the ratio of repeat gains versus losses. Our finding of three Y-STR mutations in one father-son pair (and two pairs with two mutations each) has consequences for determining the threshold of allelic differences to conclude exclusion constellations in future applications of Y-STRs in paternity testing and pedigree analyses.

Transcriptome characterization elucidates signaling networks that control human ES cell growth and differentiation.

Human embryonic stem (hES) cells hold promise for generating an unlimited supply of cells for replacement therapies. To characterize hES cells at the molecular level, we obtained 148,453 expressed sequence tags (ESTs) from undifferentiated hES cells and three differentiated derivative subpopulations. Over 32,000 different transcripts expressed in hES cells were identified, of which more than 16,000 do not match closely any gene in the UniGene public database. Queries to this EST database revealed 532 significantly upregulated and 140 significantly downregulated genes in undifferentiated hES cells. These data highlight changes in the transcriptional network that occur when hES cells differentiate. Among the differentially regulated genes are several components of signaling pathways and transcriptional regulators that likely play key roles in hES cell growth and differentiation. The genomic data presented here may facilitate the derivation of clinically useful cell types from hES cells.

Spatio-temporal patterns of intestine-specific transcription factor expression during postnatal mouse gut development.

The small intestine matures from a primitive tube into morphologically and functionally distinct regions during gut development. Maximal expression of the genes encoding the digestive enzymes lactase-phlorizin hydrolase and sucrase-isomaltase is spatially restricted to distinct segments along the anterior-posterior axis of the small intestine and is temporally regulated during postnatal maturation. Transcription factors capable of interacting with the intestinal lactase and sucrase gene promoters are candidate regulators of spatio-temporal patterning during gut development and maturation. We aimed to quantitatively examine and compare the relative expression levels of a set of intestine-specific transcription factors along the anterior-posterior gut axis during postnatal maturation. Our analysis was focused on the transcription factors capable of regulating the intestinal lactase and sucrase-isomaltase genes. A real-time PCR protocol was used to quantitatively examine and compare spatially and temporally the relative transcript abundance levels for intestine-specific factors during postnatal intestinal maturation. Distinct spatial expressions patterns were detected along the length of the small intestine for PDX-1, Cdx-2, GATA-4, GATA-5, GATA-6, HNF-1alpha, HNF-1beta and CDP transcription factor genes. There is a general decline in transcript abundance for the factor genes during postnatal maturation. Defining the spatio-temporal expression patterns for intestine-specific transcription factor genes contributes to investigation of the roles that factor gradients play in mediating gut development and differentiation.

Haplotype structure and positive selection at TLR1.

Toll-like receptor 1, when dimerized with Toll-like receptor 2, is a cell surface receptor that, upon recognition of bacterial lipoproteins, activates the innate immune system. Variants in TLR1 associate with the risk of a variety of medical conditions and diseases, including sepsis, leprosy, tuberculosis, and others. The foremost of these is rs5743618 c.2079T>G(p.(Ile602Ser)), the derived allele of which is associated with reduced risk of sepsis, leprosy, and other diseases. Interestingly, 602Ser, which shows signatures of selection, inhibits TLR1 surface trafficking and subsequent activation of NFκB upon recognition of a ligand. This suggests that reduced TLR1 activity may be beneficial for human health. To better understand TLR1 variation and its link to human health, we have typed all 7 high-frequency missense variants (>5% in at least one population) along with 17 other variants in and around TLR1 in 2548 individuals from 56 populations from around the globe. We have also found additional signatures of selection on missense variants not associated with rs5743618, suggesting that there may be multiple functional alleles under positive selection in this gene.

Progress toward an efficient panel of SNPs for ancestry inference.

Many panels of ancestry informative single nucleotide polymorphisms have been proposed in recent years for various purposes including detecting stratification in biomedical studies and determining an individual's ancestry in a forensic context. All of the panels have limitations in their generality and efficiency for routine forensic work. Some panels have used only a few populations to validate them. Some panels are based on very large numbers of SNPs thereby limiting the ability of others to test different populations. We have been working toward an efficient and globally useful panel of ancestry informative markers that is comprised of a small number of highly informative SNPs. We have developed a panel of 55 SNPs analyzed on 73 populations from around the world. We present the details of the panel and discuss its strengths and limitations.

Mini-haplotypes as lineage informative SNPs and ancestry inference SNPs.

We propose that haplotyped loci with high heterozygosity can be useful in human identification, especially within families, if recombination is very low among the sites. Three or more SNPs extending over small molecular intervals (<10 KB) can be identified in the human genome to define miniature haplotypes with moderate levels of linkage disequilibrium. Properly selected, these mini-haplotypes (or minihaps) consist of multiple haplotype lineages (alleles) that have evolved from the ancestral human haplotype but show no evidence of recurring recombination, allowing each distinct haplotype to be equated with an allele, all copies of which are essentially identical by descent. Historic recombinants, representing rare events that have drifted to common frequencies over many generations, can be identified in some cases, they do not equate to frequently recurring recombination. We have identified examples in our data collected on various projects and present eight such mini-haplotypes comprised of informative SNPs. We also discuss the ideal characteristics and advantages of minihaps for human familial identification and ancestry inference, and compare them to other types of forensic markers in use and/or that have been proposed. We expect that it is possible to carry out a systematic search and identify a useful panel of mini-haplotypes, with even better properties than the examples presented here.

Expanding data and resources for forensic use of SNPs in individual identification.

The potential value of SNPs for individual identification has been recognized by many researchers and different panels have been proposed. Here we present a new interface in the ALFRED database to access compendia of allele frequencies for several published panels of markers for forensic uses. One of those is our panel of individual identification SNPs (IISNPs) based on samples of 44 populations originating from many parts of the world. Here we also present additional data and additional statistical analyses that continue to support the value of our panel of IISNPs as a universal panel. We also describe initial developments of multiplex methods and various robustness analyses for our 45 marker IISNP panel.

Accurate, rapid and high-throughput detection of strain-specific polymorphisms in Bacillus anthracis and Yersinia pestis by next-generation sequencing.

BACKGROUND: In the event of biocrimes or infectious disease outbreaks, high-resolution genetic characterization for identifying the agent and attributing it to a specific source can be crucial for an effective response. Until recently, in-depth genetic characterization required expensive and time-consuming Sanger sequencing of a few strains, followed by genotyping of a small number of marker loci in a panel of isolates at or by gel-based approaches such as pulsed field gel electrophoresis, which by necessity ignores most of the genome. Next-generation, massively parallel sequencing (MPS) technology (specifically the Applied Biosystems sequencing by oligonucleotide ligation and detection (SOLiD™) system) is a powerful investigative tool for rapid, cost-effective and parallel microbial whole-genome characterization. RESULTS: To demonstrate the utility of MPS for whole-genome typing of monomorphic pathogens, four Bacillus anthracis and four Yersinia pestis strains were sequenced in parallel. Reads were aligned to complete reference genomes, and genomic variations were identified. Resequencing of the B. anthracis Ames ancestor strain detected no false-positive single-nucleotide polymorphisms (SNPs), and mapping of reads to the Sterne strain correctly identified 98% of the 133 SNPs that are not clustered or associated with repeats. Three geographically distinct B. anthracis strains from the A branch lineage were found to have between 352 and 471 SNPs each, relative to the Ames genome, and one strain harbored a genomic amplification. Sequencing of four Y. pestis strains from the Orientalis lineage identified between 20 and 54 SNPs per strain relative to the CO92 genome, with the single Bolivian isolate having approximately twice as many SNPs as the three more closely related North American strains. Coverage plotting also revealed a common deletion in two strains and an amplification in the Bolivian strain that appear to be due to insertion element-mediated recombination events. Most private SNPs (that is, a, variant found in only one strain in this set) selected for validation by Sanger sequencing were confirmed, although rare false-positive SNPs were associated with variable nucleotide tandem repeats. CONCLUSIONS: The high-throughput, multiplexing capability, and accuracy of this system make it suitable for rapid whole-genome typing of microbial pathogens during a forensic or epidemiological investigation. By interrogating nearly every base of the genome, rare polymorphisms can be reliably discovered, thus facilitating high-resolution strain tracking and strengthening forensic attribution.

Developmental validation of the Quantifiler Duo DNA Quantification kit for simultaneous quantification of total human and human male DNA and detection of PCR inhibitors in biological samples.

The Quantifiler Duo DNA Quantification kit enables simultaneous quantification of human DNA and human male DNA as well as detection of inhibitors of PCR in a single real-time PCR well. Pooled human male genomic DNA is used to generate standard curves for both human (ribonuclease P RNA component H1) and human male (sex determining region Y) specific targets. A shift in the cycle threshold (C(T)) values for the internal positive control monitors the presence of PCR inhibitors in a sample. The assay is human specific and exhibits a high dynamic range from 0.023 to 50 ng/microL. In addition, the multiplex assay can detect as little as 25 pg/microL of human male DNA in the presence of a 1000-fold excess of human female DNA. The multiplex assay provides assessment of the DNA extract and guidance for the selection of the appropriate AmpFlSTR Amplification Kit to obtain interpretable short tandem repeat profiles.

An efficient and high-throughput approach for experimental validation of novel human gene predictions.

A highly automated RT-PCR-based approach has been established to validate novel human gene predictions with no prior experimental evidence of mRNA splicing (ab initio predictions). Ab initio gene predictions were selected for high-throughput validation using predicted protein classification, sequence similarity to other genomes, colocalization with an MPSS tag, or microarray expression. Initial microarray prioritization followed by RT-PCR validation was the most efficient combination, resulting in approximately 35% of the ab initio predictions being validated by RT-PCR. Of the 7252 novel genes that were prioritized and processed, 796 constituted real transcripts. In addition, high-throughput RACE successfully extended the 5' and/or 3' ends of >60% of RT-PCR-validated genes. Reevaluation of these transcripts produced 574 novel transcripts using RefSeq as a reference. RT-PCR sequencing in combination with RACE on ab initio gene predictions could be used to define the transcriptome across all species.

Transcriptional regulation of the lactase-phlorizin hydrolase promoter by PDX-1.

Lactase-phlorizin hydrolase gene expression is spatially restricted along the anterior-posterior gut axis. Lactase gene transcription is maximal in the distal duodenum and jejunum in adult mammals and is barely detectable in the proximal duodenum. By contrast, pancreatic duodenal homeobox-1 (PDX-1) protein is expressed maximally in the proximal duodenum. This study aimed to determine the role of PDX-1 in regulating lactase gene promoter activity in intestinal epithelial cells. Caco-2 cells were cotransfected with lactase promoter-reporter constructs in the presence of a PDX-1 expression vector and assayed for luciferase activity. PDX-1 cotransfection results in repression of lactase promoter activity. Sequence analysis of the lactase promoter revealed a putative PDX-1 DNA binding site in the proximal 100-bp lactase gene promoter. EMSAs demonstrated that PDX-1 can interact with the lactase promoter binding site but not with a site in which the core PDX-1 binding sequence TAAT is mutated. Site-directed mutagenesis of the PDX-1 core binding site in the lactase promoter-reporter construct suggests that PDX-1 can function independently of DNA binding to its consensus binding site. Stable overexpression of PDX-1 results in repression of the endogenous human lactase gene in differentiated Caco-2 cells. Given the contrasting spatial expression pattern, PDX-1 may function to specify the anterior boundary of lactase expression in the small intestine and is thus a candidate regulator of anterior spatial restriction in the gut.