The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

Challenges in using siliceous subfossils as a tool for inferring past water level and hydroperiod in Everglades marshes

Successfully rehabilitating drained wetlands through hydrologic restoration is dependent on defining restoration targets, a process that is informed by pre-drainage conditions, as well as understanding linkages between hydrology and ecosystem structure. Paleoecological records can inform restoration goals by revealing long-term patterns of change, but are dependent on preservation of biomarkers that provide meaningful interpretations of environmental change. In the Florida Everglades, paleohydrological hind-casting could improve restoration forecasting, but frequent drying of marsh soils leads to poor preservation of many biomarkers. To determine the effectiveness of employing siliceous subfossils in paleohydrological reconstructions, we examined diatoms, plant and sponge silico-sclerids from three soil cores in the central Everglades marshes. Subfossil quality varied among cores, but the abundance of recognizable specimens was sufficient to infer 1,000–3,000 years of hydrologic change at decadal to centennial resolution. Phytolith morphotypes were linked to key marsh plant species to indirectly measure fluctuations in water depth. A modern dataset was used to derive diatom-based inferences of water depth and hydroperiod (R2 = 0.63, 0.47; RMSE = 14 cm, 120 days, respectively). Changes in subfossil quality and abundances at centennial time-scales were associated with mid-Holocene climate events including the Little Ice Age and Medieval Warm Period, while decadal-scale fluctuations in assemblage structure during the twentieth century suggested co-regulation of hydrology by cyclical climate drivers (particularly the Atlantic Multidecadal Oscillation) and water management changes. The successful reconstructions based on siliceous subfossils shown here at a coarse temporal scale (i.e., decadal to centennial) advocate for their application in more highly resolved (i.e., subdecadal) records, which should improve the ability of water managers to target the quantity and variability of water flows appropriate for hydrologic restoration.

Identification of Conserved Splicing Motifs in Mutually Exclusive Exons of 15 Insect Species

Background: During alternative splicing, the inclusion of an exon in the final mRNA molecule is determined by nuclear proteins that bind cis-regulatory sequences in a target pre-mRNA molecule. A recent study suggested that the regulatory codes of individual RNA-binding proteins may be nearly immutable between very diverse species such as mammals and insects. The model system Drosophila melanogaster therefore presents an excellent opportunity for the study of alternative splicing due to the availability of quality EST annotations in FlyBase. Methods: In this paper, we describe an in silico analysis pipeline to extract putative exonic splicing regulatory sequences from a multiple alignment of 15 species of insects. Our method, ESTs-to-ESRs (E2E), uses graph analysis of EST splicing graphs to identify mutually exclusive (ME) exons and combines phylogenetic measures, a sliding window approach along the multiple alignment and the Welch’s t statistic to extract conserved ESR motifs. Results: The most frequent 100% conserved word of length 5 bp in different insect exons was “ATGGA”. We identified 799 statistically significant “spike” hexamers, 218 motifs with either a left or right FDR corrected spike magnitude p-value < 0.05 and 83 with both left and right uncorrected p < 0.01. 11 genes were identified with highly significant motifs in one ME exon but not in the other, suggesting regulation of ME exon splicing through these highly conserved hexamers. The majority of these genes have been shown to have regulated spatiotemporal expression. 10 elements were found to match three mammalian splicing regulator databases. A putative ESR motif, GATGCAG, was identified in the ME-13b but not in the ME-13a of Drosophila N-Cadherin, a gene that has been shown to have a distinct spatiotemporal expression pattern of spliced isoforms in a recent study. Conclusions: Analysis of phylogenetic relationships and variability of sequence conservation as implemented in the E2E spikes method may lead to improved identification of ESRs. We found that approximately half of the putative ESRs in common between insects and mammals have a high statistical support (p < 0.01). Several Drosophila genes with spatiotemporal expression patterns were identified to contain putative ESRs located in one exon of the ME exon pairs but not in the other.

Better primer design for metagenomics applications by increasing taxonomic distinguishability

