Sequence analysis of the genome of an oil-bearing tree, Jatropha curcas L

Differential amplification of the subtelomeric satellite DNA JcSAT1 in the genus Jatropha L. (Euphorbiaceae)

Satellite DNAs (satDNAs) are highly repetitive sequences that occur in virtually all eukaryotic genomes and can undergo rapid copy number and nucleotide sequence variation among relatives. After chromosomal mapping of the satDNA JcSAT1, it was found a large accumulation at subtelomeres of Jatropha curcas (subgenus Curcas), but an absence of these monomers in J. integerrima (subgenus Jatropha). This fact suggests a dynamic scenario for this satellite repeat in Jatropha genomes. Here, we used a multitasking approach (sequence analysis, DNA blotting and chromosomal mapping) to investigate the molecular organization and chromosomal abundance and distribution of JcSAT1 in a broader group of species from the subgenus Jatropha (J. gossypiifolia, J. mollissima, J. podagrica, and J. multifida) in addition to J. curcas, with the aiming of understanding the evolution of this satDNA. Based on the analysis of BAC clone sequences of J. curcas, a large array (~ 30 kb) of 80 homogeneous monomers of JcSAT1 was identified in BAC 23J11. The monomer size was conserved (~ 358 bp) and contained a telomeric motif at the 5' end. PCR amplification coupled with a Southern blot revealed the presence of JcSAT1-like sequences in all species examined. However, a large set of genome copies was identified only in J. curcas, where a ladder-like pattern with multimers of different sizes was observed. In situ hybridization of BAC 23J11 confirmed the subtelomeric pattern for J. curcas, but showed no signals on chromosomes of species from the subgenus Jatropha. Our data indicate that JcSAT1 is a highly homogeneous satDNA that originated from a region near the telomeres and spread throughout the chromosomal subtermini, possibly due to frequent ectopic recombination between these regions. The abundance of JcSAT1 in the genome of J. curcas suggests that an amplification event occurred either at the base of the subgenus Curcas or at least in this species, although the repeat is shared by all species of the genus studied so far.

Genomic basis of the giga-chromosomes and giga-genome of tree peony Paeonia ostii

Tree peony (Paeonia ostii) is an economically important ornamental plant native to China. It is also notable for its seed oil, which is abundant in unsaturated fatty acids such as α-linolenic acid (ALA). Here, we report chromosome-level genome assembly (12.28 Gb) of P. ostii. In contrast to monocots with giant genomes, tree peony does not appear to have undergone lineage-specific whole-genome duplication. Instead, explosive LTR expansion in the intergenic regions within a short period (~ two million years) may have contributed to the formation of its giga-genome. In addition, expansion of five types of histone encoding genes may have helped maintain the giga-chromosomes. Further, we conduct genome-wide association studies (GWAS) on 448 accessions and show expansion and high expression of several genes in the key nodes of fatty acid biosynthetic pathway, including SAD, FAD2 and FAD3, may function in high level of ALAs synthesis in tree peony seeds. Moreover, by comparing with cultivated tree peony (P. suffruticosa), we show that ectopic expression of class A gene AP1 and reduced expression of class C gene AG may contribute to the formation of petaloid stamens. Genomic resources reported in this study will be valuable for studying chromosome/genome evolution and tree peony breeding.

Genome-Wide Identification, Expression Patterns and Sugar Transport of the Physic Nut SWEET Gene Family and a Functional Analysis of JcSWEET16 in Arabidopsis

The Sugars Will Eventually be Exported Transporters (SWEET) family is a class of sugar transporters that play key roles in phloem loading, seed filling, pollen development and the stress response in plants. Here, a total of 18 JcSWEET genes were identified in physic nut (Jatropha curcas L.) and classified into four clades by phylogenetic analysis. These JcSWEET genes share similar gene structures, and alternative splicing of messenger RNAs was observed for five of the JcSWEET genes. Three (JcSWEET1/4/5) of the JcSWEETs were found to possess transport activity for hexose molecules in yeast. Real-time quantitative PCR analysis of JcSWEETs in different tissues under normal growth conditions and abiotic stresses revealed that most are tissue-specifically expressed, and 12 JcSWEETs responded to either drought or salinity. The JcSWEET16 gene responded to drought and salinity stress in leaves, and the protein it encodes is localized in both the plasma membrane and the vacuolar membrane. The overexpression of JcSWEET16 in Arabidopsis thaliana modified the flowering time and saline tolerance levels but not the drought tolerance of the transgenic plants. Together, these results provide insights into the characteristics of SWEET genes in physic nut and could serve as a basis for cloning and further functional analysis of these genes.

Genome-Wide Characterization of DGATs and Their Expression Diversity Analysis in Response to Abiotic Stresses in Brassica napus

Triacylglycerol (TAG) is the most important storage lipid for oil plant seeds. Diacylglycerol acyltransferases (DGATs) are a key group of rate-limiting enzymes in the pathway of TAG biosynthesis. In plants, there are three types of DGATs, namely, DGAT1, DGAT2 and DGAT3. Brassica napus, an allotetraploid plant, is one of the most important oil plants in the world. Previous studies of Brassica napus DGATs (BnaDGATs) have mainly focused on BnaDGAT1s. In this study, four DGAT1s, four DGAT2s and two DGAT3s were identified and cloned from B. napus ZS11. The analyses of sequence identity, chromosomal location and collinearity, phylogenetic tree, exon/intron gene structures, conserved domains and motifs, and transmembrane domain (TMD) revealed that BnaDGAT1, BnaDGAT2 and BnaDGAT3 were derived from three different ancestors and shared little similarity in gene and protein structures. Overexpressing BnaDGATs showed that only four BnaDGAT1s can restore TAG synthesis in yeast H1246 and promote the accumulation of fatty acids in yeast H1246 and INVSc1, suggesting that the three BnaDGAT subfamilies had greater differentiation in function. Transcriptional analysis showed that the expression levels of BnaDGAT1s, BnaDGAT2s and BnaDGAT3s were different during plant development and under different stresses. In addition, analysis of fatty acid contents in roots, stems and leaves under abiotic stresses revealed that P starvation can promote the accumulation of fatty acids, but no obvious relationship was shown between the accumulation of fatty acids with the expression of BnaDGATs under P starvation. This study provides an extensive evaluation of BnaDGATs and a useful foundation for dissecting the functions of BnaDGATs in biochemical and physiological processes.

Heterologous Expression of Jatropha curcas Fatty Acyl-ACP Thioesterase A (JcFATA) and B (JcFATB) Affects Fatty Acid Accumulation and Promotes Plant Growth and Development in Arabidopsis

Plant fatty acyl-acyl carrier protein (ACP) thioesterases terminate the process of de novo fatty acid biosynthesis in plastids by hydrolyzing the acyl-ACP intermediates, and determine the chain length and levels of free fatty acids. They are of interest due to their roles in fatty acid synthesis and their potential to modify plant seed oils through biotechnology. Fatty acyl-ACP thioesterases (FAT) are divided into two families, i.e., FATA and FATB, according to their amino acid sequence and substrate specificity. The high oil content in Jatropha curcas L. seed has attracted global attention due to its potential for the production of biodiesel. However, the detailed effects of JcFATA and JcFATB on fatty acid biosynthesis and plant growth and development are still unclear. In this study, we found that JcFATB transcripts were detected in all tissues and organs examined, with especially high accumulation in the roots, leaves, flowers, and some stages of developing seeds, and JcFATA showed a very similar expression pattern. Subcellular localization of the JcFATA-GFP and JcFATB-GFP fusion protein in Arabidopsis leaf protoplasts showed that both JcFATA and JcFATB localized in chloroplasts. Heterologous expression of JcFATA and JcFATB in Arabidopsis thaliana individually generated transgenic plants with longer roots, stems and siliques, larger rosette leaves, and bigger seeds compared with those of the wild type, indicating the overall promotion effects of JcFATA and JcFATB on plant growth and development while JcFATB had a larger impact. Compositional analysis of seed oil revealed that all fatty acids except 22:0 were significantly increased in the mature seeds of JcFATA-transgenic Arabidopsis lines, especially unsaturated fatty acids, such as the predominant fatty acids of seed oil, 18:1, 18:2, and 18:3. In the mature seeds of the JcFATB-transgenic Arabidopsis lines, most fatty acids were increased compared with those in wild type too, especially saturated fatty acids, such as 16:0, 18:0, 20:0, and 22:0. Our results demonstrated the promotion effect of JcFATA and JcFATB on plant growth and development, and their possible utilization to modify the seed oil composition and content in higher plants.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae

