The biogeographic differentiation of algal microbiomes in the upper ocean from pole to pole

Eukaryotic phytoplankton are responsible for at least 20% of annual global carbon fixation. Their diversity and activity are shaped by interactions with prokaryotes as part of complex microbiomes. Although differences in their local species diversity have been estimated, we still have a limited understanding of environmental conditions responsible for compositional differences between local species communities on a large scale from pole to pole. Here, we show, based on pole-to-pole phytoplankton metatranscriptomes and microbial rDNA sequencing, that environmental differences between polar and non-polar upper oceans most strongly impact the large-scale spatial pattern of biodiversity and gene activity in algal microbiomes. The geographic differentiation of co-occurring microbes in algal microbiomes can be well explained by the latitudinal temperature gradient and associated break points in their beta diversity, with an average breakpoint at 14 °C ± 4.3, separating cold and warm upper oceans. As global warming impacts upper ocean temperatures, we project that break points of beta diversity move markedly pole-wards. Hence, abrupt regime shifts in algal microbiomes could be caused by anthropogenic climate change.

Succession of physiological stages hallmarks the transcriptomic response of the fungus Aspergillus niger to lignocellulose

BACKGROUND

Understanding how fungi degrade lignocellulose is a cornerstone of improving renewables-based biotechnology, in particular for the production of hydrolytic enzymes. Considerable progress has been made in investigating fungal degradation during time-points where CAZyme expression peaks. However, a robust understanding of the fungal survival strategies over its life time on lignocellulose is thereby missed. Here we aimed to uncover the physiological responses of the biotechnological workhorse and enzyme producer Aspergillus niger over its life time to six substrates important for biofuel production.

RESULTS

We analysed the response of A. niger to the feedstock Miscanthus and compared it with our previous study on wheat straw, alone or in combination with hydrothermal or ionic liquid feedstock pretreatments. Conserved (substrate-independent) metabolic responses as well as those affected by pretreatment and feedstock were identified via multivariate analysis of genome-wide transcriptomics combined with targeted transcript and protein analyses and mapping to a metabolic model. Initial exposure to all substrates increased fatty acid beta-oxidation and lipid metabolism transcripts. In a strain carrying a deletion of the ortholog of the Aspergillus nidulans fatty acid beta-oxidation transcriptional regulator farA, there was a reduction in expression of selected lignocellulose degradative CAZyme-encoding genes suggesting that beta-oxidation contributes to adaptation to lignocellulose. Mannan degradation expression was wheat straw feedstock-dependent and pectin degradation was higher on the untreated substrates. In the later life stages, known and novel secondary metabolite gene clusters were activated, which are of high interest due to their potential to synthesize bioactive compounds.

CONCLUSION

In this study, which includes the first transcriptional response of Aspergilli to Miscanthus, we highlighted that life time as well as substrate composition and structure (via variations in pretreatment and feedstock) influence the fungal responses to lignocellulose. We also demonstrated that the fungal response contains physiological stages that are conserved across substrates and are typically found outside of the conditions with high CAZyme expression, as exemplified by the stages that are dominated by lipid and secondary metabolism.

Identification of functional single nucleotide polymorphism of Populus trichocarpa PtrEPSP-TF and determination of its transcriptional effect

In plants, the phenylpropanoid pathway is responsible for the synthesis of a diverse array of secondary metabolites that include lignin monomers, flavonoids, and coumarins, many of which are essential for plant structure, biomass recalcitrance, stress defense, and nutritional quality. Our previous studies have demonstrated that Populus trichocarpa PtrEPSP-TF, an isoform of 5-enolpyruvylshikimate 3-phosphate (EPSP) synthase, has transcriptional activity and regulates phenylpropanoid biosynthesis in Populus. In this study, we report the identification of single nucleotide polymorphism (SNP) of PtrEPSP-TF that defines its functionality. Populus natural variants carrying this SNP were shown to have reduced lignin content. Here, we demonstrated that the SNP-induced substitution of 142nd amino acid (PtrEPSP-TFD142E) dramatically impairs the DNA-binding and transcriptional activity of PtrEPSP-TF. When introduced to a monocot species rice (Oryza sativa) in which an EPSP synthase isoform with the DNA-binding helix-turn-helix (HTH) motif is absent, the PtrEPSP-TF, but not PtrEPSP-TFD142E, activated genes in the phenylpropanoid pathway. More importantly, heterologous expression of PtrEPSP-TF uncovered five new transcriptional regulators of phenylpropanoid biosynthesis in rice. Collectively, this study identifies the key amino acid required for PtrEPSP-TF functionality and provides a strategy to uncover new transcriptional regulators in phenylpropanoid biosynthesis.

Metatranscriptomic Analyses of Diel Metabolic Functions During a Microcystis Bloom in Western Lake Erie (United States)

This study examined diel shifts in metabolic functions of Microcystis spp. during a 48-h Lagrangian survey of a toxin-producing cyanobacterial bloom in western Lake Erie in the aftermath of the 2014 Toledo Water Crisis. Transcripts mapped to the genomes of recently sequenced lower Great Lakes Microcystis isolates showed distinct patterns of gene expression between samples collected across day (10:00 h, 16:00 h) and night (22:00 h, 04:00 h). Daytime transcripts were enriched in functions related to Photosystem II (e.g., psbA), nitrogen and phosphate acquisition, cell division (ftsHZ), heat shock response (dnaK, groEL), and uptake of inorganic carbon (rbc, bicA). Genes transcribed during nighttime included those involved in phycobilisome protein synthesis and Photosystem I core subunits. Hierarchical clustering and principal component analysis (PCA) showed a tightly clustered group of nighttime expressed genes, whereas daytime transcripts were separated from each other over the 48-h duration. Lack of uniform clustering within the daytime transcripts suggested that the partitioning of gene expression in Microcystis is dependent on both circadian regulation and physicochemical changes within the environment.

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.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody species.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

Insights into Land Plant Evolution Garnered from the Marchantia polymorpha Genome

The evolution of land flora transformed the terrestrial environment. Land plants evolved from an ancestral charophycean alga from which they inherited developmental, biochemical, and cell biological attributes. Additional biochemical and physiological adaptations to land, and a life cycle with an alternation between multicellular haploid and diploid generations that facilitated efficient dispersal of desiccation tolerant spores, evolved in the ancestral land plant. We analyzed the genome of the liverwort Marchantia polymorpha, a member of a basal land plant lineage. Relative to charophycean algae, land plant genomes are characterized by genes encoding novel biochemical pathways, new phytohormone signaling pathways (notably auxin), expanded repertoires of signaling pathways, and increased diversity in some transcription factor families. Compared with other sequenced land plants, M. polymorpha exhibits low genetic redundancy in most regulatory pathways, with this portion of its genome resembling that predicted for the ancestral land plant. PAPERCLIP.

Comparative genomics reveals high biological diversity and specific adaptations in the industrially and medically important fungal genus Aspergillus

BACKGROUND