Current methods of understanding microbiome composition and structure rely on accurately estimating the number of distinct species and their relative abundance. Most of these methods require an efficient PCR whose forward and reverse primers bind well to the same, large number of identifiable species, and produce amplicons that are unique. It is therefore not surprising that currently used universal primers designed many years ago are not as efficient and fail to bind to recently cataloged species. We propose an automated general method of designing PCR primer pairs that abide by primer design rules and uses current sequence database as input. Since the method is automated, primers can be designed for targeted microbial species or updated as species are added or deleted from the database. In silico experiments and laboratory experiments confirm the efficacy of the newly designed primers for metagenomics applications.

Mechanisms of Arsenic Detoxification and Resistance

Arsenic is a ubiquitous environmental toxic substance. As a consequence of continual exposure to arsenic, nearly every organism, from Escherichia coli to humans have evolved arsenic detoxification pathways. One of the pathways is extrusion of arsenic from inside the cells, thereby conferring resistance. The R773 arsRDABC operon in E. coli encodes an ArsAB efflux pump that confers resistance to arsenite. ArsA is the catalytic subunit of the pump, while ArsB forms the oxyanion conducting pathway. ArsD is an arsenite metallochaperone that binds arsenite and transfers it to ArsA. The interaction of ArsA and ArsD allows for resistance to As(III) at environmental concentrations. The interaction between ArsA ATPase and ArsD metallochaperone was examined. A quadruple mutant in the arsD gene encoding a K2A/K37A/K62A/K104A ArsD is unable to interact with ArsA. An error-prone mutagenesis approach was used to generate random mutations in the arsA gene that restored interaction with the quadruple arsD mutant in yeast two-hybrid assays. Three such mutants encoding Q56R, F120I and D137V ArsA were able to restore interaction with the quadruple ArsD mutant. Structural models generated by in silico docking suggest that an electrostatic interface favors reversible interaction between ArsA and ArsD. Mutations in ArsA that propagate changes in hydrogen bonding and salt bridges to the ArsA-ArsD interface also affect their interactions. The second objective was to examine the mechanism of arsenite resistance through methylation and subsequent volatilization. Microbial ArsM (As(III) S-adenosylmethyltransferase) catalyzes the formation of trimethylarsine as the volatile end product. The net result is loss of arsenic from cells. The gene for CrArsM from the eukaryotic green alga Chlamydomonas reinhardtii was chemically synthesized and expressed in E. coli. The purified protein catalyzed the methylation of arsenite into methyl-, dimethyl- and trimethyl products. Synthetic purified CrArsM was crystallized in an unliganded form. Biochemical and biophysical studies conducted on CrArsM sheds new light on the pathways of biomethylation. While in microbes ArsM detoxifies arsenic, the human homolog, hAS3MT, converts inorganic arsenic into more toxic and carcinogenic forms. An understanding of the enzymatic mechanism of ArsM will be critical in deciphering its parallel roles in arsenic detoxification and carcinogenesis.

Filogenia molecular das enzimas isocitrato liase e malato sintase e sua evolução em Viridiplanta

The plant metabolism consists of a complex network of physical and chemical events resulting in photosynthesis, respiration, synthesis and degradation of organic compounds. This is only possible due to the different kinds of responses to many environmental variations that a plant could be subject through evolution, leading also to conquering new surroundings. The glyoxylate cycle is a metabolic pathway found in glyoxysomes plant, which has unique role in the seedling establishment. Considered as a variation of the citric acid cycle, it uses an acetyl coenzyme A molecule, derived from lipids beta-oxidation to synthesize compounds which are used in carbohydrate synthesis. The Malate synthase (MLS) and Isocitrate lyase (ICL) enzyme of this cycle are unique and essential in regulating the biosynthesis of carbohydrates. Because of the absence of decarboxylation steps as rate-limiting steps, detailed studies of molecular phylogeny and evolution of these proteins enables the elucidation of the effects of this route presence in the evolutionary processes involved in their distribution across the genome from different plant species. Therefore, the aim of this study was to establish a relationship between the molecular evolution of the characteristics of enzymes from the glyoxylate cycle (isocitrate lyase and malate synthase) and their molecular phylogeny, among green plants (Viridiplantae). For this, amino acid and nucleotide sequences were used, from online repositories as UniProt and Genbank. Sequences were aligned and then subjected to an analysis of the best-fit substitution models. The phylogeny was rebuilt by distance methods (neighbor-joining) and discrete methods (maximum likelihood, maximum parsimony and Bayesian analysis). The identification of structural patterns in the evolution of the enzymes was made through homology modeling and structure prediction from protein sequences. Based on comparative analyzes of in silico models and from the results of phylogenetic inferences, both enzymes show significant structure conservation and their topologies in agreement with two processes of selection and specialization of the genes. Thus, confirming the relevance of new studies to elucidate the plant metabolism from an evolutionary perspective

