Expanding the RNA virome of nematodes and other soil-inhabiting organisms

In recent years, several newly discovered viruses infecting free-living nematodes, sedentary plant-parasitic nematodes, and migratory root lesion nematodes have been described. However, to the best of our knowledge, no comprehensive research focusing exclusively on metagenomic analysis of the soil nematode community virome has thus far been carried out. In this work, we have attempted to bridge this gap by investigating viral communities that are associated with soil-inhabiting organisms, particularly nematodes. This study demonstrates a remarkable diversity of RNA viruses in the natural soil environment. Over 150 viruses were identified in different soil-inhabiting hosts, of which more than 139 are potentially new virus species. Many of these viruses belong to the nematode virome, thereby enriching our understanding of the diversity and evolution of this complex part of the natural ecosystem.

Transcriptome analysis identifies genes related to the waxy coating on blueberry fruit in two northern-adapted rabbiteye breeding populations

BACKGROUND

Blueberry is of high economic value. Most blueberry varieties selected for the fresh market have an appealing light blue coating or "bloom" on the fruit due to the presence of a visible heavy epicuticular wax layer. This waxy layer also serves as natural defense against fruit desiccation and deterioration.

RESULTS

In this study, we attempted to identify gene(s) whose expression is related to the protective waxy coating on blueberry fruit utilizing two unique germplasm populations that segregate for the waxy layer. We bulked RNA from waxy and non-waxy blueberry progenies from the two northern-adapted rabbiteye hybrid breeding populations ('Nocturne' x T 300 and 'Nocturne' x US 1212), and generated 316.85 million RNA-seq reads. We de novo assembled this data set integrated with other publicly available RNA-seq data and trimmed the assembly into a 91,861 blueberry unigene collection. All unigenes were functionally annotated, resulting in 79 genes potentially related to wax accumulation. We compared the expression pattern of waxy and non-waxy progenies using edgeR and identified overall 1125 genes in the T 300 population and 2864 genes in the US 1212 population with at least a two-fold expression difference. After validating differential expression of several genes by RT-qPCR experiments, a candidate gene, FatB, which encodes acyl-[acyl-carrier-protein] hydrolase, emerged whose expression was closely linked to the segregation of the waxy coating in our populations. This gene was expressed at more than a five-fold higher level in waxy than non-waxy plants of both populations. We amplified and sequenced the cDNA for this gene from three waxy plants of each population, but were unable to amplify the cDNA from three non-waxy plants that were tested from each population. We aligned the Vaccinium deduced FATB protein sequence to FATB protein sequences from other plant species. Within the PF01643 domain, which gives FATB its catalytic function, 80.08% of the amino acids were identical or had conservative replacements between the blueberry and the Cucumis melo sequence (XP_008467164). We then amplified and sequenced a large portion of the FatB gene itself from waxy and non-waxy individuals of both populations. Alignment of the cDNA and gDNA sequences revealed that the blueberry FatB gene consists of six exons and five introns. Although we did not sequence through two very large introns, a comparison of the exon sequences found no significant sequence differences between the waxy and non-waxy plants. This suggests that another gene, which regulates or somehow affects FatB expression, must be segregating in the populations.

CONCLUSIONS

This study is helping to achieve a greater understanding of epicuticular wax biosynthesis in blueberry. In addition, the blueberry unigene collection should facilitate functional annotation of the coming chromosomal level blueberry genome.

LUX ARRHYTHMO mediates crosstalk between the circadian clock and defense in Arabidopsis

The circadian clock is known to regulate plant innate immunity but the underlying mechanism of this regulation remains largely unclear. We show here that mutations in the core clock component LUX ARRHYTHMO (LUX) disrupt circadian regulation of stomata under free running and Pseudomonas syringae challenge conditions as well as defense signaling mediated by SA and JA, leading to compromised disease resistance. RNA-seq analysis reveals that both clock- and defense-related genes are regulated by LUX. LUX binds to clock gene promoters that have not been shown before, expanding the clock gene networks that require LUX function. LUX also binds to the promoters of EDS1 and JAZ5, likely acting through these genes to affect SA- and JA-signaling. We further show that JA signaling reciprocally affects clock activity. Thus, our data support crosstalk between the circadian clock and plant innate immunity and imply an important role of LUX in this process.

