The genome of th17 cell-inducing segmented filamentous bacteria reveals extensive auxotrophy and adaptations to the intestinal environment

Perturbations of the composition of the symbiotic intestinal microbiota can have profound consequences for host metabolism and immunity. In mice, segmented filamentous bacteria (SFB) direct the accumulation of potentially proinflammatory Th17 cells in the intestinal lamina propria. We present the genome sequence of SFB isolated from monocolonized mice, which classifies SFB phylogenetically as a unique member of Clostridiales with a highly reduced genome. Annotation analysis demonstrates that SFB depend on their environment for amino acids and essential nutrients and may utilize host and dietary glycans for carbon, nitrogen, and energy. Comparative analyses reveal that SFB are functionally related to members of the genus Clostridium and several pathogenic or commensal "minimal" genera, including Finegoldia, Mycoplasma, Borrelia, and Phytoplasma. However, SFB are functionally distinct from all 1200 examined genomes, indicating a gene complement representing biology relatively unique to their role as a gut commensal closely tied to host metabolism and immunity.     

Different degree of ileal colonization by segmented, filamentous bacteria in two strains of mice.

Segmented, filamentous bacteria (SFBs) are autochthonous, apathogenic inhabitants of the ileum of various animal species. Outbred Swiss (Cpb:SE) mice have significantly higher degrees of SFB colonization than do inbred BALB/c mice. The present studies were carried out to identify determinants of this strain difference. In a cross-fostering experiment it was shown that SFB colonization of the pups is determined by the strain of the pups themselves rather than by the strain of the nursing dam. Thus, maternal effects may not be involved in SFB colonization. In a cross-infecting experiment using germ-free and SFB-positive animals of the two mouse strains, it was found that ileal SFB colonization is determined by host characteristics rather than by origin of the SFBs. Thus, SFBs that are specific for a given mouse strain may not exist in the two strains of mice. It is concluded that the mouse strain difference in SFB colonization is determined by host characteristics, which probably have a genetic basis.     

Mapping of 195 genes in cattle and updated comparative map with man, mouse, rat and pig

Our on-going goal is to improve and update the comparative genome organization between cattle and man but also among the most detailed mammalian species genomes i.e. cattle, mouse, rat and pig. In this work, we localized 195 genes in cattle and checked all human/bovine non-concordant localizations found in the literature. Next, we compiled all the genes mapped in cattle, goat, sheep and pig (2,166) for which the human ortholog with its chromosomal position is known, added corresponding data in mouse and rat, and ordered the genes relatively to the human genome sequence. We estimate that our compilation provides bovine mapping information for about 89% of the human autosomes. Thus, a near complete, overall and detailed picture of the number, distribution and extent of bovine conserved syntenies (regardless of gene order) on human R-banded autosomes is proposed as well as a comparison with mouse, rat and pig genomes.     

Comparison of the sequence organization of related retrovirus-like multigene families in three evolutionarily distant rodent genomes.

Sequences related to mouse intracisternal A-particle (IAP) genes have been isolated from rat and Syrian hamster gene libraries as recombinants in lambda phage. The sequences are moderately reiterated in both these genomes but their sequence organization in the hamster genome is different from that in the rat genome. Restriction analysis and electron microscopy indicate that the Syrian hamster IAP sequences represent a family of relatively homogeneous well-conserved units; in this, they resemble the mouse IAP genes. The rat sequences, in contrast, are heterogeneous. Both the hamster and rat IAP sequences contain regions homologous to mouse IAP genes interspersed with regions of apparent non-homology. The interspersed regions range in size from 0.5-1.0 kilobases (Kb). The regions of homology among the mouse, rat and Syrian hamster IAP sequences have been mapped to a 5-6 Kb internal region on the mouse IAP genes. Mouse IAP long terminal repeat (LTR) sequences were not detected in the rat and Syrian hamster genomes. We used the thermal stability of hybrids between cloned and genomic IAP sequences to measure family homogeneity. Mouse and Syrian hamster IAP sequences are homogeneous by this criterion, but the rat IAP sequences are heterogeneous with a Tm 6 degrees C below the self-hybrid. The contrasting organization of IAP-related elements in the genomes of these rodents indicates that amplification or homogenization of this sequence family has occurred independently and at different periods of time during their evolution.     