Utilização de peptídeos sintéticos derivados de moléculas imunodominantes de Toxoplasma gondii para diagnóstico sorológico da toxoplasmose em ovinos e galinhas caipiras

Chapter I - The obligate intracellular parasite Toxoplasma gondii is the causative agent of toxoplasmosis, a disease that affects humans and generates economic losses in farm animals. When prevention fails disease refers to the diagnosis and subsequent treatment if the individual is diagnosed as positive. Therefore, the development of new accurate diagnostic tools for detecting T. gondii infection is a need in particular to determine the environmental source of infection which can result in more appropriate public health policies against different routes of infection and prevent potential damage that toxoplasmosis can cause when animals are infected. Chapter II - The domestic chicken (Gallus gallus domesticus) are considered epidemiological sentinels, still representing a major source of recombinant strains when predated by cats, it is common to find them with multiple infections. We evaluate the diagnostic potential of six synthetic peptides (SAG2Y, MIC1, M2AP, GRA10, ROP2 and ROP7) predicted in silico from tachyzoites immunodominant markers of T. gondii in samples from naturally infected chickens, comparing synthetic peptides with antigen total soluble (STAg). In general, our results demonstrated that reactivity rates and positivity for these peptides are similar to the STAg, and the ROP7 peptide and the combination of peptides MIC1+ROP2 have significant sensitivity, confirming them as potential diagnostic tools for the diagnosis of toxoplasmosis in chickens. Chapter III - Sheep (Ovis aries) are commonly infected with Toxoplasma gondii due to his eating habits. Infection in pregnant sheep can have serious consequences such as embryonic death, fetal resorption, mummification, and neonatal death. One concern regarding the infection in these animals is that the meat can be a source of contamination to humans and other carnivores. Therefore perform accurate diagnosis in these animals is of fundamental importance. In the present study we evaluated the potential of new synthetic peptides as a diagnostic tool. Synthetic peptides (SAG2Y, SRS52A, MIC14, GRA4, GRA10 and GRA15) were predicted in silico from immunodominant proteins of T. gondii. We determine the levels of IgG antibodies using sera obtained from two farms in the city of Uberlândia. Analyzing the results together, the peptide combination GRA10+GRA15 (Accuracy = 0,82) showed better characteristics compared with the other mixtures. This preparation could be better analyzed with an antigenic mixture potential use in the diagnosis of toxoplasmosis in sheep and other species.

Problemática nutricional en escolares de la Comunidad de Madrid

"Introducción y Objetivo" En los últimos años se han producido importantes cambios en los estilos de vida de la población y con ello en la alimentación de niños y adolescentes. El mayor desarrollo económico, el avance tecnológico alimentario y culinario, el seguimiento de hábitos alimentarios inadecuados, como por ejemplo no desayunar, la incorporación de la mujer al trabajo fuera del hogar, la nueva estructura familiar, la gran influencia de la publicidad y de la televisión, la incorporación más temprana de los niños a la escuela, junto con la mayor posibilidad de elegir sus menús sin una adecuada supervisión familiar, entre otras, condiciona un consumo de alimentos con elevado aporte calórico y baja calidad nutricional. Debido a lo anterior cada vez es más frecuente la presencia de obesidad y sobrepeso en la población escolar, así como la aparición de otras conductas del comportamiento alimentario y la ingesta insuficiente de algunos micronutrientes como el caso de la vitamina D, calcio, hierro, acido fólico, etc. Por otra parte, también es destacable el aumento de tabaquismo y del sedentarismo entre los jóvenes de nuestro país. Como consecuencia de lo anterior, se ha producido un aumento de problemas de salud entre nuestros niños y adolescentes, entre los que se encuentran las enfermedades cardiovasculares y el síndrome metabólico. Por todo lo anteriormente descrito, el objetivo de este trabajo fue conocer el estado nutricional de un colectivo de escolares españoles de la ciudad de Madrid. Material y Métodos Para conseguir datos, se agruparon los 21 distritos de Madrid Capital, en función de sus características socioeconómicas similares y buscando también que todas las zonas geográficas de Madrid queden representadas. Se constituyeron 5 grupos escogiendo al azar un distrito de cada grupo resultando al final los barrios de Carabanchel, Moratalaz, Chamberí, Villaverde y Fuencarral...

