Discovery of photosynthesis genes through whole-genome sequencing of acetate-requiring mutants of Chlamydomonas reinhardtii.

Large-scale mutant libraries have been indispensable for genetic studies, and the development of next-generation genome sequencing technologies has greatly advanced efforts to analyze mutants. In this work, we sequenced the genomes of 660 Chlamydomonas reinhardtii acetate-requiring mutants, part of a larger photosynthesis mutant collection previously generated by insertional mutagenesis with a linearized plasmid. We identified 554 insertion events from 509 mutants by mapping the plasmid insertion sites through paired-end sequences, in which one end aligned to the plasmid and the other to a chromosomal location. Nearly all (96%) of the events were associated with deletions, duplications, or more complex rearrangements of genomic DNA at the sites of plasmid insertion, and together with deletions that were unassociated with a plasmid insertion, 1470 genes were identified to be affected. Functional annotations of these genes were enriched in those related to photosynthesis, signaling, and tetrapyrrole synthesis as would be expected from a library enriched for photosynthesis mutants. Systematic manual analysis of the disrupted genes for each mutant generated a list of 253 higher-confidence candidate photosynthesis genes, and we experimentally validated two genes that are essential for photoautotrophic growth, CrLPA3 and CrPSBP4. The inventory of candidate genes includes 53 genes from a phylogenomically defined set of conserved genes in green algae and plants. Altogether, 70 candidate genes encode proteins with previously characterized functions in photosynthesis in Chlamydomonas, land plants, and/or cyanobacteria; 14 genes encode proteins previously shown to have functions unrelated to photosynthesis. Among the remaining 169 uncharacterized genes, 38 genes encode proteins without any functional annotation, signifying that our results connect a function related to photosynthesis to these previously unknown proteins. This mutant library, with genome sequences that reveal the molecular extent of the chromosomal lesions and resulting higher-confidence candidate genes, will aid in advancing gene discovery and protein functional analysis in photosynthesis.

Genome sequencing of evolved aspergilli populations reveals robust genomes, transversions in A. flavus, and sexual aberrancy in non-homologous end-joining mutants.

BACKGROUND: Aspergillus spp. comprises a very diverse group of lower eukaryotes with a high relevance for industrial applications and clinical implications. These multinucleate species are often cultured for many generations in the laboratory, which can unknowingly propagate hidden genetic mutations. To assess the likelihood of such events, we studied the genome stability of aspergilli by using a combination of mutation accumulation (MA) lines and whole genome sequencing. RESULTS: We sequenced the whole genomes of 30 asexual and 10 sexual MA lines of three Aspergillus species (A. flavus, A. fumigatus and A. nidulans) and estimated that each MA line accumulated mutations for over 4000 mitoses during asexual cycles. We estimated mutation rates of 4.2 × 10-11 (A. flavus), 1.1 × 10-11 (A. fumigatus) and 4.1 × 10-11 (A. nidulans) per site per mitosis, suggesting that the genomes are very robust. Unexpectedly, we found a very high rate of GC → TA transversions only in A. flavus. In parallel, 30 asexual lines of the non-homologous end-joining (NHEJ) mutants of the three species were also allowed to accumulate mutations for the same number of mitoses. Sequencing of these NHEJ MA lines gave an estimated mutation rate of 5.1 × 10-11 (A. flavus), 2.2 × 10-11 (A. fumigatus) and 4.5 × 10-11 (A. nidulans) per base per mitosis, which is slightly higher than in the wild-type strains and some ~ 5-6 times lower than in the yeasts. Additionally, in A. nidulans, we found a NHEJ-dependent interference of the sexual cycle that is independent of the accumulation of mutations. CONCLUSIONS: We present for the first time direct counts of the mutation rate of filamentous fungal species and find that Aspergillus genomes are very robust. Deletion of the NHEJ machinery results in a slight increase in the mutation rate, but at a rate we suggest is still safe to use for biotechnology purposes. Unexpectedly, we found GC→TA transversions predominated only in the species A. flavus, which could be generated by the hepatocarcinogen secondary metabolite aflatoxin. Lastly, a strong effect of the NHEJ mutation in self-crossing was observed and an increase in the mutations of the asexual lines was quantified.

Genome sequence of the model rice variety KitaakeX.

