GEMmaker: process massive RNA-seq datasets on heterogeneous computational infrastructure.

BACKGROUND: Quantification of gene expression from RNA-seq data is a prerequisite for transcriptome analysis such as differential gene expression analysis and gene co-expression network construction. Individual RNA-seq experiments are larger and combining multiple experiments from sequence repositories can result in datasets with thousands of samples. Processing hundreds to thousands of RNA-seq data can result in challenges related to data management, access to sufficient computational resources, navigation of high-performance computing (HPC) systems, installation of required software dependencies, and reproducibility. Processing of larger and deeper RNA-seq experiments will become more common as sequencing technology matures. RESULTS: GEMmaker, is a nf-core compliant, Nextflow workflow, that quantifies gene expression from small to massive RNA-seq datasets. GEMmaker ensures results are highly reproducible through the use of versioned containerized software that can be executed on a single workstation, institutional compute cluster, Kubernetes platform or the cloud. GEMmaker supports popular alignment and quantification tools providing results in raw and normalized formats. GEMmaker is unique in that it can scale to process thousands of local or remote stored samples without exceeding available data storage. CONCLUSIONS: Workflows that quantify gene expression are not new, and many already address issues of portability, reusability, and scale in terms of access to CPUs. GEMmaker provides these benefits and adds the ability to scale despite low data storage infrastructure. This allows users to process hundreds to thousands of RNA-seq samples even when data storage resources are limited. GEMmaker is freely available and fully documented with step-by-step setup and execution instructions.

SHR4z, a novel decoy effector from the haustorium of the parasitic weed Striga gesnerioides, suppresses host plant immunity.

Cowpea (Vigna unguiculata) cultivar B301 is resistant to races SG4 and SG3 of the root parasitic weed Striga gesnerioides, developing a hypersensitive response (HR) at the site of parasite attachment. By contrast, race SG4z overcomes B301 resistance and successfully parasitises the plant. Comparative transcriptomics and in silico analysis identified a small secreted effector protein dubbed Suppressor of Host Resistance 4z (SHR4z) in the SG4z haustorium that upon transfer to the host roots causes a loss of host immunity (i.e. decreased HR and increased parasite growth). SHR4z has significant homology to the short leucine-rich repeat (LRR) domain of SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) family proteins and functions by binding to VuPOB1, a host BTB-BACK domain-containing ubiquitin E3 ligase homologue, leading to its rapid turnover. VuPOB1 is shown to be a positive regulator of HR since silencing of VuPOB1 expression in transgenic B301 roots lowers the frequency of HR and increases the levels of successful SG4 parasitism and overexpression decreases parasitism by SG4z. These findings provide new insights into how parasitic weeds overcome host defences and could potentially contribute to the development of novel strategies for controlling Striga and other parasitic weeds thereby enhancing crop productivity and food security globally.

Risk versus reward: host dependent parasite mortality rates and phenotypes in the facultative generalist Triphysaria versicolor.

BACKGROUND: Parasitic plants engage in a complex molecular dialog with potential host plants to identify a host and overcome host defenses to initiate development of the parasitic feeding organ, the haustorium, invade host tissues, and withdraw water and nutrients. While one of two critical signaling events in the parasitic plant life cycle (germination via stimulant chemicals) has been relatively well-studied, the signaling event that triggers haustorium formation remains elusive. Elucidation of this poorly understood molecular dialogue will shed light on plant-plant communication, parasitic plant physiology, and the evolution of parasitism in plants. RESULTS: Here we present an experimental framework that develops easily quantifiable contrasts for the facultative generalist parasitic plant, Triphysaria, as it feeds across a broad range of diverse flowering plants. The contrasts, including variable parasite growth form and mortality when grown with different hosts, suggest a dynamic and host-dependent molecular dialogue between the parasite and host. Finally, by comparing transcriptome datasets from attached versus unattached parasites we gain insight into some of the physiological processes that are altered during parasitic behavior including shifts in photosynthesis-related and stress response genes. CONCLUSIONS: This work sheds light on Triphysaria's parasitic life habit and is an important step towards understanding the mechanisms of haustorium initiation factor perception, a unique form of plant-plant communication.