BACKGROUND

Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species.

RESULTS

To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree.

CONCLUSION

Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Genotyping by sequencing-based linkage map construction and identification of quantitative trait loci for yield-related traits and oil content in Jatropha (Jatropha curcas L.)

BACKGROUND

Jatropha (Jatropha curcas L.) has been considered as a potential bioenergy crop and its genetic improvement is essential for higher seed yield and oil content which has been hampered due to lack of desirable molecular markers.

METHODS AND RESULTS

An F₂ population was created using an intraspecific cross involving a Central American line RJCA9 and an Asiatic species RJCS-9 to develop a dense genetic map and for Quantitative trait loci (QTL) identification. The genotyping-by-sequencing (GBS) approach was used to genotype the mapping population of 136 F₂ individuals along with the two parental lines for classification of the genotypes based on single nucleotide polymorphism (SNPs). NextSeq 2500 sequencing technology provided a total of 517.23 million clean reads, with an average of ~ 3.8 million reads per sample. We analysed 411 SNP markers and developed 11 linkage groups. The total length of the genetic map was 4092.3 cM with an average marker interval of 10.04 cM. We have identified a total of 83 QTLs for various yield and oil content governing traits. The percentage of phenotypic variation (PV) was found to be in the range of 8.81 to 65.31%, and a QTL showed the maximum PV of 65.3% for a total seed number on the 6th linkage group (LG).

CONCLUSIONS

The QTLs detected in this study for various phenotypic traits will lay down the path for marker-assisted breeding in the future and cloning of genes that are responsible for phenotypic variation.

Ectopic Expression of JcCPL1, 2, and 4 Affects Epidermal Cell Differentiation, Anthocyanin Biosynthesis and Leaf Senescence in Arabidopsis thaliana

The CAPRICE (CPC)-like (CPL) genes belong to a single-repeat R3 MYB family, whose roles in physic nut (Jatropha curcas L.), an important energy plant, remain unclear. In this study, we identified a total of six CPL genes (JcCPL1–6) in physic nut. The JcCPL3, 4, and 6 proteins were localized mainly in the nucleus, while proteins JcCPL1, 2, and 5 were localized in both the nucleus and the cytoplasm. Ectopic overexpression of JcCPL1, 2, and 4 in Arabidopsis thaliana resulted in an increase in root hair number and decrease in trichome number. Consistent with the phenotype of reduced anthocyanin in shoots, the expression levels of anthocyanin biosynthesis genes were down-regulated in the shoots of these three transgenic A. thaliana lines. Moreover, we observed that OeJcCPL1, 2, 4 plants attained earlier leaf senescence, especially at the late developmental stage. Consistent with this, the expression levels of several senescence-associated and photosynthesis-related genes were, respectively, up-regulated and down-regulated in leaves. Taken together, our results indicate functional divergence of the six CPL proteins in physic nut. These findings also provide insight into the underlying roles of CPL transcription factors in leaf senescence.

Pathogen-driven coevolution across the CBP60 plant immune regulator subfamilies confers resilience on the regulator module

Components of the plant immune signaling network need mechanisms that confer resilience against fast-evolving pathogen effectors that target them. Among eight Arabidopsis CaM-Binding Protein (CBP) 60 family members, AtCBP60g and AtSARD1 are partially functionally redundant, major positive immune regulators, and AtCBP60a is a negative immune regulator. We investigated possible resilience-conferring evolutionary mechanisms among the CBP60a, CBP60g and SARD1 immune regulatory subfamilies. Phylogenetic analysis was used to investigate the times of CBP60 subfamily neofunctionalization. Then, using the pairwise distance rank based on the newly developed analytical platform Protein Evolution Analysis in a Euclidean Space (PEAES), hypotheses of specific coevolutionary mechanisms that could confer resilience on the regulator module were tested. The immune regulator subfamilies diversified around the time of angiosperm divergence and have been evolving very quickly. We detected significant coevolutionary interactions across the immune regulator subfamilies in all of 12 diverse core eudicot species lineages tested. The coevolutionary interactions were consistent with the hypothesized coevolution mechanisms. Despite their unusually fast evolution, members across the CBP60 immune regulator subfamilies have influenced the evolution of each other long after their diversification in a way that could confer resilience on the immune regulator module against fast-evolving pathogen effectors.

Genome-wide characterization of bZIP transcription factors and their expression patterns in response to drought and salinity stress in Jatropha curcas

The basic leucine zipper (bZIP) family is one of the largest families of transcription factors (TFs) in plants and is responsible for various functions, including regulating development and responses to abiotic/biotic stresses. However, the roles of bZIPs in the regulation of responses to drought stress and salinity stress remain poorly understood in Jatropha curcas L., a biodiesel crop. In the present study, 50 JcbZIP genes were identified and classified into ten groups. Cis-element analysis indicated that JcbZIP genes are associated with abiotic stress. Gene expression patterns and quantitative real-time PCR (qRT-PCR) showed that four JcbZIP genes (JcbZIPs 34, 36, 49 and 50) are key resistance-related genes under both drought and salinity stress conditions. On the basis of the results of cis-element and phylogenetic analyses, JcbZIP49 and JcbZIP50 are likely involved in responses to drought and salinity stress; moreover, JcbZIP34 and JcbZIP36 might also play important roles in seed development and response to abiotic stress. These findings advance our understanding of the comprehensive characteristics of JcbZIP genes and provide new insights for functional validation in the further.

Monitoring casbene synthase in Jatropha curcas tissues using targeted proteomics

BACKGROUND

Casbene synthase (CS) is responsible for the first committed step in the biosynthesis of phorbol esters (PE) in the Euphorbiaceae. PE are abundant in the seeds of the biofuel crop Jatropha curcas and its toxicity precludes the use of the protein-rich cake obtained after oil extraction as an animal feed and the toxicity of the fumes derived from burning PE containing biofuel is also a matter of concern. This toxicity is a major hindrance to exploit the potential of this crop as a source of raw material to produce biodiesel. For this reason, the current research on J. curcas is mainly focused on the understanding of the biosynthesis and site of synthesis of PE, as an avenue for the development of genotypes unable to synthesize PE in its seeds.

RESULTS

Here, we present targeted proteomics assays (SRM and PRM) to detect and quantify CS in leaves, endosperm, and roots of two J. curcas genotypes with contrasting levels of PE. These assays were based on the use of reference isotopic labeled synthetic peptides (ILSP) predicted from 12 gene models of CS from the J. curcas genome.

CONCLUSION

Our targeted proteomics methods were able to detect and quantify, for the first time, CS gene products and demonstrate the distribution of CS isoforms only in roots from J. curcas genotypes with a high and low concentration of PE. These methods can be expanded to monitor CS, at the protein level, in different tissues and genotypes of J. curcas.

Identification and characterization of long non-coding RNA (lncRNA) in the developing seeds of Jatropha curcas

Long non-coding RNAs (lncRNAs) play critical roles in plant development. However, the information of lncRNAs in Jatropha curcas remains largely unexplored. Thus, an attempt has been made in J. curcas to identify 1,850 lncRNAs based on deep sequencing of developing seeds at three typical stages. About ten percent lncRNAs (196 lncRNAs) were differentially expressed lncRNAs during seed developing process. Together with reverse transcription quantitative real-time PCR, the lncRNA expression analyses revealed the stage-specific expression patterns of some novel lncRNAs in J. curcas. The target genes of lncRNAs were annotated for their roles in various biological processes such as gene expression, metabolism, and cell growth. Besides, 10 lncRNAs were identified as the precursors of microRNAs and 26 lncRNAs were predicted to be the targets of Jatropha miRNAs. A total of 31 key lncRNAs play critical roles in the seed developing process in the context of cell growth and development, lipid metabolism, and seed maturation. Our study provides the first systematic study of lncRNAs in the developing seeds of J. curcas and facilitates the functional research of plant lncRNAs and the regulation of seed development.