BACKGROUND: The availability of thousands of complete rice genome sequences from diverse varieties and accessions has laid the foundation for in-depth exploration of the rice genome. One drawback to these collections is that most of these rice varieties have long life cycles, and/or low transformation efficiencies, which limits their usefulness as model organisms for functional genomics studies. In contrast, the rice variety Kitaake has a rapid life cycle (9 weeks seed to seed) and is easy to transform and propagate. For these reasons, Kitaake has emerged as a model for studies of diverse monocotyledonous species. RESULTS: Here, we report the de novo genome sequencing and analysis of Oryza sativa ssp. japonica variety KitaakeX, a Kitaake plant carrying the rice XA21 immune receptor. Our KitaakeX sequence assembly contains 377.6 Mb, consisting of 33 scaffolds (476 contigs) with a contig N50 of 1.4 Mb. Complementing the assembly are detailed gene annotations of 35,594 protein coding genes. We identified 331,335 genomic variations between KitaakeX and Nipponbare (ssp. japonica), and 2,785,991 variations between KitaakeX and Zhenshan97 (ssp. indica). We also compared Kitaake resequencing reads to the KitaakeX assembly and identified 219 small variations. The high-quality genome of the model rice plant KitaakeX will accelerate rice functional genomics. CONCLUSIONS: The high quality, de novo assembly of the KitaakeX genome will serve as a useful reference genome for rice and will accelerate functional genomics studies of rice and other species.

Whole-Genome Sequencing Identifies a Rice Grain Shape Mutant, gs9-1.

BACKGROUND: Breeding for genes controlling key agronomic traits is an important goal of rice genetic improvement. To gain insight into genes controlling grain morphology, we screened M3 plants derived from 1,000 whole-genome sequenced (WGS) M2 Kitaake mutants to identify lines with altered grain size. RESULTS: In this study, we isolated a mutant, named fast-neutron (FN) 60-4, which exhibits a significant reduction in grain size. We crossed FN60-4 with the parental line Kitaake and analyzed the resulting backcross population. Segregation analysis of 113 lines from the BC2F2 population revealed that the mutant phenotype is controlled by a single semi-dominant locus. Mutant FN60-4 is reduced 20% in plant height and 8.8% in 1000-grain weight compared with Kitaake. FN60-4 also exhibits an 8% reduction in cell number and a 9% reduction in cell length along the vertical axis of the glume. We carried out whole-genome sequencing of DNA pools extracted from segregants with long grains or short grains, and revealed that one gene, LOC_Os09g02650, cosegregated with the grain size phenotype in the BC1F2 and BC2F2 populations. This mutant allele was named grain shape 9-1 (gs9-1). gs9-1 carries a 3-bp deletion that affects two amino acids. This locus is a new allele of the BC12/GDD1/MTD1 gene that encodes a kinesin-like protein involved in cell-cycle progression, cellulose microfibril deposition and gibberellic acid (GA) biosynthesis. The GA biosynthesis-related gene KO2 is down-regulated in gs9-1. The dwarf phenotype of gs9-1 can be rescued by adding exogenous GA3. In contrast to the phenotypes for the other alleles, the gs9-1 is less severe, consistent with the nature of the mutation, which does not disrupt the open reading frame as observed for the other alleles. CONCLUSIONS: In this study, we isolated a mutant, which exhibits altered grain shape and identified the mutated gene, gs9-1. Our study reveals that gs9-1 is a semi-dominant gene that carries a two-amino acid mutation. gs9-1 is allelic to the BC12/GDD1/MTD1 gene involved in GA biosynthesis. These results demonstrate the efficiency and convenience of cloning genes from the whole-genome sequenced Kitaake mutant population to advance investigations into genes controlling key agronomic traits in rice.

Mediation of plant-mycorrhizal interaction by a lectin receptor-like kinase.

The molecular mechanisms underlying mycorrhizal symbioses, the most ubiquitous and impactful mutualistic plant-microbial interaction in nature, are largely unknown. Through genetic mapping, resequencing and molecular validation, we demonstrate that a G-type lectin receptor-like kinase (lecRLK) mediates the symbiotic interaction between Populus and the ectomycorrhizal fungus Laccaria bicolor. This finding uncovers an important molecular step in the establishment of symbiotic plant-fungal associations and provides a molecular target for engineering beneficial mycorrhizal relationships.

Experimental Evolution of Extreme Resistance to Ionizing Radiation in Escherichia coli after 50 Cycles of Selection.

