Snf2 Family Gene Distribution in Higher Plant Genomes Reveals DRD1 Expansion and Diversification in the Tomato Genome

As part of large protein complexes, Snf2 family ATPases are responsible for energy supply during chromatin remodeling, but the precise mechanism of action of many of these proteins is largely unknown. They influence many processes in plants, such as the response to environmental stress. This analysis is the first comprehensive study of Snf2 family ATPases in plants. We here present a comparative analysis of 1159 candidate plant Snf2 genes in 33 complete and annotated plant genomes, including two green algae. The number of Snf2 ATPases shows considerable variation across plant genomes (17-63 genes). The DRD1, Rad5/16 and Snf2 subfamily members occur most often. Detailed analysis of the plant-specific DRD1 subfamily in related plant genomes shows the occurrence of a complex series of evolutionary events. Notably tomato carries unexpected gene expansions of DRD1 gene members. Most of these genes are expressed in tomato, although at low levels and with distinct tissue or organ specificity. In contrast, the Snf2 subfamily genes tend to be expressed constitutively in tomato. The results underpin and extend the Snf2 subfamily classification, which could help to determine the various functional roles of Snf2 ATPases and to target environmental stress tolerance and yield in future breeding.     

Evolution of Genome Size and Complexity in Pinus

Background

Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood.

Methodology/Principal Findings

Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA.

Conclusions/Significance

Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.     

Phylobiochemical Characterization of Class-Ib Aspartate/Prephenate Aminotransferases Reveals Evolution of the Plant Arogenate Phenylalanine Pathway

The aromatic amino acid Phe is required for protein synthesis and serves as the precursor of abundant phenylpropanoid plant natural products. While Phe is synthesized from prephenate exclusively via a phenylpyruvate intermediate in model microbes, the alternative pathway via arogenate is predominant in plant Phe biosynthesis. However, the molecular and biochemical evolution of the plant arogenate pathway is currently unknown. Here, we conducted phylogenetically informed biochemical characterization of prephenate aminotransferases (PPA-ATs) that belong to class-Ib aspartate aminotransferases (AspAT Ibs) and catalyze the first committed step of the arogenate pathway in plants. Plant PPA-ATs and succeeding arogenate dehydratases (ADTs) were found to be most closely related to homologs from Chlorobi/Bacteroidetes bacteria. The Chlorobium tepidum PPA-AT and ADT homologs indeed efficiently converted prephenate and arogenate into arogenate and Phe, respectively. A subset of AspAT Ib enzymes exhibiting PPA-AT activity was further identified from both Plantae and prokaryotes and, together with site-directed mutagenesis, showed that Thr-84 and Lys-169 play key roles in specific recognition of dicarboxylic keto (prephenate) and amino (aspartate) acid substrates. The results suggest that, along with ADT, a gene encoding prephenate-specific PPA-AT was transferred from a Chlorobi/Bacteroidetes ancestor to a eukaryotic ancestor of Plantae, allowing efficient Phe and phenylpropanoid production via arogenate in plants today.     

A Family of Differentially Amplified Repetitive DNA Sequences in the Genus Beta Reveals Genetic Variation in Beta vulgaris Subspecies and Cultivars

Members of a highly abundant restriction satellite family have been isolated from the wild beet species Beta nana. The satellite DNA sequence is characterized by a conserved RsaI restriction site and is present in three of four sections of the genus Beta, namely Nanae, Corollinae, and Beta. It was not detected in species of the evolutionary old section Procumbentes, suggesting its amplification after separation of this section. Sequences of eight monomers were aligned revealing a size variation from 209 to 233 bp and an AT content ranging from 56.5% to 60.5%. The similarity between monomers in B. nana varied from 77.7% to 92.2%. Diverged subfamilies were identified by sequence analysis and Southern hybridization. A comparative study of this repetitive DNA element by fluorescent in situ hybridization and Southern analyses in three representative species was performed showing a variable genomic organization and heterogeneous localizations along metaphase chromosomes both within and between species. In B. nana the copy number of this satellite, with some 30,000 per haploid genome, is more than tenfold higher than in Beta lomatogona and up to 200 times higher than in Beta vulgaris, indicating different levels of sequence amplification during evolution in the genus Beta. In sugar beet (B. vulgaris), the large-scale organization of this tandem repeat was examined by pulsed-field gel electrophoresis. Southern hybridization to genomic DNA digested with DraI demonstrated that satellite arrays are located in AT-rich regions and the tandem repeat is a useful probe for the detection of genetic variation in closely related B. vulgaris cultivars, accessions, and subspecies. Received: 24 May 1996 / Accepted: 13 September 1996     

