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

Evidence of neutral transcriptome evolution in plants

The transcriptome of an organism is its set of gene transcripts (mRNAs) at a defined spatial and temporal locus. Because gene expression is affected markedly by environmental and developmental perturbations, it is widely assumed that transcriptome divergence among taxa represents adaptive phenotypic selection. This assumption has been challenged by neutral theories which propose that stochastic processes drive transcriptome evolution. To test for evidence of neutral transcriptome evolution in plants, we quantified 18 494 gene transcripts in nonsenescent leaves of 14 taxa of Brassicaceae using robust cross-species transcriptomics which includes a two-step physical and in silicobased normalization procedure based on DNA similarity among taxa. Transcriptome divergence correlates positively with evolutionary distance between taxa and with variation in gene expression among samples. Results are similar for pseudogenes and chloroplast genes evolving at different rates. Remarkably, variation in transcript abundance among root-cell samples correlates positively with transcriptome divergence among root tissues and among taxa. Because neutral processes affect transcriptome evolution in plants, many differences in gene expression among or within taxa may be nonfunctional, reflecting ancestral plasticity and founder effects. Appropriate null models are required when comparing transcriptomes in space and time

Analysis of ripening-related gene expression in papaya using an Arabidopsis-based microarray

Background Papaya (Carica papaya L.) is a commercially important crop that produces climacteric fruits with a soft and sweet pulp that contain a wide range of health promoting phytochemicals. Despite its importance, little is known about transcriptional modifications during papaya fruit ripening and their control. In this study we report the analysis of ripe papaya transcriptome by using a cross-species (XSpecies) microarray technique based on the phylogenetic proximity between papaya and Arabidopsis thaliana. Results Papaya transcriptome analyses resulted in the identification of 414 ripening-related genes with some having their expression validated by qPCR. The transcription profile was compared with that from ripening tomato and grape. There were many similarities between papaya and tomato especially with respect to the expression of genes encoding proteins involved in primary metabolism, regulation of transcription, biotic and abiotic stress and cell wall metabolism. XSpecies microarray data indicated that transcription factors (TFs) of the MADS-box, NAC and AP2/ERF gene families were involved in the control of papaya ripening and revealed that cell wall-related gene expression in papaya had similarities to the expression profiles seen in Arabidopsis during hypocotyl development. Conclusion The cross-species array experiment identified a ripening-related set of genes in papaya allowing the comparison of transcription control between papaya and other fruit bearing taxa during the ripening process

Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits

The environmental and financial costs of using inorganic phosphate fertilizers to maintain crop yield and quality are high. Breeding crops that acquire and use phosphorus (P) more efficiently could reduce these costs. The variation in shoot P concentration (shoot-P) and various measures of P use efficiency (PUE) were quantified among 355 Brassica oleracea L. accessions, 74 current commercial cultivars, and 90 doubled haploid (DH) mapping lines from a reference genetic mapping population. Accessions were grown at two or more external P concentrations in glasshouse experiments; commercial and DH accessions were also grown in replicated field experiments. Within the substantial species-wide diversity observed for shoot-P and various measures of PUE in B. oleracea, current commercial cultivars have greater PUE than would be expected by chance. This may be a consequence of breeding for increased yield, which is a significant component of most measures of PUE, or early establishment. Root development and architecture correlate with PUE; in particular, lateral root number, length, and growth rate. Significant quantitative trait loci associated with shoot-P and PUE occur on chromosomes C3 and C7. These data provide information to initiate breeding programmes to improve PUE in B. oleracea

High Resolution Melt (HRM) analysis is an efficient tool to genotype EMS mutants in complex crop genomes

Background: Targeted Induced Loci Lesions IN Genomes (TILLING) is increasingly being used to generate and identify mutations in target genes of crop genomes. TILLING populations of several thousand lines have been generated in a number of crop species including Brassica rapa. Genetic analysis of mutants identified by TILLING requires an efficient, high-throughput and cost effective genotyping method to track the mutations through numerous generations. High resolution melt (HRM) analysis has been used in a number of systems to identify single nucleotide polymorphisms (SNPs) and insertion/deletions (IN/DELs) enabling the genotyping of different types of samples. HRM is ideally suited to high-throughput genotyping of multiple TILLING mutants in complex crop genomes. To date it has been used to identify mutants and genotype single mutations. The aim of this study was to determine if HRM can facilitate downstream analysis of multiple mutant lines identified by TILLING in order to characterise allelic series of EMS induced mutations in target genes across a number of generations in complex crop genomes. Results: We demonstrate that HRM can be used to genotype allelic series of mutations in two genes, BraA.CAX1a and BraA.MET1.a in Brassica rapa. We analysed 12 mutations in BraA.CAX1.a and five in BraA.MET1.a over two generations including a back-cross to the wild-type. Using a commercially available HRM kit and the Lightscanner™ system we were able to detect mutations in heterozygous and homozygous states for both genes. Conclusions: Using HRM genotyping on TILLING derived mutants, it is possible to generate an allelic series of mutations within multiple target genes rapidly. Lines suitable for phenotypic analysis can be isolated approximately 8-9 months (3 generations) from receiving M3 seed of Brassica rapa from the RevGenUK TILLING service

