A novel Ca2+-binding protein that can rapidly transduce auxin responses during root growth.

Signaling cross talks between auxin, a regulator of plant development, and Ca2+, a universal second messenger, have been proposed to modulate developmental plasticity in plants. However, the underlying molecular mechanisms are largely unknown. Here, we report that in Arabidopsis roots, auxin elicits specific Ca2+ signaling patterns that spatially coincide with the expression pattern of auxin-regulated genes. We have identified the single EF-hand Ca2+-binding protein Ca2+-dependent modulator of ICR1 (CMI1) as an interactor of the Rho of plants (ROP) effector interactor of constitutively active ROP (ICR1). CMI1 expression is directly up-regulated by auxin, whereas the loss of function of CMI1 associates with the repression of auxin-induced Ca2+ increases in the lateral root cap and vasculature, indicating that CMI1 represses early auxin responses. In agreement, cmi1 mutants display an increased auxin response including shorter primary roots, longer root hairs, longer hypocotyls, and altered lateral root formation. Binding to ICR1 affects subcellular localization of CMI1 and its function. The interaction between CMI1 and ICR1 is Ca2+-dependent and involves a conserved hydrophobic pocket in CMI1 and calmodulin binding-like domain in ICR1. Remarkably, CMI1 is monomeric in solution and in vitro changes its secondary structure at cellular resting Ca2+ concentrations ranging between 10-9 and 10-8 M. Hence, CMI1 is a Ca2+-dependent transducer of auxin-regulated gene expression, which can function in a cell-specific fashion at steady-state as well as at elevated cellular Ca2+ levels to regulate auxin responses.

A two-step process for epigenetic inheritance in Arabidopsis.

In Arabidopsis, multisubunit RNA polymerases IV and V orchestrate RNA-directed DNA methylation (RdDM) and transcriptional silencing, but what identifies the loci to be silenced is unclear. We show that heritable silent locus identity at a specific subset of RdDM targets requires HISTONE DEACETYLASE 6 (HDA6) acting upstream of Pol IV recruitment and siRNA biogenesis. At these loci, epigenetic memory conferring silent locus identity is erased in hda6 mutants such that restoration of HDA6 activity cannot restore siRNA biogenesis or silencing. Silent locus identity is similarly lost in mutants for the cytosine maintenance methyltransferase, MET1. By contrast, pol IV or pol V mutants disrupt silencing without erasing silent locus identity, allowing restoration of Pol IV or Pol V function to restore silencing. Collectively, these observations indicate that silent locus specification and silencing are separable steps that together account for epigenetic inheritance of the silenced state.

Diversity and dynamics of the Drosophila transcriptome.

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Stilbenoid prenyltransferases define key steps in the diversification of peanut phytoalexins.

Defense responses of peanut (Arachis hypogaea) to biotic and abiotic stresses include the synthesis of prenylated stilbenoids. Members of this compound class show several protective activities in human disease studies, and the list of potential therapeutic targets continues to expand. Despite their medical and biological importance, the biosynthetic pathways of prenylated stilbenoids remain to be elucidated, and the genes encoding stilbenoid-specific prenyltransferases have yet to be identified in any plant species. In this study, we combined targeted transcriptomic and metabolomic analyses to discover prenyltransferase genes in elicitor-treated peanut hairy root cultures. Transcripts encoding five enzymes were identified, and two of these were functionally characterized in a transient expression system consisting of Agrobacterium-infiltrated leaves of Nicotiana benthamiana We observed that one of these prenyltransferases, AhR4DT-1, catalyzes a key reaction in the biosynthesis of prenylated stilbenoids, in which resveratrol is prenylated at its C-4 position to form arachidin-2, whereas another, AhR3'DT-1, added the prenyl group to C-3' of resveratrol. Each of these prenyltransferases was highly specific for stilbenoid substrates, and we confirmed their subcellular location in the plastid by fluorescence microscopy. Structural analysis of the prenylated stilbenoids suggested that these two prenyltransferase activities represent the first committed steps in the biosynthesis of a large number of prenylated stilbenoids and their derivatives in peanut. In summary, we have identified five candidate prenyltransferases in peanut and confirmed that two of them are stilbenoid-specific, advancing our understanding of this specialized enzyme family and shedding critical light onto the biosynthesis of bioactive stilbenoids.

