Alteration of oral microbiome composition in children living with pesticide-exposed farm workers.

Our prior work shows that azinphos-methyl pesticide exposure is associated with altered oral microbiomes in exposed farmworkers. Here we extend this analysis to show the same association pattern is also evident in their children. Oral buccal swab samples were analyzed at two time points, the apple thinning season in spring-summer 2005 for 78 children and 101 adults and the non-spray season in winter 2006 for 62 children and 82 adults. The pesticide exposure for the children were defined by the farmworker occupation of the cohabitating household adult and the blood azinphos-methyl detection of the cohabitating adult. Oral buccal swab 16S rRNA sequencing determined taxonomic microbiota proportional composition from concurrent samples from both adults and children. Analysis of the identified bacteria showed significant proportional changes for 12 of 23 common oral microbiome genera in association with azinphos-methyl detection and farmworker occupation. The most common significantly altered genera had reductions in the abundance of Streptococcus, suggesting an anti-microbial effect of the pesticide. Principal component analysis of the microbiome identified two primary clusters, with association of principal component 1 to azinphos-methyl blood detection and farmworker occupational status of the household. The children's buccal microbiota composition clustered with their household adult in ∼95% of the households. Household adult farmworker occupation and household pesticide exposure is associated with significant alterations in their children's oral microbiome composition. This suggests that parental occupational exposure and pesticide take-home exposure pathways elicit alteration of their children's microbiomes.

Variability in metagenomic samples from the Puget Sound: Relationship to temporal and anthropogenic impacts.

Whole-metagenome sequencing (WMS) has emerged as a powerful tool to assess potential public health risks in marine environments by measuring changes in microbial community structure and function in uncultured bacteria. In addition to monitoring public health risks such as antibiotic resistance determinants, it is essential to measure predictors of microbial variation in order to identify natural versus anthropogenic factors as well as to evaluate reproducibility of metagenomic measurements.This study expands our previous metagenomic characterization of Puget Sound by sampling new nearshore environments including the Duwamish River, an EPA superfund site, and the Hood Canal, an area characterized by highly variable oxygen levels. We also resampled a wastewater treatment plant, nearshore and open ocean sites introducing a longitudinal component measuring seasonal and locational variations and establishing metagenomics sampling reproducibility. Microbial composition from samples collected in the open sound were highly similar within the same season and location across different years, while nearshore samples revealed multi-fold seasonal variation in microbial composition and diversity. Comparisons with recently sequenced predominant marine bacterial genomes helped provide much greater species level taxonomic detail compared to our previous study. Antibiotic resistance determinants and pollution and detoxification indicators largely grouped by location showing minor seasonal differences. Metal resistance, oxidative stress and detoxification systems showed no increase in samples proximal to an EPA superfund site indicating a lack of ecosystem adaptation to anthropogenic impacts. Taxonomic analysis of common sewage influent families showed a surprising similarity between wastewater treatment plant and open sound samples suggesting a low-level but pervasive sewage influent signature in Puget Sound surface waters. Our study shows reproducibility of metagenomic data sampling in multiple Puget Sound locations while establishing baseline measurements of antibiotic resistance determinants, pollution and detoxification systems. Combining seasonal and longitudinal data across these locations provides a foundation for evaluating variation in future studies.

FARME DB: a functional antibiotic resistance element database.

Antibiotic resistance (AR) is a major global public health threat but few resources exist that catalog AR genes outside of a clinical context. Current AR sequence databases are assembled almost exclusively from genomic sequences derived from clinical bacterial isolates and thus do not include many microbial sequences derived from environmental samples that confer resistance in functional metagenomic studies. These environmental metagenomic sequences often show little or no similarity to AR sequences from clinical isolates using standard classification criteria. In addition, existing AR databases provide no information about flanking sequences containing regulatory or mobile genetic elements. To help address this issue, we created an annotated database of DNA and protein sequences derived exclusively from environmental metagenomic sequences showing AR in laboratory experiments. Our Functional Antibiotic Resistant Metagenomic Element (FARME) database is a compilation of publically available DNA sequences and predicted protein sequences conferring AR as well as regulatory elements, mobile genetic elements and predicted proteins flanking antibiotic resistant genes. FARME is the first database to focus on functional metagenomic AR gene elements and provides a resource to better understand AR in the 99% of bacteria which cannot be cultured and the relationship between environmental AR sequences and antibiotic resistant genes derived from cultured isolates.Database URL: http://staff.washington.edu/jwallace/farme.