Horizontal gene transfer is more frequent with increased heterotrophy and contributes to parasite adaptation.

Horizontal gene transfer (HGT) is the transfer of genetic material across species boundaries and has been a driving force in prokaryotic evolution. HGT involving eukaryotes appears to be much less frequent, and the functional implications of HGT in eukaryotes are poorly understood. We test the hypothesis that parasitic plants, because of their intimate feeding contacts with host plant tissues, are especially prone to horizontal gene acquisition. We sought evidence of HGTs in transcriptomes of three parasitic members of Orobanchaceae, a plant family containing species spanning the full spectrum of parasitic capabilities, plus the free-living Lindenbergia Following initial phylogenetic detection and an extensive validation procedure, 52 high-confidence horizontal transfer events were detected, often from lineages of known host plants and with an increasing number of HGT events in species with the greatest parasitic dependence. Analyses of intron sequences in putative donor and recipient lineages provide evidence for integration of genomic fragments far more often than retro-processed RNA sequences. Purifying selection predominates in functionally transferred sequences, with a small fraction of adaptively evolving sites. HGT-acquired genes are preferentially expressed in the haustorium-the organ of parasitic plants-and are strongly biased in predicted gene functions, suggesting that expression products of horizontally acquired genes are contributing to the unique adaptive feeding structure of parasitic plants.

Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty.

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative "parasitism genes." Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.

Phylogenetic analysis of pectin-related gene families in Physcomitrella patens and nine other plant species yields evolutionary insights into cell walls.

BACKGROUND: Pectins are acidic sugar-containing polysaccharides that are universally conserved components of the primary cell walls of plants and modulate both tip and diffuse cell growth. However, many of their specific functions and the evolution of the genes responsible for producing and modifying them are incompletely understood. The moss Physcomitrella patens is emerging as a powerful model system for the study of plant cell walls. To identify deeply conserved pectin-related genes in Physcomitrella, we generated phylogenetic trees for 16 pectin-related gene families using sequences from ten plant genomes and analyzed the evolutionary relationships within these families. RESULTS: Contrary to our initial hypothesis that a single ancestral gene was present for each pectin-related gene family in the common ancestor of land plants, five of the 16 gene families, including homogalacturonan galacturonosyltransferases, polygalacturonases, pectin methylesterases, homogalacturonan methyltransferases, and pectate lyase-like proteins, show evidence of multiple members in the early land plant that gave rise to the mosses and vascular plants. Seven of the gene families, the UDP-rhamnose synthases, UDP-glucuronic acid epimerases, homogalacturonan galacturonosyltransferase-like proteins, ß-1,4-galactan ß-1,4-galactosyltransferases, rhamnogalacturonan II xylosyltransferases, and pectin acetylesterases appear to have had a single member in the common ancestor of land plants. We detected no Physcomitrella members in the xylogalacturonan xylosyltransferase, rhamnogalacturonan I arabinosyltransferase, pectin methylesterase inhibitor, or polygalacturonase inhibitor protein families. CONCLUSIONS: Several gene families related to the production and modification of pectins in plants appear to have multiple members that are conserved as far back as the common ancestor of mosses and vascular plants. The presence of multiple members of these families even before the divergence of other important cell wall-related genes, such as cellulose synthases, suggests a more complex role than previously suspected for pectins in the evolution of land plants. The presence of relatively small pectin-related gene families in Physcomitrella as compared to Arabidopsis makes it an attractive target for analysis of the functions of pectins in cell walls. In contrast, the absence of genes in Physcomitrella for some families suggests that certain pectin modifications, such as homogalacturonan xylosylation, arose later during land plant evolution.

Towards identification of postharvest fruit quality transcriptomic markers in Malus domestica.

