PATRIC: the comprehensive bacterial bioinformatics resource with a focus on human pathogenic species.

Funded by the National Institute of Allergy and Infectious Diseases, the Pathosystems Resource Integration Center (PATRIC) is a genomics-centric relational database and bioinformatics resource designed to assist scientists in infectious-disease research. Specifically, PATRIC provides scientists with (i) a comprehensive bacterial genomics database, (ii) a plethora of associated data relevant to genomic analysis, and (iii) an extensive suite of computational tools and platforms for bioinformatics analysis. While the primary aim of PATRIC is to advance the knowledge underlying the biology of human pathogens, all publicly available genome-scale data for bacteria are compiled and continually updated, thereby enabling comparative analyses to reveal the basis for differences between infectious free-living and commensal species. Herein we summarize the major features available at PATRIC, dividing the resources into two major categories: (i) organisms, genomes, and comparative genomics and (ii) recurrent integration of community-derived associated data. Additionally, we present two experimental designs typical of bacterial genomics research and report on the execution of both projects using only PATRIC data and tools. These applications encompass a broad range of the data and analysis tools available, illustrating practical uses of PATRIC for the biologist. Finally, a summary of PATRIC's outreach activities, collaborative endeavors, and future research directions is provided.

Phylogeny of gammaproteobacteria.

The phylogeny of the large bacterial class Gammaproteobacteria has been difficult to resolve. Here we apply a telescoping multiprotein approach to the problem for 104 diverse gammaproteobacterial genomes, based on a set of 356 protein families for the whole class and even larger sets for each of four cohesive subregions of the tree. Although the deepest divergences were resistant to full resolution, some surprising patterns were strongly supported. A representative of the Acidithiobacillales routinely appeared among the outgroup members, suggesting that in conflict with rRNA-based phylogenies this order does not belong to Gammaproteobacteria; instead, it (and, independently, "Mariprofundus") diverged after the establishment of the Alphaproteobacteria yet before the betaproteobacteria/gammaproteobacteria split. None of the orders Alteromonadales, Pseudomonadales, or Oceanospirillales were monophyletic; we obtained strong support for clades that contain some but exclude other members of all three orders. Extreme amino acid bias in the highly A+T-rich genome of Candidatus Carsonella prevented its reliable placement within Gammaproteobacteria, and high bias caused artifacts that limited the resolution of the relationships of other insect endosymbionts, which appear to have had multiple origins, although the unbiased genome of the endosymbiont Sodalis acted as an attractor for them. Instability was observed for the root of the Enterobacteriales, with nearly equal subsets of the protein families favoring one or the other of two alternative root positions; the nematode symbiont Photorhabdus was identified as a disruptor whose omission helped stabilize the Enterobacteriales root.

Analysis of ten Brucella genomes reveals evidence for horizontal gene transfer despite a preferred intracellular lifestyle.

The facultative intracellular bacterial pathogen Brucella infects a wide range of warm-blooded land and marine vertebrates and causes brucellosis. Currently, there are nine recognized Brucella species based on host preferences and phenotypic differences. The availability of 10 different genomes consisting of two chromosomes and representing six of the species allowed for a detailed comparison among themselves and relatives in the order Rhizobiales. Phylogenomic analysis of ortholog families shows limited divergence but distinct radiations, producing four clades as follows: Brucella abortus-Brucella melitensis, Brucella suis-Brucella canis, Brucella ovis, and Brucella ceti. In addition, Brucella phylogeny does not appear to reflect the phylogeny of Brucella species' preferred hosts. About 4.6% of protein-coding genes seem to be pseudogenes, which is a relatively large fraction. Only B. suis 1330 appears to have an intact beta-ketoadipate pathway, responsible for utilization of plant-derived compounds. In contrast, this pathway in the other species is highly pseudogenized and consistent with the "domino theory" of gene death. There are distinct shared anomalous regions (SARs) found in both chromosomes as the result of horizontal gene transfer unique to Brucella and not shared with its closest relative Ochrobactrum, a soil bacterium, suggesting their acquisition occurred in spite of a predominantly intracellular lifestyle. In particular, SAR 2-5 appears to have been acquired by Brucella after it became intracellular. The SARs contain many genes, including those involved in O-polysaccharide synthesis and type IV secretion, which if mutated or absent significantly affect the ability of Brucella to survive intracellularly in the infected host.