The fungal genus Aspergillus is of critical importance to humankind. Species include those with industrial applications, important pathogens of humans, animals and crops, a source of potent carcinogenic contaminants of food, and an important genetic model. The genome sequences of eight aspergilli have already been explored to investigate aspects of fungal biology, raising questions about evolution and specialization within this genus.

RESULTS

We have generated genome sequences for ten novel, highly diverse Aspergillus species and compared these in detail to sister and more distant genera. Comparative studies of key aspects of fungal biology, including primary and secondary metabolism, stress response, biomass degradation, and signal transduction, revealed both conservation and diversity among the species. Observed genomic differences were validated with experimental studies. This revealed several highlights, such as the potential for sex in asexual species, organic acid production genes being a key feature of black aspergilli, alternative approaches for degrading plant biomass, and indications for the genetic basis of stress response. A genome-wide phylogenetic analysis demonstrated in detail the relationship of the newly genome sequenced species with other aspergilli.

CONCLUSIONS

Many aspects of biological differences between fungal species cannot be explained by current knowledge obtained from genome sequences. The comparative genomics and experimental study, presented here, allows for the first time a genus-wide view of the biological diversity of the aspergilli and in many, but not all, cases linked genome differences to phenotype. Insights gained could be exploited for biotechnological and medical applications of fungi.

Expression of Aspergillus niger CAZymes is determined by compositional changes in wheat straw generated by hydrothermal or ionic liquid pretreatments

BACKGROUND

The capacity of fungi, such as Aspergillus niger, to degrade lignocellulose is harnessed in biotechnology to generate biofuels and high-value compounds from renewable feedstocks. Most feedstocks are currently pretreated to increase enzymatic digestibility: improving our understanding of the transcriptomic responses of fungi to pretreated lignocellulosic substrates could help to improve the mix of activities and reduce the production costs of commercial lignocellulose saccharifying cocktails.

RESULTS

We investigated the responses of A. niger to untreated, ionic liquid and hydrothermally pretreated wheat straw over a 5-day time course using RNA-seq and targeted proteomics. The ionic liquid pretreatment altered the cellulose crystallinity while retaining more of the hemicellulosic sugars than the hydrothermal pretreatment. Ionic liquid pretreatment of straw led to a dynamic induction and repression of genes, which was correlated with the higher levels of pentose sugars saccharified from the ionic liquid-pretreated straw. Hydrothermal pretreatment of straw led to reduced levels of transcripts of genes encoding carbohydrate-active enzymes as well as the derived proteins and enzyme activities. Both pretreatments abolished the expression of a large set of genes encoding pectinolytic enzymes. These reduced levels could be explained by the removal of parts of the lignocellulose by the hydrothermal pretreatment. The time course also facilitated identification of temporally limited gene induction patterns.

CONCLUSIONS

The presented transcriptomic and biochemical datasets demonstrate that pretreatments caused modifications of the lignocellulose, to both specific structural features as well as the organisation of the overall lignocellulosic structure, that determined A. niger transcript levels. The experimental setup allowed reliable detection of substrate-specific gene expression patterns as well as hitherto non-expressed genes. Our data suggest beneficial effects of using untreated and IL-pretreated straw, but not HT-pretreated straw, as feedstock for CAZyme production.

Ectomycorrhizal ecology is imprinted in the genome of the dominant symbiotic fungus Cenococcum geophilum

The most frequently encountered symbiont on tree roots is the ascomycete Cenococcum geophilum, the only mycorrhizal species within the largest fungal class Dothideomycetes, a class known for devastating plant pathogens. Here we show that the symbiotic genomic idiosyncrasies of ectomycorrhizal basidiomycetes are also present in C. geophilum with symbiosis-induced, taxon-specific genes of unknown function and reduced numbers of plant cell wall-degrading enzymes. C. geophilum still holds a significant set of genes in categories known to be involved in pathogenesis and shows an increased genome size due to transposable elements proliferation. Transcript profiling revealed a striking upregulation of membrane transporters, including aquaporin water channels and sugar transporters, and mycorrhiza-induced small secreted proteins (MiSSPs) in ectomycorrhiza compared with free-living mycelium. The frequency with which this symbiont is found on tree roots and its possible role in water and nutrient transport in symbiosis calls for further studies on mechanisms of host and environmental adaptation.

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.

Comparative genomics, proteomics and transcriptomics give new insight into the exoproteome of the basidiomycete Hebeloma cylindrosporum and its involvement in ectomycorrhizal symbiosis

Extracellular proteins play crucial roles in the interaction between mycorrhizal fungi and their environment. Computational prediction and experimental detection allowed identification of 869 proteins constituting the exoproteome of Hebeloma cylindrosporum. Small secreted proteins (SSPs) and carbohydrate-active enzymes (CAZymes) were the two major classes of extracellular proteins. Twenty-eight per cent of the SSPs were secreted by free-living mycelia and five of the 10 most abundant extracellular proteins were SSPs. By contrast, 63-75% of enzymes involved in nutrient acquisition were secreted. A total of 150 extracellular protein-coding genes were differentially expressed between mycorrhizas and free-living mycelia. SSPs were the most affected. External environmental conditions also affected expression of 199 exoproteome genes in mycorrhizas. SSPs displayed different patterns of regulation in response to presence of a host plant or other environmental signals. Several of the genes most overexpressed in the presence of organic matter encoded oxidoreductases. Hebeloma cylindrosporum has not fully lost its ancestral saprotrophic capacities but rather adapted them not to harm its hosts and to use soil organic nitrogen. The complex and divergent patterns of regulation of SSPs in response to a symbiotic partner and/or organic matter suggest various roles in the biology of mycorrhizal fungi.

Identification of a Novel Nematotoxic Protein by Challenging the Model Mushroom Coprinopsis cinerea with a Fungivorous Nematode

The dung of herbivores, the natural habitat of the model mushroom Coprinopsis cinerea, is a nutrient-rich but also very competitive environment for a saprophytic fungus. We showed previously that C. cinerea expresses constitutive, tissue-specific armories against antagonists such as animal predators and bacterial competitors. In order to dissect the inducible armories against such antagonists, we sequenced the poly(A)-positive transcriptome of C. cinerea vegetative mycelium upon challenge with fungivorous and bacterivorous nematodes, Gram-negative and Gram-positive bacteria and mechanical damage. As a response to the fungivorous nematode Aphelenchus avenae, C. cinerea was found to specifically induce the transcription of several genes encoding previously characterized nematotoxic lectins. In addition, a previously not characterized gene encoding a cytoplasmic protein with several predicted Ricin B-fold domains, was found to be strongly upregulated under this condition. Functional analysis of the recombinant protein revealed a high toxicity toward the bacterivorous nematode Caenorhabditis elegans. Challenge of the mycelium with A. avenae also lead to the induction of several genes encoding putative antibacterial proteins. Some of these genes were also induced upon challenge of the mycelium with the bacteria Escherichia coli and Bacillus subtilis. These results suggest that fungi have the ability to induce specific innate defense responses similar to plants and animals.

Analysis of the Phlebiopsis gigantea genome, transcriptome and secretome provides insight into its pioneer colonization strategies of wood