The Zinc-Finger Protein SOP1 Is Required for a Subset of the Nuclear Exosome Functions in Arabidopsis.

Correct gene expression requires tight RNA quality control both at transcriptional and post-transcriptional levels. Using a splicing-defective allele of PASTICCINO2 (PAS2), a gene essential for plant development, we isolated suppressor mutations modifying pas2-1 mRNA profiles and restoring wild-type growth. Three suppressor of pas2 (sop) mutations modified the degradation of mis-spliced pas2-1 mRNA species, allowing the synthesis of a functional protein. Cloning of the suppressor mutations identified the core subunit of the exosome SOP2/RRP4, the exosome nucleoplasmic cofactor SOP3/HEN2 and a novel zinc-finger protein SOP1 that colocalizes with HEN2 in nucleoplasmic foci. The three SOP proteins counteract post-transcriptional (trans)gene silencing (PTGS), which suggests that they all act in RNA quality control. In addition, sop1 mutants accumulate some, but not all of the misprocessed mRNAs and other types of RNAs that are observed in exosome mutants. Taken together, our data show that SOP1 is a new component of nuclear RNA surveillance that is required for the degradation of a specific subset of nuclear exosome targets.

A Specialized Diacylglycerol Acyltransferase Contributes to the Extreme Medium-Chain Fatty Acid Content of Cuphea Seed Oil.

Seed oils of many Cuphea sp. contain >90% of medium-chain fatty acids, such as decanoic acid (10:0). These seed oils, which are among the most compositionally variant in the plant kingdom, arise from specialized fatty acid biosynthetic enzymes and specialized acyltransferases. These include lysophosphatidic acid acyltransferases (LPAT) and diacylglycerol acyltransferases (DGAT) that are required for successive acylation of medium-chain fatty acids in the sn-2 and sn-3 positions of seed triacylglycerols (TAGs). Here we report the identification of a cDNA for a DGAT1-type enzyme, designated CpuDGAT1, from the transcriptome of C. avigera var pulcherrima developing seeds. Microsomes of camelina (Camelina sativa) seeds engineered for CpuDGAT1 expression displayed DGAT activity with 10:0-CoA and the diacylglycerol didecanoyl, that was approximately 4-fold higher than that in camelina seed microsomes lacking CpuDGAT1. In addition, coexpression in camelina seeds of CpuDGAT1 with a C. viscosissima FatB thioesterase (CvFatB1) that generates 10:0 resulted in TAGs with nearly 15 mol % of 10:0. More strikingly, expression of CpuDGAT1 and CvFatB1 with the previously described CvLPAT2, a 10:0-CoA-specific Cuphea LPAT, increased 10:0 amounts to 25 mol % in camelina seed TAG. These TAGs contained up to 40 mol % 10:0 in the sn-2 position, nearly double the amounts obtained from coexpression of CvFatB1 and CvLPAT2 alone. Although enriched in diacylglycerol, 10:0 was not detected in phosphatidylcholine in these seeds. These findings are consistent with channeling of 10:0 into TAG through the combined activities of specialized LPAT and DGAT activities and demonstrate the biotechnological use of these enzymes to generate 10:0-rich seed oils.

A Stilbenoid-Specific Prenyltransferase Utilizes Dimethylallyl Pyrophosphate from the Plastidic Terpenoid Pathway.

