Genomic signatures of cold adaptation in the family Colwelliaceae.

Psychrophily is a phenotype describing microbial growth at low temperatures; elucidating the biomolecular and genomic adaptations necessary for survival in the cold is important for understanding life in extreme environments on Earth and in outer space. We used comparative genomics and temperature growth experiments of bacteria from the family Colwelliaceae to identify genomic factors correlated with optimal growth temperature (OGT). A phylogenomic analysis of 67 public and 39 newly sequenced strains revealed three main clades of Colwelliaceae. Temperature growth experiments revealed significant differences in mean OGT by clade, wherein strains of Colwelliaceae had similar growth rates at -1 °C but varied in their ability to tolerate 17 °C. Using amino acid compositional indices, a multiple linear regression model was constructed to predict the OGT of these organisms (RMSE 5.2 °C). Investigation of Colwelliaceae functional genes revealed a putative cold-adaptive gene cassette that was present in psychrophilic strains but absent in a closely related strain with a significantly higher OGT. This study also presents genomic evidence suggesting that the clade of Colwelliaceae containing Colwellia hornerae should be investigated as a new genus. These contributions offer key insights into the psychrophily phenotype and its underlying genomic foundation in the family Colwelliaceae.

Identification and Characterization of a Dominant Sulfolane-Degrading Rhodoferax sp. via Stable Isotope Probing Combined with Metagenomics.

Sulfolane is an industrial solvent and emerging organic contaminant affecting groundwater around the world, but little is known about microbes capable of biodegrading sulfolane or the pathways involved. We combined DNA-based stable isotope probing (SIP) with genome-resolved metagenomics to identify microorganisms associated with sulfolane biodegradation in a contaminated subarctic aquifer. In addition to 16S rRNA gene amplicon sequencing, we performed shotgun metagenomics on the 13C-labeled DNA to obtain functional and taxonomic information about the active sulfolane-degrading community. We identified the primary sulfolane degrader, comprising ~85% of the labeled community in the amplicon sequencing dataset, as closely related to Rhodoferax ferrireducens strain T118. We obtained a 99.8%-complete metagenome-assembled genome for this strain, allowing us to identify putative pathways of sulfolane biodegradation. Although the 4S dibenzothiophene desulfurization pathway has been proposed as an analog for sulfolane biodegradation, we found only a subset of the required genes, suggesting a novel pathway specific to sulfolane. DszA, the enzyme likely responsible for opening the sulfolane ring structure, was encoded on both the chromosome and a plasmid. This study demonstrates the power of integrating DNA-SIP with metagenomics to characterize emerging organic contaminant degraders without culture bias and expands the known taxonomic distribution of sulfolane biodegradation.

An inter-order horizontal gene transfer event enables the catabolism of compatible solutes by Colwellia psychrerythraea 34H.

Colwellia is a genus of mostly psychrophilic halophilic Gammaproteobacteria frequently isolated from polar marine sediments and sea ice. In exploring the capacity of Colwellia psychrerythraea 34H to survive and grow in the liquid brines of sea ice, we detected a duplicated 37 kbp genomic island in its genome based on the abnormally high G + C content. This island contains an operon encoding for heterotetrameric sarcosine oxidase and is located adjacent to several genes used in the serial demethylation of glycine betaine, a compatible solute commonly used for osmoregulation, to dimethylglycine, sarcosine, and glycine. Molecular clock inferences of important events in the adaptation of C. psychrerythraea 34H to compatible solute utilization reflect the geological evolution of the polar regions. Validating genomic predictions, C. psychrerythraea 34H was shown to grow on defined media containing either choline or glycine betaine, and on a medium with sarcosine as the sole organic source of carbon and nitrogen. Growth by 8 of 9 tested Colwellia species on a newly developed sarcosine-based defined medium suggested that the ability to catabolize glycine betaine (the catabolic precursor of sarcosine) is likely widespread in the genus Colwellia. This capacity likely provides a selective advantage to Colwellia species in cold, salty environments like sea ice, and may have contributed to the ability of Colwellia to invade these extreme niches.

Testing the infinitely many genes model for the evolution of the bacterial core genome and pangenome.

When groups of related bacterial genomes are compared, the number of core genes found in all genomes is usually much less than the mean genome size, whereas the size of the pangenome (the set of genes found on at least one of the genomes) is much larger than the mean size of one genome. We analyze 172 complete genomes of Bacilli and compare the properties of the pangenomes and core genomes of monophyletic subsets taken from this group. We then assess the capabilities of several evolutionary models to predict these properties. The infinitely many genes (IMG) model is based on the assumption that each new gene can arise only once. The predictions of the model depend on the shape of the evolutionary tree that underlies the divergence of the genomes. We calculate results for coalescent trees, star trees, and arbitrary phylogenetic trees of predefined fixed branch length. On a star tree, the pangenome size increases linearly with the number of genomes, as has been suggested in some previous studies, whereas on a coalescent tree, it increases logarithmically. The coalescent tree gives a better fit to the data, for all the examples we consider. In some cases, a fixed phylogenetic tree proved better than the coalescent tree at reproducing structure in the gene frequency spectrum, but little improvement was gained in predictions of the core and pangenome sizes. Most of the data are well explained by a model with three classes of gene: an essential class that is found in all genomes, a slow class whose rate of origination and deletion is slow compared with the time of divergence of the genomes, and a fast class showing rapid origination and deletion. Although the majority of genes originating in a genome are in the fast class, these genes are not retained for long periods, and the majority of genes present in a genome are in the slow or essential classes. In general, we show that the IMG model is useful for comparison with experimental genome data both for species level and widely divergent taxonomic groups. Software implementing the described formulae is provided at http://github.com/rec3141/pangenome.

Origin and evolution of gene families in Bacteria and Archaea.