Collectively classified as white-rot fungi, certain basidiomycetes efficiently degrade the major structural polymers of wood cell walls. A small subset of these Agaricomycetes, exemplified by Phlebiopsis gigantea, is capable of colonizing freshly exposed conifer sapwood despite its high content of extractives, which retards the establishment of other fungal species. The mechanism(s) by which P. gigantea tolerates and metabolizes resinous compounds have not been explored. Here, we report the annotated P. gigantea genome and compare profiles of its transcriptome and secretome when cultured on fresh-cut versus solvent-extracted loblolly pine wood. The P. gigantea genome contains a conventional repertoire of hydrolase genes involved in cellulose/hemicellulose degradation, whose patterns of expression were relatively unperturbed by the absence of extractives. The expression of genes typically ascribed to lignin degradation was also largely unaffected. In contrast, genes likely involved in the transformation and detoxification of wood extractives were highly induced in its presence. Their products included an ABC transporter, lipases, cytochrome P450s, glutathione S-transferase and aldehyde dehydrogenase. Other regulated genes of unknown function and several constitutively expressed genes are also likely involved in P. gigantea's extractives metabolism. These results contribute to our fundamental understanding of pioneer colonization of conifer wood and provide insight into the diverse chemistries employed by fungi in carbon cycling processes.

The wood decay fungus Phlebiopsis gigantea degrades all components of plant cell walls and is uniquely able to rapidly colonize freshly exposed conifer sapwood. However, mechanisms underlying its conversion of lignocellulose and resinous extractives have not been explored. We report here analyses of the genetic repertoire, transcriptome and secretome of P. gigantea. Numerous highly expressed hydrolases, together with lytic polysaccharide monooxygenases were implicated in P. gigantea's attack on cellulose, and an array of ligninolytic peroxidases and auxiliary enzymes were also identified. Comparisons of woody substrates with and without extractives revealed differentially expressed genes predicted to be involved in the transformation of resin. These expression patterns are likely key to the pioneer colonization of conifers by P. gigantea.

A near complete snapshot of the Zea mays seedling transcriptome revealed from ultra-deep sequencing

RNA-sequencing (RNA-seq) enables in-depth exploration of transcriptomes, but typical sequencing depth often limits its comprehensiveness. In this study, we generated nearly 3 billion RNA-Seq reads, totaling 341 Gb of sequence, from a Zea mays seedling sample. At this depth, a near complete snapshot of the transcriptome was observed consisting of over 90% of the annotated transcripts, including lowly expressed transcription factors. A novel hybrid strategy combining de novo and reference-based assemblies yielded a transcriptome consisting of 126,708 transcripts with 88% of expressed known genes assembled to full-length. We improved current annotations by adding 4,842 previously unannotated transcript variants and many new features, including 212 maize transcripts, 201 genes, 10 genes with undocumented potential roles in seedlings as well as maize lineage specific gene fusion events. We demonstrated the power of deep sequencing for large transcriptome studies by generating a high quality transcriptome, which provides a rich resource for the research community.

Insights into the maize pan-genome and pan-transcriptome

Genomes at the species level are dynamic, with genes present in every individual (core) and genes in a subset of individuals (dispensable) that collectively constitute the pan-genome. Using transcriptome sequencing of seedling RNA from 503 maize (Zea mays) inbred lines to characterize the maize pan-genome, we identified 8681 representative transcript assemblies (RTAs) with 16.4% expressed in all lines and 82.7% expressed in subsets of the lines. Interestingly, with linkage disequilibrium mapping, 76.7% of the RTAs with at least one single nucleotide polymorphism (SNP) could be mapped to a single genetic position, distributed primarily throughout the nonpericentromeric portion of the genome. Stepwise iterative clustering of RTAs suggests, within the context of the genotypes used in this study, that the maize genome is restricted and further sampling of seedling RNA within this germplasm base will result in minimal discovery. Genome-wide association studies based on SNPs and transcript abundance in the pan-genome revealed loci associated with the timing of the juvenile-to-adult vegetative and vegetative-to-reproductive developmental transitions, two traits important for fitness and adaptation. This study revealed the dynamic nature of the maize pan-genome and demonstrated that a substantial portion of variation may lie outside the single reference genome for a species.

Development of an integrated transcript sequence database and a gene expression atlas for gene discovery and analysis in switchgrass (Panicum virgatum L.)

Switchgrass (Panicum virgatum L.) is a perennial C4 grass with the potential to become a major bioenergy crop. To help realize this potential, a set of RNA-based resources were developed. Expressed sequence tags (ESTs) were generated from two tetraploid switchgrass genotypes, Alamo AP13 and Summer VS16. Over 11.5 million high-quality ESTs were generated with 454 sequencing technology, and an additional 169 079 Sanger sequences were obtained from the 5' and 3' ends of 93 312 clones from normalized, full-length-enriched cDNA libraries. AP13 and VS16 ESTs were assembled into 77 854 and 30 524 unique transcripts (unitranscripts), respectively, using the Newbler and pave programs. Published Sanger-ESTs (544 225) from Alamo, Kanlow, and 15 other cultivars were integrated with the AP13 and VS16 assemblies to create a universal switchgrass gene index (PviUT1.2) with 128 058 unitranscripts, which were annotated for function. An Affymetrix cDNA microarray chip (Pvi_cDNAa520831) containing 122 973 probe sets was designed from PviUT1.2 sequences, and used to develop a Gene Expression Atlas for switchgrass (PviGEA). The PviGEA contains quantitative transcript data for all major organ systems of switchgrass throughout development. We developed a web server that enables flexible, multifaceted analyses of PviGEA transcript data. The PviGEA was used to identify representatives of all known genes in the phenylpropanoid-monolignol biosynthesis pathway.

Comparative genome structure, secondary metabolite, and effector coding capacity across Cochliobolus pathogens

The genomes of five Cochliobolus heterostrophus strains, two Cochliobolus sativus strains, three additional Cochliobolus species (Cochliobolus victoriae, Cochliobolus carbonum, Cochliobolus miyabeanus), and closely related Setosphaeria turcica were sequenced at the Joint Genome Institute (JGI). The datasets were used to identify SNPs between strains and species, unique genomic regions, core secondary metabolism genes, and small secreted protein (SSP) candidate effector encoding genes with a view towards pinpointing structural elements and gene content associated with specificity of these closely related fungi to different cereal hosts. Whole-genome alignment shows that three to five percent of each genome differs between strains of the same species, while a quarter of each genome differs between species. On average, SNP counts among field isolates of the same C. heterostrophus species are more than 25× higher than those between inbred lines and 50× lower than SNPs between Cochliobolus species. The suites of nonribosomal peptide synthetase (NRPS), polyketide synthase (PKS), and SSP-encoding genes are astoundingly diverse among species but remarkably conserved among isolates of the same species, whether inbred or field strains, except for defining examples that map to unique genomic regions. Functional analysis of several strain-unique PKSs and NRPSs reveal a strong correlation with a role in virulence.

Porphyra (Bangiophyceae) Transcriptomes Provide Insights Into Red Algal Development And Metabolism

