Transcript profiling for early stages during embryo development in Scots pine

Background: Characterization of the expression and function of genes regulating embryo development in conifers is interesting from an evolutionary point of view. However, our knowledge about the regulation of embryo development in conifers is limited. During early embryo development in Pinus species the proembyo goes through a cleavage process, named cleavage polyembryony, giving rise to four embryos. One of these embryos develops to a dominant embryo, which will develop further into a mature, cotyledonary embryo, while the other embryos, the subordinate embryos, are degraded. The main goal of this study has been to identify processes that might be important for regulating the cleavage process and for the development of a dominant embryo. Results: RNA samples from embryos and megagametophytes at four early developmental stages during seed development in Pinus sylvestris were subjected to high-throughput sequencing. A total of 6.6 million raw reads was generated, resulting in 121,938 transcripts, out of which 36.106 contained ORFs. 18,638 transcripts were differentially expressed (DETs) in embryos and megagametophytes. GO enrichment analysis of transcripts up-regulated in embryos showed enrichment for different cellular processes, while those up-regulated in megagametophytes were enriched for accumulation of storage material and responses to stress. The highest number of DETs was detected during the initiation of the cleavage process. Transcripts related to embryogenic competence, cell wall modifications, cell division pattern, axis specification and response to hormones and stress were highly abundant and differentially expressed during early embryo development. The abundance of representative DETs was confirmed by qRT-PCR analyses. Conclusion: Based on the processes identified in the GO enrichment analyses and the expression of the selected transcripts we suggest that (i) processes related to embryogenic competence and cell wall loosening are involved in activating the cleavage process; (ii) apical-basal polarization is strictly regulated in dominant embryos but not in the subordinate embryos; (iii) the transition from the morphogenic phase to the maturation phase is not completed in subordinate embryos. This is the first genome-wide transcript expression profiling of the earliest stages during embryo development in a Pinus species. Our results can serve as a framework for future studies to reveal the functions of identified genes

Evidence for an ancient whole genome duplication in the cycad lineage

Contrary to the many whole genome duplication events recorded for angiosperms (flowering plants), whole genome duplications in gymnosperms (non-flowering seed plants) seem to be much rarer. Although ancient whole genome duplications have been reported for most gymnosperm lineages as well, some are still contested and need to be confirmed. For instance, data for ginkgo, but particularly cycads have remained inconclusive so far, likely due to the quality of the data available and flaws in the analysis. We extracted and sequenced RNA from both the cycad Encephalartos natalensis and Ginkgo biloba. This was followed by transcriptome assembly, after which these data were used to build paralog age distributions. Based on these distributions, we identified remnants of an ancient whole genome duplication in both cycads and ginkgo. The most parsimonious explanation would be that this whole genome duplication event was shared between both species and had occurred prior to their divergence, about 300 million years ago

De novo assembly of transcriptomes from a B73 maize line introgressed with a QTL for resistance to gray leaf spot disease reveals a candidate allele of a lectin receptor-like kinase

Gray leaf spot (GLS) disease in maize, caused by the fungus Cercospora zeina, is a threat to maize production globally. Understanding the molecular basis for quantitative resistance to GLS is therefore important for food security. We developed a de novo assembly pipeline to identify candidate maize resistance genes. Near-isogenic maize lines with and without a QTL for GLS resistance on chromosome 10 from inbred CML444 were produced in the inbred B73 background. The B73-QTL line showed a 20% reduction in GLS disease symptoms compared to B73 in the field (p = 0.01). B73-QTL leaf samples from this field experiment conducted under GLS disease pressure were RNA sequenced. The reads that did not map to the B73 or C. zeina genomes were expected to contain novel defense genes and were de novo assembled. A total of 141 protein-coding sequences with B73-like or plant annotations were identified from the B73-QTL plants exposed to C. zeina. To determine whether candidate gene expression was induced by C. zeina, the RNAseq reads from C. zeina-challenged and control leaves were mapped to a master assembly of all of the B73-QTL reads, and differential gene expression analysis was conducted. Combining results from both bioinformatics approaches led to the identification of a likely candidate gene, which was a novel allele of a lectin receptor-like kinase named L-RLK-CML that (i) was induced by C. zeina, (ii) was positioned in the QTL region, and (iii) had functional domains for pathogen perception and defense signal transduction. The 817AA L-RLK-CML protein had 53 amino acid differences from its 818AA counterpart in B73. A second "B73-like" allele of L-RLK was expressed at a low level in B73-QTL. Gene copy-specific RT-qPCR confirmed that the l-rlk-cml transcript was the major product induced four-fold by C. zeina. Several other expressed defense-related candidates were identified, including a wall-associated kinase, two glutathione s-transferases, a chitinase, a glucan beta-glucosidase, a plasmodesmata callose-binding protein, several other receptor-like kinases, and components of calcium signaling, vesicular trafficking, and ethylene biosynthesis. This work presents a bioinformatics protocol for gene discovery from de novo assembled transcriptomes and identifies candidate quantitative resistance genes