BACKGROUND: Comparison of complete genomes of Bacteria and Archaea shows that gene content varies considerably and that genomes evolve quite rapidly via gene duplication and deletion and horizontal gene transfer. We analyze a diverse set of 92 Bacteria and 79 Archaea in order to investigate the processes governing the origin and evolution of families of related genes within genomes. RESULTS: Genes were clustered into related groups using similarity criteria derived from BLAST. Most clusters contained genes from only one or a small number of genomes, and relatively few core clusters were found that spanned all genomes. Gene clusters found in larger numbers of genomes tended to have larger numbers of genes per genome; however, clusters with unusually large numbers of genes per genome were found among both narrowly and widely distributed clusters. Larger genomes were found to have larger mean gene family sizes and a greater proportion of families of very large size. We used a model of birth, death, and innovation to predict the distribution of gene family sizes. The key parameter is r, the ratio of duplications to deletions. It was found that the model can give a good fit to the observed distribution only if there are several classes of genes with different values of r. The preferred model in most cases had three classes of genes. CONCLUSIONS: There appears to be a rapid rate of origination of new gene families within individual genomes. Most of these gene families are deleted before they spread to large numbers of genomes, which suggests that they may not be generally beneficial to the organisms. The family size distribution is best described by a large fraction of families that tend to have only one or two genes and a small fraction of families of multi-copy genes that are highly prone to duplication. Larger families occur more frequently in larger genomes, indicating higher r in these genomes, possibly due to a greater tolerance for non-beneficial gene duplicates. The smallest genomes contain very few multi-copy families, suggesting a high rate of deletion of all but the most beneficial genes in these genomes.

Persistence of bacterial and archaeal communities in sea ice through an Arctic winter.

The structure of bacterial communities in first-year spring and summer sea ice differs from that in source seawaters, suggesting selection during ice formation in autumn or taxon-specific mortality in the ice during winter. We tested these hypotheses by weekly sampling (January-March 2004) of first-year winter sea ice (Franklin Bay, Western Arctic) that experienced temperatures from -9 degrees C to -26 degrees C, generating community fingerprints and clone libraries for Bacteria and Archaea. Despite severe conditions and significant decreases in microbial abundance, no significant changes in richness or community structure were detected in the ice. Communities of Bacteria and Archaea in the ice, as in under-ice seawater, were dominated by SAR11 clade Alphaproteobacteria and Marine Group I Crenarchaeota, neither of which is known from later season sea ice. The bacterial ice library contained clones of Gammaproteobacteria from oligotrophic seawater clades (e.g. OM60, OM182) but no clones from gammaproteobacterial genera commonly detected in later season sea ice by similar methods (e.g. Colwellia, Psychrobacter). The only common sea ice bacterial genus detected in winter ice was Polaribacter. Overall, selection during ice formation and mortality during winter appear to play minor roles in the process of microbial succession that leads to distinctive spring and summer sea ice communities.

Humans to Mars: a feasibility and cost-benefit analysis.

Mars is a compelling astrobiological target, and a human mission would provide an opportunity to collect immense amounts of scientific data. Exploration alone, however, cannot justify the increased risk. Instead, three factors drive a human mission: economics, education, and exploration. A human mission has a unique potential to inspire the next generation of young people to enter critically needed science and engineering disciplines. A mission is economically feasible, and the research and development program put in place for a human mission would propel growth in related high-technology industries. The main hurdles are human physiological responses to 1-2 years of radiation and microgravity exposure. However, enabling technologies are sufficiently mature in these areas that they can be developed within a few decade timescale. Hence, the decision of whether or not to undertake a human mission to Mars is a political decision, and thus, educational and economic benefits are the crucial factors.

An in silico assessment of gene function and organization of the phenylpropanoid pathway metabolic networks in Arabidopsis thaliana and limitations thereof.

The Arabidopsis genome sequencing in 2000 gave to science the first blueprint of a vascular plant. Its successful completion also prompted the US National Science Foundation to launch the Arabidopsis 2010 initiative, the goal of which is to identify the function of each gene by 2010. In this study, an exhaustive analysis of The Institute for Genomic Research (TIGR) and The Arabidopsis Information Resource (TAIR) databases, together with all currently compiled EST sequence data, was carried out in order to determine to what extent the various metabolic networks from phenylalanine ammonia lyase (PAL) to the monolignols were organized and/or could be predicted. In these databases, there are some 65 genes which have been annotated as encoding putative enzymatic steps in monolignol biosynthesis, although many of them have only very low homology to monolignol pathway genes of known function in other plant systems. Our detailed analysis revealed that presently only 13 genes (two PALs, a cinnamate-4-hydroxylase, a p-coumarate-3-hydroxylase, a ferulate-5-hydroxylase, three 4-coumarate-CoA ligases, a cinnamic acid O-methyl transferase, two cinnamoyl-CoA reductases) and two cinnamyl alcohol dehydrogenases can be classified as having a bona fide (definitive) function; the remaining 52 genes currently have undetermined physiological roles. The EST database entries for this particular set of genes also provided little new insight into how the monolignol pathway was organized in the different tissues and organs, this being perhaps a consequence of both limitations in how tissue samples were collected and in the incomplete nature of the EST collections. This analysis thus underscores the fact that even with genomic sequencing, presumed to provide the entire suite of putative genes in the monolignol-forming pathway, a very large effort needs to be conducted to establish actual catalytic roles (including enzyme versatility), as well as the physiological function(s) for each member of the (multi)gene families present and the metabolic networks that are operative. Additionally, one key to identifying physiological functions for many of these (and other) unknown genes, and their corresponding metabolic networks, awaits the development of technologies to comprehensively study molecular processes at the single cell level in particular tissues and organs, in order to establish the actual metabolic context.