Gene expression is highly impacted by the environment and can be reflective of past events that affected developmental processes. It is therefore expected that gene expression can serve as a signal of a current or future phenotypic traits. In this paper we identify sets of genes, which we call Prognostic Transcriptomic Biomarkers (PTBs), that can predict firmness in Malus domestica (apple) fruits. In apples, all individuals of a cultivar are clones, and differences in fruit quality are due to the environment. The apples transcriptome responds to these differences in environment, which makes PTBs an attractive predictor of future fruit quality. PTBs have the potential to enhance supply chain efficiency, reduce crop loss, and provide higher and more consistent quality for consumers. However, several questions must be addressed. In this paper we answer the question of which of two common modeling approaches, Random Forest or ElasticNet, outperforms the other. We answer if PTBs with few genes are efficient at predicting traits. This is important because we need few genes to perform qPCR, and we answer the question if qPCR is a cost-effective assay as input for PTBs modeled using high-throughput RNA-seq. To do this, we conducted a pilot study using fruit texture in the 'Gala' variety of apples across several postharvest storage regiments. Fruit texture in 'Gala' apples is highly controllable by post-harvest treatments and is therefore a good candidate to explore the use of PTBs. We find that the RandomForest model is more consistent than an ElasticNet model and is predictive of firmness (r2 = 0.78) with as few as 15 genes. We also show that qPCR is reasonably consistent with RNA-seq in a follow up experiment. Results are promising for PTBs, yet more work is needed to ensure that PTBs are robust across various environmental conditions and storage treatments.

Correction: Towards identification of postharvest fruit quality transcriptomic markers in Malus domestica.

[This corrects the article DOI: 10.1371/journal.pone.0297015.].

A chromosome-scale assembly for 'd'Anjou' pear.

Cultivated pear consists of several Pyrus species with Pyrus communis (European pear) representing a large fraction of worldwide production. As a relatively recently domesticated crop and perennial tree, pear can benefit from genome-assisted breeding. Additionally, comparative genomics within Rosaceae promises greater understanding of evolution within this economically important family. Here, we generate a fully phased chromosome-scale genome assembly of P. communis 'd'Anjou.' Using PacBio HiFi and Dovetail Omni-C reads, the genome is resolved into the expected 17 chromosomes, with each haplotype totaling nearly 540 Megabases and a contig N50 of nearly 14 Mb. Both haplotypes are highly syntenic to each other and to the Malus domestica 'Honeycrisp' apple genome. Nearly 45,000 genes were annotated in each haplotype, over 90% of which have direct RNA-seq expression evidence. We detect signatures of the known whole-genome duplication shared between apple and pear, and we estimate 57% of d'Anjou genes are retained in duplicate derived from this event. This genome highlights the value of generating phased diploid assemblies for recovering the full allelic complement in highly heterozygous crop species.

Evolution of a horizontally acquired legume gene, albumin 1, in the parasitic plant Phelipanche aegyptiaca and related species.

BACKGROUND: Parasitic plants, represented by several thousand species of angiosperms, use modified structures known as haustoria to tap into photosynthetic host plants and extract nutrients and water. As a result of their direct plant-plant connections with their host plant, parasitic plants have special opportunities for horizontal gene transfer, the nonsexual transmission of genetic material across species boundaries. There is increasing evidence that parasitic plants have served as recipients and donors of horizontal gene transfer (HGT), but the long-term impacts of eukaryotic HGT in parasitic plants are largely unknown. RESULTS: Here we show that a gene encoding albumin 1 KNOTTIN-like protein, closely related to the albumin 1 genes only known from papilionoid legumes, where they serve dual roles as food storage and insect toxin, was found in Phelipanche aegyptiaca and related parasitic species of family Orobanchaceae, and was likely acquired by a Phelipanche ancestor via HGT from a legume host based on phylogenetic analyses. The KNOTTINs are well known for their unique "disulfide through disulfide knot" structure and have been extensively studied in various contexts, including drug design. Genomic sequences from nine related parasite species were obtained, and 3D protein structure simulation tests and evolutionary constraint analyses were performed. The parasite gene we identified here retains the intron structure, six highly conserved cysteine residues necessary to form a KNOTTIN protein, and displays levels of purifying selection like those seen in legumes. The albumin 1 xenogene has evolved through >150 speciation events over ca. 16 million years, forming a small family of differentially expressed genes that may confer novel functions in the parasites. Moreover, further data show that a distantly related parasitic plant, Cuscuta, obtained two copies of albumin 1 KNOTTIN-like genes from legumes through a separate HGT event, suggesting that legume KNOTTIN structures have been repeatedly co-opted by parasitic plants. CONCLUSIONS: The HGT-derived albumins in Phelipanche represent a novel example of how plants can acquire genes from other plants via HGT that then go on to duplicate, evolve, and retain the specialized features required to perform a unique host-derived function.

