Julia Weiss

Validation of reference genes for quantitative real-time PCR during leaf and flower development in Petunia hybrida

BackgroundIdentification of genes with invariant levels of gene expression is a prerequisite for validating transcriptomic changes accompanying development. Ideally expression of these genes should be independent of the morphogenetic process or environmental condition tested as well as the methods used for RNA purification and analysis.ResultsIn an effort to identify endogenous genes meeting these criteria nine reference genes (RG) were tested in two Petunia lines (Mitchell and V30). Growth conditions differed in Mitchell and V30, and different methods were used for RNA isolation and analysis. Four different software tools were employed to analyze the data. We merged the four outputs by means of a non-weighted unsupervised rank aggregation method. The genes identified as optimal for transcriptomic analysis of Mitchell and V30 were EF1α in Mitchell and CYP in V30, whereas the least suitable gene was GAPDH in both lines.ConclusionsThe least adequate gene turned out to be GAPDH indicating that it should be rejected as reference gene in Petunia. The absence of correspondence of the best-suited genes suggests that assessing reference gene stability is needed when performing normalization of data from transcriptomic analysis of flower and leaf development.

Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida

Petunia hybrida is a popular bedding plant that has a long history as a genetic model system. We report the whole-genome sequencing and assembly of inbred derivatives of its two wild parents, P. axillaris N and P. inflata S6. The assemblies include 91.3% and 90.2% coverage of their diploid genomes (1.4 Gb; 2n = 14) containing 32,928 and 36,697 protein-coding genes, respectively. The genomes reveal that the Petunia lineage has experienced at least two rounds of hexaploidization: the older gamma event, which is shared with most Eudicots, and a more recent Solanaceae event that is shared with tomato and other solanaceous species. Transcription factors involved in the shift from bee to moth pollination reside in particularly dynamic regions of the genome, which may have been key to the remarkable diversity of floral colour patterns and pollination systems. The high-quality genome sequences will enhance the value of Petunia as a model system for research on unique biological phenomena such as small RNAs, symbiosis, self-incompatibility and circadian rhythms.

Quantitative evaluation of bias in PCR amplification and next-generation sequencing derived from metabarcoding samples

Unbiased identification of organisms by PCR reactions using universal primers followed by DNA sequencing assumes positive amplification. We used six universal loci spanning 48 plant species and quantified the bias at each step of the identification process from end point PCR to next-generation sequencing. End point amplification was significantly different for single loci and between species. Quantitative PCR revealed that Cq threshold for various loci, even within a single DNA extraction, showed 2,000-fold differences in DNA quantity after amplification. Next-generation sequencing (NGS) experiments in nine species showed significant biases towards species and specific loci using adaptor-specific primers. NGS sequencing bias may be predicted to some extent by the Cq values of qPCR amplification.

pcrEfficiency: a Web tool for PCR amplification efficiency prediction.

BackgroundRelative calculation of differential gene expression in quantitative PCR reactions requires comparison between amplification experiments that include reference genes and genes under study. Ignoring the differences between their efficiencies may lead to miscalculation of gene expression even with the same starting amount of template. Although there are several tools performing PCR primer design, there is no tool available that predicts PCR efficiency for a given amplicon and primer pair.ResultsWe have used a statistical approach based on 90 primer pair combinations amplifying templates from bacteria, yeast, plants and humans, ranging in size between 74 and 907 bp to identify the parameters that affect PCR efficiency. We developed a generalized additive model fitting the data and constructed an open source Web interface that allows the obtention of oligonucleotides optimized for PCR with predicted amplification efficiencies starting from a given sequence.ConclusionspcrEfficiency provides an easy-to-use web interface allowing the prediction of PCR efficiencies prior to web lab experiments thus easing quantitative real-time PCR set-up. A web-based service as well the source code are provided freely at http://srvgen.upct.es/efficiency.html under the GPL v2 license.

Conserved and Divergent Rhythms of Crassulacean Acid Metabolism-Related and Core Clock Gene Expression in the Cactus Opuntia ficus-indica

