Jeffrey L. Blanchard

PERSPECTIVE:EVOLUTION AND DETECTION OF GENETIC ROBUSTNESS

Abstract Robustness is the invariance of phenotypes in the face of perturbation. The robustness of phenotypes appears at various levels of biological organization, including gene expression, protein folding, metabolic flux, physiological homeostasis, development, and even organismal fitness. The mechanisms underlying robustness are diverse, ranging from thermodynamic stability at the RNA and protein level to behavior at the organismal level. Phenotypes can be robust either against heritable perturbations (e.g., mutations) or nonheritable perturbations (e.g., the weather). Here we primarily focus on the first kind of robustness—genetic robustness—and survey three growing avenues of research: (1) measuring genetic robustness in nature and in the laboratory; (2) understanding the evolution of genetic robustness; and (3) exploring the implications of genetic robustness for future evolution.

Evolutionary Origins of the Eukaryotic Shikimate Pathway: Gene Fusions, Horizontal Gene Transfer, and Endosymbiotic Replacements

ABSTRACT Currently the shikimate pathway is reported as a metabolic feature of prokaryotes, ascomycete fungi, apicomplexans, and plants. The plant shikimate pathway enzymes have similarities to prokaryote homologues and are largely active in chloroplasts, suggesting ancestry from the plastid progenitor genome. Toxoplasma gondii, which also possesses an alga-derived plastid organelle, encodes a shikimate pathway with similarities to ascomycete genes, including a five-enzyme pentafunctional arom. These data suggests that the shikimate pathway and the pentafunctional arom either had an ancient origin in the eukaryotes or was conveyed by eukaryote-to-eukaryote horizontal gene transfer (HGT). We expand sampling and analyses of the shikimate pathway genes to include the oomycetes, ciliates, diatoms, basidiomycetes, zygomycetes, and the green and red algae. Sequencing of cDNA from Tetrahymena thermophila confirmed the presence of a pentafused arom, as in fungi and T. gondii. Phylogenies and taxon distribution suggest that the arom gene fusion event may be an ancient eukaryotic innovation. Conversely, the Plantae lineage (represented here by both Viridaeplantae and the red algae) acquired different prokaryotic genes for all seven steps of the shikimate pathway. Two of the phylogenies suggest a derivation of the Plantae genes from the cyanobacterial plastid progenitor genome, but if the full Plantae pathway was originally of cyanobacterial origin, then the five other shikimate pathway genes were obtained from a minimum of two other eubacterial genomes. Thus, the phylogenies demonstrate both separate HGTs and shared derived HGTs within the Plantae clade either by primary HGT transfer or secondarily via the plastid progenitor genome. The shared derived characters support the holophyly of the Plantae lineage and a single ancestral primary plastid endosymbiosis. Our analyses also pinpoints a minimum of 50 gene/domain loss events, demonstrating that loss and replacement events have been an important process in eukaryote genome evolution.

Organellar genes: why do they end up in the nucleus?

Many mitochondrial and plastid proteins are derived from their bacterial endosymbiotic ancestors, but their genes now reside on nuclear chromosomes instead of remaining within the organelle. To become an active nuclear gene and return to the organelle as a functional protein, an organellar gene must first be assimilated into the nuclear genome. The gene must then be transcribed and acquire a transit sequence for targeting the protein back to the organelle. On reaching the organelle, the protein must be properly folded and modified, and in many cases assembled in an orderly manner into a larger protein complex. Finally, the nuclear copy must be properly regulated to achieve a fitness level comparable with the organellar gene. Given the complexity in establishing a nuclear copy, why do organellar genes end up in the nucleus? Recent data suggest that these genes are worse off than their nuclear and free-living counterparts because of a reduction in the efficiency of natural selection, but do these population-genetic processes drive the movement of genes to the nucleus? We are now at a stage where we can begin to discriminate between competing hypotheses using a combination of experimental, natural population, bioinformatic and theoretical approaches.

Deleterious mutation accumulation in organelle genomes

It is well established on theoretical grounds that the accumulation of mildly deleterious mutations in nonrecombining genomes is a major extinction risk in obligately asexual populations. Sexual populations can also incur mutational deterioration in genomic regions that experience little or no recombination, i.e., autosomal regions near centromeres, Y chromosomes, and organelle genomes. Our results suggest, for a wide array of genes (transfer RNAs, ribosomal RNAs, and proteins) in a diverse collection of species (animals, plants, and fungi), an almost universal increase in the fixation probabilities of mildly deleterious mutations arising in mitochondrial and chloroplast genomes relative to those arising in the recombining nuclear genome. This enhanced width of the selective sieve in organelle genomes does not appear to be a consequence of relaxed selection, but can be explained by the decline in the efficiency of selection that results from the reduction of effective population size induced by uniparental inheritance. Because of the very low mutation rates of organelle genomes (on the order of 10-4 per genome per year), the reduction in fitness resulting from mutation accumulation in such genomes is a very long-term process, not likely to imperil many species on time scales of less than a million years, but perhaps playing some role in phylogenetic lineage sorting on time scales of 10 to 100 million years.