Functional genomics of a generalist parasitic plant: laser microdissection of host-parasite interface reveals host-specific patterns of parasite gene expression.

BACKGROUND: Orobanchaceae is the only plant family with members representing the full range of parasitic lifestyles plus a free-living lineage sister to all parasitic lineages, Lindenbergia. A generalist member of this family, and an important parasitic plant model, Triphysaria versicolor regularly feeds upon a wide range of host plants. Here, we compare de novo assembled transcriptomes generated from laser micro-dissected tissues at the host-parasite interface to uncover details of the largely uncharacterized interaction between parasitic plants and their hosts. RESULTS: The interaction of Triphysaria with the distantly related hosts Zea mays and Medicago truncatula reveals dramatic host-specific gene expression patterns. Relative to above ground tissues, gene families are disproportionally represented at the interface including enrichment for transcription factors and genes of unknown function. Quantitative Real-Time PCR of a T. versicolor ß-expansin shows strong differential (120x) upregulation in response to the monocot host Z. mays; a result that is concordant with our read count estimates. Pathogenesis-related proteins, other cell wall modifying enzymes, and orthologs of genes with unknown function (annotated as such in sequenced plant genomes) are among the parasite genes highly expressed by T. versicolor at the parasite-host interface. CONCLUSIONS: Laser capture microdissection makes it possible to sample the small region of cells at the epicenter of parasite host interactions. The results of our analysis suggest that T. versicolor's generalist strategy involves a reliance on overlapping but distinct gene sets, depending upon the host plant it is parasitizing. The massive upregulation of a T. versicolor ß-expansin is suggestive of a mechanism for parasite success on grass hosts. In this preliminary study of the interface transcriptomes, we have shown that T. versicolor, and the Orobanchaceae in general, provide excellent opportunities for the characterization of plant genes with unknown functions.

Transcriptomic approach to uncover dynamic events in the development of mid-season sunburn in apple fruit.

Apples grown in high heat, high light, and low humidity environments are at risk for sun injury disorders like sunburn and associated crop losses. Understanding the physiological and molecular mechanisms underlying sunburn will support improvement of mitigation strategies and breeding for more resilient varieties. Numerous studies have highlighted key biochemical processes involved in sun injury, such as the phenylpropanoid and reactive oxygen species (ROS) pathways, demonstrating both enzyme activities and expression of related genes in response to sunburn conditions. Most previous studies have focused on at-harvest activity of a small number of genes in response to heat stress. Thus, it remains unclear how stress events earlier in the season affect physiology and gene expression. Here, we applied heat stress to mid-season apples in the field and collected tissue along a time course-24, 48, and 72 h following a heat stimulus-to investigate dynamic gene expression changes using a transcriptomic lens. We found a relatively small number of differentially expressed genes (DEGs) and enriched functional terms in response to heat treatments. Only a few of these belonged to pathways previously described to be involved in sunburn, such as the AsA-GSH pathway, while most DEGs had not yet been implicated in sunburn or heat stress in pome fruit.

PlantTribes2: Tools for comparative gene family analysis in plant genomics.

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

A phased, chromosome-scale genome of 'Honeycrisp' apple (Malus domestica).