The cactus Opuntia ficus-indica is a constitutive Crassulacean acid metabolism (CAM) species. Current knowledge of CAM metabolism suggests that the enzyme phosphoenolpyruvate carboxylase kinase (PPCK) is circadian regulated at the transcriptional level, whereas phosphoenolpyruvate carboxylase (PEPC), malate dehydrogenase (MDH), NADP-malic enzyme (NADP-ME), and pyruvate phosphate dikinase (PPDK) are posttranslationally controlled. As little transcriptomic data are available from obligate CAM plants, we created an expressed sequence tag database derived from different organs and developmental stages. Sequences were assembled, compared with sequences in the National Center for Biotechnology Information nonredundant database for identification of putative orthologs, and mapped using Kyoto Encyclopedia of Genes and Genomes Orthology and Gene Ontology. We identified genes involved in circadian regulation and CAM metabolism for transcriptomic analysis in plants grown in long days. We identified stable reference genes for quantitative polymerase chain reaction and found that OfiSAND, like its counterpart in Arabidopsis (Arabidopsis thaliana), and OfiTUB are generally appropriate standards for use in the quantification of gene expression in O. ficus-indica. Three kinds of expression profiles were found: transcripts of OfiPPCK oscillated with a 24-h periodicity; transcripts of the light-active OfiNADP-ME and OfiPPDK genes adapted to 12-h cycles, while transcript accumulation patterns of OfiPEPC and OfiMDH were arrhythmic. Expression of the circadian clock gene OfiTOC1, similar to Arabidopsis, oscillated with a 24-h periodicity, peaking at night. Expression of OfiCCA1 and OfiPRR9, unlike in Arabidopsis, adapted best to a 12-h rhythm, suggesting that circadian clock gene interactions differ from those of Arabidopsis. Our results indicate that the evolution of CAM metabolism could be the result of modified circadian regulation at both the transcriptional and posttranscriptional levels.

FORMOSA controls cell division and expansion during floral development in Antirrhinum majus.

Control of organ size is the product of coordinated cell division and expansion. In plants where one of these pathways is perturbed, organ size is often unaffected as compensation mechanisms are brought into play. The number of founder cells in organ primordia, dividing cells, and the period of cell proliferation determine cell number in lateral organs. We have identified the Antirrhinum FORMOSA (FO) gene as a specific regulator of floral size. Analysis of cell size and number in the fo mutant, which has increased flower size, indicates that FO is an organ-specific inhibitor of cell division and activator of cell expansion. Increased cell number in fo floral organs correlated with upregulation of genes involved in the cell cycle. In Arabidopsis the AINTEGUMENTA (ANT) gene promotes cell division. In the fo mutant increased cell number also correlates with upregulation of an Antirrhinum ANT-like gene (Am-ANT) in inflorescences that is very closely related to ANT and shares a similar expression pattern, suggesting that they may be functional equivalents. Increased cell proliferation is thought to be compensated for by reduced cell expansion to maintain organ size. In Arabidopsis petal cell expansion is inhibited by the BIGPETAL (BPE) gene, and in the fo mutant reduced cell size corresponded to upregulation of an Antirrhinum BPE-like gene (Am-BPE). Our data suggest that FO inhibits cell proliferation by negatively regulating Am-ANT, and acts upstream of Am-BPE to coordinate floral organ size. This demonstrates that organ size is modulated by the organ-specific control of both general and local gene networks.

Application of whole genome amplification and quantitative PCR for detection and quantification of spoilage yeasts in orange juice

Small cell numbers in complex food matrices and undefined PCR inhibitors often limit detection and identification of DNA species by molecular techniques. Thus in many industrial situations enrichment growths are performed. However, growth speed of different species in complex microbial mixtures in defined media is in most cases different, thus final results do not always reflect the starting situation. We tested DNA-strand displacement whole genome amplification as a possible substitute of enrichment growth. Using whole genome amplification on orange juice contaminated with Saccharomyces cerevisiae, we lowered detection level from 10(6) down to 10(2) cfu/ml. Whole genome amplification showed to be linear (R=0.992) and the relative yeast DNA copy number compared to other DNA templates did not change thus allowing quantitative estimation of initial contamination by quantitative PCR. Using melting point analysis, we were able to distinguish between the PCR products of the 5.8S-ITS region, obtained with universal primers from pure cultures of S. cerevisiae and Hanseniaspora uvarum, two major spoilage yeasts in orange juice and forming part of wine microbiota during fermentation. However, in mixed-contaminated samples, identification of both species was hampered by preferential appearance of the melting peak coinciding with H. uvarum, except when S. cerevisiae was the dominating species. Application of whole genome amplification did not prevent the preferential detection of H. uvarum. This handicap was resolved by applying an enrichment procedure up to saturation after which the melting peak of both species could clearly be identified.