LUX ARRHYTHMO mediates crosstalk between the circadian clock and defense in Arabidopsis

The circadian clock is known to regulate plant innate immunity but the underlying mechanism of this regulation remains largely unclear. We show here that mutations in the core clock component LUX ARRHYTHMO (LUX) disrupt circadian regulation of stomata under free running and Pseudomonassyringae challenge conditions as well as defense signaling mediated by SA and JA, leading to compromised disease resistance. RNA-seq analysis reveals that both clock- and defense-related genes are regulated by LUX. LUX binds to clock gene promoters that have not been shown before, expanding the clock gene networks that require LUX function. LUX also binds to the promoters of EDS1 and JAZ5, likely acting through these genes to affect SA- and JA-signaling. We further show that JA signaling reciprocally affects clock activity. Thus, our data support crosstalk between the circadian clock and plant innate immunity and imply an important role of LUX in this process.

Circadian control of plant defence likely reflects plants’ ability to coordinate development and defense. Here, Zhang et al. show that LUX regulates stomatal defense  and SA/JA signaling, leading to broad-spectrum disease resistance, and that JA signaling can, in turn, regulate clock activity.

An eFP browser for visualizing strawberry fruit and flower transcriptomes

Wild strawberry Fragaria vesca is emerging as an important model system for the cultivated strawberry due to its diploid genome and availability of extensive transcriptome data and a range of molecular genetic tools. Being able to better utilize these tools, especially the transcriptome data, will greatly facilitate research progress in strawberry and other Rosaceae fruit crops. The electronic fluorescent pictograph (eFP) software is a useful and popular tool to display transcriptome data visually, and is widely used in other model organisms including Arabidopsis and mouse. Here we applied eFP to display wild strawberry RNA sequencing (RNA-seq) data from 42 different tissues and stages, including various flower and fruit developmental stages. In addition, we generated eight additional RNA-seq data sets to represent tissues from ripening-stage receptacle fruit from yellow-colored and red-colored wild strawberry varieties. Differential gene expression analysis between these eight data sets provides additional information for understanding fruit-quality traits. Together, this work greatly facilitates the utility of the extensive transcriptome data for investigating strawberry flower and fruit development as well as fruit-quality traits.

Genome-wide functional annotation of Phomopsis longicolla isolate MSPL 10-6

Phomopsis seed decay of soybean is caused primarily by the seed-borne fungal pathogen Phomopsis longicolla (syn. Diaporthe longicolla). This disease severely decreases soybean seed quality, reduces seedling vigor and stand establishment, and suppresses yield. It is one of the most economically important soybean diseases. In this study we annotated the entire genome of P. longicolla isolate MSPL 10-6, which was isolated from field-grown soybean seed in Mississippi, USA. This study represents the first reported genome-wide functional annotation of a seed borne fungal pathogen in the Diaporthe-Phomopsis complex. The P. longicolla genome annotation will enable research into the genetic basis of fungal infection of soybean seed and provide information for the study of soybean-fungal interactions. The genome annotation will also be a valuable resource for the research and agricultural communities. It will aid in the development of new control strategies for this pathogen. The annotations can be found from: http://bioinformatics.towson.edu/phomopsis_longicolla/download.html. NCBI accession number is: AYRD00000000.

Development of Dictyostelium discoideum as a model to identify evolutionarily-conserved virulence factors in Escherichia coli