ORCAE-AOCC : a centralized portal for the annotation of African orphan crop genomes

ORCAE (Online Resource for Community Annotation of Eukaryotes) is a public genome annotation curation resource. ORCAE-AOCC is a branch that is dedicated to the genomes published as part of the African Orphan Crops Consortium (AOCC). The motivation behind the development of the ORCAE platform was to create a knowledge-based website where the research-community can make contributions to improve genome annotations. All changes to any given gene-model or gene description are stored, and the entire annotation history can be retrieved. Genomes can either be set to “public” or “restricted” mode; anonymous users can browse public genomes but cannot make any changes. Aside from providing a user- friendly interface to view genome annotations, the platform also includes tools and information (such as gene expression evidence) that enables authorized users to edit and validate genome annotations. The ORCAE-AOCC platform will enable various stakeholders from around the world to coordinate their efforts to annotate and study underutilized crops

Exploration of alternative splicing events in ten different grapevine cultivars

Background: The complex dynamics of gene regulation in plants are still far from being fully understood. Among many factors involved, alternative splicing (AS) in particular is one of the least well documented. For many years, AS has been considered of less relevant in plants, especially when compared to animals, however, since the introduction of next generation sequencing techniques the number of plant genes believed to be alternatively spliced has increased exponentially. Results: Here, we performed a comprehensive high-throughput transcript sequencing of ten different grapevine cultivars, which resulted in the first high coverage atlas of the grape berry transcriptome. We also developed findAS, a software tool for the analysis of alternatively spliced junctions. We demonstrate that at least 44 % of multi-exonic genes undergo AS and a large number of low abundance splice variants is present within the 131.622 splice junctions we have annotated from Pinot noir. Conclusions: Our analysis shows that similar to 70 % of AS events have relatively low expression levels, furthermore alternative splice sites seem to be enriched near the constitutive ones in some extent showing the noise of the splicing mechanisms. However, AS seems to be extensively conserved among the 10 cultivars

Horsetails are ancient polyploids : evidence from Equisetum giganteum

Horsetails represent an enigmatic clade within the land plants. Despite consisting only of one genus (Equisetum) that contains 15 species, they are thought to represent the oldest extant genus within the vascular plants dating back possibly as far as the Triassic. Horsetails have retained several ancient features and are also characterized by a particularly high chromosome count (n = 108). Whole-genome duplications (WGDs) have been uncovered in many angiosperm clades and have been associated with the success of angiosperms, both in terms of species richness and biomass dominance, but remain understudied in nonangiosperm clades. Here, we report unambiguous evidence of an ancient WGD in the fern linage, based on sequencing and de novo assembly of an expressed gene catalog (transcriptome) from the giant horsetail (Equisetum giganteum). We demonstrate that horsetails underwent an independent paleopolyploidy during the Late Cretaceous prior to the diversification of the genus but did not experience any recent polyploidizations that could account for their high chromosome number. We also discuss the specific retention of genes following the WGD and how this may be linked to their long-term survival

The year out

