Todd P. Michael

The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus)

Papaya, a fruit crop cultivated in tropical and subtropical regions, is known for its nutritional benefits and medicinal applications. Here we report a 3× draft genome sequence of ‘SunUp’ papaya, the first commercial virus-resistant transgenic fruit tree to be sequenced. The papaya genome is three times the size of the Arabidopsis genome, but contains fewer genes, including significantly fewer disease-resistance gene analogues. Comparison of the five sequenced genomes suggests a minimal angiosperm gene set of 13,311. A lack of recent genome duplication, atypical of other angiosperm genomes sequenced so far, may account for the smaller papaya gene number in most functional groups. Nonetheless, striking amplifications in gene number within particular functional groups suggest roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical daylengths. Transgenesis at three locations is closely associated with chloroplast insertions into the nuclear genome, and with topoisomerase I recognition sites. Papaya offers numerous advantages as a system for fruit-tree functional genomics, and this draft genome sequence provides the foundation for revealing the basis of Carica’s distinguishing morpho-physiological, medicinal and nutritional properties.

The genome of woodland strawberry ( Fragaria vesca )

The woodland strawberry, Fragaria vesca (2n = 2x = 14), is a versatile experimental plant system. This diminutive herbaceous perennial has a small genome (240 Mb), is amenable to genetic transformation and shares substantial sequence identity with the cultivated strawberry (Fragaria × ananassa) and other economically important rosaceous plants. Here we report the draft F. vesca genome, which was sequenced to ×39 coverage using second-generation technology, assembled de novo and then anchored to the genetic linkage map into seven pseudochromosomes. This diploid strawberry sequence lacks the large genome duplications seen in other rosids. Gene prediction modeling identified 34,809 genes, with most being supported by transcriptome mapping. Genes critical to valuable horticultural traits including flavor, nutritional value and flowering time were identified. Macrosyntenic relationships between Fragaria and Prunus predict a hypothetical ancestral Rosaceae genome that had nine chromosomes. New phylogenetic analysis of 154 protein-coding genes suggests that assignment of Populus to Malvidae, rather than Fabidae, is warranted.

Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.)

BackgroundSacred lotus is a basal eudicot with agricultural, medicinal, cultural and religious importance. It was domesticated in Asia about 7,000 years ago, and cultivated for its rhizomes and seeds as a food crop. It is particularly noted for its 1,300-year seed longevity and exceptional water repellency, known as the lotus effect. The latter property is due to the nanoscopic closely packed protuberances of its self-cleaning leaf surface, which have been adapted for the manufacture of a self-cleaning industrial paint, Lotusan.ResultsThe genome of the China Antique variety of the sacred lotus was sequenced with Illumina and 454 technologies, at respective depths of 101× and 5.2×. The final assembly has a contig N50 of 38.8 kbp and a scaffold N50 of 3.4 Mbp, and covers 86.5% of the estimated 929 Mbp total genome size. The genome notably lacks the paleo-triplication observed in other eudicots, but reveals a lineage-specific duplication. The genome has evidence of slow evolution, with a 30% slower nucleotide mutation rate than observed in grape. Comparisons of the available sequenced genomes suggest a minimum gene set for vascular plants of 4,223 genes. Strikingly, the sacred lotus has 16 COG2132 multi-copper oxidase family proteins with root-specific expression; these are involved in root meristem phosphate starvation, reflecting adaptation to limited nutrient availability in an aquatic environment.ConclusionsThe slow nucleotide substitution rate makes the sacred lotus a better resource than the current standard, grape, for reconstructing the pan-eudicot genome, and should therefore accelerate comparative analysis between eudicots and monocots.

SOPRA: Scaffolding algorithm for paired reads via statistical optimization

BackgroundHigh throughput sequencing (HTS) platforms produce gigabases of short read (<100 bp) data per run. While these short reads are adequate for resequencing applications, de novo assembly of moderate size genomes from such reads remains a significant challenge. These limitations could be partially overcome by utilizing mate pair technology, which provides pairs of short reads separated by a known distance along the genome.ResultsWe have developed SOPRA, a tool designed to exploit the mate pair/paired-end information for assembly of short reads. The main focus of the algorithm is selecting a sufficiently large subset of simultaneously satisfiable mate pair constraints to achieve a balance between the size and the quality of the output scaffolds. Scaffold assembly is presented as an optimization problem for variables associated with vertices and with edges of the contig connectivity graph. Vertices of this graph are individual contigs with edges drawn between contigs connected by mate pairs. Similar graph problems have been invoked in the context of shotgun sequencing and scaffold building for previous generation of sequencing projects. However, given the error-prone nature of HTS data and the fundamental limitations from the shortness of the reads, the ad hoc greedy algorithms used in the earlier studies are likely to lead to poor quality results in the current context. SOPRA circumvents this problem by treating all the constraints on equal footing for solving the optimization problem, the solution itself indicating the problematic constraints (chimeric/repetitive contigs, etc.) to be removed. The process of solving and removing of constraints is iterated till one reaches a core set of consistent constraints. For SOLiD sequencer data, SOPRA uses a dynamic programming approach to robustly translate the color-space assembly to base-space. For assessing the quality of an assembly, we report the no-match/mismatch error rate as well as the rates of various rearrangement errors.ConclusionsApplying SOPRA to real data from bacterial genomes, we were able to assemble contigs into scaffolds of significant length (N50 up to 200 Kb) with very few errors introduced in the process. In general, the methodology presented here will allow better scaffold assemblies of any type of mate pair sequencing data.