Development of genome wide transposable elements based repeat junction markers in Jatropha (Jatropha curcas L.)

Jatropha curcas is a potential biodiesel crop and a highly adaptable species to various agro-climatic conditions. In this study, we have utilized transposable elements' (TE) repeat junctions (RJs) which are an important constituent of the genome, used to form a genome-wide molecular marker platform owing to its use in genomic studies of plants. We screened our previously generated Jatropha hybrid genome assembly of size 265 Mbp using RJPrimers pipeline software and identified a total of 1274 TE junctions. For the predicted RJs, we designed 2868 polymerase chain reaction (PCR) based RJ markers (RJMs) flanking the junction regions. In addition to marker design, the identified RJs were utilized to detect 225,517 TEs across the genome. The different types of transposable repeat elements mainly were scattered into Retro, LTR, Copia and Gypsy categories. The efficacy of the designed markers was tested by utilizing a subset of RJMs selected randomly. We have validated 96 randomly selected RJ primers in a group of 32 J. curcas genotypes and more than 90% of the markers effectively intensified as amplicons. Of these, 10 primers were shown to be polymorphic in estimating genetic diversity among the 32 Jatropha lines. UPGMA cluster analysis revealed the formation of two clusters such as A and B exhibiting 85.5% and 87% similarity coefficient respectively. The various RJMs identified in this study could be utilized as a significant asset in Jatropha functional genomics including genome determination, mapping and marker-assisted selection.

Identification of potential candidate genes controlling pea aphid tolerance in a Pisum fulvum high-density integrated DArTseq SNP-based genetic map

BACKGROUND

Pea (Pisum sativum) is one of the most important temperate grain legumes in the world, and its production is severely constrained by the pea aphid (Acyrthosiphon pisum). Wild relatives, such as P. fulvum, are valuable sources of allelic diversity to improve the genetic resistance of cultivated pea species against A. pisum attack. To unravel the genetic control underlying resistance to the pea aphid attack, a quantitative trait loci (QTL) analysis was performed using the previously developed high density integrated genetic linkage map originated from an intraspecific recombinant inbred line (RIL) population (P. fulvum: IFPI3260 × IFPI3251).

RESULTS

We accurately evaluated specific resistance responses to pea aphid that allowed the identification, for the first time, of genomic regions that control plant damage and aphid reproduction. Eight QTLs associated with tolerance to pea aphid were identified in LGs I, II, III, IV and V, which individually explained from 17.0% to 51.2% of the phenotypic variation depending on the trait scored, and as a whole from 17.0% to 88.6%. The high density integrated genetic linkage map also allowed the identification of potential candidate genes co-located with the QTLs identified.

CONCLUSIONS

Our work shows how the survival of P. fulvum after the pea aphid attack depends on the triggering of a multi-component protection strategy that implies a quantitative tolerance. The genomic regions associated with the tolerance responses of P. fulvum during A. pisum infestation have provided six potential candidate genes that could be useful in marker-assisted selection (MAS) and genomic assisted breeding (GAB) after functional validation in the future. © 2019 Society of Chemical Industry.

Gene expression profiling identifies pathways involved in seed maturation of Jatropha curcas

Background

Jatropha curcas, a tropical shrub, is a promising biofuel crop, which produces seeds with high content of oil and protein. To better understand the maturation process of J. curcas seeds and to improve its agronomic performance, a two-step approach was performed in six different maturation stages of seeds: 1) generation of the entire transcriptome of J. curcas seeds using 454-Roche sequencing of a cDNA library, 2) comparison of transcriptional expression levels using a custom Agilent 8x60K oligonucleotide microarray.

Results

A total of 793,875 high-quality reads were assembled into 19,382 unique full-length contigs, of which 13,507 could be annotated with Gene Ontology (GO) terms. Microarray data analysis identified 9111 probes (out of 57,842 probes), which were differentially expressed between the six maturation stages. The expression results were validated for 75 selected transcripts based on expression levels, predicted function, pathway, and length.

Result from cluster analyses showed that transcripts associated with fatty acid, flavonoid, and phenylpropanoid biosynthesis were over-represented in the early stages, while those of lipid storage were over-represented in the late stages. Expression analyses of different maturation stages of J. curcas seed showed that most changes in transcript abundance occurred between the two last stages, suggesting that the timing of metabolic pathways during seed maturation in J. curcas occurs in late stages. The co-expression results showed that the hubs (CB5-D, CDR1, TT8, DFR, HVA22) with the highest number of edges, associated with fatty acid and flavonoid biosynthesis, are showing a decrease in their expression during seed maturation. Furthermore, seed development and hormone pathways are significantly well connected.

Conclusion

The obtained results revealed differentially expressed sequences (DESs) regulating important pathways related to seed maturation, which could contribute to the understanding of the complex regulatory network during seed maturation with the focus on lipid, flavonoid and phenylpropanoid biosynthesis. This study provides detailed information on transcriptional changes during J. curcas seed maturation and provides a starting point for a genomic survey of seed quality traits. The results highlighted specific genes and processes relevant to the molecular mechanisms involved in Jatropha seed maturation. These data can also be utilized regarding other Euphorbiaceae species.

Exploitation of Hi-C sequencing for improvement of genome assembly and in-vitro validation of differentially expressing genes in Jatropha curcas L

Jatropha curcas is one of the major sources of renewable energy due to potential use of its oil as a biofuel. The genome of this crop is constituted by the high content of repetitive elements. We employed the Hi-C proximity ligation technique to re-scaffold our existing hybrid genome assembly of an elite genotype (RJC1) developed using Illumina and Pacbio technologies. We assembled 99.81% of non-truncated reads to achieve 266.80 Mbp of the genome with an N50 value of 1.58 Mb. Furthermore, we compared the efficiency of Hi-C-augmented genome assembly with the hybrid genome assembly and observed a ~ 50% reduction in scaffolds and a tenfold increase in the N50 value. The gene ontology analysis revealed the identification of terms for molecular function (45.52%), cellular component (33.47%), and biological function (20.99%). Comparative genomic analysis of 13-plant species showed the conservation of 414 lipid metabolizing genes identified in the KEGG pathway analysis. Differential gene expression (DGE) studies were conducted in the healthy and Jatropha mosaic virus-infected leaves via RNA-seq analysis and observed gene expression changes for 2185 genes. Out of these, we observed 546 genes having more than two-fold change of transcript level and among these 259 genes were down-regulated and 287 genes were up-regulated. To validate RNA-seq data, two DEGs were selected for gene expression analysis using qRT-PCR and the data was in correlation with in silico results. RNA-seq analysis further shows the identification of some of the candidate genes and may be useful to develop JMV resistant plants after functional validation. This Hi-C genome assembly provides a detailed accurate reference genome which could be utilized to improve Jatropha and other economically important Euphorbiaceae family members.

The Chromosome-Based Rubber Tree Genome Provides New Insights into Spurge Genome Evolution and Rubber Biosynthesis

The rubber tree, Hevea brasiliensis, produces natural rubber that serves as an essential industrial raw material. Here, we present a high-quality reference genome for a rubber tree cultivar GT1 using single-molecule real-time sequencing (SMRT) and Hi-C technologies to anchor the ∼1.47-Gb genome assembly into 18 pseudochromosomes. The chromosome-based genome analysis enabled us to establish a model of spurge chromosome evolution, since the common paleopolyploid event occurred before the split of Hevea and Manihot. We show recent and rapid bursts of the three Hevea-specific LTR-retrotransposon families during the last 10 million years, leading to the massive expansion by ∼65.88% (∼970 Mbp) of the whole rubber tree genome since the divergence from Manihot. We identify large-scale expansion of genes associated with whole rubber biosynthesis processes, such as basal metabolic processes, ethylene biosynthesis, and the activation of polysaccharide and glycoprotein lectin, which are important properties for latex production. A map of genomic variation between the cultivated and wild rubber trees was obtained, which contains ∼15.7 million high-quality single-nucleotide polymorphisms. We identified hundreds of candidate domestication genes with drastically lowered genomic diversity in the cultivated but not wild rubber trees despite a relatively short domestication history of rubber tree, some of which are involved in rubber biosynthesis. This genome assembly represents key resources for future rubber tree research and breeding, providing novel targets for improving plant biotic and abiotic tolerance and rubber production.