Figure S5. Distribution of up-regulated TF family members in embryos and megagametophytes. Presented data are based on TF family members differentially accumulated (FC > 2) during seed development in any of the pairwise comparisons between embryos and megagametophytes. Orange bars show the number of TFs belonging to each family in embryos and green bars in megagametophytes. TFs were classified into TF families by using the publicly available PlantTFDB v 3.0 database. Figure S6. Abundance of the ten largest TF families differentially expressed between embryos and megagametophytes during seed development shown in Fig. 5. Number of members in each TF family detected at different developmental stages in (A) embryos and (B) megagametophytes. Subordinate embryos were excluded from this analysis. (PDF 312 kb

Genomics of clinal local adaptation in Pinus sylvestris under continuous environmental and spatial genetic setting

Understanding the consequences of local adaptation at the genomic diversity is a central goal in evolutionary genetics of natural populations. In species with large continuous geographical distributions the phenotypic signal of local adaptation is frequently clear, but the genetic basis often remains elusive. We examined the patterns of genetic diversity inPinus sylvestris, a keystone species in many Eurasian ecosystems with a huge distribution range and decades of forestry research showing that it is locally adapted to the vast range of environmental conditions. MakingP. sylvestrisan even more attractive subject of local adaptation study, population structure has been shown to be weak previously and in this study. However, little is known about the molecular genetic basis of adaptation, as the massive size of gymnosperm genomes has prevented large scale genomic surveys. We generated a both geographically and genomically extensive dataset using a targeted sequencing approach. By applying divergence-based and landscape genomics methods we identified several loci contributing to local adaptation, but only few with large allele frequency changes across latitude. We also discovered a very large (ca. 300 Mbp) putative inversion potentially under selection, which to our knowledge is the first such discovery in conifers. Our results call for more detailed analysis of structural variation in relation to genomic basis of local adaptation, emphasize the lack of large effect loci contributing to local adaptation in the coding regions and thus point out the need for more attention toward multi-locus analysis of polygenic adaptation

In Silico Survey of the Mitochondrial Protein Uptake and Maturation Systems in the Brown Alga Ectocarpus siliculosus

The acquisition of mitochondria was a key event in eukaryote evolution. The aim of this study was to identify homologues of the components of the mitochondrial protein import machinery in the brown alga Ectocarpus and to use this information to investigate the evolutionary history of this fundamental cellular process. Detailed searches were carried out both for components of the protein import system and for related peptidases. Comparative and phylogenetic analyses were used to investigate the evolution of mitochondrial proteins during eukaryote diversification. Key observations include phylogenetic evidence for very ancient origins for many protein import components (Tim21, Tim50, for example) and indications of differences between the outer membrane receptors that recognize the mitochondrial targeting signals, suggesting replacement, rearrangement and/or emergence of new components across the major eukaryotic lineages. Overall, the mitochondrial protein import components analysed in this study confirmed a high level of conservation during evolution, indicating that most are derived from very ancient, ancestral proteins. Several of the protein import components identified in Ectocarpus, such as Tim21, Tim50 and metaxin, have also been found in other stramenopiles and this study suggests an early origin during the evolution of the eukaryotes

Evidence for an ancient whole genome duplication in the cycad lineage

Contrary to the many whole genome duplication events recorded for angiosperms (flowering plants), whole genome duplications in gymnosperms (non-flowering seed plants) seem to be much rarer. Although ancient whole genome duplications have been reported for most gymnosperm lineages as well, some are still contested and need to be confirmed. For instance, data for ginkgo, but particularly cycads have remained inconclusive so far, likely due to the quality of the data available and flaws in the analysis. We extracted and sequenced RNA from both the cycad Encephalartos natalensis and Ginkgo biloba. This was followed by transcriptome assembly, after which these data were used to build paralog age distributions. Based on these distributions, we identified remnants of an ancient whole genome duplication in both cycads and ginkgo. The most parsimonious explanation would be that this whole genome duplication event was shared between both species and had occurred prior to their divergence, about 300 million years ago.The UP Research Development Programme (RDP), the Genomics Research Institute (GRI), Multidisciplinary Research Partnership `Bioinformatics: from nucleotides to networks' Project (no. 01MR0310W) of Ghent University, and funding from the European Union Seventh Framework Programme (FP7/2007-2013) under European Research Council Advanced Grant Agreement 322739 ± DOUBLEUP.http://www.plosone.orgam2017Forestry and Agricultural Biotechnology Institute (FABI)Genetic