Développement d'approches quantitatives de type "structure-activité" pour la modélisation pharmacocinétique

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Identification de gènes ciblès par ETV6-AML1, un facteur de transcription chimérique retrouvé dans la leucémie de l'enfant

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Variable mutation rates as an adaptive strategy in replicator populations

For evolving populations of replicators, there is much evidence that the effect of mutations on fitness depends on the degree of adaptation to the selective pressures at play. In optimized populations, most mutations have deleterious effects, such that low mutation rates are favoured. In contrast to this, in populations thriving in changing environments a larger fraction of mutations have beneficial effects, providing the diversity necessary to adapt to new conditions. What is more, non-adapted populations occasionally benefit from an increase in the mutation rate. Therefore, there is no optimal universal value of the mutation rate and species attempt to adjust it to their momentary adaptive needs. In this work we have used stationary populations of RNA molecules evolving in silico to investigate the relationship between the degree of adaptation of an optimized population and the value of the mutation rate promoting maximal adaptation in a short time to a new selective pressure. Our results show that this value can significantly differ from the optimal value at mutation-selection equilibrium, being strongly influenced by the structure of the population when the adaptive process begins. In the short-term, highly optimized populations containing little variability respond better to environmental changes upon an increase of the mutation rate, whereas populations with a lower degree of optimization but higher variability benefit from reducing the mutation rate to adapt rapidly. These findings show a good agreement with the behaviour exhibited by actual organisms that replicate their genomes under broadly different mutation rates. © 2010 Stich et al.

An investigation of the molecular interactions between statins and bacterial pathogens, and their combined impact on the human immune system