The red seaweed Porphyra (Bangiophyceae) and related Bangiales have global economic importance. Here, we report the analysis of a comprehensive transcriptome comprising ca. 4.7 million expressed sequence tag (EST) reads from P. umbilicalis (L.) J. Agardh and P. purpurea (Roth) C. Agardh (ca. 980 Mbp of data generated using 454 FLX pyrosequencing). These ESTs were isolated from the haploid gametophyte (blades from both species) and diploid conchocelis stage (from P. purpurea). In a bioinformatic analysis, only 20% of the contigs were found to encode proteins of known biological function. Comparative analysis of predicted protein functions in mesophilic (including Porphyra) and extremophilic red algae suggest that the former has more putative functions related to signaling, membrane transport processes, and establishment of protein complexes. These enhanced functions may reflect general mesophilic adaptations. A near-complete repertoire of genes encoding histones and ribosomal proteins was identified, with some differentially regulated between the blade and conchocelis stage in P. purpurea. This finding may reflect specific regulatory processes associated with these distinct phases of the life history. Fatty acid desaturation patterns, in combination with gene expression profiles, demonstrate differences from seed plants with respect to the transport of fatty acid/lipid among subcellular compartments and the molecular machinery of lipid assembly. We also recovered a near-complete gene repertoire for enzymes involved in the formation of sterols and carotenoids, including candidate genes for the biosynthesis of lutein. Our findings provide key insights into the evolution, development, and biology of Porphyra, an important lineage of red algae.

The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes

Wood is a major pool of organic carbon that is highly resistant to decay, owing largely to the presence of lignin. The only organisms capable of substantial lignin decay are white rot fungi in the Agaricomycetes, which also contains non-lignin-degrading brown rot and ectomycorrhizal species. Comparative analyses of 31 fungal genomes (12 generated for this study) suggest that lignin-degrading peroxidases expanded in the lineage leading to the ancestor of the Agaricomycetes, which is reconstructed as a white rot species, and then contracted in parallel lineages leading to brown rot and mycorrhizal species. Molecular clock analyses suggest that the origin of lignin degradation might have coincided with the sharp decrease in the rate of organic carbon burial around the end of the Carboniferous period.

Comparative genomics of the white-rot fungi, Phanerochaete carnosa and P. chrysosporium, to elucidate the genetic basis of the distinct wood types they colonize

Background

Softwood is the predominant form of land plant biomass in the Northern hemisphere, and is among the most recalcitrant biomass resources to bioprocess technologies. The white rot fungus, Phanerochaete carnosa, has been isolated almost exclusively from softwoods, while most other known white-rot species, including Phanerochaete chrysosporium, were mainly isolated from hardwoods. Accordingly, it is anticipated that P. carnosa encodes a distinct set of enzymes and proteins that promote softwood decomposition. To elucidate the genetic basis of softwood bioconversion by a white-rot fungus, the present study reports the P. carnosa genome sequence and its comparative analysis with the previously reported P. chrysosporium genome.

Results

P. carnosa encodes a complete set of lignocellulose-active enzymes. Comparative genomic analysis revealed that P. carnosa is enriched with genes encoding manganese peroxidase, and that the most divergent glycoside hydrolase families were predicted to encode hemicellulases and glycoprotein degrading enzymes. Most remarkably, P. carnosa possesses one of the largest P450 contingents (266 P450s) among the sequenced and annotated wood-rotting basidiomycetes, nearly double that of P. chrysosporium. Along with metabolic pathway modeling, comparative growth studies on model compounds and chemical analyses of decomposed wood components showed greater tolerance of P. carnosa to various substrates including coniferous heartwood.

Conclusions

The P. carnosa genome is enriched with genes that encode P450 monooxygenases that can participate in extractives degradation, and manganese peroxidases involved in lignin degradation. The significant expansion of P450s in P. carnosa, along with differences in carbohydrate- and lignin-degrading enzymes, could be correlated to the utilization of heartwood and sapwood preparations from both coniferous and hardwood species.

MAJOR DEVELOPMENTAL REGULATORS AND THEIR EXPRESSION IN TWO CLOSELY RELATED SPECIES OF PORPHYRA (RHODOPHYTA)(1)

Little is known about the genetic and biochemical mechanisms that underlie red algal development, for example, why the group failed to evolve complex parenchyma and tissue differentiation. Here we examined expressed sequence tag (EST) data from two closely related species, Porphyra umbilicalis (L.) J. Agardh and P. purpurea (Roth) C. Agardh, for conserved developmental regulators known from model eukaryotes, and their expression levels in several developmental stages. Genes for most major developmental families were present, including MADS-box and homeodomain (HD) proteins, SNF2 chromatin-remodelers, and proteins involved in sRNA biogenesis. Some of these genes displayed altered expression correlating with different life history stages or cell types. Notably, two ESTs encoding HD proteins showed eightfold higher expression in the P. purpurea sporophyte (conchocelis) than in the gametophyte (blade), whereas two MADS domain-containing paralogs showed significantly different patterns of expression in the conchocelis and blade respectively. These developmental gene families do not appear to have undergone the kinds of dramatic expansions in copy number found in multicellular land plants and animals, which are important for regulating developmental processes in those groups. Analyses of small RNAs did not validate the presence of miRNAs, but homologs of Argonaute were present. In general, it appears that red algae began with a similar molecular toolkit for directing development as did other multicellular eukaryotes, but probably evolved altered roles for many key proteins, as well as novel mechanisms yet to be discovered.

Insight into trade-off between wood decay and parasitism from the genome of a fungal forest pathogen

Parasitism and saprotrophic wood decay are two fungal strategies fundamental for succession and nutrient cycling in forest ecosystems. An opportunity to assess the trade-off between these strategies is provided by the forest pathogen and wood decayer Heterobasidion annosum sensu lato. We report the annotated genome sequence and transcript profiling, as well as the quantitative trait loci mapping, of one member of the species complex: H. irregulare. Quantitative trait loci critical for pathogenicity, and rich in transposable elements, orphan and secreted genes, were identified. A wide range of cellulose-degrading enzymes are expressed during wood decay. By contrast, pathogenic interaction between H. irregulare and pine engages fewer carbohydrate-active enzymes, but involves an increase in pectinolytic enzymes, transcription modules for oxidative stress and secondary metabolite production. Our results show a trade-off in terms of constrained carbohydrate decomposition and membrane transport capacity during interaction with living hosts. Our findings establish that saprotrophic wood decay and necrotrophic parasitism involve two distinct, yet overlapping, processes.

The ectomycorrhizal fungus Paxillus involutus converts organic matter in plant litter using a trimmed brown-rot mechanism involving Fenton chemistry