Integrative analysis of chromatin states in Arabidopsis identified potential regulatory mechanisms for natural antisense transcript production

Genome-wide analyses of epigenomic and transcriptomic profiles provide extensive resources for discovering epigenetic regulatory mechanisms. However, the construction of functionally relevant hypotheses from correlative patterns and the rigorous testing of these hypotheses may be challenging. We combined bioinformatics-driven hypothesis building with mutant analyses to identify potential epigenetic mechanisms using the model plant Arabidopsis thaliana. Genome-wide maps of nine histone modifications produced by ChIP-seq were used together with a strand-specific RNA-seq dataset to profile the epigenome and transcriptome of Arabidopsis. Combinatorial chromatin patterns were described by 42 major chromatin states with selected states validated using the re-ChIP assay. The functional relevance of chromatin modifications was analyzed using the ANchored CORrelative Pattern (ANCORP) method and a newly developed state-specific effects analysis (SSEA) method, which interrogates individual chromatin marks in the context of combinatorial chromatin states. Based on results from these approaches, we propose the hypothesis that cytosine methylation (5mC) and histone methylation H3K36me may synergistically repress production of natural antisense transcripts (NATs) in the context of actively expressed genes. Mutant analyses supported this proposed model at a significant proportion of the tested loci. We further identified polymerase-associated factor as a potential repressor for NAT abundance. Although the majority of tested NATs were found to localize to the nucleus, we also found evidence for cytoplasmically partitioned NATs. The significance of the subcellular localization of NATs and their biological functions remain to be defined.

Global Profiling of Rice and Poplar Transcriptomes Highlights Key Conserved Circadian-Controlled Pathways and cis-Regulatory Modules

BACKGROUND Circadian clocks provide an adaptive advantage through anticipation of daily and seasonal environmental changes. In plants, the central clock oscillator is regulated by several interlocking feedback loops. It was shown that a substantial proportion of the Arabidopsis genome cycles with phases of peak expression covering the entire day. Synchronized transcriptome cycling is driven through an extensive network of diurnal and clock-regulated transcription factors and their target cis-regulatory elements. Study of the cycling transcriptome in other plant species could thus help elucidate the similarities and differences and identify hubs of regulation common to monocot and dicot plants. METHODOLOGY/PRINCIPAL FINDINGS Using a combination of oligonucleotide microarrays and data mining pipelines, we examined daily rhythms in gene expression in one monocotyledonous and one dicotyledonous plant, rice and poplar, respectively. Cycling transcriptomes were interrogated under different diurnal (driven) and circadian (free running) light and temperature conditions. Collectively, photocycles and thermocycles regulated about 60% of the expressed nuclear genes in rice and poplar. Depending on the condition tested, up to one third of oscillating Arabidopsis-poplar-rice orthologs were phased within three hours of each other suggesting a high degree of conservation in terms of rhythmic gene expression. We identified clusters of rhythmically co-expressed genes and searched their promoter sequences to identify phase-specific cis-elements, including elements that were conserved in the promoters of Arabidopsis, poplar, and rice. CONCLUSIONS/SIGNIFICANCE Our results show that the cycling patterns of many circadian clock genes are highly conserved across poplar, rice, and Arabidopsis. The expression of many orthologous genes in key metabolic and regulatory pathways is diurnal and/or circadian regulated and phased to similar times of day. Our results confirm previous findings in Arabidopsis of three major classes of cis-regulatory modules within the plant circadian network: the morning (ME, GBOX), evening (EE, GATA), and midnight (PBX/TBX/SBX) modules. Identification of identical overrepresented motifs in the promoters of cycling genes from different species suggests that the core diurnal/circadian cis-regulatory network is deeply conserved between mono- and dicotyledonous species.

Cytochrome P450 Monooxygenases as Reporters for Circadian-Regulated Pathways