Genomics of clostridial pathogens: implication of extrachromosomal elements in pathogenicity

The recently decoded genomes of the major clostridial toxin-producing pathogens Clostridium perfringens, Clostridium tetani, Clostridium botulinum and Clostridium difficile have provided a huge amount of new sequence data. Recent studies have focused on the identification and investigation of pathogenic determinants and the regulatory events governing their expression. The sequence data revealed also the genomic background of virulence genes, as well as the contribution of extrachromosomal elements to a pathogenic phenotype. This has generated new insights in clostridial pathogenesis - and will continue to do so in the future - and has deepened our understanding of the anaerobic lifestyle of clostridial species.     

Protein coding potential of retroviruses and other transposable elements in vertebrate genomes

We suggest an annotation strategy for genes encoded by retroviruses and transposable elements (RETRA genes) based on a set of marker protein domains. Usually RETRA genes are masked in vertebrate genomes prior to the application of automated gene prediction pipelines under the assumption that they provide no selective advantage to the host. Yet, we show that about 1000 genes in four vertebrate gene sets analyzed contain at least one RETRA gene marker domain. Using the conservation of genomic neighborhood (synteny), we were able to discriminate between RETRA genes with putative functionality in the vertebrates and those that probably function only in the context of mobile elements. We identified 35 such genes in human, along with their corresponding mouse and rat orthologs; which included almost all known human genes with similarity to mobile elements. The results also imply that the vast majority of the remaining RETRA genes in current gene sets are unlikely to encode vertebrate functions. To automatically annotate RETRA genes in other vertebrate genomes, we provide as a tool a set of marker protein domains and a manually refined list of domesticated or ancestral RETRA genes for rescuing genes with vertebrate functions.     

Comprehensive discovery of CRISPR-targeted terminally redundant sequences in the human gut metagenome: Viruses,plasmids, and more

Viruses are the most numerous biological entity, existing in all environments and infecting all cellular organisms. Compared with cellular life, the evolution and origin of viruses are poorly understood; viruses are enormously diverse, and most lack sequence similarity to cellular genes. To uncover viral sequences without relying on either reference viral sequences from databases or marker genes that characterize specific viral taxa, we developed an analysis pipeline for virus inference based on clustered regularly interspaced short palindromic repeats (CRISPR). CRISPR is a prokaryotic nucleic acid restriction system that stores the memory of previous exposure. Our protocol can infer CRISPR-targeted sequences, including viruses, plasmids, and previously uncharacterized elements, and predict their hosts using unassembled short-read metagenomic sequencing data. By analyzing human gut metagenomic data, we extracted 11,391 terminally redundant CRISPR-targeted sequences, which are likely complete circular genomes. The sequences included 2,154 tailed-phage genomes, together with 257 complete crAssphage genomes, 11 genomes larger than 200 kilobases, 766 genomes of Microviridae species, 56 genomes of Inoviridae species, and 95 previously uncharacterized circular small genomes that have no reliably predicted protein-coding gene. We predicted the host(s) of approximately 70% of the discovered genomes at the taxonomic level of phylum by linking protospacers to taxonomically assigned CRISPR direct repeats. These results demonstrate that our protocol is efficient for de novo inference of CRISPR-targeted sequences and their host prediction.     

The family of toxin-related ecto-ADP-ribosyltransferases in humans and the mouse