Soils in boreal forests contain large stocks of carbon. Plants are the main source of this carbon through tissue residues and root exudates. A major part of the exudates are allocated to symbiotic ectomycorrhizal fungi. In return, the plant receives nutrients, in particular nitrogen from the mycorrhizal fungi. To capture the nitrogen, the fungi must at least partly disrupt the recalcitrant organic matter-protein complexes within which the nitrogen is embedded. This disruption process is poorly characterized. We used spectroscopic analyses and transcriptome profiling to examine the mechanism by which the ectomycorrhizal fungus Paxillus involutus degrades organic matter when acquiring nitrogen from plant litter. The fungus partially degraded polysaccharides and modified the structure of polyphenols. The observed chemical changes were consistent with a hydroxyl radical attack, involving Fenton chemistry similar to that of brown-rot fungi. The set of enzymes expressed by Pa. involutus during the degradation of the organic matter was similar to the set of enzymes involved in the oxidative degradation of wood by brown-rot fungi. However, Pa. involutus lacked transcripts encoding extracellular enzymes needed for metabolizing the released carbon. The saprotrophic activity has been reduced to a radical-based biodegradation system that can efficiently disrupt the organic matter-protein complexes and thereby mobilize the entrapped nutrients. We suggest that the released carbon then becomes available for further degradation and assimilation by commensal microbes, and that these activities have been lost in ectomycorrhizal fungi as an adaptation to symbiotic growth on host photosynthate. The interdependence of ectomycorrhizal symbionts and saprophytic microbes would provide a key link in the turnover of nutrients and carbon in forest ecosystems.

The genome of the polar eukaryotic microalga Coccomyxa subellipsoidea reveals traits of cold adaptation

Background

Little is known about the mechanisms of adaptation of life to the extreme environmental conditions encountered in polar regions. Here we present the genome sequence of a unicellular green alga from the division chlorophyta, Coccomyxa subellipsoidea C-169, which we will hereafter refer to as C-169. This is the first eukaryotic microorganism from a polar environment to have its genome sequenced.

Results

The 48.8 Mb genome contained in 20 chromosomes exhibits significant synteny conservation with the chromosomes of its relatives Chlorella variabilis and Chlamydomonas reinhardtii. The order of the genes is highly reshuffled within synteny blocks, suggesting that intra-chromosomal rearrangements were more prevalent than inter-chromosomal rearrangements. Remarkably, Zepp retrotransposons occur in clusters of nested elements with strictly one cluster per chromosome probably residing at the centromere. Several protein families overrepresented in C. subellipsoidae include proteins involved in lipid metabolism, transporters, cellulose synthases and short alcohol dehydrogenases. Conversely, C-169 lacks proteins that exist in all other sequenced chlorophytes, including components of the glycosyl phosphatidyl inositol anchoring system, pyruvate phosphate dikinase and the photosystem 1 reaction center subunit N (PsaN).

Conclusions

We suggest that some of these gene losses and gains could have contributed to adaptation to low temperatures. Comparison of these genomic features with the adaptive strategies of psychrophilic microbes suggests that prokaryotes and eukaryotes followed comparable evolutionary routes to adapt to cold environments.

The genome of wine yeast Dekkera bruxellensis provides a tool to explore its food-related properties

The yeast Dekkera/Brettanomyces bruxellensis can cause enormous economic losses in wine industry due to production of phenolic off-flavor compounds. D. bruxellensis is a distant relative of baker's yeast Saccharomyces cerevisiae. Nevertheless, these two yeasts are often found in the same habitats and share several food-related traits, such as production of high ethanol levels and ability to grow without oxygen. In some food products, like lambic beer, D. bruxellensis can importantly contribute to flavor development. We determined the 13.4 Mb genome sequence of the D. bruxellensis strain Y879 (CBS2499) and deduced the genetic background of several "food-relevant" properties and evolutionary history of this yeast. Surprisingly, we find that this yeast is phylogenetically distant to other food-related yeasts and most related to Pichia (Komagataella) pastoris, which is an aerobic poor ethanol producer. We further show that the D. bruxellensis genome does not contain an excess of lineage specific duplicated genes nor a horizontally transferred URA1 gene, two crucial events that promoted the evolution of the food relevant traits in the S. cerevisiae lineage. However, D. bruxellensis has several independently duplicated ADH and ADH-like genes, which are likely responsible for metabolism of alcohols, including ethanol, and also a range of aromatic compounds.

The genome of the xerotolerant mold Wallemia sebi reveals adaptations to osmotic stress and suggests cryptic sexual reproduction

Wallemia (Wallemiales, Wallemiomycetes) is a genus of xerophilic Fungi of uncertain phylogenetic position within Basidiomycota. Most commonly found as food contaminants, species of Wallemia have also been isolated from hypersaline environments. The ability to tolerate environments with reduced water activity is rare in Basidiomycota. We sequenced the genome of W. sebi in order to understand its adaptations for surviving in osmotically challenging environments, and we performed phylogenomic and ultrastructural analyses to address its systematic placement and reproductive biology. W. sebi has a compact genome (9.8 Mb), with few repeats and the largest fraction of genes with functional domains compared with other Basidiomycota. We applied several approaches to searching for osmotic stress-related proteins. In silico analyses identified 93 putative osmotic stress proteins; homology searches showed the HOG (High Osmolarity Glycerol) pathway to be mostly conserved. Despite the seemingly reduced genome, several gene family expansions and a high number of transporters (549) were found that also provide clues to the ability of W. sebi to colonize harsh environments. Phylogenetic analyses of a 71-protein dataset support the position of Wallemia as the earliest diverging lineage of Agaricomycotina, which is confirmed by septal pore ultrastructure that shows the septal pore apparatus as a variant of the Tremella-type. Mating type gene homologs were identified although we found no evidence of meiosis during conidiogenesis, suggesting there may be aspects of the life cycle of W. sebi that remain cryptic.

The transcriptome of the arbuscular mycorrhizal fungus Glomus intraradices (DAOM 197198) reveals functional tradeoffs in an obligate symbiont

• The arbuscular mycorrhizal symbiosis is arguably the most ecologically important eukaryotic symbiosis, yet it is poorly understood at the molecular level. To provide novel insights into the molecular basis of symbiosis-associated traits, we report the first genome-wide analysis of the transcriptome from Glomus intraradices DAOM 197198. • We generated a set of 25,906 nonredundant virtual transcripts (NRVTs) transcribed in germinated spores, extraradical mycelium and symbiotic roots using Sanger and 454 sequencing. NRVTs were used to construct an oligoarray for investigating gene expression. • We identified transcripts coding for the meiotic recombination machinery, as well as meiosis-specific proteins, suggesting that the lack of a known sexual cycle in G. intraradices is not a result of major deletions of genes essential for sexual reproduction and meiosis. Induced expression of genes encoding membrane transporters and small secreted proteins in intraradical mycelium, together with the lack of expression of hydrolytic enzymes acting on plant cell wall polysaccharides, are all features of G. intraradices that are shared with ectomycorrhizal symbionts and obligate biotrophic pathogens. • Our results illuminate the genetic basis of symbiosis-related traits of the most ancient lineage of plant biotrophs, advancing future research on these agriculturally and ecologically important symbionts.

Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris

Thermostable enzymes and thermophilic cell factories may afford economic advantages in the production of many chemicals and biomass-based fuels. Here we describe and compare the genomes of two thermophilic fungi, Myceliophthora thermophila and Thielavia terrestris. To our knowledge, these genomes are the first described for thermophilic eukaryotes and the first complete telomere-to-telomere genomes for filamentous fungi. Genome analyses and experimental data suggest that both thermophiles are capable of hydrolyzing all major polysaccharides found in biomass. Examination of transcriptome data and secreted proteins suggests that the two fungi use shared approaches in the hydrolysis of cellulose and xylan but distinct mechanisms in pectin degradation. Characterization of the biomass-hydrolyzing activity of recombinant enzymes suggests that these organisms are highly efficient in biomass decomposition at both moderate and high temperatures. Furthermore, we present evidence suggesting that aside from representing a potential reservoir of thermostable enzymes, thermophilic fungi are amenable to manipulation using classical and molecular genetics.

The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi

Brown rot decay removes cellulose and hemicellulose from wood--residual lignin contributing up to 30% of forest soil carbon--and is derived from an ancestral white rot saprotrophy in which both lignin and cellulose are decomposed. Comparative and functional genomics of the "dry rot" fungus Serpula lacrymans, derived from forest ancestors, demonstrated that the evolution of both ectomycorrhizal biotrophy and brown rot saprotrophy were accompanied by reductions and losses in specific protein families, suggesting adaptation to an intercellular interaction with plant tissue. Transcriptome and proteome analysis also identified differences in wood decomposition in S. lacrymans relative to the brown rot Postia placenta. Furthermore, fungal nutritional mode diversification suggests that the boreal forest biome originated via genetic coevolution of above- and below-ground biota.

Rapid evolution of coral proteins responsible for interaction with the environment

Background

Corals worldwide are in decline due to climate change effects (e.g., rising seawater temperatures), pollution, and exploitation. The ability of corals to cope with these stressors in the long run depends on the evolvability of the underlying genetic networks and proteins, which remain largely unknown. A genome-wide scan for positively selected genes between related coral species can help to narrow down the search space considerably.

Methodology/Principal Findings

We screened a set of 2,604 putative orthologs from EST-based sequence datasets of the coral species Acropora millepora and Acropora palmata to determine the fraction and identity of proteins that may experience adaptive evolution. 7% of the orthologs show elevated rates of evolution. Taxonomically-restricted (i.e. lineage-specific) genes show a positive selection signature more frequently than genes that are found across many animal phyla. The class of proteins that displayed elevated evolutionary rates was significantly enriched for proteins involved in immunity and defense, reproduction, and sensory perception. We also found elevated rates of evolution in several other functional groups such as management of membrane vesicles, transmembrane transport of ions and organic molecules, cell adhesion, and oxidative stress response. Proteins in these processes might be related to the endosymbiotic relationship corals maintain with dinoflagellates in the genus Symbiodinium.

Conclusion/Relevance

This study provides a birds-eye view of the processes potentially underlying coral adaptation, which will serve as a foundation for future work to elucidate the rates, patterns, and mechanisms of corals' evolutionary response to global climate change.

The Selaginella genome identifies genetic changes associated with the evolution of vascular plants

Vascular plants appeared ~410 million years ago, then diverged into several lineages of which only two survive: the euphyllophytes (ferns and seed plants) and the lycophytes. We report here the genome sequence of the lycophyte Selaginella moellendorffii (Selaginella), the first nonseed vascular plant genome reported. By comparing gene content in evolutionarily diverse taxa, we found that the transition from a gametophyte- to a sporophyte-dominated life cycle required far fewer new genes than the transition from a nonseed vascular to a flowering plant, whereas secondary metabolic genes expanded extensively and in parallel in the lycophyte and angiosperm lineages. Selaginella differs in posttranscriptional gene regulation, including small RNA regulation of repetitive elements, an absence of the trans-acting small interfering RNA pathway, and extensive RNA editing of organellar genes.

Comparative genomics of citric-acid-producing Aspergillus niger ATCC 1015 versus enzyme-producing CBS 513.88

The filamentous fungus Aspergillus niger exhibits great diversity in its phenotype. It is found globally, both as marine and terrestrial strains, produces both organic acids and hydrolytic enzymes in high amounts, and some isolates exhibit pathogenicity. Although the genome of an industrial enzyme-producing A. niger strain (CBS 513.88) has already been sequenced, the versatility and diversity of this species compel additional exploration. We therefore undertook whole-genome sequencing of the acidogenic A. niger wild-type strain (ATCC 1015) and produced a genome sequence of very high quality. Only 15 gaps are present in the sequence, and half the telomeric regions have been elucidated. Moreover, sequence information from ATCC 1015 was used to improve the genome sequence of CBS 513.88. Chromosome-level comparisons uncovered several genome rearrangements, deletions, a clear case of strain-specific horizontal gene transfer, and identification of 0.8 Mb of novel sequence. Single nucleotide polymorphisms per kilobase (SNPs/kb) between the two strains were found to be exceptionally high (average: 7.8, maximum: 160 SNPs/kb). High variation within the species was confirmed with exo-metabolite profiling and phylogenetics. Detailed lists of alleles were generated, and genotypic differences were observed to accumulate in metabolic pathways essential to acid production and protein synthesis. A transcriptome analysis supported up-regulation of genes associated with biosynthesis of amino acids that are abundant in glucoamylase A, tRNA-synthases, and protein transporters in the protein producing CBS 513.88 strain. Our results and data sets from this integrative systems biology analysis resulted in a snapshot of fungal evolution and will support further optimization of cell factories based on filamentous fungi.

Comparative genome sequence analysis underscores mycoparasitism as the ancestral life style of Trichoderma

BACKGROUND

Mycoparasitism, a lifestyle where one fungus is parasitic on another fungus, has special relevance when the prey is a plant pathogen, providing a strategy for biological control of pests for plant protection. Probably, the most studied biocontrol agents are species of the genus Hypocrea/Trichoderma.

RESULTS

Here we report an analysis of the genome sequences of the two biocontrol species Trichoderma atroviride (teleomorph Hypocrea atroviridis) and Trichoderma virens (formerly Gliocladium virens, teleomorph Hypocrea virens), and a comparison with Trichoderma reesei (teleomorph Hypocrea jecorina). These three Trichoderma species display a remarkable conservation of gene order (78 to 96%), and a lack of active mobile elements probably due to repeat-induced point mutation. Several gene families are expanded in the two mycoparasitic species relative to T. reesei or other ascomycetes, and are overrepresented in non-syntenic genome regions. A phylogenetic analysis shows that T. reesei and T. virens are derived relative to T. atroviride. The mycoparasitism-specific genes thus arose in a common Trichoderma ancestor but were subsequently lost in T. reesei.

CONCLUSIONS

The data offer a better understanding of mycoparasitism, and thus enforce the development of improved biocontrol strains for efficient and environmentally friendly protection of plants.

Discovery and annotation of small proteins using genomics, proteomics, and computational approaches

