Cell walls of developing wheat starchy endosperm: comparison of composition and RNA-Seq transcriptome

The transcriptome of the developing starchy endosperm of hexaploid wheat (Triticum aestivum) was determined using RNA-Seq isolated at five stages during grain fill. This resource represents an excellent way to identify candidate genes responsible for the starchy endosperm cell wall, which is dominated by arabinoxylan (AX), accounting for 70% of the cell wall polysaccharides, with 20% (1,3; 1,4)-beta-D-glucan, 7% glucomannan, and 4% cellulose. A complete inventory of transcripts of 124 glycosyltransferase (GT) and 72 glycosylhydrolase (GH) genes associated with cell walls is presented. The most highly expressed GT transcript (excluding those known to be involved in starch synthesis) was a GT47 family transcript similar to Arabidopsis (Arabidopsis thaliana) IRX10 involved in xylan extension, and the second most abundant was a GT61. Profiles for GT43 IRX9 and IRX14 putative orthologs were consistent with roles in AX synthesis. Low abundances were found for transcripts from genes in the acyl-coA transferase BAHD family, for which a role in AX feruloylation has been postulated. The relative expression of these was much greater in whole grain compared with starchy endosperm, correlating with the levels of bound ferulate. Transcripts associated with callose (GSL), cellulose (CESA), pectin (GAUT), and glucomannan (CSLA) synthesis were also abundant in starchy endosperm, while the corresponding cell wall polysaccharides were confirmed as low abundance (glucomannan and callose) or undetectable (pectin) in these samples. Abundant transcripts from GH families associated with the hydrolysis of these polysaccharides were also present, suggesting that they may be rapidly turned over. Abundant transcripts in the GT31 family may be responsible for the addition of Gal residues to arabinogalactan peptide.

Mass spectroscopic characterization of the coronavirus infectious bronchitis virus nucleoprotein and elucidation of the role of phosphorylation in RNA binding by using surface plasmon resonance

Phosphorylation of the coronavirus nucleoprotein (N protein) has been predicted to play a role in RNA binding. To investigate this hypothesis, we examined the kinetics of RNA binding between nonphosphorylated and phosphorylated infectious bronchitis virus N protein with nonviral and viral RNA by surface plasmon resonance (Biacore). Mass spectroscopic analysis of N protein identified phosphorylation sites that were proximal to RNA binding domains. Kinetic analysis, by surface plasmon resonance, indicated that nonphospborylated N protein bound with the same affinity to viral RNA as phosphorylated N protein. However, phosphorylated N protein bound to viral RNA with a higher binding affinity than nonviral RNA, suggesting that phosphorylation of N protein determined the recognition of virus RNA. The data also indicated that a known N protein binding site (involved in transcriptional regulation) consisting of a conserved core sequence present near the 5' end of the genome (in the leader sequence) functioned by promoting high association rates of N protein binding. Further analysis of the leader sequence indicated that the core element was not the only binding site for N protein and that other regions functioned to promote high-affinity binding.

Transcription profiling in human platelets reveals LRRFIP1 as a novel protein regulating platelet function

Within the healthy population, there is substantial, heritable, and interindividual variability in the platelet response. We explored whether a proportion of this variability could be accounted for by interindividual variation in gene expression. Through a correlative analysis of genome-wide platelet RNA expression data from 37 subjects representing the normal range of platelet responsiveness within a cohort of 500 subjects, we identified 63 genes in which transcript levels correlated with variation in the platelet response to adenosine diphosphate and/or the collagen-mimetic peptide, cross-linked collagen-related peptide. Many of these encode proteins with no reported function in platelets. An association study of 6 of the 63 genes in 4235 cases and 6379 controls showed a putative association with myocardial infarction for COMMD7 (COMM domain-containing protein 7) and a major deviation from the null hypo thesis for LRRFIP1 [leucine-rich repeat (in FLII) interacting protein 1]. Morpholino-based silencing in Danio rerio identified a modest role for commd7 and a significant effect for lrrfip1 as positive regulators of thrombus formation. Proteomic analysis of human platelet LRRFIP1-interacting proteins indicated that LRRFIP1 functions as a component of the platelet cytoskeleton, where it interacts with the actin-remodeling proteins Flightless-1 and Drebrin. Taken together, these data reveal novel proteins regulating the platelet response.