The apple cultivar 'Honeycrisp' has superior fruit quality traits, cold hardiness, and disease resistance, making it a popular breeding parent. However, it suffers from several physiological disorders, production, and postharvest issues. Despite several available apple genome sequences, understanding of the genetic mechanisms underlying cultivar-specific traits remains lacking. Here, we present a highly contiguous, fully phased, chromosome-level genome of 'Honeycrisp' apples, using PacBio HiFi, Omni-C, and Illumina sequencing platforms, with two assembled haplomes of 674 Mbp and 660 Mbp, and contig N50 values of 32.8 Mbp and 31.6 Mbp, respectively. Overall, 47,563 and 48,655 protein-coding genes were annotated from each haplome, capturing 96.8-97.4% complete BUSCOs in the eudicot database. Gene family analysis reveals most 'Honeycrisp' genes are assigned into orthogroups shared with other genomes, with 121 'Honeycrisp'-specific orthogroups. This resource is valuable for understanding the genetic basis of important traits in apples and related Rosaceae species to enhance breeding efforts.

Building a foundation for gene family analysis in Rosaceae genomes with a novel workflow: A case study in Pyrus architecture genes.

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, 'd'Anjou.' Our comparative gene family approach revealed significant issues with the most recent 'Bartlett' genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the 'Bartlett' genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

Comparative Analysis of the Apple Root Transcriptome as Affected by Rootstock Genotype and Brassicaceae Seed Meal Soil Amendment: Implications for Plant Health.

Brassicaceae seed meal (SM) soil amendment has been utilized as an effective strategy to control the biological complex of organisms, which includes oomycetes, fungi, and parasitic nematodes, that incites the phenomenon termed apple replant disease. Soil-borne disease control attained in response to Brassicaceae SM amendment is reliant on multiple chemical and biological attributes, including specific SM-generated modifications to the soil/rhizosphere microbiome. In this study, we conducted a comparative analyses of apple root gene expression as influenced by rootstock genotype combined with a seed meal (SM) soil amendment. Apple replant disease (ARD) susceptible (M.26) and tolerant (G.210) rootstocks cultivated in SM-amended soil exhibited differential gene expression relative to corresponding non-treated control (NTC) orchard soil. The temporal dynamics of gene expression indicated that the SM-amended soil system altered the trajectory of the root transcriptome in a genotype-specific manner. In both genotypes, the expression of genes related to plant defense and hormone signaling were altered in SM-amended soil, suggesting SM-responsive phytohormone regulation. Altered gene expression was temporally associated with changes in rhizosphere microbiome density and composition in the SM-treated soil. Gene expression analysis across the two rootstocks cultivated in the pathogen-infested NTC soil showed genotype-specific responses indicative of different defensive strategies. These results are consistent with previously described resistance mechanisms of ARD "tolerant" rootstock cultivars and also add to our understanding of the multiple mechanisms by which SM soil amendment and the resulting rhizosphere microbiome affect apple rootstock physiology. Future studies which assess transcriptomic and metagenomic data in parallel will be important for illuminating important connections between specific rhizosphere microbiota, gene-regulation, and plant health.

Transcriptomics of Differential Ripening in 'd'Anjou' Pear (Pyrus communis L.).

Estimating maturity in pome fruits is a critical task that directs virtually all postharvest supply chain decisions. This is especially important for European pear (Pyrus communis) cultivars because losses due to spoilage and senescence must be minimized while ensuring proper ripening capacity is achieved (in part by satisfying a fruit chilling requirement). Reliable methods are lacking for accurate estimation of pear fruit maturity, and because ripening is maturity dependent it makes predicting ripening capacity a challenge. In this study of the European pear cultivar 'd'Anjou', we sorted fruit at harvest based upon on-tree fruit position to build contrasts of maturity. Our sorting scheme showed clear contrasts of maturity between canopy positions, yet there was substantial overlap in the distribution of values for the index of absorbance difference (I AD ), a non-destructive spectroscopic measurement that has been used as a proxy for pome fruit maturity. This presented an opportunity to explore a contrast of maturity that was more subtle than I AD could differentiate, and thus guided our subsequent transcriptome analysis of tissue samples taken at harvest and during storage. Using a novel approach that tests for condition-specific differences of co-expressed genes, we discovered genes with a phased character that mirrored our sorting scheme. The expression patterns of these genes are associated with fruit quality and ripening differences across the experiment. Functional profiles of these co-expressed genes are concordant with previous findings, and also offer new clues, and thus hypotheses, about genes involved in pear fruit quality, maturity, and ripening. This work may lead to new tools for enhanced postharvest management based on activity of gene co-expression modules, rather than individual genes. Further, our results indicate that modules may have utility within specific windows of time during postharvest management of 'd'Anjou' pear.