In previous work (D. R. Harris et al., J Bacteriol 191:5240-5252, 2009, https://doi.org/10.1128/JB.00502-09; B. T. Byrne et al., Elife 3:e01322, 2014, https://doi.org/10.7554/eLife.01322), we demonstrated that Escherichia coli could acquire substantial levels of resistance to ionizing radiation (IR) via directed evolution. Major phenotypic contributions involved adaptation of organic systems for DNA repair. We have now undertaken an extended effort to generate E. coli populations that are as resistant to IR as Deinococcus radiodurans After an initial 50 cycles of selection using high-energy electron beam IR, four replicate populations exhibit major increases in IR resistance but have not yet reached IR resistance equivalent to D. radiodurans Regular deep sequencing reveals complex evolutionary patterns with abundant clonal interference. Prominent IR resistance mechanisms involve novel adaptations to DNA repair systems and alterations in RNA polymerase. Adaptation is highly specialized to resist IR exposure, since isolates from the evolved populations exhibit highly variable patterns of resistance to other forms of DNA damage. Sequenced isolates from the populations possess between 184 and 280 mutations. IR resistance in one isolate, IR9-50-1, is derived largely from four novel mutations affecting DNA and RNA metabolism: RecD A90E, RecN K429Q, and RpoB S72N/RpoC K1172I. Additional mechanisms of IR resistance are evident.IMPORTANCE Some bacterial species exhibit astonishing resistance to ionizing radiation, with Deinococcus radiodurans being the archetype. As natural IR sources rarely exceed mGy levels, the capacity of Deinococcus to survive 5,000 Gy has been attributed to desiccation resistance. To understand the molecular basis of true extreme IR resistance, we are using experimental evolution to generate strains of Escherichia coli with IR resistance levels comparable to Deinococcus Experimental evolution has previously generated moderate radioresistance for multiple bacterial species. However, these efforts could not take advantage of modern genomic sequencing technologies. In this report, we examine four replicate bacterial populations after 50 selection cycles. Genomic sequencing allows us to follow the genesis of mutations in populations throughout selection. Novel mutations affecting genes encoding DNA repair proteins and RNA polymerase enhance radioresistance. However, more contributors are apparent.

Corrigendum: Regulation of multiple carbon monoxide consumption pathways in anaerobic bacteria.

[This corrects the article on p. 147 in vol. 2, PMID: 21808633.].

Genomic Sequence Variation Analysis by Resequencing.

Whole-genome resequencing is a method for determining the differences between individuals and a reference genome. The experiments are performed by sequencing the individuals, aligning generated reads to a common reference and discovering variation within the data set by analysis of the alignment with software tools. When correlated with phenotypic information, sites of causative genomic variation may be putatively assigned.While the analysis is generally straightforward, there are many nuances, and we aim to help you understand how to generate an initial result, sift through it to identify likely candidates for a phenotype of interest, and flag false positive calls.

Extensive gene content variation in the Brachypodium distachyon pan-genome correlates with population structure.

While prokaryotic pan-genomes have been shown to contain many more genes than any individual organism, the prevalence and functional significance of differentially present genes in eukaryotes remains poorly understood. Whole-genome de novo assembly and annotation of 54 lines of the grass Brachypodium distachyon yield a pan-genome containing nearly twice the number of genes found in any individual genome. Genes present in all lines are enriched for essential biological functions, while genes present in only some lines are enriched for conditionally beneficial functions (e.g., defense and development), display faster evolutionary rates, lie closer to transposable elements and are less likely to be syntenic with orthologous genes in other grasses. Our data suggest that differentially present genes contribute substantially to phenotypic variation within a eukaryote species, these genes have a major influence in population genetics, and transposable elements play a key role in pan-genome evolution.

Population genomics of picophytoplankton unveils novel chromosome hypervariability.

Tiny photosynthetic microorganisms that form the picoplankton (between 0.3 and 3 µm in diameter) are at the base of the food web in many marine ecosystems, and their adaptability to environmental change hinges on standing genetic variation. Although the genomic and phenotypic diversity of the bacterial component of the oceans has been intensively studied, little is known about the genomic and phenotypic diversity within each of the diverse eukaryotic species present. We report the level of genomic diversity in a natural population of Ostreococcus tauri (Chlorophyta, Mamiellophyceae), the smallest photosynthetic eukaryote. Contrary to the expectations of clonal evolution or cryptic species, the spectrum of genomic polymorphism observed suggests a large panmictic population (an effective population size of 1.2 × 107) with pervasive evidence of sexual reproduction. De novo assemblies of low-coverage chromosomes reveal two large candidate mating-type loci with suppressed recombination, whose origin may pre-date the speciation events in the class Mamiellophyceae. This high genetic diversity is associated with large phenotypic differences between strains. Strikingly, resistance of isolates to large double-stranded DNA viruses, which abound in their natural environment, is positively correlated with the size of a single hypervariable chromosome, which contains 44 to 156 kb of strain-specific sequences. Our findings highlight the role of viruses in shaping genome diversity in marine picoeukaryotes.

The Sequences of 1504 Mutants in the Model Rice Variety Kitaake Facilitate Rapid Functional Genomic Studies.

The availability of a whole-genome sequenced mutant population and the cataloging of mutations of each line at a single-nucleotide resolution facilitate functional genomic analysis. To this end, we generated and sequenced a fast-neutron-induced mutant population in the model rice cultivar Kitaake (Oryza sativa ssp japonica), which completes its life cycle in 9 weeks. We sequenced 1504 mutant lines at 45-fold coverage and identified 91,513 mutations affecting 32,307 genes, i.e., 58% of all rice genes. We detected an average of 61 mutations per line. Mutation types include single-base substitutions, deletions, insertions, inversions, translocations, and tandem duplications. We observed a high proportion of loss-of-function mutations. We identified an inversion affecting a single gene as the causative mutation for the short-grain phenotype in one mutant line. This result reveals the usefulness of the resource for efficient, cost-effective identification of genes conferring specific phenotypes. To facilitate public access to this genetic resource, we established an open access database called KitBase that provides access to sequence data and seed stocks. This population complements other available mutant collections and gene-editing technologies. This work demonstrates how inexpensive next-generation sequencing can be applied to generate a high-density catalog of mutations.

Expansion of Signal Transduction Pathways in Fungi by Extensive Genome Duplication.

Plants and fungi use light and other signals to regulate development, growth, and metabolism. The fruiting bodies of the fungus Phycomyces blakesleeanus are single cells that react to environmental cues, including light, but the mechanisms are largely unknown [1]. The related fungus Mucor circinelloides is an opportunistic human pathogen that changes its mode of growth upon receipt of signals from the environment to facilitate pathogenesis [2]. Understanding how these organisms respond to environmental cues should provide insights into the mechanisms of sensory perception and signal transduction by a single eukaryotic cell, and their role in pathogenesis. We sequenced the genomes of P. blakesleeanus and M. circinelloides and show that they have been shaped by an extensive genome duplication or, most likely, a whole-genome duplication (WGD), which is rarely observed in fungi [3-6]. We show that the genome duplication has expanded gene families, including those involved in signal transduction, and that duplicated genes have specialized, as evidenced by differences in their regulation by light. The transcriptional response to light varies with the developmental stage and is still observed in a photoreceptor mutant of P. blakesleeanus. A phototropic mutant of P. blakesleeanus with a heterozygous mutation in the photoreceptor gene madA demonstrates that photosensor dosage is important for the magnitude of signal transduction. We conclude that the genome duplication provided the means to improve signal transduction for enhanced perception of environmental signals. Our results will help to understand the role of genome dynamics in the evolution of sensory perception in eukaryotes.

Genome-wide selective sweeps and gene-specific sweeps in natural bacterial populations.

Multiple models describe the formation and evolution of distinct microbial phylogenetic groups. These evolutionary models make different predictions regarding how adaptive alleles spread through populations and how genetic diversity is maintained. Processes predicted by competing evolutionary models, for example, genome-wide selective sweeps vs gene-specific sweeps, could be captured in natural populations using time-series metagenomics if the approach were applied over a sufficiently long time frame. Direct observations of either process would help resolve how distinct microbial groups evolve. Here, from a 9-year metagenomic study of a freshwater lake (2005-2013), we explore changes in single-nucleotide polymorphism (SNP) frequencies and patterns of gene gain and loss in 30 bacterial populations. SNP analyses revealed substantial genetic heterogeneity within these populations, although the degree of heterogeneity varied by >1000-fold among populations. SNP allele frequencies also changed dramatically over time within some populations. Interestingly, nearly all SNP variants were slowly purged over several years from one population of green sulfur bacteria, while at the same time multiple genes either swept through or were lost from this population. These patterns were consistent with a genome-wide selective sweep in progress, a process predicted by the 'ecotype model' of speciation but not previously observed in nature. In contrast, other populations contained large, SNP-free genomic regions that appear to have swept independently through the populations prior to the study without purging diversity elsewhere in the genome. Evidence for both genome-wide and gene-specific sweeps suggests that different models of bacterial speciation may apply to different populations coexisting in the same environment.

Expanding the Limits of Thermoacidophily in the Archaeon Sulfolobus solfataricus by Adaptive Evolution.

Extremely thermoacidophilic Crenarchaeota belonging to the order Sulfolobales flourish in hot acidic habitats that are strongly oxidizing. The pH extremes of these habitats, however, often exceed the acid tolerance of type species and strains. Here, adaptive laboratory evolution was used over a 3-year period to test whether such organisms harbor additional thermoacidophilic capacity. Three distinct cell lines derived from a single type species were subjected to high-temperature serial passage while culture acidity was gradually increased. A 178-fold increase in thermoacidophily was achieved after 29 increments of shifted culture pH resulting in growth at pH 0.8 and 80°C. These strains were named super-acid-resistant Crenarchaeota (SARC). Mathematical modeling using growth parameters predicted the limits of acid resistance, while genome resequencing and transcriptome resequencing were conducted for insight into mechanisms responsible for the evolved trait. Among the mutations that were detected, a set of eight nonsynonymous changes may explain the heritability of increased acid resistance despite an unexpected lack of transposition. Four multigene components of the SARC transcriptome implicated oxidative stress as a primary challenge accompanying growth at acid extremes. These components included accelerated membrane biogenesis, induction of the mer operon, and an increased capacity for the generation of energy and reductant.

Complete Genome Sequences of Evolved Arsenate-Resistant Metallosphaera sedula Strains.

Metallosphaera sedula is a thermoacidophilic crenarchaeote with a 2.19-Mb genome. Here, we report the genome sequences of several evolved derivatives of M. sedula generated through adaptive laboratory evolution for enhanced arsenate resistance.

The molecular dimension of microbial species: 3. Comparative genomics of Synechococcus strains with different light responses and in situ diel transcription patterns of associated putative ecotypes in the Mushroom Spring microbial mat.

Genomes were obtained for three closely related strains of Synechococcus that are representative of putative ecotypes (PEs) that predominate at different depths in the 1 mm-thick, upper-green layer in the 60°C mat of Mushroom Spring, Yellowstone National Park, and exhibit different light adaptation and acclimation responses. The genomes were compared to the published genome of a previously obtained, closely related strain from a neighboring spring, and differences in both gene content and orthologous gene alleles between high-light-adapted and low-light-adapted strains were identified. Evidence of genetic differences that relate to adaptation to light intensity and/or quality, CO2uptake, nitrogen metabolism, organic carbon metabolism, and uptake of other nutrients were found between strains of the different putative ecotypes. In situ diel transcription patterns of genes, including genes unique to either low-light-adapted or high-light-adapted strains and different alleles of an orthologous photosystem gene, revealed that expression is fine-tuned to the different light environments experienced by ecotypes prevalent at various depths in the mat. This study suggests that strains of closely related PEs have different genomic adaptations that enable them to inhabit distinct ecological niches while living in close proximity within a microbial community.

Acute Limonene Toxicity in Escherichia coli Is Caused by Limonene Hydroperoxide and Alleviated by a Point Mutation in Alkyl Hydroperoxidase AhpC.

Limonene, a major component of citrus peel oil, has a number of applications related to microbiology. The antimicrobial properties of limonene make it a popular disinfectant and food preservative, while its potential as a biofuel component has made it the target of renewable production efforts through microbial metabolic engineering. For both applications, an understanding of microbial sensitivity or tolerance to limonene is crucial, but the mechanism of limonene toxicity remains enigmatic. In this study, we characterized a limonene-tolerant strain of Escherichia coli and found a mutation in ahpC, encoding alkyl hydroperoxidase, which alleviated limonene toxicity. We show that the acute toxicity previously attributed to limonene is largely due to the common oxidation product limonene hydroperoxide, which forms spontaneously in aerobic environments. The mutant AhpC protein with an L-to-Q change at position 177 (AhpC(L177Q)) was able to alleviate this toxicity by reducing the hydroperoxide to a more benign compound. We show that the degree of limonene toxicity is a function of its oxidation level and that nonoxidized limonene has relatively little toxicity to wild-type E. coli cells. Our results have implications for both the renewable production of limonene and the applications of limonene as an antimicrobial.

Complete Genome Sequence of Sulfolobus solfataricus Strain 98/2 and Evolved Derivatives.

Sulfolobus solfataricus is a thermoacidophilic crenarcheote with a 3.0-Mb genome. Here, we report the genome sequence of S. solfataricus strain 98/2, along with several evolved derivatives generated through experimental microbial evolution for enhanced thermoacidophily.