ADP-ribosyltransferases including toxins secreted by Vibrio cholera, Pseudomonas aerurginosa, and other pathogenic bacteria inactivate the function of human target proteins by attaching ADP-ribose onto a critical amino acid residue. Cross-species polymerase chain reaction (PCR) and database mining identified the orthologs of these ADP-ribosylating toxins in humans and the mouse. The human genome contains four functional toxin-related ADP-ribosyltransferase genes (ARTs) and two related intron-containing pseudogenes; the mouse has six functional orthologs. The human and mouse ART genes map to chromosomal regions with conserved linkage synteny. The individual ART genes reveal highly restricted expression patterns, which are largely conserved in humans and the mouse. We confirmed the predicted extracellular location of the ART proteins by expressing recombinant ARTs in insect cells. Two human and four mouse ARTs contain the active site motif (R-S-EXE) typical of arginine-specific ADP-ribosyltransferases and exhibit the predicted enzyme activities. Two other human ARTs and their murine orthologues deviate in the active site motif and lack detectable enzyme activity. Conceivably, these ARTs may have acquired a new specificity or function. The position-sensitive iterative database search program PSI-BLAST connected the mammalian ARTs with most known bacterial ADP-ribosylating toxins. In contrast, no related open reading frames occur in the four completed genomes of lower eucaryotes (yeast, worm, fly, and mustard weed). Interestingly, these organisms also lack genes for ADP-ribosylhydrolases, the enzymes that reverse protein ADP-ribosylation. This suggests that the two enzyme families that catalyze reversible mono-ADP-ribosylation either were lost from the genomes of these nonchordata eucaryotes or were subject to horizontal gene transfer between kingdoms.     

关于黄曲条跳甲的寄主范围

通过非选择性试验方法,测定黄曲条跳甲Ｐｈｙｌｌｏｔｒｅｔａ ｓｔｒｉｏｌａｔａ（Ｆａｂｒｕｃｉｕｓ）在不同植物上的存活率,及番茄、茄子枝叶提取物对曲条跳甲的拒食作用结果表明,在供试的中,仅十字花科植物能使之正常存活,其他科的植物均全部或大量死亡,不可能作为黄曲条跳甲的寄主,甚至茄科的番茄、茄子提取物对黄曲条跳甲在显著的拒食作用,拒食率高达３０．９％和６０．５２％。从化学生态学的角度讨论了典曲条跳甲     

Comprehensive analysis of keratin gene clusters in humans and rodents

Here, we present the comparative analysis of the two keratin (K) gene clusters in the genomes of man, mouse and rat. Overall, there is a remarkable but not perfect synteny among the clusters of the three mammalian species. The human type I keratin gene cluster consists of 27 genes and 4 pseudogenes, all in the same orientation. It is interrupted by a domain of multiple genes encoding keratin-associated proteins (KAPs). Cytokeratin, hair and inner root sheath keratin genes are grouped together in small subclusters, indicating that evolution occurred by duplication events. At the end of the rodent type I gene cluster, a novel gene related to K14 and K17 was identified, which is converted to a pseudogene in humans. The human type II cluster consists of 27 genes and 5 pseudogenes, most of which are arranged in the same orientation. Of the 26 type II murine keratin genes now known, the expression of two new genes was identified by RT-PCR. Kb20, the first gene in the cluster, was detected in lung tissue. Kb39, a new ortholog of K1, is expressed in certain stratified epithelia. It represents a candidate gene for those hyperkeratotic skin syndromes in which no K1 mutations were identified so far. Most remarkably, the human K3 gene which causes Meesmann's corneal dystrophy when mutated, lacks a counterpart in the mouse genome. While the human genome has 138 pseudogenes related to K8 and K18, the mouse and rat genomes contain only 4 and 6 such pseudogenes. Our results also provide the basis for a unified keratin nomenclature and for future functional studies.     

Comparative sequence analysis of a single-gene conserved segment in mouse and human

Comparative mapping and sequencing of the mouse and human genomes have defined large, conserved chromosomal segments in which gene content and order are highly conserved. These regions span megabase-sized intervals and together comprise the vast majority of both genomes. However, the evolutionary relationships among the small remaining portions of these genomes are not as well characterized. Here we describe the sequencing and annotation of a 341-kb region of mouse Chr 2 containing nine genes, including biliverdin reductase A (Blvra), and its comparison with the orthologous regions of the human and rat genomes. These analyses reveal that the known conserved synteny between mouse Chromosome (Chr) 2 and human Chr 7 reflects an interval containing one gene (Blvra/BLVRA) that is, at most, just 34 kb in the mouse genome. In the mouse, this segment is flanked proximally by genes orthologous to human chromosome 15q21 and distally by genes orthologous to human Chr 2q11. The observed differences between the human and mouse genomes likely resulted from one or more rearrangements in the rodent lineage. In addition to the resulting changes in gene order and location, these rearrangements also appear to have included genomic deletions that led to the loss of at least one gene in the rodent lineage. Finally, we also have identified a recent mouse-specific segmental duplication. These finding illustrate that small genomic regions outside the large mouse–human conserved segments can contain a single gene as well as sequences that are apparently unique to one genome. The nucleotide sequence data reported in this paper have been submitted to GenBank and assigned the accession numbers AC074224 and AC074041.     