Identification of candidate genes for calcium and magnesium accumulation in Brassica napus L. by association genetics

Calcium (Ca) and magnesium (Mg) are essential plant nutrients and vital for human and animal nutrition. Biofortification of crops has previously been suggested to alleviate widespread human Ca and Mg deficiencies. In this study, new candidate genes influencing the leaf accumulation of Ca and Mg were identified in young Brassica napus plants using associative transcriptomics of bionomics datasets. A total of 247 and 166 SNP markers were associated with leaf Ca and Mg concentration, respectively, after false discovery rate correction and removal of SNPs with low second allele frequency. Gene expression markers at similar positions were also associated with leaf Ca and Mg concentration, including loci on chromosomes A10 and C2, within which lie previously identified transporter genes ACA8 and MGT7. Further candidate genes were selected from seven loci and the mineral composition of whole Arabidopsis thaliana shoots were characterized from lines mutated in orthologous genes. Four and two mutant lines had reduced shoot Ca and Mg concentration, respectively, compared to wild type plants. Three of these mutations were found to have tissue specific effects; notably reduced silique Ca in all three such mutant lines. This knowledge could be applied in targeted breeding, with the possibility of increasing Ca and Mg in plant tissue for improving human and livestock nutrition

Recycled Components in Mantle Plumes Deduced From Variations in Halogens (Cl, Br, and I), Trace Elements, and He-3/He-4 Along the Hawaiian-Emperor Seamount Chain

Halogens are primarily located within surface reservoirs of the Earth; as such they have proven to be effective tracers for the identification of subducted volatiles within the mantle. Subducting lithologies exhibit a wide variety of halogen compositions, yet the mantle maintains a fairly uniform signature, suggesting halogens may be homogenized during subduction to the mantle or during eruption. Here we present halogen (Cl, Br, and I), K, noble gas, and major and trace element data on olivines from three seamounts along the Hawaiian-Emperor seamount chain to determine if the deep mantle source has retained evidence of halogen heterogeneities introduced through subduction. High Ni contents indicate that the Hawaiian-Emperor mantle source contains a recycled oceanic crust component in the form of pyroxenite, which increases from the 46% in the oldest (Detroit) to 70% in the younger seamount (Koko). Detroit seamount retains mid-ocean ridge basalts (MORB)-like Br/Cl and I/Cl, while the Br/Cl and I/Cl of Suiko and Koko seamounts are higher than MORB and similar to altered oceanic crust and dehydrated serpentinite. Helium isotopes show a similar evolution, from MORB-like values at Detroit seamount toward higher values at Suiko and Koko seamounts. The correlation between pyroxenite contributions, Br/Cl, I/Cl, and He-3/He-4 indicates that subducted material has been incorporated into the primordial undegassed Hawaiian mantle plume source. The identification of recycled oceanic crustal signatures in both the trace elements and halogens indicates that subduction and dehydration of altered oceanic crust may exert control on the cycling of volatile elements to the deep mantle

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca2+ transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 x R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca2+ transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization

Tandem quadruplication of HMA4 in the zinc (Zn) and cadmium (Cd) hyperaccumulator noccaea caerulescens

Zinc (Zn) and cadmium (Cd) hyperaccumulation may have evolved twice in the Brassicaceae, in Arabidopsis halleri and in the Noccaea genus. Tandem gene duplication and deregulated expression of the Zn transporter, HMA4, has previously been linked to Zn/Cd hyperaccumulation in A. halleri. Here, we tested the hypothesis that tandem duplication and deregulation of HMA4 expression also occurs in Noccaea. A Noccaea caerulescens genomic library was generated, containing 36,864 fosmid pCC1FOSTM clones with insert sizes ~20–40 kbp, and screened with a PCR-generated HMA4 genomic probe. Gene copy number within the genome was estimated through DNA fingerprinting and pooled fosmid pyrosequencing. Gene copy numbers within individual clones was determined by PCR analyses with novel locus specific primers. Entire fosmids were then sequenced individually and reads equivalent to 20-fold coverage were assembled to generate complete whole contigs. Four tandem HMA4 repeats were identified in a contiguous sequence of 101,480 bp based on sequence overlap identities. These were flanked by regions syntenous with up and downstream regions of AtHMA4 in Arabidopsis thaliana. Promoter-reporter b-glucuronidase (GUS) fusion analysis of a NcHMA4 in A. thaliana revealed deregulated expression in roots and shoots, analogous to AhHMA4 promoters, but distinct from AtHMA4 expression which localised to the root vascular tissue. This remarkable consistency in tandem duplication and deregulated expression of metal transport genes between N. caerulescens and A. halleri, which last shared a common ancestor >40 mya, provides intriguing evidence that parallel evolutionary pathways may underlie Zn/Cd hyperaccumulation in Brassicaceae