Tempo and mode of genome evolution in a 50,000-generation experiment

Adaptation by natural selection depends on the rates, effects and interactions of many mutations, making it difficult to determine what proportion of mutations in an evolving lineage are beneficial. Here we analysed 264 complete genomes from 12 Escherichia coli populations to characterize their dynamics over 50,000 generations. The populations that retained the ancestral mutation rate support a model in which most fixed mutations are beneficial, the fraction of beneficial mutations declines as fitness rises, and neutral mutations accumulate at a constant rate. We also compared these populations to mutation-accumulation lines evolved under a bottlenecking regime that minimizes selection. Nonsynonymous mutations, intergenic mutations, insertions and deletions are overrepresented in the long-term populations, further supporting the inference that most mutations that reached high frequency were favoured by selection. These results illuminate the shifting balance of forces that govern genome evolution in populations adapting to a new environment.

Pervasive migration of organellar DNA to the nucleus in plants

A surprisingly large number of plant nuclear DNA sequences inferred to be remnants of chloroplast and mitochondrial DNA migration events were detected through computer-assisted database searches. Nineteen independent organellar DNA insertions, with a median size of 117 by (range of 38 to >785 bp), occur in the proximity of 15 nuclear genes. One fragment appears to have been passed through a RNA intermediate, based on the presence of an edited version of the mitochondrial gene in the nucleus. Tandemly arranged fragments from disparate regions of organellar genomes and from different organellar genomes indicate that the fragments joined together from an intracellular pool of RNA and/or DNA before they integrated into the nuclear genome. Comparisons of integrated sequences to genes lacking the insertions, as well as the occurrence of coligated fragments, support a model of random integration by end joining. All transferred sequences were found in noncoding regions, but the positioning of organellar-derived DNA in introns, as well as regions 5′ and 3′ to nuclear genes, suggests that the random integration of organellar DNA has the potential to influence gene expression patterns. A semiquantitative estimate was performed on the amount of organellar DNA being transferred and assimilated into the nucleus. Based on this database survey, we estimate that 3–7% of the plant nuclear genomic sequence files contain organellar-derived DNA. The timing and the magnitude of genetic flux to the nuclear genome suggest that random integration is a substantial and ongoing process for creating sequence variation.

Rapid Changes in Gene Expression Dynamics in Response to Superoxide Reveal SoxRS-Dependent and Independent Transcriptional Networks

Background SoxR and SoxS constitute an intracellular signal response system that rapidly detects changes in superoxide levels and modulates gene expression in E. coli. A time series microarray design was used to identify co-regulated SoxRS-dependent and independent genes modulated by superoxide minutes after exposure to stress. Methodology/Principal Findings soxS mRNA levels surged to near maximal levels within the first few minutes of exposure to paraquat, a superoxide-producing compound, followed by a rise in mRNA levels of known SoxS-regulated genes. Based on a new method for determining the biological significance of clustering results, a total of 138 genic regions, including several transcription factors and putative sRNAs were identified as being regulated through the SoxRS signaling pathway within 10 minutes of paraquat treatment. A statistically significant two-block SoxS motif was identified through analysis of the SoxS-regulated genes. The SoxRS-independent response included members of the OxyR, CysB, IscR, BirA and Fur regulons. Finally, the relative sensitivity to superoxide was measured in 94 strains carrying deletions in individual, superoxide-regulated genes. Conclusions/Significance By integrating our microarray time series results with other microarray data, E. coli databases and the primary literature, we propose a model of the primary transcriptional response containing 226 protein-coding and sRNA sequences. From the SoxS dependent network the first statistically significant SoxS-related motif was identified.

The non-photosynthetic plastid in malarial parasites and other apicomplexans is derived from outside the green plastid lineage.