Quantitative levels of Deficiens and Globosa during late petal development show a complex transcriptional network topology of B function

The transcriptional network topology of B function in Antirrhinum, required for petal and stamen development, is thought to rely on initial activation of transcription of DEFICIENS (DEF) and GLOBOSA (GLO), followed by a positive autoregulatory loop maintaining gene expression levels. Here, we show that the mutant compacta (co), whose vegetative growth and petal size are affected, plays a role in B function. Late events in petal morphogenesis such as development of conical cell area and scent emissions were reduced in co and def (nicotianoides) (def (nic) ), and absent in co def (nic) double mutants, suggesting a role for CO in petal identity. Expression of DEF was down-regulated in co but surprisingly GLO was not affected. We investigated the levels of DEF and GLO at late stages of petal development in the co, def (nic) and glo-1 mutants, and established a reliable transformation protocol that yielded RNAi-DEF lines. We show that the threshold levels of DEF or GLO required to obtain petal tissue are approximately 11% of wild-type. The relationship between DEF and GLO transcripts is not equal or constant and changes during development. Furthermore, down-regulation of DEF or GLO does not cause parallel down-regulation of the partner. Our results demonstrate that, at late stages of petal development, the B function transcriptional network topology is not based on positive autoregulation, and has additional components of transcriptional maintenance. Our results suggest changes in network topology that may allow changes in protein complexes that would explain the fact that not all petal traits appear early in development.

The mutants compacta ahnlich, Nitida and Grandiflora define developmental compartments and a compensation mechanism in floral development in Antirrhinum majus

In order to improve our understanding of floral size control we characterised three mutants of Antirrhinummajus with different macroscopic floral phenotypes. The recessive mutant compacta ähnlich has smaller flowers affected mainly in petal lobe expansion, the dominant mutant Grandiflora has overall larger organs, whilst the semidominant mutation Nitida exhibits smaller flowers in a dose-dependent manner. We developed a cell map in order to establish the cellular phenotypes of the mutants. Changes in organ size were both organ- and region-specific. Nitida and compacta ähnlich affected cell expansion in proximal and distal petal regions, respectively, suggesting differential regulation between petal lobe regions. Although petal size was smaller in compacta ähnlich than in wild type, conical cells were significantly bigger, suggesting a compensation mechanism involved in petal development. Grandiflora had larger cells in petals and increased cell division in stamens and styles, suggesting a relationship between genes controlling organ size and organ identity. The level of ploidy in petals of Grandiflora and coan was found to be equivalent to wild type petals and leaves, ruling out an excess of growth via endoreduplication. We discuss our results in terms of current models about control of lateral organ size.

A molecular recombination map of Antirrhinum majus.

BackgroundGenetic recombination maps provide important frameworks for comparative genomics, identifying gene functions, assembling genome sequences and for breeding. The molecular recombination map currently available for the model eudicot Antirrhinum majus is the result of a cross with Antirrhinum molle, limiting its usefulness within A. majus.ResultsWe created a molecular linkage map of A. majus based on segregation of markers in the F2 population of two inbred lab strains of A. majus. The resulting map consisted of over 300 markers in eight linkage groups, which could be aligned with a classical recombination map and the A. majus karyotype. The distribution of recombination frequencies and distorted transmission of parental alleles differed from those of a previous inter-species hybrid. The differences varied in magnitude and direction between chromosomes, suggesting that they had multiple causes. The map, which covered an estimated of 95% of the genome with an average interval of 2 cM, was used to analyze the distribution of a newly discovered family of MITE transposons and tested for its utility in positioning seven mutations that affect aspects of plant size.ConclusionsThe current map has an estimated interval of 1.28 Mb between markers. It shows a lower level of transmission ratio distortion and a longer length than the previous inter-species map, making it potentially more useful. The molecular recombination map further indicates that the IDLE MITE transposons are distributed throughout the genome and are relatively stable. The map proved effective in mapping classical morphological mutations of A. majus.