Small proteins (10-200 amino acids [aa] in length) encoded by short open reading frames (sORF) play important regulatory roles in various biological processes, including tumor progression, stress response, flowering, and hormone signaling. However, ab initio discovery of small proteins has been relatively overlooked. Recent advances in deep transcriptome sequencing make it possible to efficiently identify sORFs at the genome level. In this study, we obtained ~2.6 million expressed sequence tag (EST) reads from Populus deltoides leaf transcriptome and reconstructed full-length transcripts from the EST sequences. We identified an initial set of 12,852 sORFs encoding proteins of 10-200 aa in length. Three computational approaches were then used to enrich for bona fide protein-coding sORFs from the initial sORF set: (1) coding-potential prediction, (2) evolutionary conservation between P. deltoides and other plant species, and (3) gene family clustering within P. deltoides. As a result, a high-confidence sORF candidate set containing 1469 genes was obtained. Analysis of the protein domains, non-protein-coding RNA motifs, sequence length distribution, and protein mass spectrometry data supported this high-confidence sORF set. In the high-confidence sORF candidate set, known protein domains were identified in 1282 genes (higher-confidence sORF candidate set), out of which 611 genes, designated as highest-confidence candidate sORF set, were supported by proteomics data. Of the 611 highest-confidence candidate sORF genes, 56 were new to the current Populus genome annotation. This study not only demonstrates that there are potential sORF candidates to be annotated in sequenced genomes, but also presents an efficient strategy for discovery of sORFs in species with no genome annotation yet available.

Comparative genomics of the social amoebae Dictyostelium discoideum and Dictyostelium purpureum

Background

The social amoebae (Dictyostelia) are a diverse group of Amoebozoa that achieve multicellularity by aggregation and undergo morphogenesis into fruiting bodies with terminally differentiated spores and stalk cells. There are four groups of dictyostelids, with the most derived being a group that contains the model species Dictyostelium discoideum.

Results

We have produced a draft genome sequence of another group dictyostelid, Dictyostelium purpureum, and compare it to the D. discoideum genome. The assembly (8.41 × coverage) comprises 799 scaffolds totaling 33.0 Mb, comparable to the D. discoideum genome size. Sequence comparisons suggest that these two dictyostelids shared a common ancestor approximately 400 million years ago. In spite of this divergence, most orthologs reside in small clusters of conserved synteny. Comparative analyses revealed a core set of orthologous genes that illuminate dictyostelid physiology, as well as differences in gene family content. Interesting patterns of gene conservation and divergence are also evident, suggesting function differences; some protein families, such as the histidine kinases, have undergone little functional change, whereas others, such as the polyketide synthases, have undergone extensive diversification. The abundant amino acid homopolymers encoded in both genomes are generally not found in homologous positions within proteins, so they are unlikely to derive from ancestral DNA triplet repeats. Genes involved in the social stage evolved more rapidly than others, consistent with either relaxed selection or accelerated evolution due to social conflict.

Conclusions

The findings from this new genome sequence and comparative analysis shed light on the biology and evolution of the Dictyostelia.

The ecoresponsive genome of Daphnia pulex

We describe the draft genome of the microcrustacean Daphnia pulex, which is only 200 megabases and contains at least 30,907 genes. The high gene count is a consequence of an elevated rate of gene duplication resulting in tandem gene clusters. More than a third of Daphnia's genes have no detectable homologs in any other available proteome, and the most amplified gene families are specific to the Daphnia lineage. The coexpansion of gene families interacting within metabolic pathways suggests that the maintenance of duplicated genes is not random, and the analysis of gene expression under different environmental conditions reveals that numerous paralogs acquire divergent expression patterns soon after duplication. Daphnia-specific genes, including many additional loci within sequenced regions that are otherwise devoid of annotations, are the most responsive genes to ecological challenges.

The ecoresponsive genome of Daphnia pulex

We describe the draft genome of the microcrustacean Daphnia pulex, which is only 200 megabases and contains at least 30,907 genes. The high gene count is a consequence of an elevated rate of gene duplication resulting in tandem gene clusters. More than a third of Daphnia's genes have no detectable homologs in any other available proteome, and the most amplified gene families are specific to the Daphnia lineage. The coexpansion of gene families interacting within metabolic pathways suggests that the maintenance of duplicated genes is not random, and the analysis of gene expression under different environmental conditions reveals that numerous paralogs acquire divergent expression patterns soon after duplication. Daphnia-specific genes, including many additional loci within sequenced regions that are otherwise devoid of annotations, are the most responsive genes to ecological challenges.

The Chlorella variabilis NC64A genome reveals adaptation to photosymbiosis, coevolution with viruses, and cryptic sex

Chlorella variabilis NC64A, a unicellular photosynthetic green alga (Trebouxiophyceae), is an intracellular photobiont of Paramecium bursaria and a model system for studying virus/algal interactions. We sequenced its 46-Mb nuclear genome, revealing an expansion of protein families that could have participated in adaptation to symbiosis. NC64A exhibits variations in GC content across its genome that correlate with global expression level, average intron size, and codon usage bias. Although Chlorella species have been assumed to be asexual and nonmotile, the NC64A genome encodes all the known meiosis-specific proteins and a subset of proteins found in flagella. We hypothesize that Chlorella might have retained a flagella-derived structure that could be involved in sexual reproduction. Furthermore, a survey of phytohormone pathways in chlorophyte algae identified algal orthologs of Arabidopsis thaliana genes involved in hormone biosynthesis and signaling, suggesting that these functions were established prior to the evolution of land plants. We show that the ability of Chlorella to produce chitinous cell walls likely resulted from the capture of metabolic genes by horizontal gene transfer from algal viruses, prokaryotes, or fungi. Analysis of the NC64A genome substantially advances our understanding of the green lineage evolution, including the genomic interplay with viruses and symbiosis between eukaryotes.

Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri

The multicellular green alga Volvox carteri and its morphologically diverse close relatives (the volvocine algae) are well suited for the investigation of the evolution of multicellularity and development. We sequenced the 138-mega-base pair genome of V. carteri and compared its approximately 14,500 predicted proteins to those of its unicellular relative Chlamydomonas reinhardtii. Despite fundamental differences in organismal complexity and life history, the two species have similar protein-coding potentials and few species-specific protein-coding gene predictions. Volvox is enriched in volvocine-algal-specific proteins, including those associated with an expanded and highly compartmentalized extracellular matrix. Our analysis shows that increases in organismal complexity can be associated with modifications of lineage-specific proteins rather than large-scale invention of protein-coding capacity.

The genome of the Western clawed frog Xenopus tropicalis

The western clawed frog Xenopus tropicalis is an important model for vertebrate development that combines experimental advantages of the African clawed frog Xenopus laevis with more tractable genetics. Here we present a draft genome sequence assembly of X. tropicalis. This genome encodes more than 20,000 protein-coding genes, including orthologs of at least 1700 human disease genes. Over 1 million expressed sequence tags validated the annotation. More than one-third of the genome consists of transposable elements, with unusually prevalent DNA transposons. Like that of other tetrapods, the genome of X. tropicalis contains gene deserts enriched for conserved noncoding elements. The genome exhibits substantial shared synteny with human and chicken over major parts of large chromosomes, broken by lineage-specific chromosome fusions and fissions, mainly in the mammalian lineage.