Temporal Dynamics of the Soil Metabolome and Microbiome During Simulated Anaerobic Soil Disinfestation.

Significant interest exists in engineering the soil microbiome to attain suppression of soil-borne plant diseases. Anaerobic soil disinfestation (ASD) has potential as a biologically regulated disease control method; however, the role of specific metabolites and microbial community dynamics contributing to ASD mediated disease control is mostly uncharacterized. Understanding the trajectory of co-evolutionary processes leading to syntrophic generation of functional metabolites during ASD is a necessary prelude to the predictive utilization of this disease management approach. Consequently, metabolic and microbial community profiling were used to generate highly dimensional datasets and network analysis to identify sequential transformations through aerobic, facultatively anaerobic, and anaerobic soil phases of the ASD process and distinct groups of metabolites and microorganisms linked with those stages. Transient alterations in abundance of specific microbial groups, not consistently accounted for in previous studies of the ASD process, were documented in this time-course study. Such events initially were associated with increases and subsequent diminution in highly labile metabolites conferred by the carbon input. Proliferation and dynamic compositional changes in the Firmicutes community continued throughout the anaerobic phase and was linked to temporal changes in metabolite abundance including accumulation of small chain organic acids, methyl sulfide compounds, hydrocarbons, and p-cresol with antimicrobial properties. Novel potential modes of disease control during ASD were identified and the importance of the amendment and "community metabolism" for temporally supplying specific classes of labile compounds were revealed.

Genome Sequence of Striga asiatica Provides Insight into the Evolution of Plant Parasitism.

Parasitic plants in the genus Striga, commonly known as witchweeds, cause major crop losses in sub-Saharan Africa and pose a threat to agriculture worldwide. An understanding of Striga parasite biology, which could lead to agricultural solutions, has been hampered by the lack of genome information. Here, we report the draft genome sequence of Striga asiatica with 34,577 predicted protein-coding genes, which reflects gene family contractions and expansions that are consistent with a three-phase model of parasitic plant genome evolution. Striga seeds germinate in response to host-derived strigolactones (SLs) and then develop a specialized penetration structure, the haustorium, to invade the host root. A family of SL receptors has undergone a striking expansion, suggesting a molecular basis for the evolution of broad host range among Striga spp. We found that genes involved in lateral root development in non-parasitic model species are coordinately induced during haustorium development in Striga, suggesting a pathway that was partly co-opted during the evolution of the haustorium. In addition, we found evidence for horizontal transfer of host genes as well as retrotransposons, indicating gene flow to S. asiatica from hosts. Our results provide valuable insights into the evolution of parasitism and a key resource for the future development of Striga control strategies.

A practical examination of RNA isolation methods for European pear (Pyrus communis).

OBJECTIVE: With the goal of identifying fast, reliable, and broadly applicable RNA isolation methods in European pear fruit for downstream transcriptome analysis, we evaluated several commercially available kit-based RNA isolation methods, plus our modified version of a published cetyl trimethyl ammonium bromide (CTAB)-based method. RESULTS: We confirmed previous work indicating chaotropic agent-based kits produced sufficient, high-quality RNA in freshly harvested, mature 'Bartlett' fruit. However, RNA isolation from 'd'Anjou' pear peel and especially cortical tissues of fruit stored for 11 months proved challenging to all but our modified CTAB-based method. Generally, more RNA was recovered from peel tissues than cortical tissues. Less toxic dithiothreitol was confirmed to be an acceptable reducing agent as its substitution for 2-mercaptoethanol often yielded high quality RNA. Finally, we present evidence that erroneous signals in the 5S region of Bioanalyzer RNA size plot histograms, that interfered with RNA integrity number calculation, were small RNA fragments that are reduced by simple cleanup procedures, not artifacts as previously reported.