Human Oral Buccal Microbiomes Are Associated with Farmworker Status and Azinphos-Methyl Agricultural Pesticide Exposure.

In a longitudinal agricultural community cohort sampling of 65 adult farmworkers and 52 adult nonfarmworkers, we investigated agricultural pesticide exposure-associated changes in the oral buccal microbiota. We found a seasonally persistent association between the detected blood concentration of the insecticide azinphos-methyl and the taxonomic composition of the buccal swab oral microbiome. Blood and buccal samples were collected concurrently from individual subjects in two seasons, spring/summer 2005 and winter 2006. Mass spectrometry quantified blood concentrations of the organophosphate insecticide azinphos-methyl. Buccal oral microbiome samples were 16S rRNA gene DNA sequenced, assigned to the bacterial taxonomy, and analyzed after "centered-log-ratio" transformation to handle the compositional nature of the proportional abundances of bacteria per sample. Nonparametric analysis of the transformed microbiome data for individuals with and without azinphos-methyl blood detection showed significant perturbations in seven common bacterial taxa (>0.5% of sample mean read depth), including significant reductions in members of the common oral bacterial genus Streptococcus Diversity in centered-log-ratio composition between individuals' microbiomes was also investigated using principal-component analysis (PCA) to reveal two primary PCA clusters of microbiome types. The spring/summer "exposed" microbiome cluster with significantly less bacterial diversity was enriched for farmworkers and contained 27 of the 30 individuals who also had azinphos-methyl agricultural pesticide exposure detected in the blood. IMPORTANCE: In this study, we show in human subjects that organophosphate pesticide exposure is associated with large-scale significant alterations of the oral buccal microbiota composition, with extinctions of whole taxa suggested in some individuals. The persistence of this association from the spring/summer to the winter also suggests that long-lasting effects on the commensal microbiota have occurred. The important health-related outcomes of these agricultural community individuals' pesticide-associated microbiome perturbations are not understood at this time. Future investigations should index medical and dental records for common and chronic diseases that may be interactively caused by this association between pesticide exposure and microbiome alteration.

Metagenomic frameworks for monitoring antibiotic resistance in aquatic environments.

BACKGROUND: High-throughput genomic technologies offer new approaches for environmental health monitoring, including metagenomic surveillance of antibiotic resistance determinants (ARDs). Although natural environments serve as reservoirs for antibiotic resistance genes that can be transferred to pathogenic and human commensal bacteria, monitoring of these determinants has been infrequent and incomplete. Furthermore, surveillance efforts have not been integrated into public health decision making. OBJECTIVES: We used a metagenomic epidemiology-based approach to develop an ARD index that quantifies antibiotic resistance potential, and we analyzed this index for common modal patterns across environmental samples. We also explored how metagenomic data such as this index could be conceptually framed within an early risk management context. METHODS: We analyzed 25 published data sets from shotgun pyrosequencing projects. The samples consisted of microbial community DNA collected from marine and freshwater environments across a gradient of human impact. We used principal component analysis to identify index patterns across samples. RESULTS: We observed significant differences in the overall index and index subcategory levels when comparing ecosystems more proximal versus distal to human impact. The selection of different sequence similarity thresholds strongly influenced the index measurements. Unique index subcategory modes distinguished the different metagenomes. CONCLUSIONS: Broad-scale screening of ARD potential using this index revealed utility for framing environmental health monitoring and surveillance. This approach holds promise as a screening tool for establishing baseline ARD levels that can be used to inform and prioritize decision making regarding management of ARD sources and human exposure routes. CITATION: Port JA, Cullen AC, Wallace JC, Smith MN, Faustman EM. 2014. Metagenomic frameworks for monitoring antibiotic resistance in aquatic environments. Environ Health Perspect 122:222­228; http://dx.doi.org/10.1289/ehp.1307009

Identification of G protein-coupled receptor signaling pathway proteins in marine diatoms using comparative genomics.