De novo genome assembly and Hi-C analysis reveal an association between chromatin architecture alterations and sex differentiation in the woody plant Jatropha curcas

BACKGROUND

Chromatin architecture is an essential factor regulating gene transcription in different cell types and developmental phases. However, studies on chromatin architecture in perennial woody plants and on the function of chromatin organization in sex determination have not been reported.

RESULTS

Here, we produced a chromosome-scale de novo genome assembly of the woody plant Jatropha curcas with a total length of 379.5 Mb and a scaffold N50 of 30.7 Mb using Pacific Biosciences long reads combined with genome-wide chromosome conformation capture (Hi-C) technology. Based on this high-quality reference genome, we detected chromatin architecture differences between monoecious and gynoecious inflorescence buds of Jatropha. Differentially expressed genes were significantly enriched in the changed A/B compartments and topologically associated domain regions and occurred preferentially in differential contact regions between monoecious and gynoecious inflorescence buds. Twelve differentially expressed genes related to flower development or hormone synthesis displayed significantly different genomic interaction patterns in monoecious and gynoecious inflorescence buds. These results demonstrate that chromatin organization participates in the regulation of gene transcription during the process of sex differentiation in Jatropha.

CONCLUSIONS

We have revealed the features of chromatin architecture in perennial woody plants and investigated the possible function of chromatin organization in Jatropha sex differentiation. These findings will facilitate understanding of the regulatory mechanisms of sex determination in higher plants.

JCDB: a comprehensive knowledge base for Jatropha curcas, an emerging model for woody energy plants

Background

Jatropha curcas is an oil-bearing plant, and has seeds with high oil content (~ 40%). Several advantages, such as easy genetic transformation and short generation duration, have led to the emergence of J. curcas as a model for woody energy plants. With the development of high-throughput sequencing, the genome of Jatropha curcas has been sequenced by different groups and a mass of transcriptome data was released. How to integrate and analyze these omics data is crucial for functional genomics research on J. curcas.

Results

By establishing pipelines for processing novel gene identification, gene function annotation, and gene network construction, we systematically integrated and analyzed a series of J. curcas transcriptome data. Based on these data, we constructed a J. curcas database (JCDB), which not only includes general gene information, gene functional annotation, gene interaction networks, and gene expression matrices but also provides tools for browsing, searching, and downloading data, as well as online BLAST, the JBrowse genome browser, ID conversion, heatmaps, and gene network analysis tools.

Conclusions

JCDB is the most comprehensive and well annotated knowledge base for J. curcas. We believe it will make a valuable contribution to the functional genomics study of J. curcas. The database is accessible at http://jcdb.xtbg.ac.cn.

Novel findings to the biosynthetic pathway of magnoflorine and taspine through transcriptomic and metabolomic analysis of Croton draco (Euphorbiaceae)

BACKGROUND

Croton draco is an arboreal species and its latex as well as some other parts of the plant, are traditionally used in the treatment of a wide range of ailments and diseases. Alkaloids, such as magnoflorine, prevent early atherosclerosis progression while taspine, an abundant constituent of latex, has been described as a wound-healer and antitumor-agent. Despite the great interest for these and other secondary metabolites, no omics resources existed for the species and the biosynthetic pathways of these alkaloids remain largely unknown.

RESULTS

To gain insights into the pathways involved in magnoflorine and taspine biosynthesis by C. draco and identify the key enzymes in these processes, we performed an integrated analysis of the transcriptome and metabolome in the major organs (roots, stem, leaves, inflorescences, and flowers) of this species. Transcript profiles were generated through high-throughput RNA-sequencing analysis while targeted and high resolution untargeted metabolomic profiling was also performed. The biosynthesis of these compounds appears to occur in the plant organs examined, but intermediaries may be translocated from the cells in which they are produced to other cells in which they accumulate.

CONCLUSIONS

Our results provide a framework to better understand magnoflorine and taspine biosynthesis in C. draco. In addition, we demonstrate the potential of multi-omics approaches to identify candidate genes involved in the biosynthetic pathways of interest.

Polysaccharide biosynthetic pathway profiling and putative gene mining of Dendrobium moniliforme using RNA-Seq in different tissues

BACKGROUND

Dendrobium moniliforme (Linnaeus) Swartz is a well-known plant used in traditional Chinese medicine due to bioactive constituents. Polysaccharides are the main medicinal ingredients, yet no studies have been published on polysaccharide biosynthesis in D. moniliforme. To comprehensively investigate the polysaccharide at the transcription level, we performed de novo transcriptome sequencing for the first time to produce a comprehensive transcriptome of D. moniliforme.

RESULTS

In our study, a database of 562,580 unigenes (average length = 1115.67 bases) was generated by performing transcriptome sequencing. Based on the gene annotation of the transcriptome, we identified 1204 carbohydrate-active related unigenes against CAZy database, including 417 glycosyltransferase genes (GTs), 780 glycoside hydrolases (GHs), 19 carbohydrate esterases (CEs), 75 carbohydrate-binding modules (CBMs), and 44 polysaccharide lyases (PLs). In the cellulose synthase family, 21 differential expression genes (DEGs) related to polysaccharide were identified. Subsequently, the tissue-specific expression patterns of the genes involved in polysaccharide pathway were investigated, which provide understanding of the biosynthesis and regulation of DMP at the molecular level. The two key enzyme genes (Susy and SPS) involved in the polysaccharide pathway were identified, and their expression patterns in different tissues were further analyzed using quantitative real-time PCR.

CONCLUSIONS

We determined the content of polysaccharides from Dendrobium moniliforme under different tissues, and we obtained a large number of differential genes by transcriptome sequencing. This database provides a pool of candidate genes involved in biosynthesis of polysaccharides in D. moniliforme. Furthermore, the comprehensive analysis and characterization of the significant pathways are expected to give a better insight regarding the diversity of chemical composition, synthetic characteristics, and the regulatory mechanism which operate in this medical herb.

Hybrid de novo transcriptome assembly of poinsettia (Euphorbia pulcherrima Willd. Ex Klotsch) bracts

BACKGROUND

Poinsettia is a popular and important ornamental crop, mostly during the Christmas season. Its bract coloration ranges from pink/red to creamy/white shades. Despite its ornamental value, there is a lack of knowledge about the genetics and molecular biology of poinsettia, especially on the mechanisms of color formation. We performed an RNA-Seq analysis in order to shed light on the transcriptome of poinsettia bracts. Moreover, we analyzed the transcriptome differences of red- and white-bracted poinsettia varieties during bract development and coloration. For the assembly of a bract transcriptome, two paired-end cDNA libraries from a red and white poinsettia pair were sequenced with the Illumina technology, and one library from a red-bracted variety was used for PacBio sequencing. Both short and long reads were assembled using a hybrid de novo strategy. Samples of red- and white-bracted poinsettias were sequenced and comparatively analyzed in three color developmental stages in order to understand the mechanisms of color formation and accumulation in the species.

RESULTS

The final transcriptome contains 288,524 contigs, with 33% showing confident protein annotation against the TAIR10 database. The BUSCO pipeline, which is based on near-universal orthologous gene groups, was applied to assess the transcriptome completeness. From a total of 1440 BUSCO groups searched, 77% were categorized as complete (41% as single-copy and 36% as duplicated), 10% as fragmented and 13% as missing BUSCOs. The gene expression comparison between red and white varieties of poinsettia showed a differential regulation of the flavonoid biosynthesis pathway only at particular stages of bract development. An initial impairment of the flavonoid pathway early in the color accumulation process for the white poinsettia variety was observed, but these differences were no longer present in the subsequent stages of bract development. Nonetheless, GSTF11 and UGT79B10 showed a lower expression in the last stage of bract development for the white variety and, therefore, are potential candidates for further studies on poinsettia coloration.

CONCLUSIONS

In summary, this transcriptome analysis provides a valuable foundation for further studies on poinsettia, such as plant breeding and genetics, and highlights crucial information on the molecular mechanism of color formation.

Transcriptome analysis of two inflorescence branching mutants reveals cytokinin is an important regulator in controlling inflorescence architecture in the woody plant Jatropha curcas