Prenylated stilbenoids synthesized in some legumes exhibit plant pathogen defense properties and pharmacological activities with potential benefits to human health. Despite their importance, the biosynthetic pathways of these compounds remain to be elucidated. Peanut (Arachis hypogaea) hairy root cultures produce a diverse array of prenylated stilbenoids upon treatment with elicitors. Using metabolic inhibitors of the plastidic and cytosolic isoprenoid biosynthetic pathways, we demonstrated that the prenyl moiety on the prenylated stilbenoids derives from a plastidic pathway. We further characterized, to our knowledge for the first time, a membrane-bound stilbenoid-specific prenyltransferase activity from the microsomal fraction of peanut hairy roots. This microsomal fraction-derived resveratrol 4-dimethylallyl transferase utilizes 3,3-dimethylallyl pyrophosphate as a prenyl donor and prenylates resveratrol to form arachidin-2. It also prenylates pinosylvin to chiricanine A and piceatannol to arachidin-5, a prenylated stilbenoid identified, to our knowledge, for the first time in this study. This prenyltransferase exhibits strict substrate specificity for stilbenoids and does not prenylate flavanone, flavone, or isoflavone backbones, even though it shares several common features with flavonoid-specific prenyltransferases.

Structurally divergent lysophosphatidic acid acyltransferases with high selectivity for saturated medium chain fatty acids from Cuphea seeds.

Lysophosphatidic acid acyltransferase (LPAT) catalyzes acylation of the sn-2 position on lysophosphatidic acid by an acyl CoA substrate to produce the phosphatidic acid precursor of polar glycerolipids and triacylglycerols (TAGs). In the case of TAGs, this reaction is typically catalyzed by an LPAT2 from microsomal LPAT class A that has high specificity for C18 fatty acids containing Δ9 unsaturation. Because of this specificity, the occurrence of saturated fatty acids in the TAG sn-2 position is infrequent in seed oils. To identify LPATs with variant substrate specificities, deep transcriptomic mining was performed on seeds of two Cuphea species producing TAGs that are highly enriched in saturated C8 and C10 fatty acids. From these analyses, cDNAs for seven previously unreported LPATs were identified, including cDNAs from Cuphea viscosissima (CvLPAT2) and Cuphea avigera var. pulcherrima (CpuLPAT2a) encoding microsomal, seed-specific class A LPAT2s and a cDNA from C. avigera var. pulcherrima (CpuLPATB) encoding a microsomal, seed-specific LPAT from the bacterial-type class B. The activities of these enzymes were characterized in Camelina sativa by seed-specific co-expression with cDNAs for various Cuphea FatB acyl-acyl carrier protein thioesterases (FatB) that produce a variety of saturated medium-chain fatty acids. CvLPAT2 and CpuLPAT2a expression resulted in accumulation of 10:0 fatty acids in the Camelina sativa TAG sn-2 position, indicating a 10:0 CoA specificity that has not been previously described for plant LPATs. CpuLPATB expression generated TAGs with 14:0 at the sn-2 position, but not 10:0. Identification of these LPATs provides tools for understanding the structural basis of LPAT substrate specificity and for generating altered oil functionalities.

FEATnotator: A tool for integrated annotation of sequence features and variation, facilitating interpretation in genomics experiments.

As approaches are sought for more efficient and democratized uses of non-model and expanded model genomics references, ease of integration of genomic feature datasets is especially desirable in multidisciplinary research communities. Valuable conclusions are often missed or slowed when researchers refer experimental results to a single reference sequence that lacks integrated pan-genomic and multi-experiment data in accessible formats. Association of genomic positional information, such as results from an expansive variety of next-generation sequencing experiments, with annotated reference features such as genes or predicted protein binding sites, provides the context essential for conclusions and ongoing research. When the experimental system includes polymorphic genomic inputs, rapid calculation of gene structural and protein translational effects of sequence variation from the reference can be invaluable. Here we present FEATnotator, a lightweight, fast and easy to use open source software program that integrates and reports overlap and proximity in genomic information from any user-defined datasets including those from next generation sequencing applications. We illustrate use of the tool by summarizing whole genome sequence variation of a widely used natural isolate of Arabidopsis thaliana in the context of gene models of the reference accession. Previous discovery of a protein coding deletion influencing root development is replicated rapidly. Appropriate even in investigations of a single gene or genic regions such as QTL, comprehensive reports provided by FEATnotator better prepare researchers for interpretation of their experimental results. The tool is available for download at http://featnotator.sourceforge.net.