Escherichia coli are a typical component of mammalian gut microbiota; however, a number of pathogenic strains remain a health concern for developing countries. Identification and disruption of virulence factors that promote infection could lead to reduction in infection. To identify novel bacterial virulence factors, we are using the soil amoeba, Dictyostelium discoideum, which preys upon bacteria, including E. coli, in the soil, as a model for mammalian phagocytic cells. Using random mutagenesis, we have created mutant stocks of MG1655 E. coli. These stocks were fed to D. discoideum, enriching the culture and selecting for mutants that were resistant to killing upon phagocytosis. After several rounds of selection, resistant mutant strains were isolated, and we found that some strains showed up to 100-fold increase in survival over a 2 h period. Once isolated, the genomes of the resistant mutant strains were sequenced and the sequences compared to that of the initial WT strain using Illumina’s MiSeq technology. While analyzing the sequences from these mutant strains, we found mutations in several potential genes of interest and are working to characterize their roles in the cell. We are testing E. coli strains reconstituted with specific mutations to identify particular mutations that are associated with phagocytosis-resistant phenotypes. These studies are aimed at identifying particular bacterial genes that may be involved in regulation of virulence. In addition, given that the development of mechanisms to resist D. discoideum predation may have contributed to the selection and maintenance of bacterial virulence factors against mammalian hosts, these studies may provide insight on the evolution of host-pathogen interactions.

Proteomic Investigation of Rhizoctonia solani AG 4 Identifies Secretome and Mycelial Proteins with Roles in Plant Cell Wall Degradation and Virulence

Rhizoctonia solani AG 4 is a soilborne necrotrophic fungal plant pathogen that causes economically important diseases on agronomic crops worldwide. This study used a proteomics approach to characterize both intracellular proteins and the secretome of R. solani AG 4 isolate Rs23A under several growth conditions, the secretome being highly important in pathogenesis. From over 500 total secretome and soluble intracellular protein spots from 2-D gels, 457 protein spots were analyzed and 318 proteins positively matched with fungal proteins of known function by comparison with available R. solani genome databases specific for anastomosis groups 1-IA, 1-IB, and 3. These proteins were categorized to possible cellular locations and functional groups and for some proteins their putative roles in plant cell wall degradation and virulence. The majority of the secreted proteins were grouped to extracellular regions and contain hydrolase activity.

Draft genome sequence of Phomopsis longicolla isolate MSPL 10-6

Phomopsis longicolla is the primary cause of Phomopsis seed decay in soybean. This disease severely affects soybean seed quality by reducing seed viability and oil content, altering seed composition, and increasing frequencies of moldy and/or split beans. It is one of the most economically important soybean diseases. Here, we report the de novo assembled draft genome sequence of the P. longicolla isolate MSPL10-6, which was isolated from field-grown soybean seed in Mississippi, USA. This study represents the first reported genome sequence of a seedborne fungal pathogen in the Diaporthe-Phomopsis complex. The P. longicolla genome sequence will enable research into the genetic basis of fungal infection of soybean seed and provide information for the study of soybean-fungal interactions. The genome sequence will also be valuable for molecular genetic marker development, manipulation of pathogenicity-related genes and development of new control strategies for this pathogen.

Floral Transcriptomes in Woodland Strawberry Uncover Developing Receptacle and Anther Gene Networks

Flowers are reproductive organs and precursors to fruits and seeds. While the basic tenets of the ABCE model of flower development are conserved in angiosperms, different flowering plants exhibit different and sometimes unique characteristics. A distinct feature of strawberry (Fragaria spp.) flowers is the development of several hundreds of individual apocarpous (unfused) carpels. These individual carpels are arranged in a spiral pattern on the subtending stem tip, the receptacle. Therefore, the receptacle is an integral part of the strawberry flower and is of significant agronomic importance, being the precursor to strawberry fruit. Taking advantage of next-generation sequencing and laser capture microdissection, we generated different tissue- and stage-transcriptomic profiling of woodland strawberry (Fragaria vesca) flower development. Using pairwise comparisons and weighted gene coexpression network analysis, we identified modules of coexpressed genes and hub genes of tissue-specific networks. Of particular importance is the discovery of a developing receptacle-specific module exhibiting similar molecular features to those of young floral meristems. The strawberry homologs of a number of meristem regulators, including LOST MERISTEM and WUSCHEL, are identified as hub genes operating in the developing receptacle network. Furthermore, almost 25% of the F-box genes in the genome are transiently induced in developing anthers at the meiosis stage, indicating active protein degradation. Together, this work provides important insights into the molecular networks underlying strawberry's unique reproductive developmental processes. This extensive floral transcriptome data set is publicly available and can be readily queried at the project Web site, serving as an important genomic resource for the plant biology research community.