BACKGROUND

In higher plants, inflorescence architecture is an important agronomic trait directly determining seed yield. However, little information is available on the regulatory mechanism of inflorescence development in perennial woody plants. Based on two inflorescence branching mutants, we investigated the transcriptome differences in inflorescence buds between two mutants and wild-type (WT) plants by RNA-Seq to identify the genes and regulatory networks controlling inflorescence architecture in Jatropha curcas L., a perennial woody plant belonging to Euphorbiaceae.

RESULTS

Two inflorescence branching mutants were identified in germplasm collection of Jatropha. The duo xiao hua (dxh) mutant has a seven-order branch inflorescence, and the gynoecy (g) mutant has a three-order branch inflorescence, while WT Jatropha has predominantly four-order branch inflorescence, occasionally the three- or five-order branch inflorescences in fields. Using weighted gene correlation network analysis (WGCNA), we identified several hub genes involved in the cytokinin metabolic pathway from modules highly associated with inflorescence phenotypes. Among them, Jatropha ADENOSINE KINASE 2 (JcADK2), ADENINE PHOSPHORIBOSYL TRANSFERASE 1 (JcAPT1), CYTOKININ OXIDASE 3 (JcCKX3), ISOPENTENYLTRANSFERASE 5 (JcIPT5), LONELY GUY 3 (JcLOG3) and JcLOG5 may participate in cytokinin metabolic pathway in Jatropha. Consistently, exogenous application of cytokinin (6-benzyladenine, 6-BA) on inflorescence buds induced high-branch inflorescence phenotype in both low-branch inflorescence mutant (g) and WT plants. These results suggested that cytokinin is an important regulator in controlling inflorescence branching in Jatropha. In addition, comparative transcriptome analysis showed that Arabidopsis homologous genes Jatropha AGAMOUS-LIKE 6 (JcAGL6), JcAGL24, FRUITFUL (JcFUL), LEAFY (JcLFY), SEPALLATAs (JcSEPs), TERMINAL FLOWER 1 (JcTFL1), and WUSCHEL-RELATED HOMEOBOX 3 (JcWOX3), were differentially expressed in inflorescence buds between dxh and g mutants and WT plants, indicating that they may participate in inflorescence development in Jatropha. The expression of JcTFL1 was downregulated, while the expression of JcLFY and JcAP1 were upregulated in inflorescences in low-branch g mutant.

CONCLUSIONS

Cytokinin is an important regulator in controlling inflorescence branching in Jatropha. The regulation of inflorescence architecture by the genes involved in floral development, including TFL1, LFY and AP1, may be conservative in Jatropha and Arabidopsis. Our results provide helpful information for elucidating the regulatory mechanism of inflorescence architecture in Jatropha.

Identification and validation of SNP markers linked to seed toxicity in Jatropha curcas L

Edible/non-toxic varieties of Jatropha curcas L. are gaining increasing attention, providing both oil as biofuel feedstock or even as edible oil and the seed kernel meal as animal feed ingredient. They are a viable alternative to the limitation posed by the presence of phorbol esters in toxic varieties. Accurate genotyping of toxic/non-toxic accessions is critical to breeding management. The aim of this study was to identify SNP markers linked to seed toxicity in J. curcas. For SNP discovery, NGS technology was used to sequence the whole genomes of a toxic and non-toxic parent along with a bulk of 51 toxic and 30 non-toxic F₂ plants. To ascertain the association between SNP markers and seed toxicity trait, candidate SNPs were genotyped on 672 individuals segregating for seed toxicity and two collections of J. curcas composed of 96 individuals each. In silico SNP discovery approaches led to the identification of 64 candidate SNPs discriminating non-toxic and toxic samples. These SNPs were mapped on Chromosome 8 within the Linkage Group 8 previously identified as a genomic region important for phorbol ester biosynthesis. The association study identified two new SNPs, SNP_J22 and SNP_J24 significantly linked to low toxicity with R² values of 0.75 and 0.54, respectively. Our study released two valuable SNP markers for high-throughput, marker-assisted breeding of seed toxicity in J. curcas.

Early Sex-Chromosome Evolution in the Diploid Dioecious Plant Mercurialis annua

Suppressed recombination allows divergence between homologous sex chromosomes and the functionality of their genes. Here, we reveal patterns of the earliest stages of sex-chromosome evolution in the diploid dioecious herb Mercurialis annua on the basis of cytological analysis, de novo genome assembly and annotation, genetic mapping, exome resequencing of natural populations, and transcriptome analysis. The genome assembly contained 34,105 expressed genes, of which 10,076 were assigned to linkage groups. Genetic mapping and exome resequencing of individuals across the species range both identified the largest linkage group, LG1, as the sex chromosome. Although the sex chromosomes of M. annua are karyotypically homomorphic, we estimate that about one-third of the Y chromosome, containing 568 transcripts and spanning 22.3 cM in the corresponding female map, has ceased recombining. Nevertheless, we found limited evidence for Y-chromosome degeneration in terms of gene loss and pseudogenization, and most X- and Y-linked genes appear to have diverged in the period subsequent to speciation between M. annua and its sister species M. huetii, which shares the same sex-determining region. Taken together, our results suggest that the M. annua Y chromosome has at least two evolutionary strata: a small old stratum shared with M. huetii, and a more recent larger stratum that is probably unique to M. annua and that stopped recombining ∼1 MYA. Patterns of gene expression within the nonrecombining region are consistent with the idea that sexually antagonistic selection may have played a role in favoring suppressed recombination.

Identifying High Confidence microRNAs in the Developing Seeds of Jatropha curcas

MicroRNAs (miRNAs) are endogenously short noncoding regulatory RNAs implicated in plant development and physiology. Nine small RNA (sRNA) libraries from three typical seed developmental stages (young, intermediate, and mature) were generated by deep sequencing to identify the miRNAs of J. curcas, a potential oilseed crop for the production of renewable oil. Strict criteria were adopted to identify 93 high confidence miRNAs including 48 conserved miRNAs and 45 novel miRNAs. Target genes of these miRNAs were involved in a broad range of physiological functions, including gene expression regulation, primary & secondary metabolism, growth & development, signal transduction, and stress response. About one third (29 out of 93) miRNAs showed significant changes in expression levels during the seed developmental process, indicating that the miRNAs might regulate its targets by their changes of transcription levels in seed development. However, most miRNAs were found differentially expressed in the late stage of seed development, suggesting that miRNAs play more important roles in the stage when seed accumulating organic matters and suffering dehydration stress. This study presents the first large scale identification of high confidence miRNAs in the developing seeds of J. curcas.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

Genome sequence of Jatropha curcas L., a non-edible biodiesel plant, provides a resource to improve seed-related traits

Jatropha curcas (physic nut), a non-edible oilseed crop, represents one of the most promising alternative energy sources due to its high seed oil content, rapid growth and adaptability to various environments. We report ~339 Mbp draft whole genome sequence of J. curcas var. Chai Nat using both the PacBio and Illumina sequencing platforms. We identified and categorized differentially expressed genes related to biosynthesis of lipid and toxic compound among four stages of seed development. Triacylglycerol (TAG), the major component of seed storage oil, is mainly synthesized by phospholipid:diacylglycerol acyltransferase in Jatropha, and continuous high expression of homologs of oleosin over seed development contributes to accumulation of high level of oil in kernels by preventing the breakdown of TAG. A physical cluster of genes for diterpenoid biosynthetic enzymes, including casbene synthases highly responsible for a toxic compound, phorbol ester, in seed cake, was syntenically highly conserved between Jatropha and castor bean. Transcriptomic analysis of female and male flowers revealed the up-regulation of a dozen family of TFs in female flower. Additionally, we constructed a robust species tree enabling estimation of divergence times among nine Jatropha species and five commercial crops in Malpighiales order. Our results will help researchers and breeders increase energy efficiency of this important oil seed crop by improving yield and oil content, and eliminating toxic compound in seed cake for animal feed.

Seed size evolution and biogeography of Plukenetia (Euphorbiaceae), a pantropical genus with traditionally cultivated oilseed species

BACKGROUND