Statins are a class of drug that inhibits cholesterol biosynthesis, and are used to treat patients with high serum cholesterol levels. They exert this function by competitively binding to the enzyme 3-hydroxy-3-methylglutaryl-CoenzymeA reductase (HMGR), which catalyses the formation of mevalonate, a rate-limiting step in cholesterol biosynthesis. In addition, statins have what are called “pleiotropic effects”, which include the reduction of inflammation, immunomodulation, and antimicrobial effects. Statins can also improve survival of patients with sepsis and pneumonia. Cystic fibrosis (CF) is the most common recessive inherited disease in the Caucasian population, which is characterised by factors including, but not limited to, excessive lung inflammation and increased susceptibility to infection. Therefore, the overall objective of this study was to examine the effects of statins on CFassociated bacterial pathogens and the host response. In this work, the prevalence of HMGR was examined in respiratory pathogens, and several CF-associated pathogens were found to possess homologues of this enzyme. HMGR homology was analysed in Staphylococcus aureus, Burkholderia cenocepacia and Streptococcus pneumoniae, and the HMGR of B. cenocepacia was found to have significant conservation to that of Pseudomonas mevalonii, which is the most widely-characterised bacterial HMGR. However, in silico analysis revealed that, unlike S. aureus and S. pneumoniae, B. cenocepacia did not possess homologues of other mevalonate pathway proteins, and that the HMGR of B. cenocepacia appeared to be involved in an alternative metabolic pathway. The effect of simvastatin was subsequently tested on the growth and virulence of S. aureus, B. cenocepacia and S. pneumoniae. Simvastatin inhibited the growth of all 3 species in a dose-dependent manner. In addition, statin treatment also attenuated biofilm formation of all 3 species, and reduced in vitro motility of S. aureus. Interestingly, simvastatin also increased the potency of the aminoglycoside antibiotic gentamicin against B. cenocepacia. The impact of statins was subsequently tested on the predominant CF-associated pathogen Pseudomonas aeruginosa, which does not possess a HMGR homologue. Mevastatin, lovastatin and simvastatin did not influence the growth of this species. However, sub-inhibitory statin concentrations reduced the swarming motility and biofilm formation of P. aeruginosa. The influence of statins was also examined on Type 3 toxin secretion, quorum sensing and chemotaxis, and no statin effect was observed on any of these phenotypes. Statins did not appear to have a characteristic effect on the P. aeruginosa transcriptome. However, a mutant library screen revealed that the effect of statins on P. aeruginosa biofilm was mediated through the PvrR regulator and the Cup fimbrial biosynthesis genes. Furthermore, proteomic analysis demonstrated that 6 proteins were reproducibly induced by simvastatin in the P. aeruginosa swarming cells. The effect of statins on the regulation of the host-P. aeruginosa immune response was also investigated. Statin treatment increased expression of the pro-inflammatory cytokines IL-8 and CCL20 in lung epithelial cells, but did not attenuate P. aeruginosa-mediated inflammatory gene induction. In fact, simvastatin and P. aeruginosa caused a synergistic effect on CCL20 expression. The expression of the transcriptional regulators KLF2 and KLF6 was also increased by statins and P. aeruginosa, with the induction of KLF6 by simvastatin proving to be a novel effect. Interestingly, both statins and P. aeruginosa were capable of inducing alternative splicing of KLF6. P. aeruginosa was found to induce KLF6 alternative splicing by way of the type 3 secreted toxin ExoS. In addition, a mechanistic role was elucidated for KLF6 in the lung, as it was determined that statin-mediated induction of this protein was responsible for the induction of the host response genes CCL20 and iNOS. Moreover, statin treatment caused a slight increase in infection-related cytotoxicity, and increased bacterial adhesion to cells. Taken together, these data demonstrate that statins can reduce the virulence of CFassociated bacterial pathogens and alter host response effectors. Furthermore, novel statin effectors were identified in both bacterial and host cells.

Neural Network Evidence for the Coupling of Presaccadic Visual Remapping to Predictive Eye Position Updating.

As we look around a scene, we perceive it as continuous and stable even though each saccadic eye movement changes the visual input to the retinas. How the brain achieves this perceptual stabilization is unknown, but a major hypothesis is that it relies on presaccadic remapping, a process in which neurons shift their visual sensitivity to a new location in the scene just before each saccade. This hypothesis is difficult to test in vivo because complete, selective inactivation of remapping is currently intractable. We tested it in silico with a hierarchical, sheet-based neural network model of the visual and oculomotor system. The model generated saccadic commands to move a video camera abruptly. Visual input from the camera and internal copies of the saccadic movement commands, or corollary discharge, converged at a map-level simulation of the frontal eye field (FEF), a primate brain area known to receive such inputs. FEF output was combined with eye position signals to yield a suitable coordinate frame for guiding arm movements of a robot. Our operational definition of perceptual stability was "useful stability," quantified as continuously accurate pointing to a visual object despite camera saccades. During training, the emergence of useful stability was correlated tightly with the emergence of presaccadic remapping in the FEF. Remapping depended on corollary discharge but its timing was synchronized to the updating of eye position. When coupled to predictive eye position signals, remapping served to stabilize the target representation for continuously accurate pointing. Graded inactivations of pathways in the model replicated, and helped to interpret, previous in vivo experiments. The results support the hypothesis that visual stability requires presaccadic remapping, provide explanations for the function and timing of remapping, and offer testable hypotheses for in vivo studies. We conclude that remapping allows for seamless coordinate frame transformations and quick actions despite visual afferent lags. With visual remapping in place for behavior, it may be exploited for perceptual continuity.

Interpreting Human Genetic Variation through Genomics and In Vivo Models