De novo Genome Assembly of the Fungal Plant Pathogen Pyrenophora semeniperda

Pyrenophora semeniperda (anamorph Drechslera campulata) is a necrotrophic fungal seed pathogen that has a wide host range within the Poaceae. One of its hosts is cheatgrass (Bromus tectorum), a species exotic to the United States that has invaded natural ecosystems of the Intermountain West. As a natural pathogen of cheatgrass, P. semeniperda has potential as a biocontrol agent due to its effectiveness at killing seeds within the seed bank; however, few genetic resources exist for the fungus. Here, the genome of P. semeniperda isolate assembled from sequence reads of 454 pyrosequencing is presented. The total assembly is 32.5 Mb and includes 11,453 gene models encoding putative proteins larger than 24 amino acids. The models represent a variety of putative genes that are involved in pathogenic pathways typically found in necrotrophic fungi. In addition, extensive rearrangements, including inter- and intrachromosomal rearrangements, were found when the P. semeniperda genome was compared to P. tritici-repentis, a related fungal species.     

In Vivo Mutational Characterization of DndE Involved in DNA Phosphorothioate Modification

DNA phosphorothioate (PT) modification is a recently identified epigenetic modification that occurs in the sugar-phosphate backbone of prokaryotic DNA. Previous studies have demonstrated that DNA PT modification is governed by the five DndABCDE proteins in a sequence-selective and R _P stereo-specific manner. Bacteria may have acquired this physiological modification along with dndFGH as a restriction-modification system. However, little is known about the biological function of Dnd proteins, especially the smallest protein, DndE, in the PT modification pathway. DndE was reported to be a DNA-binding protein with a preference for nicked dsDNA in vitro; the binding of DndE to DNA occurs via six positively charged lysine residues on its surface. The substitution of these key lysine residues significantly decreased the DNA binding affinities of DndE proteins to undetectable levels. In this study, we conducted site-directed mutagenesis of dndE on a plasmid and measured DNA PT modifications under physiological conditions by mass spectrometry. We observed distinctive differences from the in vitro binding assays. Several mutants with lysine residues mutated to alanine decreased the total frequency of PT modifications, but none of the mutants completely eliminated PT modification. Our results suggest that the nicked dsDNA-binding capacity of DndE may not be crucial for PT modification and/or that DndE may have other biological functions in addition to binding to dsDNA.     

In Situ Characterization of Intrahepatic Non-Parenchymal Cells in PSC Reveals Phenotypic Patterns Associated with Disease Severity

Liver-infiltrating T cells have been implicated in the pathogenesis of primary sclerosing cholangitis (PSC), however little information is available about changes in other cellular compartments in the liver during PSC. This study aimed to characterize non-parenchymal intrahepatic cells in PSC livers and to find associations between phenotypes and disease severity. Using immunohistochemistry, followed by automated image analysis and quantification and a principal component analysis, we have studied non-parenchymal intrahepatic cells in PSC-patient livers (n = 17) and controls (n = 17). We observed a significant increase of T cells in the PSC patients, localized to the fibrotic areas. MAIT cells, normally present at high numbers in the liver, were not increased to the same extent. PSC patients had lower expression of MHC class I than controls. However, the levels of NKp46+ NK cells were similar between patients and controls, nevertheless, NKp46 was identified as a phenotypic marker that distinguished PSC patients with mild from those with severe fibrosis. Beyond that, a group of PSC patients had lost expression of Caldesmon and this was associated with more extensive bile duct proliferation and higher numbers of T cells. Our data reveals phenotypic patterns in PSC patients associated with disease severity.     

Patterns of Intraspecific DNA Variation in the Daphnia Nuclear Genome