BACKGROUND: The G protein-coupled receptor (GPCR) signaling pathway plays an essential role in signal transmission and response to external stimuli in mammalian cells. Protein components of this pathway have been characterized in plants and simpler eukaryotes such as yeast, but their presence and role in unicellular photosynthetic eukaryotes have not been determined. We use a comparative genomics approach using whole genome sequences and gene expression libraries of four diatoms (Pseudo-nitzschia multiseries, Thalassiosira pseudonana, Phaeodactylum tricornutum and Fragilariopsis cylindrus) to search for evidence of GPCR signaling pathway proteins that share sequence conservation to known GPCR pathway proteins. RESULTS: The majority of the core components of GPCR signaling were well conserved in all four diatoms, with protein sequence similarity to GPCRs, human G protein α- and ß-subunits and downstream effectors. There was evidence for the Gγ-subunit and thus a full heterotrimeric G protein only in T. pseudonana. Phylogenetic analysis of putative diatom GPCRs indicated similarity but deep divergence to the class C GPCRs, with branches basal to the GABAB receptor subfamily. The extracellular and intracellular regions of these putative diatom GPCR sequences exhibited large variation in sequence length, and seven of these sequences contained the necessary ligand binding domain for class C GPCR activation. Transcriptional data indicated that a number of the putative GPCR sequences are expressed in diatoms under various stress conditions in culture, and that many of the GPCR-activated signaling proteins, including the G protein, are also expressed. CONCLUSIONS: The presence of sequences in all four diatoms that code for the proteins required for a functional mammalian GPCR pathway highlights the highly conserved nature of this pathway and suggests a complex signaling machinery related to environmental perception and response in these unicellular organisms. The lack of evidence for some GPCR pathway proteins in one or more of the diatoms, such as the Gγ-subunit, may be due to differences in genome completeness and genome coverage for the four diatoms. The high divergence of putative diatom GPCR sequences to known class C GPCRs suggests these sequences may represent another, potentially ancestral, subfamily of class C GPCRs.

A resource management tool for public health continuity of operations during disasters.

OBJECTIVE: We developed and validated a user-centered information system to support the local planning of public health continuity of operations for the Community Health Services Division, Public Health - Seattle & King County, Washington. METHODS: The Continuity of Operations Data Analysis (CODA) system was designed as a prototype developed using requirements identified through participatory design. CODA uses open-source software that links personnel contact and licensing information with needed skills and clinic locations for 821 employees at 14 public health clinics in Seattle and King County. Using a web-based interface, CODA can visualize locations of personnel in relationship to clinics to assist clinic managers in allocating public health personnel and resources under dynamic conditions. RESULTS: Based on user input, the CODA prototype was designed as a low-cost, user-friendly system to inventory and manage public health resources. In emergency conditions, the system can run on a stand-alone battery-powered laptop computer. A formative evaluation by managers of multiple public health centers confirmed the prototype design's usefulness. Emergency management administrators also provided positive feedback about the system during a separate demonstration. CONCLUSIONS: Validation of the CODA information design prototype by public health managers and emergency management administrators demonstrates the potential usefulness of building a resource management system using open-source technologies and participatory design principles.

Metagenomic profiling of microbial composition and antibiotic resistance determinants in Puget Sound.

Human-health relevant impacts on marine ecosystems are increasing on both spatial and temporal scales. Traditional indicators for environmental health monitoring and microbial risk assessment have relied primarily on single species analyses and have provided only limited spatial and temporal information. More high-throughput, broad-scale approaches to evaluate these impacts are therefore needed to provide a platform for informing public health. This study uses shotgun metagenomics to survey the taxonomic composition and antibiotic resistance determinant content of surface water bacterial communities in the Puget Sound estuary. Metagenomic DNA was collected at six sites in Puget Sound in addition to one wastewater treatment plant (WWTP) that discharges into the Sound and pyrosequenced. A total of ~550 Mbp (1.4 million reads) were obtained, 22 Mbp of which could be assembled into contigs. While the taxonomic and resistance determinant profiles across the open Sound samples were similar, unique signatures were identified when comparing these profiles across the open Sound, a nearshore marina and WWTP effluent. The open Sound was dominated by α-Proteobacteria (in particular Rhodobacterales sp.), γ-Proteobacteria and Bacteroidetes while the marina and effluent had increased abundances of Actinobacteria, ß-Proteobacteria and Firmicutes. There was a significant increase in the antibiotic resistance gene signal from the open Sound to marina to WWTP effluent, suggestive of a potential link to human impacts. Mobile genetic elements associated with environmental and pathogenic bacteria were also differentially abundant across the samples. This study is the first comparative metagenomic survey of Puget Sound and provides baseline data for further assessments of community composition and antibiotic resistance determinants in the environment using next generation sequencing technologies. In addition, these genomic signals of potential human impact can be used to guide initial public health monitoring as well as more targeted and functionally-based investigations.

Long-term skeletal and dental stability of mandibular symphyseal distraction osteogenesis with a hybrid distractor.