Cytochrome P450 monooxygenases (P450s) play important roles in the synthesis of diverse secondary compounds in Arabidopsis (Arabidopsis thaliana). Comparison of four data sets analyzing seedlings harvested over a 2-d period of constant conditions after growth with varying photoperiods and thermocycles recorded a total of 98 P450 loci as circadian regulated for at least one of the four conditions. Here, we further describe the circadian-regulated pathways using, as reporters, individual P450 loci that are likely to be rate limiting in secondary metabolic pathways. Reverse transcription-polymerase chain reaction gel blot analyses have confirmed circadian regulation of P450s in phenylpropanoid, carotenoid, oxylipin, glucosinolate, and brassinosteroid biosyntheses and have shown that both P450 and non-P450 genes in the many branches of the phenylpropanoid pathway have similar circadian patterns of expression. In silico analyses of the subsets of coregulated promoters have identified overrepresented promoter elements in various biosynthetic pathway genes, including MYB and MYB4 elements that are significantly more abundant in promoters for the core and lignin sections of phenylpropanoid metabolism. Interactions with these elements important for circadian regulation do not involve the MYB transcription factor PAP1, as previously proposed, since the expression patterns of circadian-regulated P450s are the same in pap1-D mutant seedlings as in wild-type seedlings. Further analysis of circadian-regulated promoters in other biochemical pathways provides us with the opportunity to identify novel promoter motifs that might be important in P450 circadian regulation.

Analysis of global gene expression in Brachypodium distachyon reveals extensive network plasticity in response to abiotic stress.

Brachypodium distachyon is a close relative of many important cereal crops. Abiotic stress tolerance has a significant impact on productivity of agriculturally important food and feedstock crops. Analysis of the transcriptome of Brachypodium after chilling, high-salinity, drought, and heat stresses revealed diverse differential expression of many transcripts. Weighted Gene Co-Expression Network Analysis revealed 22 distinct gene modules with specific profiles of expression under each stress. Promoter analysis implicated short DNA sequences directly upstream of module members in the regulation of 21 of 22 modules. Functional analysis of module members revealed enrichment in functional terms for 10 of 22 network modules. Analysis of condition-specific correlations between differentially expressed gene pairs revealed extensive plasticity in the expression relationships of gene pairs. Photosynthesis, cell cycle, and cell wall expression modules were down-regulated by all abiotic stresses. Modules which were up-regulated by each abiotic stress fell into diverse and unique gene ontology GO categories. This study provides genomics resources and improves our understanding of abiotic stress responses of Brachypodium.

Isolation and analysis of high quality nuclear DNA with reduced organellar DNA for plant genome sequencing and resequencing.

BackgroundHigh throughput sequencing (HTS) technologies have revolutionized the field of genomics by drastically reducing the cost of sequencing, making it feasible for individual labs to sequence or resequence plant genomes. Obtaining high quality, high molecular weight DNA from plants poses significant challenges due to the high copy number of chloroplast and mitochondrial DNA, as well as high levels of phenolic compounds and polysaccharides. Multiple methods have been used to isolate DNA from plants; the CTAB method is commonly used to isolate total cellular DNA from plants that contain nuclear DNA, as well as chloroplast and mitochondrial DNA. Alternatively, DNA can be isolated from nuclei to minimize chloroplast and mitochondrial DNA contamination.ResultsWe describe optimized protocols for isolation of nuclear DNA from eight different plant species encompassing both monocot and eudicot species. These protocols use nuclei isolation to minimize chloroplast and mitochondrial DNA contamination. We also developed a protocol to determine the number of chloroplast and mitochondrial DNA copies relative to the nuclear DNA using quantitative real time PCR (qPCR). We compared DNA isolated from nuclei to total cellular DNA isolated with the CTAB method. As expected, DNA isolated from nuclei consistently yielded nuclear DNA with fewer chloroplast and mitochondrial DNA copies, as compared to the total cellular DNA prepared with the CTAB method. This protocol will allow for analysis of the quality and quantity of nuclear DNA before starting a plant whole genome sequencing or resequencing experiment.ConclusionsExtracting high quality, high molecular weight nuclear DNA in plants has the potential to be a bottleneck in the era of whole genome sequencing and resequencing. The methods that are described here provide a framework for researchers to extract and quantify nuclear DNA in multiple types of plants.

The carnivorous bladderwort (Utricularia, Lentibulariaceae): a system inflates

Carnivorous plants inhabit nutrient-poor environments, where prominent targets of prey capture are organic nitrogen and phosphorus. Some carnivorous plants also acquire carbon from their victims. A new report focusing on Utricularia, the bladderwort, demonstrates that carbon assimilated from photosynthesis is paradoxically secreted into the trapping environment, where it may help to support a mutualistic bacterial community. This bacterial community may also secrete allelochemicals that attract microcrustaceans which bear a strong overt resemblance to bladderwort traps. Furthermore, Utricularia and its sister genus Genlisea share anomalous molecular evolutionary features, such as highly increased rates of nucleotide substitution and dynamic evolution of genome size, from approximately 60-1500 megabases depending on the species or even population. A mechanistic hypothesis, based on the mutagenic action of reactive oxygen species (ROS) is proposed to underlie these phenomena, involving error-prone repair at the level of DNA bases and double-strand breaks. It is argued that these plants are prime candidates for further research on the complexities of plant physiology associated with carnivory, metagenomic surveys of trap microbial communities, novel plant nitrogen/nutrient utilization pathways, the ecology of prey attraction, whole-plant and trap comparative development, and, finally, evolution of the minimal angiosperm genome.