Genome-scale transcriptomic insights into early-stage fruit development in woodland strawberry Fragaria vesca

Fragaria vesca, a diploid woodland strawberry with a small and sequenced genome, is an excellent model for studying fruit development. The strawberry fruit is unique in that the edible flesh is actually enlarged receptacle tissue. The true fruit are the numerous dry achenes dotting the receptacle's surface. Auxin produced from the achene is essential for the receptacle fruit set, a paradigm for studying crosstalk between hormone signaling and development. To investigate the molecular mechanism underlying strawberry fruit set, next-generation sequencing was employed to profile early-stage fruit development with five fruit tissue types and five developmental stages from floral anthesis to enlarged fruits. This two-dimensional data set provides a systems-level view of molecular events with precise spatial and temporal resolution. The data suggest that the endosperm and seed coat may play a more prominent role than the embryo in auxin and gibberellin biosynthesis for fruit set. A model is proposed to illustrate how hormonal signals produced in the endosperm and seed coat coordinate seed, ovary wall, and receptacle fruit development. The comprehensive fruit transcriptome data set provides a wealth of genomic resources for the strawberry and Rosaceae communities as well as unprecedented molecular insight into fruit set and early stage fruit development.

Generation and analysis of blueberry transcriptome sequences from leaves, developing fruit, and flower buds from cold acclimation through deacclimation

BACKGROUND

There has been increased consumption of blueberries in recent years fueled in part because of their many recognized health benefits. Blueberry fruit is very high in anthocyanins, which have been linked to improved night vision, prevention of macular degeneration, anti-cancer activity, and reduced risk of heart disease. Very few genomic resources have been available for blueberry, however. Further development of genomic resources like expressed sequence tags (ESTs), molecular markers, and genetic linkage maps could lead to more rapid genetic improvement. Marker-assisted selection could be used to combine traits for climatic adaptation with fruit and nutritional quality traits.

RESULTS

Efforts to sequence the transcriptome of the commercial highbush blueberry (Vaccinium corymbosum) cultivar Bluecrop and use the sequences to identify genes associated with cold acclimation and fruit development and develop SSR markers for mapping studies are presented here. Transcriptome sequences were generated from blueberry fruit at different stages of development, flower buds at different stages of cold acclimation, and leaves by next-generation Roche 454 sequencing. Over 600,000 reads were assembled into approximately 15,000 contigs and 124,000 singletons. The assembled sequences were annotated and functionally mapped to Gene Ontology (GO) terms. Frequency of the most abundant sequences in each of the libraries was compared across all libraries to identify genes that are potentially differentially expressed during cold acclimation and fruit development. Real-time PCR was performed to confirm their differential expression patterns. Overall, 14 out of 17 of the genes examined had differential expression patterns similar to what was predicted from their reads alone. The assembled sequences were also mined for SSRs. From these sequences, 15,886 blueberry EST-SSR loci were identified. Primers were designed from 7,705 of the SSR-containing sequences with adequate flanking sequence. One hundred primer pairs were tested for amplification and polymorphism among parents of two blueberry populations currently being used for genetic linkage map construction. The tetraploid mapping population was based on a cross between the highbush cultivars Draper and Jewel (V. darrowii is also in the background of 'Jewel'). The diploid mapping population was based on a cross between an F1 hybrid of V. darrowii and diploid V. corymbosum and another diploid V. corymbosum. The overall amplification rate of the SSR primers was 68% and the polymorphism rate was 43%.

CONCLUSIONS

These results indicate that this large collection of 454 ESTs will be a valuable resource for identifying genes that are potentially differentially expressed and play important roles in flower bud development, cold acclimation, chilling unit accumulation, and fruit development in blueberry and related species. In addition, the ESTs have already proved useful for the development of SSR and EST-PCR markers, and are currently being used for construction of genetic linkage maps in blueberry.

The mRNA-destabilizing protein tristetraprolin is suppressed in many cancers, altering tumorigenic phenotypes and patient prognosis