Plukenetia is a small pantropical genus of lianas and vines with variably sized edible oil-rich seeds that presents an ideal system to investigate neotropical and pantropical diversification patterns and seed size evolution. We assessed the biogeography and seed evolution of Plukenetia through phylogenetic analyses of a 5069 character molecular dataset comprising five nuclear and two plastid markers for 86 terminals in subtribe Plukenetiinae (representing 20 of ~ 23 Plukenetia species). Two nuclear genes, KEA1 and TEB, were used for phylogenetic reconstruction for the first time. Our goals were: (1) produce a robust, time-dependent evolutionary framework for Plukenetia using BEAST; (2) reconstruct its biogeographical history with ancestral range estimation in BIOGEOBEARS; (3) define seed size categories; (4) identify patterns of seed size evolution using ancestral state estimation; and (5) conduct regression analyses with putative drivers of seed size using the threshold model.

RESULTS

Plukenetia was resolved into two major groups, which we refer to as the pinnately- and palmately-veined clades. Our analyses suggest Plukenetia originated in the Amazon or Atlantic Forest of Brazil during the Oligocene (28.7 Mya) and migrated/dispersed between those regions and Central America/Mexico throughout the Miocene. Trans-oceanic dispersals explain the pantropical distribution of Plukenetia, including from the Amazon to Africa in the Early Miocene (17.4 Mya), followed by Africa to Madagascar and Africa to Southeast Asia in the Late Miocene (9.4 Mya) and Pliocene (4.5 Mya), respectively. We infer a single origin of large seeds in the ancestor of Plukenetia. Seed size fits a Brownian motion model of trait evolution and is moderately to strongly associated with plant size, fruit type/dispersal syndrome, and seedling ecology. Biome shifts were not drivers of seed size, although there was a weak association with a transition to fire prone semi-arid savannas.

CONCLUSIONS

The major relationships among the species of Plukenetia are now well-resolved. Our biogeographical analyses support growing evidence that many pantropical distributions developed by periodic trans-oceanic dispersals throughout the Miocene and Pliocene. Selection on a combination of traits contributed to seed size variation, while movement between forest edge/light gap and canopy niches likely contributed to the seed size extremes in Plukenetia.

De Novo Sequencing and Hybrid Assembly of the Biofuel Crop Jatropha curcas L.: Identification of Quantitative Trait Loci for Geminivirus Resistance

Jatropha curcas is an important perennial, drought tolerant plant that has been identified as a potential biodiesel crop. We report here the hybrid de novo genome assembly of J. curcas generated using Illumina and PacBio sequencing technologies, and identification of quantitative loci for Jatropha Mosaic Virus (JMV) resistance. In this study, we generated scaffolds of 265.7 Mbp in length, which correspond to 84.8% of the gene space, using Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Additionally, 96.4% of predicted protein-coding genes were captured in RNA sequencing data, which reconfirms the accuracy of the assembled genome. The genome was utilized to identify 12,103 dinucleotide simple sequence repeat (SSR) markers, which were exploited in genetic diversity analysis to identify genetically distinct lines. A total of 207 polymorphic SSR markers were employed to construct a genetic linkage map for JMV resistance, using an interspecific F₂ mapping population involving susceptible J. curcas and resistant Jatropha integerrima as parents. Quantitative trait locus (QTL) analysis led to the identification of three minor QTLs for JMV resistance, and the same has been validated in an alternate F₂ mapping population. These validated QTLs were utilized in marker-assisted breeding for JMV resistance. Comparative genomics of oil-producing genes across selected oil producing species revealed 27 conserved genes and 2986 orthologous protein clusters in Jatropha. This reference genome assembly gives an insight into the understanding of the complex genetic structure of Jatropha, and serves as source for the development of agronomically improved virus-resistant and oil-producing lines.

Global Phosphoproteomic Analysis Reveals the Defense and Response Mechanisms of Jatropha Curcas Seedling under Chilling Stress

As a promising energy plant for biodiesel, Jatropha curcas is a tropical and subtropical shrub and its growth is affected by one of major abiotic stress, chilling. Therefore, we adopt the phosphoproteomic analysis, physiological measurement and ultrastructure observation to illustrate the responsive mechanism of J. curcas seedling under chilling (4 °C) stress. After chilling for 6 h, 308 significantly changed phosphoproteins were detected. Prolonged the chilling treatment for 24 h, obvious physiological injury can be observed and a total of 332 phosphoproteins were examined to be significantly changed. After recovery (28 °C) for 24 h, 291 phosphoproteins were varied at the phosphorylation level. GO analysis showed that significantly changed phosphoproteins were mainly responsible for cellular protein modification process, transport, cellular component organization and signal transduction at the chilling and recovery periods. On the basis of protein-protein interaction network analysis, phosphorylation of several protein kinases, such as SnRK2, MEKK1, EDR1, CDPK, EIN2, EIN4, PI4K and 14-3-3 were possibly responsible for cross-talk between ABA, Ca²⁺, ethylene and phosphoinositide mediated signaling pathways. We also highlighted the phosphorylation of HOS1, APX and PIP2 might be associated with response to chilling stress in J. curcas seedling. These results will be valuable for further study from the molecular breeding perspective.

Two jatropha karyotypes constructed from meiotic pachytene chromosomes: Pericentric distribution of heterochromatin and variation in repetitive DNAs

Jatropha (Jatropha curcas) is an oil-bearing plant used for biodiesel production. Construction of its standard karyotype and identification of the euchromatin/heterochromatin distribution associated with gene expression and meiotic recombination are essential to fully characterize its genome. Here, we developed a J. curcas karyotype based on meiotic pachytene chromosomes. In addition, a karyotype of J. integerrima, a useful species for jatropha breeding, was also constructed. Five out of eleven J. curcas chromosomes were metacentric, but only two were metacentric in J. integerrima. Almost all of the heterochromatin was distributed around the pericentric regions. The interstitial and distal regions were euchromatic without heterochromatic knobs, except for small heterochromatin regions associated with the subtelomeric repeat sequence JcSat1. These pericentric heterochromatin distribution patterns, together with chromosome structure data and the results of FISH probing with rDNA and JcSat1, allowed us to classify all chromosomes of both species. The two species had two 35S rDNA loci and one 5S rDNA locus; one 35S rDNA locus in J. integerrima was located on the interstitial region of the short arms. In addition, JcSat1 was found at only the heterochromatic ends of the J. curcas chromosome, not the J. integerrima chromosome. Despite the same chromosome number, the two pachytene chromosome-based karyotypes suggest variation in chromosome structure and distribution of repetitive DNAs in these two species.

miR172 Regulates both Vegetative and Reproductive Development in the Perennial Woody Plant Jatropha curcas

Jatropha curcas is a promising feedstock for biofuel production because its oil is highly suitable for processing bio-jet fuels and biodiesel. However, Jatropha exhibits a long juvenile stage in subtropical areas. miR172, a conserved small non-protein-coding RNA molecule with 21 nucleotides, regulates a wide range of developmental processes. To date, however, no studies have examined the function of miR172 in Jatropha. There are five miR172 precursors encoding two mature miR172s in Jatropha, which are expressed in all tissues, with the highest expression level in leaves, and the levels are up-regulated with age. Overexpression of JcmiR172a resulted in early flowering, abnormal flowers, and altered leaf morphology in transgenic Arabidopsis and Jatropha. The expression levels of miR172 target genes were down-regulated, and the flower identity genes were up-regulated in the JcmiR172a-overexpressing transgenic plants. Interestingly, we showed that JcmiR172 might be involved in regulation of stem vascular development through manipulating the expression of cellulose and lignin biosynthesis genes. Overexpression of JcmiR172a enhanced xylem development and reduced phloem and pith development. This study helped elucidate the functions of miR172 in perennial plants, a known age-related miRNA involved in the regulation of perennial plant phase change and organ development.

Transcriptome Analysis Reveals Dynamic Fat Accumulation in the Walnut Kernel