The discovery of a non‐photosynthetic plastid genome in Plasmodium falciparum and other apicomplexans has provided a new drug target, but the evolutionary origin of the plastid has been muddled by the lack of characters, that typically define major plastid lineages. To clarify the ancestry of the plastid, we undertook a comprehensive analysis of all genomic characters shared by completely sequenced plastid genomes. Cladistic analysis of the pattern of plastid gene loss and gene rearrangements suggests that the apicomplexan plastid is derived from an ancestor outside of the green plastid lineage. Phylogenetic analysis of primary sequence data (DNA and amino acid characters) produces results that are generally independent of the analytical method, but similar genes (i.e. rpoB and rpoC) give similar topologies. The conflicting phylogenies in primary sequence data sets make it difficult to determine the exact origin of the apicomplexan plastid and the apparent artifactual association of apicomplexan and euglenoid sequences suggests that DNA sequence data may be an inappropriate set of characters to address this phylogenetic question. At present we cannot reject our null hypothesis that the apicomplexan plastid is derived from a shared common ancestor between apicomplexans and dinoflagellates. During the analysis, we noticed that the Plasmodium tRNA‐Met is probably tRNA‐fMet and the tRNA‐fMet is probably tRNA‐Ile. We suggest that P. falciparum has lost the elongator type tRNA‐Met and that similar to metazoan mitochondria there is only one species of methionine tRNA. In P. falciparum, this has been accomplished by recruiting the fMet‐type tRNA to dually function in initiation and elongation. The tRNA‐Ile has an unusual stem‐loop in the variable region. The insertion in this region appears to have occurred after the primary origin of the plastid and further supports the monophyletic ancestory of plastids.

Long-term forest soil warming alters microbial communities in temperate forest soils

Soil microbes are major drivers of soil carbon cycling, yet we lack an understanding of how climate warming will affect microbial communities. Three ongoing field studies at the Harvard Forest Long-term Ecological Research (LTER) site (Petersham, MA) have warmed soils 5°C above ambient temperatures for 5, 8, and 20 years. We used this chronosequence to test the hypothesis that soil microbial communities have changed in response to chronic warming. Bacterial community composition was studied using Illumina sequencing of the 16S ribosomal RNA gene, and bacterial and fungal abundance were assessed using quantitative PCR. Only the 20-year warmed site exhibited significant change in bacterial community structure in the organic soil horizon, with no significant changes in the mineral soil. The dominant taxa, abundant at 0.1% or greater, represented 0.3% of the richness but nearly 50% of the observations (sequences). Individual members of the Actinobacteria, Alphaproteobacteria and Acidobacteria showed strong warming responses, with one Actinomycete decreasing from 4.5 to 1% relative abundance with warming. Ribosomal RNA copy number can obfuscate community profiles, but is also correlated with maximum growth rate or trophic strategy among bacteria. Ribosomal RNA copy number correction did not affect community profiles, but rRNA copy number was significantly decreased in warming plots compared to controls. Increased bacterial evenness, shifting beta diversity, decreased fungal abundance and increased abundance of bacteria with low rRNA operon copy number, including Alphaproteobacteria and Acidobacteria, together suggest that more or alternative niche space is being created over the course of long-term warming.

Involvement of a bacterial microcompartment in the metabolism of fucose and rhamnose by Clostridium phytofermentans.

Background Clostridium phytofermentans, an anaerobic soil bacterium, can directly convert plant biomass into biofuels. The genome of C. phytofermentans contains three loci with genes encoding shell proteins of bacterial microcompartments (BMC), organelles composed entirely of proteins. Methodology and Principal Findings One of the BMC loci has homology to a BMC-encoding locus implicated in the conversion of fucose to propanol and propionate in a human gut commensal, Roseburia inulinivorans. We hypothesized that it had a similar role in C. phytofermentans. When C. phytofermentans was grown on fucose, the major products identified were ethanol, propanol and propionate. Transmission electron microscopy of fucose- and rhamnose-grown cultures revealed polyhedral structures, presumably BMCs. Microarray analysis indicated that during growth on fucose, operons coding for the BMC locus, fucose dissimilatory enzymes, and an ATP-binding cassette transporter became the dominant transcripts. These data are consistent with fucose fermentation producing a 1,2-propanediol intermediate that is further metabolized in the microcompartment encoded in the BMC locus. Growth on another deoxyhexose sugar, rhamnose, resulted in the expression of the same BMC locus and similar fermentation products. However, a different set of dissimilatory enzymes and transport system genes were induced. Quite surprisingly, growth on fucose or rhamnose also led to the expression of a diverse array of complex plant polysaccharide-degrading enzymes. Conclusions/Significance Based on physiological, genomic, and microarray analyses, we propose a model for the fermentation of fucose and rhamnose in C. phytofermentans that includes enzymes encoded in the same BMC locus. Comparative genomic analysis suggests that this BMC may be present in other clostridial species.