Understanding nucleotide variation in natural populations has been a subject of great interest for decades. However, many taxonomic groups, especially those with atypical life history attributes remain unstudied, and Drosophila is the only arthropod genus for which DNA polymorphism data are presently abundant. As a result of the recent release of the complete genome sequence and a wide variety of new genomic resources, the Daphnia system is quickly becoming a promising new avenue for expanding our knowledge of nucleotide variation in natural populations. Here, we examine nucleotide variation in six protein-coding loci for Daphnia pulex and its congeners with particular emphasis on D. pulicaria, the closest extant relative of D. pulex. Levels of synonymous intraspecific variation, π_s, averaged 0.0136 for species in the Daphnia genus, and are slightly lower than most prior estimates in invertebrates. Tests of neutrality indicated that segregating variation conforms to neutral model expectations for the loci that we examined in most species, while K_a/K_s ratios revealed strong purifying selection. Using a full maximum-likelihood coalescent-based method, the ratio of the recombination rate to the mutation rate (c/u), averaged 0.5255 for species of the Daphnia genus. Lastly, a divergence population-genetics approach was used to investigate gene flow and divergence between D. pulex and D. pulicaria.     

The Complete Chloroplast Genome Sequence of the Medicinal Plant Salvia miltiorrhiza

Salvia miltiorrhiza is an important medicinal plant with great economic and medicinal value. The complete chloroplast (cp) genome sequence of Salvia miltiorrhiza, the first sequenced member of the Lamiaceae family, is reported here. The genome is 151,328 bp in length and exhibits a typical quadripartite structure of the large (LSC, 82,695 bp) and small (SSC, 17,555 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,539 bp). It contains 114 unique genes, including 80 protein-coding genes, 30 tRNAs and four rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Four forward, three inverted and seven tandem repeats were detected in the Salvia miltiorrhiza cp genome. Simple sequence repeat (SSR) analysis among the 30 asterid cp genomes revealed that most SSRs are AT-rich, which contribute to the overall AT richness of these cp genomes. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the non-coding regions, indicating an uneven distribution of SSRs within the cp genomes. Entire cp genome comparison of Salvia miltiorrhiza and three other Lamiales cp genomes showed a high degree of sequence similarity and a relatively high divergence of intergenic spacers. Sequence divergence analysis discovered the ten most divergent and ten most conserved genes as well as their length variation, which will be helpful for phylogenetic studies in asterids. Our analysis also supports that both regional and functional constraints affect gene sequence evolution. Further, phylogenetic analysis demonstrated a sister relationship between Salvia miltiorrhiza and Sesamum indicum. The complete cp genome sequence of Salvia miltiorrhiza reported in this paper will facilitate population, phylogenetic and cp genetic engineering studies of this medicinal plant.     

Bifidobacterium asteroides PRL2011 Genome Analysis Reveals Clues for Colonization of the Insect Gut

Bifidobacteria are known as anaerobic/microaerophilic and fermentative microorganisms, which commonly inhabit the gastrointestinal tract of various animals and insects. Analysis of the 2,167,301 bp genome of Bifidobacterium asteroides PRL2011, a strain isolated from the hindgut of Apis mellifera var. ligustica, commonly known as the honey bee, revealed its predicted capability for respiratory metabolism. Conservation of the latter gene clusters in various B. asteroides strains enforces the notion that respiration is a common metabolic feature of this ancient bifidobacterial species, which has been lost in currently known mammal-derived Bifidobacterium species. In fact, phylogenomic based analyses suggested an ancient origin of B. asteroides and indicates it as an ancestor of the genus Bifidobacterium. Furthermore, the B. asteroides PRL2011 genome encodes various enzymes for coping with toxic products that arise as a result of oxygen-mediated respiration.     

In Vivo Time-Resolved Microtomography Reveals the Mechanics of the Blowfly Flight Motor