Evolutionary divergence of core and post-translational circadian clock genes in the pitcher-plant mosquito, Wyeomyia smithii.

BACKGROUND: Internal circadian (circa, about; dies, day) clocks enable organisms to maintain adaptive timing of their daily behavioral activities and physiological functions. Eukaryotic clocks consist of core transcription-translation feedback loops that generate a cycle and post-translational modifiers that maintain that cycle at about 24 h. We use the pitcher-plant mosquito, Wyeomyia smithii (subfamily Culicini, tribe Sabethini), to test whether evolutionary divergence of the circadian clock genes in this species, relative to other insects, has involved primarily genes in the core feedback loops or the post-translational modifiers. Heretofore, there is no reference transcriptome or genome sequence for any mosquito in the tribe Sabethini, which includes over 375 mainly circumtropical species. METHODS: We sequenced, assembled and annotated the transcriptome of W. smithii containing nearly 95 % of conserved single-copy orthologs in animal genomes. We used the translated contigs and singletons to determine the average rates of circadian clock-gene divergence in W. smithii relative to three other mosquito genera, to Drosophila, to the butterfly, Danaus, and to the wasp, Nasonia. RESULTS: Over 1.08 million cDNA sequence reads were obtained consisting of 432.5 million nucleotides. Their assembly produced 25,904 contigs and 54,418 singletons of which 62 % and 28 % are annotated as protein-coding genes, respectively, sharing homology with other animal proteomes. DISCUSSION: The W. smithii transcriptome includes all nine circadian transcription-translation feedback-loop genes and all eight post-translational modifier genes we sought to identify (Fig. 1). After aligning translated W. smithii contigs and singletons from this transcriptome with other insects, we determined that there was no significant difference in the average divergence of W. smithii from the six other taxa between the core feedback-loop genes and post-translational modifiers. CONCLUSIONS: The characterized transcriptome is sufficiently complete and of sufficient quality to have uncovered all of the insect circadian clock genes we sought to identify (Fig. 1). Relative divergence does not differ between core feedback-loop genes and post-translational modifiers of those genes in a Sabethine species (W. smithii) that has experienced a continual northward dispersal into temperate regions of progressively longer summer day lengths as compared with six other insect taxa. An associated microarray platform derived from this work will enable the investigation of functional genomics of circadian rhythmicity, photoperiodic time measurement, and diapause along a photic and seasonal geographic gradient.

Deep sequencing of the Mexican avocado transcriptome, an ancient angiosperm with a high content of fatty acids.

BACKGROUND: Avocado (Persea americana) is an economically important tropical fruit considered to be a good source of fatty acids. Despite its importance, the molecular and cellular characterization of biochemical and developmental processes in avocado is limited due to the lack of transcriptome and genomic information. RESULTS: The transcriptomes of seeds, roots, stems, leaves, aerial buds and flowers were determined using different sequencing platforms. Additionally, the transcriptomes of three different stages of fruit ripening (pre-climacteric, climacteric and post-climacteric) were also analyzed. The analysis of the RNAseqatlas presented here reveals strong differences in gene expression patterns between different organs, especially between root and flower, but also reveals similarities among the gene expression patterns in other organs, such as stem, leaves and aerial buds (vegetative organs) or seed and fruit (storage organs). Important regulators, functional categories, and differentially expressed genes involved in avocado fruit ripening were identified. Additionally, to demonstrate the utility of the avocado gene expression atlas, we investigated the expression patterns of genes implicated in fatty acid metabolism and fruit ripening. CONCLUSIONS: A description of transcriptomic changes occurring during fruit ripening was obtained in Mexican avocado, contributing to a dynamic view of the expression patterns of genes involved in fatty acid biosynthesis and the fruit ripening process.

Oil biosynthesis in a basal angiosperm: transcriptome analysis of Persea Americana mesocarp.