Walnut (Juglans regia L.) is an important woody oilseed species cultivated throughout the world. In this study, comparative transcript profiling was performed using high-throughput RNA sequencing technology at the following three stages of walnut fat synthesis in the “Lvling” walnut cultivar: the initial developmental stage (L1), the fast developing stage (L2), and the last developing stage (L3). A total of 68.18 GB of data were obtained on the three developmental stages, and 92% to 94% of clean data were able to be located to the reference genome. Further comparisons of the transcripts in the three libraries revealed that 724, 2027, and 4817 genes were differentially expressed between the L2 and L1 (L2vsL1), L3 and L2 (L3vsL2), and L3 and L1 (L3vsL1) samples, respectively. Through the GO gene enrichment analysis, differentially expressed genes (DEGs) in L2vsL1, L3vsL2, and L3vsL1 were enriched into 3, 0, and 2 functional categories, respectively. According to the KEGG enrichment analysis, DEGs in L2vsL1, L3vsL2, and L3vsL1 were annotated into 77, 110, and 3717 taxonomic metabolic pathways in the KEGG database, respectively. Next, we analyzed expression levels of genes related to fat synthesis. Our results indicated that ACCase, LACS, and FAD7 were the key genes related to fat synthesis. The high-throughput transcriptome sequencing of walnut in different developmental stages has greatly enriched the current genomic available resources. The comparison of DEGs under different developmental stages identified a wealth of candidate genes involved in fat synthesis, which will facilitate further genetic improvement and molecular studies of the walnut.

Synthesis, Racemic X-ray Crystallographic, and Permeability Studies of Bioactive Orbitides from Jatropha Species

Orbitides are small cyclic peptides with a diverse range of therapeutic bioactivities. They are produced by many plant species, including those of the Jatropha genus. Here, the objective was to provide new structural information on orbitides to complement the growing knowledge base on orbitide sequences and activities by focusing on three Jatropha orbitides: ribifolin (1), pohlianin C (7), and jatrophidin (12). To determine three-dimensional structures, racemic crystallography, an emerging structural technique that enables rapid crystallization of biomolecules by combining equal amounts of the two enantiomers, was used. The high-resolution structure of ribifolin (0.99 Å) was elucidated from its racemate and showed it was identical to the structure crystallized from its l-enantiomer only (1.35 Å). Racemic crystallography was also used to elucidate high-resolution structures of pohlianin C (1.20 Å) and jatrophidin (1.03 Å), for which there was difficulty forming crystals without using racemic mixtures. The structures were used to interpret membrane permeability data in PAMPA and a Caco-2 cell assay, showing they had poor permeability. Overall, the results show racemic crystallography can be used to obtain high-resolution structures of orbitides and is useful when enantiopure samples are difficult to crystallize or solution structures from NMR are of low resolution.

Overview of Next-generation Sequencing Platforms Used in Published Draft Plant Genomes in Light of Genotypization of Immortelle Plant (Helichrysium Arenarium)

INTRODUCTION

Major advancements in DNA sequencing methods introduced in the first decade of the new millennium initiated a rapid expansion of sequencing studies, which yielded a tremendous amount of DNA sequence data, including whole sequenced genomes of various species, including plants. A set of novel sequencing platforms, often collectively named as "next-generation sequencing" (NGS) completely transformed the life sciences, by allowing extensive throughput, while greatly reducing the necessary time, labor and cost of any sequencing endeavor.

PURPOSE

of this paper is to present an overview NGS platforms used to produce the current compendium of published draft genomes of various plants, namely the Roche/454, ABI/SOLiD, and Solexa/Illumina, and to determine the most frequently used platform for the whole genome sequencing of plants in light of genotypization of immortelle plant.

MATERIALS AND METHODS

45 papers were selected (with 47 presented plant genome draft sequences), and utilized sequencing techniques and NGS platforms (Roche/454, ABI/SOLiD and Illumina/Solexa) in selected papers were determined. Subsequently, frequency of usage of each platform or combination of platforms was calculated.

RESULTS

Illumina/Solexa platforms are by used either as sole sequencing tool in 40.42% of published genomes, or in combination with other platforms - additional 48.94% of published genomes, followed by Roche/454 platforms, used in combination with traditional Sanger sequencing method (10.64%), and never as a sole tool. ABI/SOLiD was only used in combination with Illumina/Solexa and Roche/454 in 4.25% of publications.

CONCLUSIONS

Illumina/Solexa platforms are by far most preferred by researchers, most probably due to most affordable sequencing costs. Taking into consideration the current economic situation in the Balkans region, Illumina Solexa is the best (if not the only) platform choice if the sequencing of immortelle plant (Helichrysium arenarium) is to be performed by the researchers in this region.

GinMicrosatDb: a genome-wide microsatellite markers database for sesame (Sesamum indicum L.)

Molecular breeding in sesame is still at infancy due to limited number of microsatellite markers available and the low level of polymorphism exhibited by them. Therefore, whole genome sequencing was used for development of microsatellite markers so as to ensure availability of substantial number of polymorphic markers for use in marker assisted breeding programs. Whole genome sequencing of sesame variety 'Swetha' was done using Illumina paired-end sequencing and Roche 454 shotgun sequencing technologies (GCA_000975565.1 in GenBank). 'GinMicrosatDb', a genome-wide microsatellite marker database has been developed using the whole genome sequence data of sesame variety 'Swetha'. The database consists of microsatellites localized on both linkage groups and scaffolds with their genomic co-ordinates. It provides five sets of forward and reverse primers for each of the microsatellite loci along with the flanking sequences, primer GC content, product size and melting temperature etc. The distribution of microsatellites can be viewed and selected through a genome browser as well as through a physical map. The newly identified microsatellite markers are expected to help sesame breeders in developing marker tags for traits of economic importance thereby bringing about greater efficiency in marker-assisted selection programs.

Identification and expression analysis of cytokinin metabolic genes IPTs, CYP735A and CKXs in the biofuel plant Jatropha curcas

The seed oil of Jatropha curcas is considered a potential bioenergy source that could replace fossil fuels. However, the seed yield of Jatropha is low and has yet to be improved. We previously reported that exogenous cytokinin treatment increased the seed yield of Jatropha. Cytokinin levels are directly regulated by isopentenyl transferase (IPT), cytochrome P450 monooxygenase, family 735, subfamily A (CYP735A), and cytokinin oxidase/dehydrogenase (CKX). In this study, we cloned six IPT genes, one JcCYP735A gene, and seven JcCKX genes. The expression patterns of these 14 genes in various organs were determined using real-time quantitative PCR. JcIPT1 was primarily expressed in roots and seeds, JcIPT2 was expressed in roots, apical meristems, and mature leaves, JcIPT3 was expressed in stems and mature leaves, JcIPT5 was expressed in roots and mature leaves, JcIPT6 was expressed in seeds at 10 days after pollination, and JcIPT9 was expressed in mature leaves. JcCYP735A was mainly expressed in roots, flower buds, and seeds. The seven JcCKX genes also showed different expression patterns in different organs of Jatropha. In addition, CK levels were detected in flower buds and seeds at different stages of development. The concentration of N⁶-(Δ²-isopentenyl)-adenine (iP), iP-riboside, and trans-zeatin (tZ) increased with flower development, and the concentration of iP decreased with seed development, while that of tZ increased. We further analyzed the function of JcCYP735A using the CRISPR-Cas9 system, and found that the concentrations of tZ and tZ-riboside decreased significantly in the Jccyp735a mutants, which showed severely retarded growth. These findings will be helpful for further studies of the functions of cytokinin metabolic genes and understanding the roles of cytokinins in Jatropha growth and development.

Karyotype and genome size comparative analyses among six species of the oilseed-bearing genus Jatropha (Euphorbiaceae)

Jatropha is an important genus of Euphorbiaceae, with species largely used for various purposes, including the manufacturing of soaps and pharmaceutical products and applications in the bioenergetic industry. Although there have been several studies focusing J. curcas in various aspects, the karyotype features of Jatropha species are poorly known. Therefore, we analyzed six Jatropha species through fluorochrome staining (CMA/DAPI), fluorescent in situ hybridization (FISH) with 5S and 45S rDNA probes and genome size estimation by flow cytometry. Our results revealed several chromosome markers by both CMA/DAPI and FISH for the analyzed species. Five Jatropha species (J. curcas, J. gossypiifolia, J. integerrima, J. multifida and J. podagrica) showed four CMA-positive (CMA+) bands associated with the 5S and 45S rDNA sites (one and two pairs, respectively). However, J. mollissima displayed six CMA+/DAPI- bands co-localized with both 5S and 45S rDNA, which showed a FISH superposition. A gradual variation in the genome sizes was observed (2C = 0.64 to 0.86 pg), although an association between evidenced heterochromatin and genome sizes was not found among species. Except for the unique banding pattern of J. mollissima and the pericentromeric heterochromatin of J. curcas and J. podagrica, our data evidenced relatively conserved karyotypes.