Dipteran flies are amongst the smallest and most agile of flying animals. Their wings are driven indirectly by large power muscles, which cause cyclical deformations of the thorax that are amplified through the intricate wing hinge. Asymmetric flight manoeuvres are controlled by 13 pairs of steering muscles acting directly on the wing articulations. Collectively the steering muscles account for <3% of total flight muscle mass, raising the question of how they can modulate the vastly greater output of the power muscles during manoeuvres. Here we present the results of a synchrotron-based study performing micrometre-resolution, time-resolved microtomography on the 145 Hz wingbeat of blowflies. These data represent the first four-dimensional visualizations of an organism''s internal movements on sub-millisecond and micrometre scales. This technique allows us to visualize and measure the three-dimensional movements of five of the largest steering muscles, and to place these in the context of the deforming thoracic mechanism that the muscles actuate. Our visualizations show that the steering muscles operate through a diverse range of nonlinear mechanisms, revealing several unexpected features that could not have been identified using any other technique. The tendons of some steering muscles buckle on every wingbeat to accommodate high amplitude movements of the wing hinge. Other steering muscles absorb kinetic energy from an oscillating control linkage, which rotates at low wingbeat amplitude but translates at high wingbeat amplitude. Kinetic energy is distributed differently in these two modes of oscillation, which may play a role in asymmetric power management during flight control. Structural flexibility is known to be important to the aerodynamic efficiency of insect wings, and to the function of their indirect power muscles. We show that it is integral also to the operation of the steering muscles, and so to the functional flexibility of the insect flight motor.     

Meiosis Drives Extraordinary Genome Plasticity in the Haploid Fungal Plant Pathogen Mycosphaerella graminicola

Meiosis in the haploid plant-pathogenic fungus Mycosphaerella graminicola results in eight ascospores due to a mitotic division following the two meiotic divisions. The transient diploid phase allows for recombination among homologous chromosomes. However, some chromosomes of M. graminicola lack homologs and do not pair during meiosis. Because these chromosomes are not present universally in the genome of the organism they can be considered to be dispensable. To analyze the meiotic transmission of unequal chromosome numbers, two segregating populations were generated by crossing genetically unrelated parent isolates originating from Algeria and The Netherlands that had pathogenicity towards durum or bread wheat, respectively. Detailed genetic analyses of these progenies using high-density mapping (1793 DArT, 258 AFLP and 25 SSR markers) and graphical genotyping revealed that M. graminicola has up to eight dispensable chromosomes, the highest number reported in filamentous fungi. These chromosomes vary from 0.39 to 0.77 Mb in size, and represent up to 38% of the chromosomal complement. Chromosome numbers among progeny isolates varied widely, with some progeny missing up to three chromosomes, while other strains were disomic for one or more chromosomes. Between 15–20% of the progeny isolates lacked one or more chromosomes that were present in both parents. The two high-density maps showed no recombination of dispensable chromosomes and hence, their meiotic processing may require distributive disjunction, a phenomenon that is rarely observed in fungi. The maps also enabled the identification of individual twin isolates from a single ascus that shared the same missing or doubled chromosomes indicating that the chromosomal polymorphisms were mitotically stable and originated from nondisjunction during the second division and, less frequently, during the first division of fungal meiosis. High genome plasticity could be among the strategies enabling this versatile pathogen to quickly overcome adverse biotic and abiotic conditions in wheat fields.     

Extreme Recombination Frequencies Shape Genome Variation and Evolution in the Honeybee,Apis mellifera

Meiotic recombination is a fundamental cellular process, with important consequences for evolution and genome integrity. However, we know little about how recombination rates vary across the genomes of most species and the molecular and evolutionary determinants of this variation. The honeybee, Apis mellifera, has extremely high rates of meiotic recombination, although the evolutionary causes and consequences of this are unclear. Here we use patterns of linkage disequilibrium in whole genome resequencing data from 30 diploid honeybees to construct a fine-scale map of rates of crossing over in the genome. We find that, in contrast to vertebrate genomes, the recombination landscape is not strongly punctate. Crossover rates strongly correlate with levels of genetic variation, but not divergence, which indicates a pervasive impact of selection on the genome. Germ-line methylated genes have reduced crossover rate, which could indicate a role of methylation in suppressing recombination. Controlling for the effects of methylation, we do not infer a strong association between gene expression patterns and recombination. The site frequency spectrum is strongly skewed from neutral expectations in honeybees: rare variants are dominated by AT-biased mutations, whereas GC-biased mutations are found at higher frequencies, indicative of a major influence of GC-biased gene conversion (gBGC), which we infer to generate an allele fixation bias 5 – 50 times the genomic average estimated in humans. We uncover further evidence that this repair bias specifically affects transitions and favours fixation of CpG sites. Recombination, via gBGC, therefore appears to have profound consequences on genome evolution in honeybees and interferes with the process of natural selection. These findings have important implications for our understanding of the forces driving molecular evolution.     