An optimized method for the extraction of bacterial mRNA from plant roots infected with Escherichia coli O157:H7

Analysis of microbial gene expression during host colonization provides valuable information on the nature of interaction, beneficial or pathogenic, and the adaptive processes involved. Isolation of bacterial mRNA for in planta analysis can be challenging where host nucleic acid may dominate the preparation, or inhibitory compounds affect downstream analysis, e.g., quantitative reverse transcriptase PCR (qPCR), microarray, or RNA-seq. The goal of this work was to optimize the isolation of bacterial mRNA of food-borne pathogens from living plants. Reported methods for recovery of phytopathogen-infected plant material, using hot phenol extraction and high concentration of bacterial inoculation or large amounts of infected tissues, were found to be inappropriate for plant roots inoculated with Escherichia coli O157:H7. The bacterial RNA yields were too low and increased plant material resulted in a dominance of plant RNA in the sample. To improve the yield of bacterial RNA and reduce the number of plants required, an optimized method was developed which combines bead beating with directed bacterial lysis using SDS and lysozyme. Inhibitory plant compounds, such as phenolics and polysaccharides, were counteracted with the addition of high-molecular-weight polyethylene glycol and hexadecyltrimethyl ammonium bromide. The new method increased the total yield of bacterial mRNA substantially and allowed assessment of gene expression by qPCR. This method can be applied to other bacterial species associated with plant roots, and also in the wider context of food safety.

Global transcriptome analysis and identification of differentially expressed genes after infection of Fragaria vesca with powdery mildew (Podosphaera aphanis)

Background: Podosphaera aphanis, the causal agent of strawberry powdery mildew causes significant economic loss worldwide. Methods: We used the diploid strawberry species Fragaria vesca as a model to study plant pathogen interactions. RNA-seq was employed to generate a transcriptome dataset from two accessions, F. vesca ssp. vesca Hawaii 4 (HW) and F. vesca f. semperflorens Yellow Wonder 5AF7 (YW) at 1 d (1 DAI) and 8 d (8 DAI) after infection. Results: Of the total reads identified about 999 million (92%) mapped to the F. vesca genome. These transcripts were derived from a total of 23,470 and 23,464 genes in HW and YW, respectively from the three time points (control, 1 and 8 DAI). Analysis identified 1,567, 1,846 and 1,145 up-regulated genes between control and 1 DAI, control and 8 DAI, and 1 and 8 DAI, respectively in HW. Similarly, 1,336, 1,619 and 968 genes were up-regulated in YW. Also 646, 1,098 and 624 down-regulated genes were identified in HW, while 571, 754 and 627 genes were down-regulated in YW between all three time points, respectively. Conclusion: Investigation of differentially expressed genes (log2 fold changes �5) between control and 1 DAI in both HW and YW identified a large number of genes related to secondary metabolism, signal transduction; transcriptional regulation and disease resistance were highly expressed. These included flavonoid 3´-monooxygenase, peroxidase 15, glucan endo-1,3-β-glucosidase 2, receptor-like kinases, transcription factors, germin-like proteins, F-box proteins, NB-ARC and NBS-LRR proteins. This is the first application of RNA-seq to any pathogen interaction in strawberry

From insect to man: Photorhabdus sheds light on the emergence of human pathogenicity