AU-rich element-binding proteins (ARE-BP) regulate the stability and/or translational efficiency of mRNAs containing cognate binding sites. Many targeted transcripts encode factors that control processes such as cell division, apoptosis, and angiogenesis, suggesting that dysregulated ARE-BP expression could dramatically influence oncogenic phenotypes. Using several approaches, we evaluated the expression of four well-characterized ARE-BPs across a variety of human neoplastic syndromes. AUF1, TIA-1, and HuR mRNAs were not systematically dysregulated in cancers; however, tristetraprolin mRNA levels were significantly decreased across many tumor types, including advanced cancers of the breast and prostate. Restoring tristetraprolin expression in an aggressive tumor cell line suppressed three key tumorgenic phenotypes: cell proliferation, resistance to proapoptotic stimuli, and expression of vascular endothelial growth factor mRNA. However, the cellular consequences of tristetraprolin expression varied across different cell models. Analyses of gene array data sets revealed that suppression of tristetraprolin expression is a negative prognostic indicator in breast cancer, because patients with low tumor tristetraprolin mRNA levels were more likely to present increased pathologic tumor grade, vascular endothelial growth factor expression, and mortality from recurrent disease. Collectively, these data establish that tristetraprolin expression is frequently suppressed in human cancers, which in turn can alter tumorigenic phenotypes that influence patient outcomes.

Laser capture microdissection (LCM) and expression analyses of Glycine max (soybean) syncytium containing root regions formed by the plant pathogen Heterodera glycines (soybean cyst nematode)

Roots of soybean, Glycine max cv. Kent L. Merr., plants susceptible to the soybean cyst nematode (SCN), Heterodera glycines Ichinohe, were inoculated and allowed to develop feeding sites (syncytia) for 8 days. Root samples enriched in syncytial cells were collected using laser capture microdissection (LCM). RNA was extracted and used to make a cDNA library and expressed sequence tags (ESTs) were produced and used for a Gene Ontology (GO) analysis. RT-PCR results indicated enhanced expression of an aquaporin (GmPIP2,2), alpha-tubulin (GmTubA1), beta-tubulin (GmTubB4) and several other genes in syncytium-enriched samples as compared to samples extracted from whole roots. While RT-PCR data showed increased transcript levels of GmPIP2,2 from LCM tissue enriched in syncytial cells, in situ hybridization showed prominent GmPIP2,2 hybridization to RNA in the parenchymal cells tightly juxtaposed to the syncytium. Immunolocalization indicated stronger alpha-tubulin signal within the syncytium as compared to surrounding tissue. However, alpha-tubulin labeling appeared diffuse or clumped. Thus, LCM allowed for the isolation of tissue enriched for syncytial cells, providing material suitable for a variety of molecular analyses.

Microarray analysis of gene expression in soybean roots susceptible to the soybean cyst nematode two days post invasion

Soybean root cells undergo dramatic morphological and biochemical changes during the establishment of a feeding site in a compatible interaction with the soybean cyst nematode (SCN). We constructed a cDNA microarray with approximately 1,300 cDNA inserts targeted to identify differentially expressed genes during the compatible interaction of SCN with soybean roots 2 days after infection. Three independent biological replicates were grown and inoculated with SCN, and 2 days later RNA was extracted for hybridization to microarrays and compared to noninoculated controls. Statistical analysis indicated that approximately 8% of the genes monitored were induced and more than 50% of these were genes of unknown function. Notable genes that were more highly expressed 2 days after inoculation with SCN as compared to noninoculated roots included the repetitive proline-rich glycoprotein, the stress-induced gene SAM22, ss-1,3-endoglucanase, peroxidase, and those involved in carbohydrate metabolism, plant defense, and signaling.

Analysis of expressed sequence tags from roots of resistant soybean infected by the soybean cyst nematode