Host Species-specific Repetitive DNA Sequence in the Genome of Magnaporthe grisea,the Rice Blast Fungus

We cloned a repetitive sequence to show RFLPs in the genome of Magnaporthe grisea, a fungal pathogen responsible for rice blast. As the sequence was 0.8 kb in length and dispersed in the genome, it was named MGSR1 (for Magnaporthe grisea short repeat 1). MGSR1 was conserved highly in the genome of rice pathogens, but poorly in the genome of pathogens of grasses other than rice. And the RFLPs, displayed with the sequence, could distinguish between clonal lineages in rice-pathogenic isolates. The nucleotide sequence showed the presence of an internal promoter of RNA polymerase III, a 3?-poly(T), and an 8-bp direct repeat in it.     

Whole Genome Sequencing of Mutation Accumulation Lines Reveals a Low Mutation Rate in the Social Amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10⁻⁹, with a Poisson confidence interval of 4.1×10⁻⁹ − 9.5×10⁻⁹, per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10⁻¹¹, with a Poisson confidence interval ranging from 7.4×10⁻¹³ to 1.6×10⁻¹⁰, is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.     

Interactome-Wide Prediction of Protein-Protein Binding Sites Reveals Effects of Protein Sequence Variation in Arabidopsis thaliana

The specificity of protein-protein interactions is encoded in those parts of the sequence that compose the binding interface. Therefore, understanding how changes in protein sequence influence interaction specificity, and possibly the phenotype, requires knowing the location of binding sites in those sequences. However, large-scale detection of protein interfaces remains a challenge. Here, we present a sequence- and interactome-based approach to mine interaction motifs from the recently published Arabidopsis thaliana interactome. The resultant proteome-wide predictions are available via www.ab.wur.nl/sliderbio and set the stage for further investigations of protein-protein binding sites. To assess our method, we first show that, by using a priori information calculated from protein sequences, such as evolutionary conservation and residue surface accessibility, we improve the performance of interface prediction compared to using only interactome data. Next, we present evidence for the functional importance of the predicted sites, which are under stronger selective pressure than the rest of protein sequence. We also observe a tendency for compensatory mutations in the binding sites of interacting proteins. Subsequently, we interrogated the interactome data to formulate testable hypotheses for the molecular mechanisms underlying effects of protein sequence mutations. Examples include proteins relevant for various developmental processes. Finally, we observed, by analysing pairs of paralogs, a correlation between functional divergence and sequence divergence in interaction sites. This analysis suggests that large-scale prediction of binding sites can cast light on evolutionary processes that shape protein-protein interaction networks.     

Variation in the Fitness Effects of Mutations with Population Density and Size in Escherichia coli

The fitness effects of mutations are context specific and depend on both external (e.g., environment) and internal (e.g., cellular stress, genetic background) factors. The influence of population size and density on fitness effects are unknown, despite the central role population size plays in the supply and fixation of mutations. We addressed this issue by comparing the fitness of 92 Keio strains (Escherichia coli K12 single gene knockouts) at comparatively high (1.2×10⁷ CFUs/mL) and low (2.5×10² CFUs/mL) densities, which also differed in population size (high: 1.2×10⁸; low: 1.25×10³). Twenty-eight gene deletions (30%) exhibited a fitness difference, ranging from 5 to 174% (median: 35%), between the high and low densities. Our analyses suggest this variation among gene deletions in fitness responses reflected in part both gene orientation and function, of the gene properties we examined (genomic position, length, orientation, and function). Although we could not determine the relative effects of population density and size, our results suggest fitness effects of mutations vary with these two factors, and this variation is gene-specific. Besides being a mechanism for density-dependent selection (r-K selection), the dependence of fitness effects on population density and size has implications for any population that varies in size over time, including populations undergoing evolutionary rescue, species invasions into novel habitats, and cancer progression and metastasis. Further, combined with recent advances in understanding the roles of other context-specific factors in the fitness effects of mutations, our results will help address theoretical and applied biological questions more realistically.