Isolation and characterization of an AINTEGUMENTA-like gene in different coconut (Cocos nucifera L.) varieties from Sri Lanka

Development of an efficient tissue culture protocol in coconut is hampered by numerous technical constraints. Thus a greater understanding of the fundamental aspects of embryogenesis is essential. The role of AINTEGUMENTA-like genes in embryogenesis has been elucidated not only in model plants but also in economically important crops. A coconut gene, CnANT, that encodes two APETALA2 (AP2) domains and a conserved linker region similar to those of the BABY BOOM transcription factor was cloned, characterized, and its tissue specific expression was examined. The full-length cDNA of 1,780 bp contains a 1,425-bp open reading frame that encodes a putative peptide of 474 amino acids. The genomic DNA sequence includes 2,317 bp and consists of nine exons interrupted by eight introns. The exon/intron organization of CnANT is similar to that of homologous genes in other plant species. Analysis of differential tissue expression by real-time polymerase chain reaction indicated that CnANT is expressed more highly in in vitro grown tissues than in other vegetative tissues. Sequence comparison of the genomic sequence of CnANT in different coconut varieties revealed one single nucleotide polymorphism and one indel in the first exon and first intron, respectively, which differentiate the Tall group of trees from Dwarfs. The indel sequence, which can be considered a simple sequence repeats marker, was successfully used to distinguish the Tall and Dwarf groups as well as to develop a marker system, which may be of value in the identification of parental varieties that are used in coconut breeding programs in Sri Lanka.

Construction of a Festuca pratensis BAC library for map-based cloning in Festulolium substitution lines

Introgression in Festulolium is a potentially powerful tool to isolate genes for a large number of traits which differ between Festuca pratensis Huds. and Lolium perenne L. Not only are hybrids between the two species fertile, but the two genomes can be distinguished by genomic in situ hybridisation and a high frequency of recombination occurs between homoeologous chromosomes and chromosome segments. By a programme of introgression and a series of backcrosses, L. perenne lines have been produced which contain small F. pratensis substitutions. This material is a rich source of polymorphic markers targeted towards any trait carried on the F. pratensis substitution not observed in the L. perenne background. We describe here the construction of an F. pratensis BAC library, which establishes the basis of a map-based cloning strategy in L. perenne. The library contains 49,152 clones, with an average insert size of 112 kbp, providing coverage of 2.5 haploid genome equivalents. We have screened the library for eight amplified fragment length polymorphism (AFLP) derived markers known to be linked to an F. pratensis gene introgressed into L. perenne and conferring a staygreen phenotype as a consequence of a mutation in primary chlorophyll catabolism. While for four of the markers it was possible to identify bacterial artificial chromosome (BAC) clones, the other four AFLPs were too repetitive to enable reliable identification of locus-specific BACs. Moreover, when the four BACs were partially sequenced, no obvious coding regions could be identified. This contrasted to BACs identified using cDNA sequences, when multiple genes were identified on the same BAC.

Error, reproducibility and sensitivity: a pipeline for data processing of Agilent oligonucleotide expression arrays