Photorhabdus are highly effective insect pathogenic bacteria that exist in a mutualistic relationship with Heterorhabditid nematodes. Unlike other members of the genus, Photorhabdus asymbiotica can also infect humans. Most Photorhabdus cannot replicate above 34°C, limiting their host-range to poikilothermic invertebrates. In contrast, P. asymbiotica must necessarily be able to replicate at 37°C or above. Many well-studied mammalian pathogens use the elevated temperature of their host as a signal to regulate the necessary changes in gene expression required for infection. Here we use RNA-seq, proteomics and phenotype microarrays to examine temperature dependent differences in transcription, translation and phenotype of P. asymbiotica at 28°C versus 37°C, relevant to the insect or human hosts respectively. Our findings reveal relatively few temperature dependant differences in gene expression. There is however a striking difference in metabolism at 37°C, with a significant reduction in the range of carbon and nitrogen sources that otherwise support respiration at 28°C. We propose that the key adaptation that enables P. asymbiotica to infect humans is to aggressively acquire amino acids, peptides and other nutrients from the human host, employing a so called “nutritional virulence” strategy. This would simultaneously cripple the host immune response while providing nutrients sufficient for reproduction. This might explain the severity of ulcerated lesions observed in clinical cases of Photorhabdosis. Furthermore, while P. asymbiotica can invade mammalian cells they must also resist immediate killing by humoral immunity components in serum. We observed an increase in the production of the insect Phenol-oxidase inhibitor Rhabduscin normally deployed to inhibit the melanisation immune cascade. Crucially we demonstrated this molecule also facilitates protection against killing by the alternative human complement pathway.

Phylogenomics of eukaryotes: Impact of missing data on large alignments

Resolving the relationships between Metazoa and other eukaryotic groups as well as between metazoan phyla is central to the understanding of the origin and evolution of animals. The current view is based on limited data sets, either a single gene with many species (e.g., ribosomal RNA) or many genes but with only a few species. Because a reliable phylogenetic inference simultaneously requires numerous genes and numerous species, we assembled a very large data set containing 129 orthologous proteins (similar to30,000 aligned amino acid positions) for 36 eukaryotic species. Included in the alignments are data from the choanoflagellate Monosiga ovata, obtained through the sequencing of about 1,000 cDNAs. We provide conclusive support for choanoflagellates as the closest relative of animals and for fungi as the second closest. The monophyly of Plantae and chromalveolates was recovered but without strong statistical support. Within animals, in contrast to the monophyly of Coelomata observed in several recent large-scale analyses, we recovered a paraphyletic Coelamata, with nematodes and platyhelminths nested within. To include a diverse sample of organisms, data from EST projects were used for several species, resulting in a large amount of missing data in our alignment (about 25%). By using different approaches, we verify that the inferred phylogeny is not sensitive to these missing data. Therefore, this large data set provides a reliable phylogenetic framework for studying eukaryotic and animal evolution and will be easily extendable when large amounts of sequence information become available from a broader taxonomic range.

SARS coronavirus unique domain: three-domain molecular architecture in solution and RNA binding

Nonstructural protein 3 of the severe acute respiratory syndrome (SARS) coronavirus includes a "SARS-unique domain" (SUD) consisting of three globular domains separated by short linker peptide segments. This work reports NMR structure determinations of the C-terminal domain (SUD-C) and a two-domain construct (SUD-MC) containing the middle domain (SUD-M) and the C-terminal domain, and NMR data on the conformational states of the N-terminal domain (SUD-N) and the SUD-NM two-domain construct. Both SUD-N and SUD-NM are monomeric and globular in solution; in SUD-NM, there is high mobility in the two-residue interdomain linking sequence, with no preferred relative orientation of the two domains. SUD-C adopts a frataxin like fold and has structural similarity to DNA-binding domains of DNA-modifying enzymes. The structures of both SUD-M (previously determined) and SUD-C (from the present study) are maintained in SUD-MC, where the two domains are flexibly linked. Gel-shift experiments showed that both SUD-C and SUD-MC bind to single-stranded RNA and recognize purine bases more strongly than pyrimidine bases, whereby SUD-MC binds to a more restricted set of purine-containing RNA sequences than SUD-M. NMR chemical shift perturbation experiments with observations of (15)N-labeled proteins further resulted in delineation of RNA binding sites (i.e., in SUD-M, a positively charged surface area with a pronounced cavity, and in SUD-C, several residues of an anti-parallel beta-sheet). Overall, the present data provide evidence for molecular mechanisms involving the concerted actions of SUD-M and SUD-C, which result in specific RNA binding that might be unique to the SUD and, thus, to the SARS coronavirus.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.