INTRODUCTION: The purpose of this study was to examine the long-term skeletal and dental stability of mandibular symphyseal distraction osteogenesis (MSDO) with a tooth-borne and bone-borne hybrid distractor. To differentiate the effects of MSDO from the orthodontic movement and relapse, each phase of treatment was analyzed. METHODS: Twenty-five patients were included in the study, ranging in age from 12.0 to 30.9 years at the initiation of treatment (mean, 15.8 ± 4.8 years). Of this group, 16 patients were recalled at a mean of 7.5 ± 0.9 years (range, 6.3-9.6 years) after distraction for long-term analysis of skeletal and dental changes. Orthodontic records were taken at 5 times: T1, pretreatment; T2, predistraction; T3, postdistraction; T4, posttreatment, and T5, postretention. The data were statistically analyzed by using repeated-measures analysis of variance (ANOVA). RESULTS: There were significant increases in all interdental transverse measurements except the mandibular intersecond molar distance from T1 to T4. The largest overall expansion was achieved between the mandibular second premolars (4.32 ± 0.60 mm), followed by the interfirst premolar (3.44 ± 0.44 mm), the interfirst molar (2.60 ± 0.65 mm), and the intercanine (1.87 ± 0.44 mm) widths. The overall amount of transverse dental expansion was substantially less when analyzed from the time of the mandibular symphyseal osteotomy to posttreatment (T2-T4). From T3 to T4, there were significant decreases between the mandibular intersecond premolars (-3.10 ± 0.52 mm), interfirst premolars (-3.90 ± 0.35 m), intercanines (-4.47 ± 0.38 mm), and intercentral incisors (-5.60 ± 0.32 mm). There were no significant changes in bicondylar, bigonial, and biantigonial widths. At the long-term follow-up, there were no significant changes in the interdental or skeletal measurements between T4 and T5, except for interincisor apices. The irregularity index significantly decreased during the orthodontic treatment but significantly increased in the long-term follow-up period (T4-T5). After the MSDO, T3 to T5, the results indicated symphyseal basal bone skeletal stability. CONCLUSIONS: The results indicate that the expansion of the mandibular arch with MSDO and conventional orthodontic mechanics produces no statistically significant transverse changes from posttreatment to long-term follow-up. The risks of using a surgical procedure and MSDO to achieve additional expansion should be evaluated by the clinician and compared with more traditional orthodontic methods.

Integrating open-source technologies to build low-cost information systems for improved access to public health data.

Effective public health practice relies on the availability of public health data sources and assessment tools to convey information to investigators, practitioners, policy makers, and the general public. Emerging communication technologies on the Internet can deliver all components of the "who, what, when, and where" quartet more quickly than ever with a potentially higher level of quality and assurance, using new analysis and visualization tools. Open-source software provides the opportunity to build low-cost information systems allowing health departments with modest resources access to modern data analysis and visualization tools. In this paper, we integrate open-source technologies and public health data to create a web information system which is accessible to a wide audience through the Internet. Our web application, "EpiVue," was tested using two public health datasets from the Washington State Cancer Registry and Washington State Center for Health Statistics. A third dataset shows the extensibility and scalability of EpiVue in displaying gender-based longevity statistics over a twenty-year interval for 3,143 United States counties. In addition to providing an integrated visualization framework, EpiVue's highly interactive web environment empowers users by allowing them to upload their own geospatial public health data in either comma-separated text files or MS Excel spreadsheet files and visualize the geospatial datasets with Google Maps.

High-density rhesus macaque oligonucleotide microarray design using early-stage rhesus genome sequence information and human genome annotations.

BACKGROUND: Until recently, few genomic reagents specific for non-human primate research have been available. To address this need, we have constructed a macaque-specific high-density oligonucleotide microarray by using highly fragmented low-pass sequence contigs from the rhesus genome project together with the detailed sequence and exon structure of the human genome. Using this method, we designed oligonucleotide probes to over 17,000 distinct rhesus/human gene orthologs and increased by four-fold the number of available genes relative to our first-generation expressed sequence tag (EST)-derived array. RESULTS: We constructed a database containing 248,000 exon sequences from 23,000 human RefSeq genes and compared each human exon with its best matching sequence in the January 2005 version of the rhesus genome project list of 486,000 DNA contigs. Best matching rhesus exon sequences for each of the 23,000 human genes were then concatenated in the proper order and orientation to produce a rhesus "virtual transcriptome." Microarray probes were designed, one per gene, to the region closest to the 3' untranslated region (UTR) of each rhesus virtual transcript. Each probe was compared to a composite rhesus/human transcript database to test for cross-hybridization potential yielding a final probe set representing 18,296 rhesus/human gene orthologs, including transcript variants, and over 17,000 distinct genes. We hybridized mRNA from rhesus brain and spleen to both the EST- and genome-derived microarrays. Besides four-fold greater gene coverage, the genome-derived array also showed greater mean signal intensities for genes present on both arrays. Genome-derived probes showed 99.4% identity when compared to 4,767 rhesus GenBank sequence tag site (STS) sequences indicating that early stage low-pass versions of complex genomes are of sufficient quality to yield valuable functional genomic information when combined with finished genome information from a closely related species. CONCLUSION: The number of different genes represented on microarrays for unfinished genomes can be greatly increased by matching known gene transcript annotations from a closely related species with sequence data from the unfinished genome. Signal intensity on both EST- and genome-derived arrays was highly correlated with probe distance from the 3' UTR, information often missing from ESTs yet present in early-stage genome projects.