The soybean cyst nematode (SCN) Heterodera glycines is the most devastating pest of soybean in the U.S.A. The resistance response elicited by SCN in soybean is complex, and genes involved in the response to a large extent are unknown and not well characterized. We constructed cDNA libraries made from mRNA extracted from roots of the resistant soybean Glycine max L. Merr. 'Peking' at 12 h, 2 to 4 days, and 6 to 8 days post inoculation with the soybean cyst nematode, population NL1-RHp, similar to race 3. Expressed sequence tag analysis of the libraries provides rapid discovery of genes involved in the response of soybean to the nematode. A total of 3454 cDNA clones were examined from the three libraries, of which 25 cDNAs were derived from nematode RNA. The levels of certain stress-induced genes such as SAM22 and glutathione S-transferase (GST8) were elevated in the SCN-infected roots relative to uninoculated roots. Early defense response genes, particularly ascorbate peroxidase and lipoxygenase, were abundant in the 12-h library. By 6-8 days, the expression of most of those genes was not as abundant, whereas genes coding for unknown proteins and stress-induced proteins continued to be highly expressed. These ESTs and associated information will be useful to scientists examining gene and protein interactions between nematodes and plants.

Profiling local gene expression changes associated with Eimeria maxima and Eimeria acervulina using cDNA microarray

Eimeria parasites show preferential sites of invasion in the avian intestine and produce a species-specific host immune response. Two economically important species, Eimeria acervulina and Eimeria maxima, preferentially invade and develop in the avian duodenum and jejunum/ileum, respectively. To investigate local host immune responses induced by parasite infection, global transcriptional changes in intestinal intraepithelial lymphocytes (IELs) induced by oral inoculation of chickens with E. acervulina or E. maxima were monitored using cDNA microarrays containing 400 unique chicken genes. Multiple gene transcripts were significantly up- or down-regulated following primary or secondary infection with E. acervulina or E. maxima. In general, infection by either parasite resulted in the expression changes of more genes following primary infection than following secondary infection, and E. acervulina caused more changes than did E. maxima. Although different regions of the small intestine were infected, similar changes in the levels of several cytokine mRNAs were observed in both Eimeria species following primary infection. Also identified was a set of transcripts whose expression was commonly enhanced or repressed in intestinal IELs of chickens infected with either parasite. Microarray analysis of chicken genes induced or repressed following Eimeria infection offers a powerful tool to enhance our understanding of host-parasite interactions leading to protective immunity.

Ascaris suum: cDNA microarray analysis of 4th stage larvae (L4) during self-cure from the intestine

There is spontaneous cure of a large portion of Ascaris suum 4th-stage larvae (L4) from the jejunum of infected pigs between 14 and 21 days after inoculation (DAI). Those L4 that remain in the jejunum continue to develop while those that have moved to the ileum are eventually expelled from the intestines. Although increases in intestinal mucosal mast cells and changes in local host immunity are coincidental with spontaneous cure, the population of L4 that continue to develop in the jejunum may counteract host protective mechanisms by the differential production of factors related to parasitism. To this end, a cDNA library was constructed from L4 isolated from pig jejunum at 21 DAI, and 93% of 1920 original clones containing a single amplicon in the range 400-1500 bp were verified by gel electrophoresis and printed onto glass slides for microarray analysis. Fluorescent probes were prepared from total RNA isolated from: (1) 3rd stage-larvae from lung at 7 DAI, (L3); (2) L4 from jejunum at 14 DAI (L4-14-J); (3) L4 from jejunum at 21 DAI (L4-21-J); (4) L4 from ileum at 21 DAI (L4-21-I, and; (5) adults (L5). Cy3-labeled L3, L4-14-J, L4-21-I and L5 cDNA, and Cy5-labeled L4-21-J cDNA were simultaneously used to screen the printed arrays containing the L4-21-J-derived cDNA library. Several clones showed consistent differential gene expression over two separate experiments and were grouped into 3 distinct transcription patterns. The data showed that sequences from muscle actin and myosin, ribosomal protein L11, glyceraldehyde-3-phosphate dehydrogenase and the flavoprotein subunit of succinate dehydrogenase were highly expressed in L4-21-J, but not in L4-21-I; as were a collection of un-annotated genes derived from a worm body wall-hypodermis library, and a testes germinal zone tissue library. These results suggest that only actively developing A. suum L4 are destined to parasitize the host and successfully neutralize host protective responses.