BACKGROUND: The mechanism by which plants synthesize and store high amounts of triacylglycerols (TAG) in tissues other than seeds is not well understood. The comprehension of controls for carbon partitioning and oil accumulation in nonseed tissues is essential to generate oil-rich biomass in perennial bioenergy crops. Persea americana (avocado), a basal angiosperm with unique features that are ancestral to most flowering plants, stores ~ 70 % TAG per dry weight in its mesocarp, a nonseed tissue. Transcriptome analyses of select pathways, from generation of pyruvate and leading up to TAG accumulation, in mesocarp tissues of avocado was conducted and compared with that of oil-rich monocot (oil palm) and dicot (rapeseed and castor) tissues to identify tissue- and species-specific regulation and biosynthesis of TAG in plants. RESULTS: RNA-Seq analyses of select lipid metabolic pathways of avocado mesocarp revealed patterns similar to that of other oil-rich species. However, only some predominant orthologs of the fatty acid biosynthetic pathway genes in this basal angiosperm were similar to those of monocots and dicots. The accumulation of TAG, rich in oleic acid, was associated with higher transcript levels for a putative stearoyl-ACP desaturase and endoplasmic reticulum (ER)-associated acyl-CoA synthetases, during fruit development. Gene expression levels for enzymes involved in terminal steps to TAG biosynthesis in the ER further indicated that both acyl-CoA-dependent and -independent mechanisms might play a role in TAG assembly, depending on the developmental stage of the fruit. Furthermore, in addition to the expression of an ortholog of WRINKLED1 (WRI1), a regulator of fatty acid biosynthesis, high transcript levels for WRI2-like and WRI3-like suggest a role for additional transcription factors in nonseed oil accumulation. Plastid pyruvate necessary for fatty acid synthesis is likely driven by the upregulation of genes involved in glycolysis and transport of its intermediates. Together, a comparative transcriptome analyses for storage oil biosynthesis in diverse plants and tissues suggested that several distinct and conserved features in this basal angiosperm species might contribute towards its rich TAG content. CONCLUSIONS: Our work represents a comprehensive transcriptome resource for a basal angiosperm species and provides insight into their lipid metabolism in mesocarp tissues. Furthermore, comparison of the transcriptome of oil-rich mesocarp of avocado, with oil-rich seed and nonseed tissues of monocot and dicot species, revealed lipid gene orthologs that are highly conserved during evolution. The orthologs that are distinctively expressed in oil-rich mesocarp tissues of this basal angiosperm, such as WRI2, ER-associated acyl-CoA synthetases, and lipid-droplet associated proteins were also identified. This study provides a foundation for future investigations to increase oil-content and has implications for metabolic engineering to enhance storage oil content in nonseed tissues of diverse species.

Toward production of jet fuel functionality in oilseeds: identification of FatB acyl-acyl carrier protein thioesterases and evaluation of combinatorial expression strategies in Camelina seeds.

Seeds of members of the genus Cuphea accumulate medium-chain fatty acids (MCFAs; 8:0-14:0). MCFA- and palmitic acid- (16:0) rich vegetable oils have received attention for jet fuel production, given their similarity in chain length to Jet A fuel hydrocarbons. Studies were conducted to test genes, including those from Cuphea, for their ability to confer jet fuel-type fatty acid accumulation in seed oil of the emerging biofuel crop Camelina sativa. Transcriptomes from Cuphea viscosissima and Cuphea pulcherrima developing seeds that accumulate >90% of C8 and C10 fatty acids revealed three FatB cDNAs (CpuFatB3, CvFatB1, and CpuFatB4) expressed predominantly in seeds and structurally divergent from typical FatB thioesterases that release 16:0 from acyl carrier protein (ACP). Expression of CpuFatB3 and CvFatB1 resulted in Camelina oil with capric acid (10:0), and CpuFatB4 expression conferred myristic acid (14:0) production and increased 16:0. Co-expression of combinations of previously characterized Cuphea and California bay FatBs produced Camelina oils with mixtures of C8-C16 fatty acids, but amounts of each fatty acid were less than obtained by expression of individual FatB cDNAs. Increases in lauric acid (12:0) and 14:0, but not 10:0, in Camelina oil and at the sn-2 position of triacylglycerols resulted from inclusion of a coconut lysophosphatidic acid acyltransferase specialized for MCFAs. RNA interference (RNAi) suppression of Camelina ß-ketoacyl-ACP synthase II, however, reduced 12:0 in seeds expressing a 12:0-ACP-specific FatB. Camelina lines presented here provide platforms for additional metabolic engineering targeting fatty acid synthase and specialized acyltransferases for achieving oils with high levels of jet fuel-type fatty acids.