Rickettsia phylogenomics: unwinding the intricacies of obligate intracellular life.

BACKGROUND: Completed genome sequences are rapidly increasing for Rickettsia, obligate intracellular alpha-proteobacteria responsible for various human diseases, including epidemic typhus and Rocky Mountain spotted fever. In light of phylogeny, the establishment of orthologous groups (OGs) of open reading frames (ORFs) will distinguish the core rickettsial genes and other group specific genes (class 1 OGs or C1OGs) from those distributed indiscriminately throughout the rickettsial tree (class 2 OG or C2OGs). METHODOLOGY/PRINCIPAL FINDINGS: We present 1823 representative (no gene duplications) and 259 non-representative (at least one gene duplication) rickettsial OGs. While the highly reductive (approximately 1.2 MB) Rickettsia genomes range in predicted ORFs from 872 to 1512, a core of 752 OGs was identified, depicting the essential Rickettsia genes. Unsurprisingly, this core lacks many metabolic genes, reflecting the dependence on host resources for growth and survival. Additionally, we bolster our recent reclassification of Rickettsia by identifying OGs that define the AG (ancestral group), TG (typhus group), TRG (transitional group), and SFG (spotted fever group) rickettsiae. OGs for insect-associated species, tick-associated species and species that harbor plasmids were also predicted. Through superimposition of all OGs over robust phylogeny estimation, we discern between C1OGs and C2OGs, the latter depicting genes either decaying from the conserved C1OGs or acquired laterally. Finally, scrutiny of non-representative OGs revealed high levels of split genes versus gene duplications, with both phenomena confounding gene orthology assignment. Interestingly, non-representative OGs, as well as OGs comprised of several gene families typically involved in microbial pathogenicity and/or the acquisition of virulence factors, fall predominantly within C2OG distributions. CONCLUSION/SIGNIFICANCE: Collectively, we determined the relative conservation and distribution of 14354 predicted ORFs from 10 rickettsial genomes across robust phylogeny estimation. The data, available at PATRIC (PathoSystems Resource Integration Center), provide novel information for unwinding the intricacies associated with Rickettsia pathogenesis, expanding the range of potential diagnostic, vaccine and therapeutic targets.

The human obesity gene map: the 2004 update.

This paper presents the eleventh update of the human obesity gene map, which incorporates published results up to the end of October 2004. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, transgenic and knockout murine models relevant to obesity, quantitative trait loci (QTLs) from animal cross-breeding experiments, association studies with candidate genes, and linkages from genome scans is reviewed. As of October 2004, 173 human obesity cases due to single-gene mutations in 10 different genes have been reported, and 49 loci related to Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. There are 166 genes which, when mutated or expressed as transgenes in the mouse, result in phenotypes that affect body weight and adiposity. The number of QTLs reported from animal models currently reaches 221. The number of human obesity QTLs derived from genome scans continues to grow, and we have now 204 QTLs for obesity-related phenotypes from 50 genome-wide scans. A total of 38 genomic regions harbor QTLs replicated among two to four studies. The number of studies reporting associations between DNA sequence variation in specific genes and obesity phenotypes has also increased considerably with 358 findings of positive associations with 113 candidate genes. Among them, 18 genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, >600 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful publications and genomic and other relevant sites can be found at http://obesitygene.pbrc.edu.

The human obesity gene map: the 2003 update.

This is the tenth update of the human obesity gene map, incorporating published results up to the end of October 2003 and continuing the previous format. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) from human genome-wide scans and animal crossbreeding experiments, and association and linkage studies with candidate genes and other markers is reviewed. Transgenic and knockout murine models relevant to obesity are also incorporated (N = 55). As of October 2003, 41 Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. QTLs reported from animal models currently number 183. There are 208 human QTLs for obesity phenotypes from genome-wide scans and candidate regions in targeted studies. A total of 35 genomic regions harbor QTLs replicated among two to five studies. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 272 studies reporting positive associations with 90 candidate genes. Fifteen such candidate genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, more than 430 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful sites can be found at http://obesitygene.pbrc.edu.

Serial analysis of gene expression in the skeletal muscle of endurance athletes compared to sedentary men.