The genome of Naegleria gruberi illuminates early eukaryotic versatility

Genome sequences of diverse free-living protists are essential for understanding eukaryotic evolution and molecular and cell biology. The free-living amoeboflagellate Naegleria gruberi belongs to a varied and ubiquitous protist clade (Heterolobosea) that diverged from other eukaryotic lineages over a billion years ago. Analysis of the 15,727 protein-coding genes encoded by Naegleria's 41 Mb nuclear genome indicates a capacity for both aerobic respiration and anaerobic metabolism with concomitant hydrogen production, with fundamental implications for the evolution of organelle metabolism. The Naegleria genome facilitates substantially broader phylogenomic comparisons of free-living eukaryotes than previously possible, allowing us to identify thousands of genes likely present in the pan-eukaryotic ancestor, with 40% likely eukaryotic inventions. Moreover, we construct a comprehensive catalog of amoeboid-motility genes. The Naegleria genome, analyzed in the context of other protists, reveals a remarkably complex ancestral eukaryote with a rich repertoire of cytoskeletal, sexual, signaling, and metabolic modules.

The porcelain crab transcriptome and PCAD, the porcelain crab microarray and sequence database

Background

With the emergence of a completed genome sequence of the freshwater crustacean Daphnia pulex, construction of genomic-scale sequence databases for additional crustacean sequences are important for comparative genomics and annotation. Porcelain crabs, genus Petrolisthes, have been powerful crustacean models for environmental and evolutionary physiology with respect to thermal adaptation and understanding responses of marine organisms to climate change. Here, we present a large-scale EST sequencing and cDNA microarray database project for the porcelain crab Petrolisthes cinctipes.

Methodology/Principal Findings

A set of ∼30K unique sequences (UniSeqs) representing ∼19K clusters were generated from ∼98K high quality ESTs from a set of tissue specific non-normalized and mixed-tissue normalized cDNA libraries from the porcelain crab Petrolisthes cinctipes. Homology for each UniSeq was assessed using BLAST, InterProScan, GO and KEGG database searches. Approximately 66% of the UniSeqs had homology in at least one of the databases. All EST and UniSeq sequences along with annotation results and coordinated cDNA microarray datasets have been made publicly accessible at the Porcelain Crab Array Database (PCAD), a feature-enriched version of the Stanford and Longhorn Array Databases.

Conclusions/Significance

The EST project presented here represents the third largest sequencing effort for any crustacean, and the largest effort for any crab species. Our assembly and clustering results suggest that our porcelain crab EST data set is equally diverse to the much larger EST set generated in the Daphnia pulex genome sequencing project, and thus will be an important resource to the Daphnia research community. Our homology results support the pancrustacea hypothesis and suggest that Malacostraca may be ancestral to Branchiopoda and Hexapoda. Our results also suggest that our cDNA microarrays cover as much of the transcriptome as can reasonably be captured in EST library sequencing approaches, and thus represent a rich resource for studies of environmental genomics.

Assembly of 500,000 inter-specific catfish expressed sequence tags and large scale gene-associated marker development for whole genome association studies

Twelve cDNA libraries from two species of catfish have been sequenced, resulting in the generation of nearly 500,000 ESTs.

Background

Through the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energy's Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification.

Results

A total of 438,321 quality ESTs were generated from 8 channel catfish (Ictalurus punctatus) and 4 blue catfish (Ictalurus furcatus) libraries, bringing the number of catfish ESTs to nearly 500,000. Assembly of all catfish ESTs resulted in 45,306 contigs and 66,272 singletons. Over 35% of the unique sequences had significant similarities to known genes, allowing the identification of 14,776 unique genes in catfish. Over 300,000 putative SNPs have been identified, of which approximately 48,000 are high-quality SNPs identified from contigs with at least four sequences and the minor allele presence of at least two sequences in the contig. The EST resource should be valuable for identification of microsatellites, genome annotation, large-scale expression analysis, and comparative genome analysis.

Conclusions

This project generated a large EST resource for catfish that captured the majority of the catfish transcriptome. The parallel analysis of ESTs from two closely related Ictalurid catfishes should also provide powerful means for the evaluation of ancient and recent gene duplications, and for the development of high-density microarrays in catfish. The inter- and intra-specific SNPs identified from all catfish EST dataset assembly will greatly benefit the catfish introgression breeding program and whole genome association studies.

Genome sequence of the palaeopolyploid soybean

Soybean (Glycine max) is one of the most important crop plants for seed protein and oil content, and for its capacity to fix atmospheric nitrogen through symbioses with soil-borne microorganisms. We sequenced the 1.1-gigabase genome by a whole-genome shotgun approach and integrated it with physical and high-density genetic maps to create a chromosome-scale draft sequence assembly. We predict 46,430 protein-coding genes, 70% more than Arabidopsis and similar to the poplar genome which, like soybean, is an ancient polyploid (palaeopolyploid). About 78% of the predicted genes occur in chromosome ends, which comprise less than one-half of the genome but account for nearly all of the genetic recombination. Genome duplications occurred at approximately 59 and 13 million years ago, resulting in a highly duplicated genome with nearly 75% of the genes present in multiple copies. The two duplication events were followed by gene diversification and loss, and numerous chromosome rearrangements. An accurate soybean genome sequence will facilitate the identification of the genetic basis of many soybean traits, and accelerate the creation of improved soybean varieties.

Transcriptomic response of the mycoparasitic fungus Trichoderma atroviride to the presence of a fungal prey

Background

Combating the action of plant pathogenic microorganisms by mycoparasitic fungi has been announced as an attractive biological alternative to the use of chemical fungicides since two decades. The fungal genus Trichoderma includes a high number of taxa which are able to recognize, combat and finally besiege and kill their prey. Only fragments of the biochemical processes related to this ability have been uncovered so far, however.

Results

We analyzed genome-wide gene expression changes during the begin of physical contact between Trichoderma atroviride and two plant pathogens Botrytis cinerea and Rhizoctonia solani, and compared with gene expression patterns of mycelial and conidiating cultures, respectively. About 3000 ESTs, representing about 900 genes, were obtained from each of these three growth conditions. 66 genes, represented by 442 ESTs, were specifically and significantly overexpressed during onset of mycoparasitism, and the expression of a subset thereof was verified by expression analysis. The upregulated genes comprised 18 KOG groups, but were most abundant from the groups representing posttranslational processing, and amino acid metabolism, and included components of the stress response, reaction to nitrogen shortage, signal transduction and lipid catabolism. Metabolic network analysis confirmed the upregulation of the genes for amino acid biosynthesis and of those involved in the catabolism of lipids and aminosugars.

Conclusion

The analysis of the genes overexpressed during the onset of mycoparasitism in T. atroviride has revealed that the fungus reacts to this condition with several previously undetected physiological reactions. These data enable a new and more comprehensive interpretation of the physiology of mycoparasitism, and will aid in the selection of traits for improvement of biocontrol strains by recombinant techniques.