Camelina seed transcriptome: a tool for meal and oil improvement and translational research.

Camelina (Camelina sativa), a Brassicaceae oilseed, has received recent interest as a biofuel crop and production platform for industrial oils. Limiting wider production of camelina for these uses is the need to improve the quality and content of the seed protein-rich meal and oil, which is enriched in oxidatively unstable polyunsaturated fatty acids that are deleterious for biodiesel. To identify candidate genes for meal and oil quality improvement, a transcriptome reference was built from 2047 Sanger ESTs and more than 2 million 454-derived sequence reads, representing genes expressed in developing camelina seeds. The transcriptome of approximately 60K transcripts from 22 597 putative genes includes camelina homologues of nearly all known seed-expressed genes, suggesting a high level of completeness and usefulness of the reference. These sequences included candidates for 12S (cruciferins) and 2S (napins) seed storage proteins (SSPs) and nearly all known lipid genes, which have been compiled into an accessible database. To demonstrate the utility of the transcriptome for seed quality modification, seed-specific RNAi lines deficient in napins were generated by targeting 2S SSP genes, and high oleic acid oil lines were obtained by targeting FATTY ACID DESATURASE 2 (FAD2) and FATTY ACID ELONGASE 1 (FAE1). The high sequence identity between Arabidopsis thaliana and camelina genes was also exploited to engineer high oleic lines by RNAi with Arabidopsis FAD2 and FAE1 sequences. It is expected that these transcriptomic data will be useful for breeding and engineering of additional camelina seed traits and for translating findings from the model Arabidopsis to an oilseed crop.

Identifying viral integration sites using SeqMap 2.0.

UNLABELLED: Retroviral integration has been implicated in several biomedical applications, including identification of cancer-associated genes and malignant transformation in gene therapy clinical trials. We introduce an efficient and scalable method for fast identification of viral vector integration sites from long read high-throughput sequencing. Individual sequence reads are masked to remove non-genomic sequence, aligned to the host genome and assembled into contiguous fragments used to pinpoint the position of integration. AVAILABILITY AND IMPLEMENTATION: The method is implemented in a publicly accessible web server platform, SeqMap 2.0, containing analysis tools and both private and shared lab workspaces that facilitate collaboration among researchers. Available at http://seqmap.compbio.iupui.edu/.

Sequencing of a QTL-rich region of the Theobroma cacao genome using pooled BACs and the identification of trait specific candidate genes.

BACKGROUND: BAC-based physical maps provide for sequencing across an entire genome or a selected sub-genomic region of biological interest. Such a region can be approached with next-generation whole-genome sequencing and assembly as if it were an independent small genome. Using the minimum tiling path as a guide, specific BAC clones representing the prioritized genomic interval are selected, pooled, and used to prepare a sequencing library. RESULTS: This pooled BAC approach was taken to sequence and assemble a QTL-rich region, of ~3 Mbp and represented by twenty-seven BACs, on linkage group 5 of the Theobroma cacao cv. Matina 1-6 genome. Using various mixtures of read coverages from paired-end and linear 454 libraries, multiple assemblies of varied quality were generated. Quality was assessed by comparing the assembly of 454 reads with a subset of ten BACs individually sequenced and assembled using Sanger reads. A mixture of reads optimal for assembly was identified. We found, furthermore, that a quality assembly suitable for serving as a reference genome template could be obtained even with a reduced depth of sequencing coverage. Annotation of the resulting assembly revealed several genes potentially responsible for three T. cacao traits: black pod disease resistance, bean shape index, and pod weight. CONCLUSIONS: Our results, as with other pooled BAC sequencing reports, suggest that pooling portions of a minimum tiling path derived from a BAC-based physical map is an effective method to target sub-genomic regions for sequencing. While we focused on a single QTL region, other QTL regions of importance could be similarly sequenced allowing for biological discovery to take place before a high quality whole-genome assembly is completed.