Improvements in genomic technology, both in the increased speed and reduced cost of sequencing, have expanded the appreciation of the abundance of human genetic variation. However the sheer amount of variation, as well as the varying type and genomic content of variation, poses a challenge in understanding the clinical consequence of a single mutation. This work uses several methodologies to interpret the observed variation in the human genome, and presents novel strategies for the prediction of allele pathogenicity.

Using the zebrafish model system as an in vivo assay of allele function, we identified a novel driver of Bardet-Biedl Syndrome (BBS) in CEP76. A combination of targeted sequencing of 785 cilia-associated genes in a cohort of BBS patients and subsequent in vivo functional assays recapitulating the human phenotype gave strong evidence for the role of CEP76 mutations in the pathology of an affected family. This portion of the work demonstrated the necessity of functional testing in validating disease-associated mutations, and added to the catalogue of known BBS disease genes.

Further study into the role of copy-number variations (CNVs) in a cohort of BBS patients showed the significant contribution of CNVs to disease pathology. Using high-density array comparative genomic hybridization (aCGH) we were able to identify pathogenic CNVs as small as several hundred bp. Dissection of constituent gene and in vivo experiments investigating epistatic interactions between affected genes allowed for an appreciation of several paradigms by which CNVs can contribute to disease. This study revealed that the contribution of CNVs to disease in BBS patients is much higher than previously expected, and demonstrated the necessity of consideration of CNV contribution in future (and retrospective) investigations of human genetic disease.

Finally, we used a combination of comparative genomics and in vivo complementation assays to identify second-site compensatory modification of pathogenic alleles. These pathogenic alleles, which are found compensated in other species (termed compensated pathogenic deviations [CPDs]), represent a significant fraction (from 3 – 10%) of human disease-associated alleles. In silico pathogenicity prediction algorithms, a valuable method of allele prioritization, often misrepresent these alleles as benign, leading to omission of possibly informative variants in studies of human genetic disease. We created a mathematical model that was able to predict CPDs and putative compensatory sites, and functionally showed in vivo that second-site mutation can mitigate the pathogenicity of disease alleles. Additionally, we made publically available an in silico module for the prediction of CPDs and modifier sites.

These studies have advanced the ability to interpret the pathogenicity of multiple types of human variation, as well as made available tools for others to do so as well.

Pro-inflammatory flagellin proteins of prevalent motile commensal bacteria are variably abundant in the intestinal microbiome of elderly humans

Some Eubacterium and Roseburia species are among the most prevalent motile bacteria present in the intestinal microbiota of healthy adults. These flagellate species contribute "cell motility" category genes to the intestinal microbiome and flagellin proteins to the intestinal proteome. We reviewed and revised the annotation of motility genes in the genomes of six Eubacterium and Roseburia species that occur in the human intestinal microbiota and examined their respective locus organization by comparative genomics. Motility gene order was generally conserved across these loci. Five of these species harbored multiple genes for predicted flagellins. Flagellin proteins were isolated from R. inulinivorans strain A2-194 and from E. rectale strains A1-86 and M104/1. The amino-termini sequences of the R. inulinivorans and E. rectale A1-86 proteins were almost identical. These protein preparations stimulated secretion of interleukin-8 (IL-8) from human intestinal epithelial cell lines, suggesting that these flagellins were pro-inflammatory. Flagellins from the other four species were predicted to be pro-inflammatory on the basis of alignment to the consensus sequence of pro-inflammatory flagellins from the beta- and gamma-proteobacteria. Many fliC genes were deduced to be under the control of sigma(28). The relative abundance of the target Eubacterium and Roseburia species varied across shotgun metagenomes from 27 elderly individuals. Genes involved in the flagellum biogenesis pathways of these species were variably abundant in these metagenomes, suggesting that the current depth of coverage used for metagenomic sequencing (3.13-4.79 Gb total sequence in our study) insufficiently captures the functional diversity of genomes present at low (<= 1%) relative abundance. E. rectale and R. inulinivorans thus appear to synthesize complex flagella composed of flagellin proteins that stimulate IL-8 production. A greater depth of sequencing, improved evenness of sequencing and improved metagenome assembly from short reads will be required to facilitate in silico analyses of complete complex biochemical pathways for low-abundance target species from shotgun metagenomes.