Comparative genomics of odorant and pheromone receptor genes in rodents

We applied a comprehensive data-mining strategy to examine the repertoires of rat and mouse odorant receptors (ORs) and type 1 pheromone receptors (V1Rs) using the mm5 (mouse) and rn3 (rat) genomes. We identified 1576 rat OR genes, including 292 pseudogenes. The rat V1R repertoire is composed of 115 intact genes and 72 pseudogenes. The mouse OR and V1R databases were updated using the new assembly mm5, from which 1375 mouse ORs and 308 V1Rs were identified, with more than 100 putative pseudogenes from mm2 now identified as intact because of the higher sequence quality. With these new data we have conducted a series of genomic analyses of the OR and V1R genes from mouse and rat. Orthologous OR clusters were identified in mouse and rat and comparison analysis was performed at three incremental levels: families, coding sequences, and motifs. At the family level, we found that V1R genes have more species-specific families than OR genes. About 20% of intact V1R genes have no orthologous counterpart in the same family, whereas less than 1% of intact ORs are similarly isolated. At the coding sequence level, OR genes are more conserved between mouse and rat than V1R genes. OR genes share greater similarity with their orthologous counterparts than with their closest neighbor, whereas V1R genes show the opposite tendency. Motifs were identified to obtain biological insights. Motifs specific for species or families were found in OR and V1R genes, which may result in the differential pheromone-dependent behaviors and perception of odors between mouse and rat.     

A genome-wide survey on basic helix-loop-helix transcription factors in rat and mouse

The basic helix-loop-helix (bHLH) proteins play essential roles in a wide range of developmental processes in higher organisms. bHLH family members have been identified in over 20 organisms, including nematode, fruit fly, and human. Our study identified 114 rat and 14 additional mouse bHLH members in rat and mouse genomes, respectively. Phylogenetic analyses revealed that both rat and mouse had 49, 26, 15, 4, 12, and 4 bHLH members in groups A, B, C, D, E, and F, respectively. Only the rat Mxi1 gene has two copies in the genome. All other rat bHLH genes and all mouse bHLH genes are single-copy genes. The chromosomal distribution pattern of mouse, rat, and human bHLH genes suggests the emergence of some bHLH genes through gene duplication, which probably happened at least before the divergence of vertebrates from invertebrates. The present study provides useful information for future studies using rat as a model animal for mammalian development. X. Zheng and Y. Wang are jointly first authors. An erratum to this article can be found at     

The Mammalian 2′-5′ Oligoadenylate Synthetase Gene Family: Evidence for Concerted Evolution of Paralogous Oas1 Genes in Rodentia and Artiodactyla

Multiple 2′-5′ oligoadenylate (2-5A) synthetases are important components of innate immunity in mammals. Gene families encoding these proteins have previously been studied mainly in humans and mice. To reconstruct the evolution of this gene family in mammals, a search for additional 2-5A synthetase genes was performed in rat, cattle, pig, and dog. Twelve 2′-5′ oligoadenylate synthetase (Oas) genes were identified in the rat genome, including eight Oas1 genes, two Oas1 pseudogenes, single copies of Oas2 and Oas3, and two Oas-like genes, Oasl1 and Oasl2. Four OAS genes were detected in the pig genome and five OAS genes were found in both the cattle and dog genomes. An OAS3 gene was not found in either the cattle or the pig genome. While two tandemly duplicated OAS-like (OASL) genes were identified in the dog genome, only a single OASL orthologue was found in both the cattle and the pig genomes. The bovine and porcine OASL genes contain premature stop codons and encode truncated proteins, which lack the typical C-terminal double ubiquitin domains. The cDNA sequences of the rat, cattle, pig, and dog OAS genes were amplified, sequenced and compared with each other and with those in the human, mouse, horse, and chicken genomes. Evidence of concerted evolution of paralogous 2′-5′ oligoadenylate synthetase 1 genes was obtained in rodents (Rodentia) and even-toed ungulates (Artiodactyla). Calculations using the nonparametric Kolmogorov-Smirnov test suggested that the homogenization of paralogous OAS1 sequences was due to gene conversion rather than stabilizing selection. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. Martin Kreitman     