A garter snake transcriptome: pyrosequencing, de novo assembly, and sex-specific differences.

BACKGROUND: The reptiles, characterized by both diversity and unique evolutionary adaptations, provide a comprehensive system for comparative studies of metabolism, physiology, and development. However, molecular resources for ectothermic reptiles are severely limited, hampering our ability to study the genetic basis for many evolutionarily important traits such as metabolic plasticity, extreme longevity, limblessness, venom, and freeze tolerance. Here we use massively parallel sequencing (454 GS-FLX Titanium) to generate a transcriptome of the western terrestrial garter snake (Thamnophis elegans) with two goals in mind. First, we develop a molecular resource for an ectothermic reptile; and second, we use these sex-specific transcriptomes to identify differences in the presence of expressed transcripts and potential genes of evolutionary interest. RESULTS: Using sex-specific pools of RNA (one pool for females, one pool for males) representing 7 tissue types and 35 diverse individuals, we produced 1.24 million sequence reads, which averaged 366 bp in length after cleaning. Assembly of the cleaned reads from both sexes with NEWBLER and MIRA resulted in 96,379 contigs containing 87% of the cleaned reads. Over 34% of these contigs and 13% of the singletons were annotated based on homology to previously identified proteins. From these homology assignments, additional clustering, and ORF predictions, we estimate that this transcriptome contains ~13,000 unique genes that were previously identified in other species and over 66,000 transcripts from unidentified protein-coding genes. Furthermore, we use a graph-clustering method to identify contigs linked by NEWBLER-split reads that represent divergent alleles, gene duplications, and alternatively spliced transcripts. Beyond gene identification, we identified 95,295 SNPs and 31,651 INDELs. From these sex-specific transcriptomes, we identified 190 genes that were only present in the mRNA sequenced from one of the sexes (84 female-specific, 106 male-specific), and many highly variable genes of evolutionary interest. CONCLUSIONS: This is the first large-scale, multi-organ transcriptome for an ectothermic reptile. This resource provides the most comprehensive set of EST sequences available for an individual ectothermic reptile species, increasing the number of snake ESTs 50-fold. We have identified genes that appear to be under evolutionary selection and those that are sex-specific. This resource will assist studies on gene expression and comparative genomics, and will facilitate the study of evolutionarily important traits at the molecular level.

Corrigendum: A Larger Chocolate Chip-Development of a 15K Theobroma cacao L. SNP Array to Create High-Density Linkage Maps.

[This corrects the article on p. 2008 in vol. 8, PMID: 29259608.].

Population genomic analyses of the chocolate tree, Theobroma cacao L., provide insights into its domestication process.

Domestication has had a strong impact on the development of modern societies. We sequenced 200 genomes of the chocolate plant Theobroma cacao L. to show for the first time to our knowledge that a single population, the Criollo population, underwent strong domestication ~3600 years ago (95% CI: 2481-13,806 years ago). We also show that during the process of domestication, there was strong selection for genes involved in the metabolism of the colored protectants anthocyanins and the stimulant theobromine, as well as disease resistance genes. Our analyses show that domesticated populations of T. cacao (Criollo) maintain a higher proportion of high-frequency deleterious mutations. We also show for the first time the negative consequences of the increased accumulation of deleterious mutations during domestication on the fitness of individuals (significant reduction in kilograms of beans per hectare per year as Criollo ancestry increases, as estimated from a GLM, P = 0.000425).