The Pattern and Distribution of Induced Mutations in J. curcas Using Reduced Representation Sequencing

Mutagenesis in combination with Genotyping by Sequencing (GBS) is a powerful tool for introducing variation, studying gene function and identifying causal mutations underlying phenotypes of interest in crop plant genomes. About 400 million paired-end reads were obtained from 82 ethylmethane sulfonate (EMS) induced mutants and 14 wild-type accessions of Jatropha curcas for the detection of Single Nucleotide Polymorphisms (SNPs) and Insertion/Deletions (InDels) by two different approaches (nGBS and ddGBS) on an Illumina HiSeq 2000 sequencer. Using bioinformatics analyses, 1,452 induced SNPs and InDels were identified in coding regions, which were distributed across 995 genes. The predominantly observed mutations were G/C to A/T transitions (64%), while transversions were observed at a lower frequency (36%). Regarding the effect of mutations on gene function, 18% of the mutations were located in intergenic regions. In fact, mutants with the highest number of heterozygous SNPs were found in samples treated with 0.8% EMS for 3 h. Reconstruction of the metabolic pathways showed that in total 16 SNPs were located in six KEGG pathways by nGBS and two pathways by ddGBS. The most highly represented pathways were ether-lipid metabolism and glycerophospholipid metabolism, followed by starch and sucrose metabolism by nGBS and triterpenoid biosynthesis as well as steroid biosynthesis by ddGBS. Furthermore, high genome methylation was observed in J. curcas, which might help to understand the plasticity of the Jatropha genome in response to environmental factors. At last, the results showed that continuously vegetatively propagated tissue is a fast, efficient and accurate method to dissolve chimeras, especially for long-lived plants like J. curcas. Obtained data showed that allelic variations and in silico analyses of gene functions (gene function prediction), which control important traits, could be identified in mutant populations using nGBS and ddGBS. However, the handling of GBS data is more difficult and more challenging than the traditional TILLING strategy in mutated plants, since the Jatropha genome sequence is incomplete, which makes alignment and variant analysis of target sequence reads challenging to perform and interpret. Therefore, providing a complete Jatropha reference genome sequence with high quality should be a priority for any breeding program.

Current Strategies for the Detoxification of Jatropha curcas Seed Cake: A Review

Jatropha curcas is an important oilseed plant, with considerable potential in the development of biodiesel. Although Jatropha seed cake, the byproduct of oil extraction, is a residue rich in nitrogen, phosphorus, potassium, and carbon, with high protein content suitable for application in animal feed, the presence of toxic phorbol esters limits its application in feed supplements and fertilizers. This review summarizes the current methods available for detoxification of this residue, based upon chemical, physical, biological, or combined processes. The advantages and disadvantages of each process are discussed, and future directions involving genomic and proteomic approaches for advancing our understanding of biodegradation processes involving microorganisms are highlighted.

Genome-wide comparative analysis of papain-like cysteine protease family genes in castor bean and physic nut

Papain-like cysteine proteases (PLCPs) are a class of proteolytic enzymes involved in many plant processes. Compared with the extensive research in Arabidopsis thaliana, little is known in castor bean (Ricinus communis) and physic nut (Jatropha curcas), two Euphorbiaceous plants without any recent whole-genome duplication. In this study, a total of 26 or 23 PLCP genes were identified from the genomes of castor bean and physic nut respectively, which can be divided into nine subfamilies based on the phylogenetic analysis: RD21, CEP, XCP, XBCP3, THI, SAG12, RD19, ALP and CTB. Although most of them harbor orthologs in Arabidopsis, several members in subfamilies RD21, CEP, XBCP3 and SAG12 form new groups or subgroups as observed in other species, suggesting specific gene loss occurred in Arabidopsis. Recent gene duplicates were also identified in these two species, but they are limited to the SAG12 subfamily and were all derived from local duplication. Expression profiling revealed diverse patterns of different family members over various tissues. Furthermore, the evolution characteristics of PLCP genes were also compared and discussed. Our findings provide a useful reference to characterize PLCP genes and investigate the family evolution in Euphorbiaceae and species beyond.

Construction of an ultrahigh-density genetic linkage map for Jatropha curcas L. and identification of QTL for fruit yield

BACKGROUND

As an important biofuel plant, the demand for higher yield Jatropha curcas L. is rapidly increasing. However, genetic analysis of Jatropha and molecular breeding for higher yield have been hampered by the limited number of molecular markers available.

RESULTS

An ultrahigh-density linkage map for a Jatropha mapping population of 153 individuals was constructed and covered 1380.58 cM of the Jatropha genome, with average marker density of 0.403 cM. The genetic linkage map consisted of 3422 SNP and indel markers, which clustered into 11 linkage groups. With this map, 13 repeatable QTLs (reQTLs) for fruit yield traits were identified. Ten reQTLs, qNF-1, qNF-2a, qNF-2b, qNF-2c, qNF-3, qNF-4, qNF-6, qNF-7a, qNF-7b and qNF-8, that control the number of fruits (NF) mapped to LGs 1, 2, 3, 4, 6, 7 and 8, whereas three reQTLs, qTWF-1, qTWF-2 and qTWF-3, that control the total weight of fruits (TWF) mapped to LGs 1, 2 and 3, respectively. It is interesting that there are two candidate critical genes, which may regulate Jatropha fruit yield. We also identified three pleiotropic reQTL pairs associated with both the NF and TWF traits.

CONCLUSION

This study is the first to report an ultrahigh-density Jatropha genetic linkage map construction, and the markers used in this study showed great potential for QTL mapping. Thirteen fruit-yield reQTLs and two important candidate genes were identified based on this linkage map. This genetic linkage map will be a useful tool for the localization of other economically important QTLs and candidate genes for Jatropha.

Current Strategies of Polyploid Plant Genome Sequence Assembly

Polyploidy or duplication of an entire genome occurs in the majority of angiosperms. The understanding of polyploid genomes is important for the improvement of those crops, which humans rely on for sustenance and basic nutrition. As climate change continues to pose a potential threat to agricultural production, there will increasingly be a demand for plant cultivars that can resist biotic and abiotic stresses and also provide needed and improved nutrition. In the past decade, Next Generation Sequencing (NGS) has fundamentally changed the genomics landscape by providing tools for the exploration of polyploid genomes. Here, we review the challenges of the assembly of polyploid plant genomes, and also present recent advances in genomic resources and functional tools in molecular genetics and breeding. As genomes of diploid and less heterozygous progenitor species are increasingly available, we discuss the lack of complexity of these currently available reference genomes as they relate to polyploid crops. Finally, we review recent approaches of haplotyping by phasing and the impact of third generation technologies on polyploid plant genome assembly.

Genome-wide analysis of SSR and ILP markers in trees: diversity profiling, alternate distribution, and applications in duplication

Molecular markers are efficient tools for breeding and genetic studies. However, despite their ecological and economic importance, their development and application have long been hampered. In this study, we identified 524,170 simple sequence repeat (SSR), 267,636 intron length polymorphism (ILP), and 11,872 potential intron polymorphism (PIP) markers from 16 tree species based on recently available genome sequences. Larger motifs, including hexamers and heptamers, accounted for most of the seven different types of SSR loci. Within these loci, A/T bases comprised a significantly larger proportion of sequence than G/C. SSR and ILP markers exhibited an alternative distribution pattern. Most SSRs were monomorphic markers, and the proportions of polymorphic markers were positively correlated with genome size. By verifying with all 16 tree species, 54 SSR, 418 ILP, and four PIP universal markers were obtained, and their efficiency was examined by PCR. A combination of five SSR and six ILP markers were used for the phylogenetic analysis of 30 willow samples, revealing a positive correlation between genetic diversity and geographic distance. We also found that SSRs can be used as tools for duplication analysis. Our findings provide important foundations for the development of breeding and genetic studies in tree species.