Discovery of functional noncoding elements by digital analysis of chromatin structure.

We developed a quantitative methodology, digital analysis of chromatin structure (DACS), for high-throughput, automated mapping of DNase I-hypersensitive sites and associated cis-regulatory sequences in the human and other complex genomes. We used 19/20-bp genomic DNA tags to localize individual DNase I cutting events in nuclear chromatin and produced approximately 257,000 tags from erythroid cells. Tags were mapped to the human genome, and a quantitative algorithm was applied to discriminate statistically significant clusters of independent DNase I cutting events. We show that such clusters identify both known regulatory sequences and previously unrecognized functional elements across the genome. We used in silico simulation to demonstrate that DACS is capable of efficient and accurate localization of the majority of DNase I-hypersensitive sites in the human genome without requiring an independent validation step. A unique feature of DACS is that it permits unbiased evaluation of the chromatin state of regulatory sequences from widely separated genomic loci. We found surprisingly large differences in the accessibility of distant regulatory sequences, suggesting the existence of a hierarchy of nuclear organization that escapes detection by conventional chromatin assays.

Genome-wide identification of DNaseI hypersensitive sites using active chromatin sequence libraries.

Comprehensive identification of sequences that regulate transcription is one of the major goals of genome biology. Focal alteration in chromatin structure in vivo, detectable through hypersensitivity to DNaseI and other nucleases, is the sine qua non of a diverse cast of transcriptional regulatory elements including enhancers, promoters, insulators, and locus control regions. We developed an approach for genome-scale identification of DNaseI hypersensitive sites (HSs) via isolation and cloning of in vivo DNaseI cleavage sites to create libraries of active chromatin sequences (ACSs). Here, we describe analysis of >61,000 ACSs derived from erythroid cells. We observed peaks in the density of ACSs at the transcriptional start sites of known genes at non-gene-associated CpG islands, and, to a lesser degree, at evolutionarily conserved noncoding sequences. Peaks in ACS density paralleled the distribution of DNaseI HSs. ACSs and DNaseI HSs were distributed between both expressed and nonexpressed genes, suggesting that a large proportion of genes reside within open chromatin domains. The results permit a quantitative approximation of the distribution of HSs and classical cis-regulatory sequences in the human genome.

Unilateral Brodie bite treated with distraction osteogenesis.

Mandibular widening by distraction osteogenesis has recently been shown to be a very predictable, successful, and stable procedure with limited morbidity. One indication for mandibular widening by distraction osteogenesis is a Brodie bite (scissors or buccal crossbite). This case report demonstrates a technique used to treat a 12-year-old girl with a true unilateral buccal crossbite (no functional shift). It details how distraction osteogenesis can successfully treat unilateral problems with a distractor, a splint, and elastics. Details on this treatment technique have not been previously published.

High-throughput localization of functional elements by quantitative chromatin profiling.

Identification of functional, noncoding elements that regulate transcription in the context of complex genomes is a major goal of modern biology. Localization of functionality to specific sequences is a requirement for genetic and computational studies. Here, we describe a generic approach, quantitative chromatin profiling, that uses quantitative analysis of in vivo chromatin structure over entire gene loci to rapidly and precisely localize cis-regulatory sequences and other functional modalities encoded by DNase I hypersensitive sites. To demonstrate the accuracy of this approach, we analyzed approximately 300 kilobases of human genome sequence from diverse gene loci and cleanly delineated functional elements corresponding to a spectrum of classical cis-regulatory activities including enhancers, promoters, locus control regions and insulators as well as novel elements. Systematic, high-throughput identification of functional elements coinciding with DNase I hypersensitive sites will substantially expand our knowledge of transcriptional regulation and should simplify the search for noncoding genetic variation with phenotypic consequences.