Activation of the SCPx promoter in mouse adrenocortical Y1 cells

Sterol carrier protein X (SCPx) is a peroxisomal protein with both lipid transfer and thiolase activity. Treatment of mouse adrenal Y1 cells with cAMP for 24h caused a significant induction of SCPx mRNA levels. Reporter gene studies demonstrated that treatment with cAMP and SF-1 was able to activate the SCPx promoter. Sequence analysis revealed the presence of three putative steroidogenic factor-1 (SF-1) binding motifs (designated SFB1, SFB2, and SFB3) and one CRE. Only SFB1 and SFB3 were able to bind recombinant SF-1 protein in electrophoretic mobility shift assays. The CRE was able to form a DNA/protein complex in the presence of Y1 nuclear extracts. Mutational analysis studies demonstrated that SFB3 is required for full activation of the SCPx promoter by cAMP treatment. Regulation of the SCPx gene by SF-1 and cAMP is similar to the regulatory mechanisms observed for other steroidogenic genes.     

Clustering of main orthologs for multiple genomes

The identification of orthologous genes shared by multiple genomes is critical for both functional and evolutionary studies in comparative genomics. While it is usually done by sequence similarity search and reconciled tree construction in practice, recently a new combinatorial approach and high-throughput system MSOAR for ortholog identification between closely related genomes based on genome rearrangement and gene duplication has been proposed in Fu et al. MSOAR assumes that orthologous genes correspond to each other in the most parsimonious evolutionary scenario, minimizing the number of genome rearrangement and (postspeciation) gene duplication events. However, the parsimony approach used by MSOAR limits it to pairwise genome comparisons. In this paper, we extend MSOAR to multiple (closely related) genomes and propose an ortholog clustering method, called MultiMSOAR, to infer main orthologs in multiple genomes. As a preliminary experiment, we apply MultiMSOAR to rat, mouse, and human genomes, and validate our results using gene annotations and gene function classifications in the public databases. We further compare our results to the ortholog clusters predicted by MultiParanoid, which is an extension of the well-known program InParanoid for pairwise genome comparisons. The comparison reveals that MultiMSOAR gives more detailed and accurate orthology information, since it can effectively distinguish main orthologs from inparalogs.     

Long-branch attraction phenomenon and the impact of among-site rate variation on rodent phylogeny

The phylogenetic relationships among major lineages of rodents is one of the issues most debated by both paleontologists and molecular biologists. In the present study, we have analyzed all complete mammalian mitochondrial genomes available in the databases, including five rodent species (rat, mouse, dormouse, squirrel and guinea-pig). Phylogenetic analyses were performed on H-strand amino acid sequences by means of maximum-likelihood and on H-strand protein-coding and ribosomal genes by means of distance methods. Also, log-likelihood ratio tests were applied to different tree topologies under the assumption of rodent monophyly, paraphyly or polyphyly. The analyses significantly rejected rodent monophyly and showed the existence of two differentiated clades, one containing non-murids (dormouse, squirrel and guinea-pig) and the other containing murids (rat and mouse). Long-branch attraction between murids and the outgroups could not be responsible for the existence of two different rodent clades, as no significant differences in evolutionary rate have been observed, except in the case of the squirrel, which shows a lower rate. The impact of among-site rate variation models on the phylogeny of rodents has been evaluated using the gamma distribution model. Results have shown that relationships among rodents remained unchanged, and the general topology of the tree was not affected, even though some branches were not properly resolved, most likely due to a lack of fit between estimated and real rate heterogeneity parameters.