Physical exercise produces several adaptive changes in skeletal muscle. However, the molecular mechanisms of these effects are poorly understood. We performed serial analysis of gene expression (SAGE) to quantify the global gene expression profile in sedentary and endurance-trained muscle. A total of 10869 SAGE tags was sequenced and represented 4727 genes. The genes most expressed in muscle are mainly involved in contraction and energy metabolism. Thirty-three genes were differentially expressed between endurance athletes and sedentary individuals. Four genes such as myosin binding protein C fast-type, glycogen phosphorylase, and pyruvate kinase were expressed less in endurance athletes, whereas eight genes coding for expressed sequence tag similar to (EST) crystallin alpha B, EST myosin light chain 2, EST surfactant pulmonary-associated protein A1, EST thrombospondin, EST fructose-bisphosphate aldolase A, EST cytochrome oxidase 1, NADH dehydrogenase 3, and G8 protein were up-regulated. Most of the up-regulated tags corresponded to novel genes. On the other hand, different isoforms of fructose-bisphosphate aldolase A were also differentially expressed. The current study underlying the most highly expressed genes allows a better understanding of global muscle characteristics in normal and endurance-trained individuals. Moreover, the current data suggest novel candidate genes that may be responsible for enhanced endurance performance.

The human obesity gene map: the 2002 update.

This is the ninth update of the human obesity gene map, incorporating published results through October 2002 and continuing the previous format. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) from human genome-wide scans and various animal crossbreeding experiments, and association and linkage studies with candidate genes and other markers is reviewed. For the first time, transgenic and knockout murine models exhibiting obesity as a phenotype are incorporated (N = 38). As of October 2002, 33 Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and the causal genes or strong candidates have been identified for 23 of these syndromes. QTLs reported from animal models currently number 168; there are 68 human QTLs for obesity phenotypes from genome-wide scans. Additionally, significant linkage peaks with candidate genes have been identified in targeted studies. Seven genomic regions harbor QTLs replicated among two to five studies. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 222 studies reporting positive associations with 71 candidate genes. Fifteen such candidate genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. More than 300 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful sites can be found at http://obesitygene.pbrc.edu.

Melanocortin 4 receptor sequence variations are seldom a cause of human obesity: the Swedish Obese Subjects, the HERITAGE Family Study, and a Memphis cohort.

The prevalence of mutations within and in the flanking regions of the gene encoding the melanocortin 4 receptor was investigated in severely obese and normal-weight subjects from the Swedish Obese Subjects study, the Health, Risk Factors, Exercise Training, and Genetics (HERITAGE) Family study, and a Memphis cohort. A total of 433 white and 95 black subjects (94% females) were screened for mutations by direct sequencing. Three previously described missense variants and nine novel (three missense, six silent) variants were detected. None of them showed significant association with obesity or related phenotypes. In addition, two novel deletions were found in two heterozygous obese women: a -65_-64delTG mutation within the 5' noncoding region and a 171delC frameshift mutation predicted to result in a truncated nonfunctional receptor. No pathogenic mutations were found among obese blacks or nonobese controls. Furthermore, none of the null mutations found in other populations was present in this sample. In conclusion, our results do not support the prevailing notion that sequence variation in the melanocortin 4 receptor gene is a frequent cause of human obesity.

The human obesity gene map: the 2001 update.

This report constitutes the eighth update of the human obesity gene map, incorporating published results up to the end of October 2001. Evidence from the rodent and human obesity cases caused by single-gene mutations, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) uncovered in human genome-wide scans and in crossbreeding experiments in various animal models, association and linkage studies with candidate genes and other markers is reviewed. The human cases of obesity related in some way to single-gene mutations in six different genes are incorporated. Twenty-five Mendelian disorders exhibiting obesity as one of their clinical manifestations have now been mapped. The number of different QTLs reported from animal models currently reaches 165. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 174 studies reporting positive associations with 58 candidate genes. Finally, 59 loci have been linked to obesity indicators in genomic scans and other linkage study designs. The obesity gene map depicted in Figure 1 reveals that putative loci affecting obesity-related phenotypes can be found on all chromosomes except chromosome Y. A total of 54 new loci have been added to the map in the past 12 months, and the number of genes, markers, and chromosomal regions that have been associated or linked with human obesity phenotypes is now above 250. Likewise, the number of negative studies, which are only partially reviewed here, is also on the rise.