Background: Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples. Results: We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2 of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log(2) units (6 of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators. Conclusions: This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells.

Changes in gene expression in Arabidopsis shoots during phosphate starvation and the potential for developing smart plants

Our aim was to generate and prove the concept of "smart" plants to monitor plant phosphorus (P) status in Arabidopsis. Smart plants can be genetically engineered by transformation with a construct containing the promoter of a gene up-regulated specifically by P starvation in an accessible tissue upstream of a marker gene such as beta-glucuronidase (GUS). First, using microarrays, we identified genes whose expression changed more than 2.5-fold in shoots of plants growing hydroponically when P, but not N or K, was withheld from the nutrient solution. The transient changes in gene expression occurring immediately (4 h) after P withdrawal were highly variable, and many nonspecific, shock-induced genes were up-regulated during this period. However, two common putative cis-regulatory elements (a PHO-like element and a TATA box-like element) were present significantly more often in the promoters of genes whose expression increased 4 h after the withdrawal of P compared with their general occurrence in the promoters of all genes represented on the microarray. Surprisingly, the expression of only four genes differed between shoots of P-starved and -replete plants 28 h after P was withdrawn. This lull in differential gene expression preceded the differential expression of a new group of 61 genes 100 h after withdrawing P. A literature survey indicated that the expression of many of these "late" genes responded specifically to P starvation. Shoots had reduced P after 100 h, but growth was unaffected. The expression of SQD1, a gene involved in the synthesis of sulfolipids, responded specifically to P starvation and was increased 100 h after withdrawing P. Leaves of Arabidopsis bearing a SQD1::GUS construct showed increased GUS activity after P withdrawal, which was detectable before P starvation limited growth. Hence, smart plants can monitor plant P status. Transferring this technology to crops would allow precision management of P fertilization, thereby maintaining yields while reducing costs, conserving natural resources, and preventing pollution.

Gene expression changes in phosphorus deficient potato (Solanum tuberosum L.) leaves and the potential for diagnostic gene expression markers

Background: There are compelling economic and environmental reasons to reduce our reliance on inorganic phosphate (Pi) fertilisers. Better management of Pi fertiliser applications is one option to improve the efficiency of Pi fertiliser use, whilst maintaining crop yields. Application rates of Pi fertilisers are traditionally determined from analyses of soil or plant tissues. Alternatively, diagnostic genes with altered expression under Pi limiting conditions that suggest a physiological requirement for Pi fertilisation, could be used to manage Pifertiliser applications, and might be more precise than indirect measurements of soil or tissue samples. Results: We grew potato (Solanum tuberosum L.) plants hydroponically, under glasshouse conditions, to control their nutrient status accurately. Samples of total leaf RNA taken periodically after Pi was removed from the nutrient solution were labelled and hybridised to potato oligonucleotide arrays. A total of 1,659 genes were significantly differentially expressed following Pi withdrawal. These included genes that encode proteins involved in lipid, protein, and carbohydrate metabolism, characteristic of Pi deficient leaves and included potential novel roles for genes encoding patatin like proteins in potatoes. The array data were analysed using a support vector machine algorithm to identify groups of genes that could predict the Pi status of the crop. These groups of diagnostic genes were tested using field grown potatoes that had either been fertilised or unfertilised. A group of 200 genes could correctly predict the Pi status of field grown potatoes. Conclusions: This paper provides a proof-of-concept demonstration for using microarrays and class prediction tools to predict the Pi status of a field grown potato crop. There is potential to develop this technology for other biotic and abiotic stresses in field grown crops. Ultimately, a better understanding of crop stresses may improve our management of the crop, improving the sustainability of agriculture.

Cardiac protein kinases: the cardiomyocyte kinome and differential kinase expression in human failing hearts

Aims. Protein kinases are potential therapeutic targets for heart failure, but most studies of cardiac protein kinases derive from other systems, an approach that fails to account for specific kinases expressed in the heart and the contractile cardiomyocytes. We aimed to define the cardiomyocyte kinome (i.e. the protein kinases expressed in cardiomyocytes) and identify kinases with altered expression in human failing hearts. Methods and Results. Expression profiling (Affymetrix microarrays) detected >400 protein kinase mRNAs in rat neonatal ventricular myocytes (NVMs) and/or adult ventricular myocytes (AVMs), 32 and 93 of which were significantly upregulated or downregulated (>2-fold), respectively, in AVMs. Data for AGC family members were validated by qPCR. Proteomics analysis identified >180 cardiomyocyte protein kinases, with high relative expression of mitogen-activated protein kinase cascades and other known cardiomyocyte kinases (e.g. CAMKs, cAMP-dependent protein kinase). Other kinases are poorly-investigated (e.g. Slk, Stk24, Oxsr1). Expression of Akt1/2/3, BRaf, ERK1/2, Map2k1, Map3k8, Map4k4, MST1/3, p38-MAPK, PKCδ, Pkn2, Ripk1/2, Tnni3k and Zak was confirmed by immunoblotting. Relative to total protein, Map3k8 and Tnni3k were upregulated in AVMs vs NVMs. Microarray data for human hearts demonstrated variation in kinome expression that may influence responses to kinase inhibitor therapies. Furthermore, some kinases were upregulated (e.g. NRK, JAK2, STK38L) or downregulated (e.g. MAP2K1, IRAK1, STK40) in human failing hearts. Conclusions. This characterization of the spectrum of kinases expressed in cardiomyocytes and the heart (cardiomyocyte and cardiac kinomes) identified novel kinases, some of which are differentially expressed in failing human hearts and could serve as potential therapeutic targets.

Arabinose and protocatechuate catabolism genes are important for growth of Rhizobium leguminosarum biovar viciae in the pea rhizosphere

Background and aims: To form nitrogen-fixing nodules on pea roots, Rhizobium leguminosarum biovar viciae must be competitive in the rhizosphere. Our aim was to identify genes important for rhizosphere fitness. Methods: Signature-tagged mutants were screened using microarrays to identify mutants reduced for growth in pea rhizospheres. Candidate mutants were assessed relative to controls for growth in minimal medium, growth in pea rhizospheres and for infection of peas in mixed inoculants. Mutated genes were identified by DNA sequencing and confirmed by transduction. Results: Of 5508 signature-tagged mutants, microarrays implicated 50 as having decreased rhizosphere fitness. Growth tests identified six mutants with rhizosphere-specific phenotypes. The mutation in one of the genes (araE) was in an arabinose catabolism operon and blocked growth on arabinose. The mutation in another gene (pcaM), encoding a predicted solute binding protein for protocatechuate and hydroxybenzoate uptake, decreased growth on protocatechuate. Both mutants were decreased for nodule infection competitiveness with mixed inoculants, but nodulated peas normally when inoculated alone. Other mutants with similar phenotypes had mutations predicted to affect secondary metabolism. Conclusions: Catabolism of arabinose and protocatechuate in the pea rhizosphere is important for competitiveness of R.l. viciae. Other genes predicted to be involved in secondary metabolism are also important.

Relationships Among Murine CD11chigh Dendritic Cell Subsets as Revealed by Baseline Gene Expression Patterns

The functional relationships and properties of different subtypes of dendritic cells (DC) remain largely undefined. To better characterize these cells, we used global gene analysis to determine gene expression patterns among murine CD11c(high) DC subsets. CD4(+), CD8alpha(+), and CD8alpha(-) CD4(-) (double negative (DN)) DC were purified from spleens of normal C57/BL6 mice and analyzed using Affymetrix microarrays. The CD4(+) and CD8alpha(+) DC subsets showed distinct basal expression profiles differing by >200 individual genes. These included known DC subset markers as well as previously unrecognized, differentially expressed CD Ags such as CD1d, CD5, CD22, and CD72. Flow cytometric analysis confirmed differential expression in nine of nine cases, thereby validating the microarray analysis. Interestingly, the microarray expression profiles for DN cells strongly resembled those of CD4(+) DC, differing from them by <25 genes. This suggests that CD4(+) and DN DC are closely related phylogenetically, whereas CD8alpha(+) DC represent a more distant lineage, supporting the historical distinction between CD8alpha(+) and CD8alpha(-) DC. However, staining patterns revealed that in contrast to CD4(+) DC, the DN subset is heterogeneous and comprises at least two subpopulations. Gene Ontology and literature mining analyses of genes expressed differentially among DC subsets indicated strong associations with immune response parameters as well as cell differentiation and signaling. Such associations offer clues to possible unique functions of the CD11c(high) DC subsets that to date have been difficult to define as rigid distinctions.

Transcriptomic analysis of Enterohaemorrhagic Escherichia coli O157:H7 in response to plant extracts

Enterohaemorrhagic Escherichia coli (EHEC) are a group of food and contact-borne pathogens responsible for haemorrhagic colitis. The bacteria can be transmitted by contaminated meat, but importantly, also by plants. The bacteria can use plants as an alternative host, where they associate with both the leaves and the roots. Colonisation in the rhizosphere of plants is thought to be the main habitat for colonisation. Four different plant species, commonly associated with EHEC outbreaks, were infected with EHEC O157:H7 isolates Sakai and TUV 93-0 over ten days to assess the colonisation potential of the bacteria in both the phyllosphere and rhizosphere of plants. The rhizosphere was found to sustain a higher population level of bacteria over time in comparison to the phyllosphere, yet both strains were unable to utilize root exudates for growth. Global gene expression changes of EHEC O157:H7 strain Sakai were measured in response to plant extracts such as leaf lysates, root exudates and leaf cell wall polysaccharides from spinach cultivar Amazon and lettuce cultivar Salinas. Microarrays analysis showed a significant change in expression of 17 % of genes on exposure to leaf lysates of spinach. A more specific response was seen to spinach leaf cell wall polysaccharides with only a 1.5 % change. In contrast, when exposed to lettuce leaf cell wall polysaccharides a higher change of 4.8 % was seen. Genes that were differentially expressed belonged to multiple functional groups, including metabolism, indicating the utilization of plant-specific polysaccharides. Several areas of further investigation have been determined from this project, including the importance of culturing bacterial strains at a relevant temperature, the proposed lack of the type III secretion system in plant colonization by EHEC O157:H7 and the utilization of plant components for growth and persistence in the plant environment.

Transcriptome analysis for discovering candidate genes involve in embryogenesis in coconut (Cocos nucifera L.) through 454 pyrosequencing

Coconut, Cocos nucifera L. is a major plantation crop, which ensures income for millions of people in the tropical region. Detailed molecular studies on zygotic embryo development would provide valuable clues for the identification of molecular markers to improve somatic embryogenesis. Since there is no ongoing genome project for this species, coconut expressed sequence tags (EST) would be an interesting technique to identify important coconut embryo specific genes as well as other functional genes in different biochemical pathways. The goal of this study was to analyse the ESTs by examining the transcriptome data of the different embryo tissue types together with one somatic tissue. Here, four cDNA libraries from immature embryo, mature embryo, microspore derived embryo and mature leaves were constructed. cDNA was sequenced by the Roche-454 GS-FLX system and assembled into 32621 putative unigenes and 155017 singletons. Of these unigenes, 18651 had significant sequence similarities to non-redundant protein database, from which 16153 were assigned to one or more gene ontology categories. Homologue genes, which are responsible for embryo development such as chitinase, beta-1,3-glucanase, ATP synthase CF0 subunit, thaumatin-like protein and metallothionein-like protein were identified among the embryo EST collection. Of the unigenes, 6694 were mapped into 139 KEGG pathways including carbohydrate metabolism, energy metabolism, lipid metabolism, amino acid metabolism and nucleotide metabolism. This collection of 454-derived EST data generated from different tissue types provides a significant resource for genome wide studies and gene discovery of coconut, a non-model species.

Using U0126 to dissect the role of the extracellular signal-regulated kinase 1/2 (ERK1/2) cascade in the regulation of gene expression by endothelin-1 in cardiac myocytes

The hypertrophic agonist endothelin-1 rapidly but transiently activates the extracellular signal-regulated kinase 1/2 (ERK1/2) cascade (and other signalling pathways) in cardiac myocytes, but the events linking this to hypertrophy are not understood. Using Affymetrix rat U34A microarrays, we identified the short-term (2-4 h) changes in gene expression induced in neonatal myocytes by endothelin-1 alone or in combination with the ERK1/2 cascade inhibitor, U0126. Expression of 15 genes was significantly changed by U0126 alone, and expression of an additional 78 genes was significantly changed by endothelin-1. Of the genes upregulated by U0126, four are classically induced through the aryl hydrocarbon receptor (AhR) by dioxins suggesting that U0126 activates the xenobiotic response element in cardiac myocytes potentially independently of effects on ERK1/2 signalling. The 78 genes showing altered expression with endothelin-1 formed five clusters: (i) three clusters showing upregulation by endothelin-1 according to time course (4 h > 2 h; 2 h > 4 h; 2 h approximately 4 h) with at least partial inhibition by U0126; (ii) a cluster of 11 genes upregulated by endothelin-1 but unaffected by U0126 suggesting regulation through signalling pathways other than ERK1/2; (iii) a cluster of six genes downregulated by endothelin-1 with attenuation by U0126. Thus, U0126 apparently activates the AhR in cardiac myocytes (which must be taken into account in protracted studies), but careful analysis allows identification of genes potentially regulated acutely via the ERK1/2 cascade. Our data suggest that the majority of changes in gene expression induced by endothelin-1 are mediated by the ERK1/2 cascade.

Cardiac myocyte gene expression profiling during H2O2-induced apoptosis

High levels of oxidative stress promote cardiac myocyte death, though lower levels are potentially cytoprotective/anabolic. We examined the changes in gene expression in rat neonatal cardiac myocytes exposed to apoptotic (0.2 mM) or nontoxic (0.04 mM) concentrations of H2O2 (2, 4, or 24 h) using Affymetrix microarrays. Using U34B arrays, we identified a ubiquitously expressed, novel H2O2-responsive gene [putative peroxide-inducible transcript 1 (Perit1)], which generates two alternatively spliced transcripts. Using 230 2.0 arrays, H2O2 (0.04 mM) promoted significant changes in expression of only 32 genes, all of which were seen with 0.2 mM H2O2. We failed to detect any increase in the rate of protein synthesis in cardiac myocytes exposed to <0.1 mM H2O2, further suggesting that global, low concentrations of H2O2 are not anabolic in this system. H2O2 (0.2 mM) promoted significant (P < 0.05, >1.75-fold) changes in expression of 649 mRNAs and 187 RNAs corresponding to no established gene. Of the mRNAs, 114 encoded transcriptional regulators including Krüppel-like factors (Klfs). Quantitative PCR independently verified the changes in Klf expression. Thus, H2O2-induced cardiac myocyte apoptosis is associated with dynamic changes in gene expression. The expression of these genes and their protein products potentially influences the progression of the apoptotic response.

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.