The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis.

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

Chromosome-level genome of Camellia lanceoleosa provides a valuable resource for understanding genome evolution and self-incompatibility.

The section Oleifera (Theaceae) has attracted attention for the high levels of unsaturated fatty acids found in its seeds. Here, we report the chromosome-scale genome of the sect. Oleifera using diploid wild Camellia lanceoleosa with a final size of 3.00 Gb and an N50 scaffold size of 186.43 Mb. Repetitive sequences accounted for 80.63% and were distributed unevenly across the genome. Camellia lanceoleosa underwent a whole-genome duplication event approximately 65 million years ago (65 Mya), prior to the divergence of C. lanceoleosa and Camellia sinensis (approx. 6-7 Mya). Syntenic comparisons of these two species elucidated the genomic rearrangement, appearing to be driven in part by the activity of transposable elements. The expanded and positively selected genes in C. lanceoleosa were significantly enriched in oil biosynthesis, and the expansion of homomeric acetyl-coenzyme A carboxylase (ACCase) genes and the seed-biased expression of genes encoding heteromeric ACCase, diacylglycerol acyltransferase, glyceraldehyde-3-phosphate dehydrogenase and stearoyl-ACP desaturase could be of primary importance for the high oil and oleic acid content found in C. lanceoleosa. Theanine and catechins were present in the leaves of C. lanceoleosa. However, caffeine can not be dectected in the leaves but was abundant in the seeds and roots. The functional and transcriptional divergence of genes encoding SAM-dependent N-methyltransferases may be associated with caffeine accumulation and distribution. Gene expression profiles, structural composition and chromosomal location suggest that the late-acting self-incompatibility of C. lanceoleosa is likely to have favoured a novel mechanism co-occurring with gametophytic self-incompatibility. This study provides valuable resources for quantitative and qualitative improvements and genome assembly of polyploid plants in sect. Oleifera.

The complex genome and adaptive evolution of polyploid Chinese pepper (Zanthoxylum armatum and Zanthoxylum bungeanum).

Zanthoxylum armatum and Zanthoxylum bungeanum, known as 'Chinese pepper', are distinguished by their extraordinary complex genomes, phenotypic innovation of adaptive evolution and species-special metabolites. Here, we report reference-grade genomes of Z. armatum and Z. bungeanum. Using high coverage sequence data and comprehensive assembly strategies, we derived 66 pseudochromosomes comprising 33 homologous phased groups of two subgenomes, including autotetraploid Z. armatum. The genomic rearrangements and two whole-genome duplications created large (~4.5 Gb) complex genomes with a high ratio of repetitive sequences (>82%) and high chromosome number (2n = 4x = 132). Further analysis of the high-quality genomes shed lights on the genomic basis of involutional reproduction, allomones biosynthesis and adaptive evolution in Chinese pepper, revealing a high consistent relationship between genomic evolution, environmental factors and phenotypic innovation. Our study provides genomic resources and new insights for investigating diversification and phenotypic innovation in Chinese pepper, with broader implications for the protection of plants under severe environmental changes.

Fine-Mapping of the Powdery Mildew Resistance Gene mlxbd in the Common Wheat Landrace Xiaobaidong.

Powdery mildew, which is caused by Blumeria graminis f. sp. tritici (Bgt), is a disease of wheat worldwide. Xiaobaidong is a Chinese wheat landrace, which still maintains good resistance against powdery mildew. To obtain more genetic markers closely linked to the powdery mildew resistance gene mlxbd and narrow the candidate region for its isolation, new simple sequence repeats and cross intron-spanning markers were designed based on the genome sequence of Triticum aestivum cultivar Chinese Spring chromosome 7BL. The flanking markers 7BLSSR49 and WGGC5746 were found to be tightly linked to mlxbd at genetic distances of 0.4 cM and 0.3 cM, respectively. The resistance locus was mapped to a 63.40 kb and 0.29 Mb region of the Chinese Spring genome and Zavitan genome, respectively. The linked markers of mlxbd could be used as diagnostic markers for mlxbd. The linked molecular markers and delineated genomic region in the sequenced Chinese Spring genome will assist the future map-based cloning of mlxbd.

iTRAQ-based comparative proteomic analysis provides insights into somatic embryogenesis in Gossypium hirsutum L.

KEY MESSAGE: iTRAQ based proteomic identified key proteins and provided new insights into the molecular mechanisms underlying somatic embryogenesis in cotton. Somatic embryogenesis, which involves cell dedifferentiation and redifferentiation, has been used as a model system for understanding molecular events of plant embryo development in vitro. In this study, we performed comparative proteomics analysis using samples of non-embryogenic callus (NEC), embryogenic callus (EC) and somatic embryo (SE) using the isobaric tags for relative and absolute quantitation (iTRAQ) technology. In total, 5892 proteins were identified amongst the three samples. The majority of these proteins (93.4%) were found to have catalytic activity, binding activity, transporter activity or structural molecular activity. Of these proteins, 1024 and 858 were differentially expressed in NEC versus EC and EC versus SE, respectively. Compared to NEC, EC had 452 and 572 down- and up-regulated proteins, respectively, and compared to EC, SE had 647 and 221 down- and up-regulated proteins, respectively. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis indicated that genetic information transmission, plant hormone transduction, glycolysis, fatty acid biosynthesis and metabolism, galactose metabolism were the top pathways involved in somatic embryogenesis. Our proteomics results not only confirmed our previous transcriptomic results on the role of the polyamine metabolic pathways and stress responses in cotton somatic embryogenesis, but identified key proteins important for cotton somatic embryogenesis and provided new insights into the molecular mechanisms underlying somatic embryogenesis in cotton.

De novo transcriptome analysis reveals insights into dynamic homeostasis regulation of somatic embryogenesis in upland cotton (G. hirsutum L.).

Plant regeneration via somatic embryogenesis (SE) is the key step for genetic improvement of cotton (Gossypium hirsutum L.) through genetic engineering mediated by Agrobacteria, but the molecular mechanisms underlying SE in cotton is still unclear. Here, RNA-Sequencing was used to analyze the genes expressed during SE and their expression dynamics using RNAs isolated from non-embryogenic callus (NEC), embryogenic callus (EC) and somatic embryos (SEs). A total of 101, 670 unigenes were de novo assembled. The genes differentially expressed (DEGs) amongst NEC, EC and SEs were identified, annotated and classified. More DEGs were found between SEs and EC than between EC and NEC. A significant number of DEGs were related to hormone homeostasis, stress and ROS responses, and metabolism of polyamines. To confirm the expression dynamics of selected DEGs involved in various pathways, experiments were set up to investigate the effects of hormones (Indole-3-butytric acid, IBA; Kinetin, KT), polyamines, H2O2 and stresses on SE. Our results showed that exogenous application of IBA and KT positively regulated the development of EC and SEs, and that polyamines and H2O2 promoted the conversion of EC into SEs. Furthermore, we found that low and moderate stress is beneficial for proliferation of EC and SEs formation. Together, our global analysis of transcriptomic dynamics reveals that hormone homeostasis, polyamines, and stress response synergistically regulating SE in cotton.

Construction of Host Plant Insect-Resistance Mutant Library by High-Throughput CRISPR/Cas9 System and Identification of A Broad-Spectrum Insect Resistance Gene.

Insects pose significant challenges in cotton-producing regions. Here, they describe a high-throughput CRISPR/Cas9-mediated large-scale mutagenesis library targeting endogenous insect-resistance-related genes in cotton. This library targeted 502 previously identified genes using 968 sgRNAs, generated ≈2000 T0 plants and achieved 97.29% genome editing with efficient heredity, reaching upto 84.78%. Several potential resistance-related mutants (10% of 200 lines) their identified that may contribute to cotton-insect molecular interaction. Among these, they selected 139 and 144 lines showing decreased resistance to pest infestation and targeting major latex-like protein 423 (GhMLP423) for in-depth study. Overexpression of GhMLP423 enhanced insect resistance by activating the plant systemic acquired resistance (SAR) of salicylic acid (SA) and pathogenesis-related (PR) genes. This activation is induced by an elevation of cytosolic calcium [Ca2+ ]cyt flux eliciting reactive oxygen species (ROS), which their demoted in GhMLP423 knockout (CR) plants. Protein-protein interaction assays revealed that GhMLP423 interacted with a human epidermal growth factor receptor substrate15 (EPS15) protein at the cell membrane. Together, they regulated the systemically propagating waves of Ca2+ and ROS, which in turn induced SAR. Collectively, this large-scale mutagenesis library provides an efficient strategy for functional genomics research of polyploid plant species and serves as a solid platform for genetic engineering of insect resistance.

Transcriptome analysis provides insights into the molecular mechanism of GhSAMDC1 involving in rapid vegetative growth and early flowering in tobacco.

In previous study, ectopic expression of GhSAMDC1 improved vegetative growth and early flowering in tobacco, which had been explained through changes of polyamine content, polyamines and flowering relate genes expression. To further disclose the transcript changes of ectopic expression of GhSAMDC1 in tobacco, the leaves from wild type and two transgenic lines at seedling (30 days old), bolting (60 days old) and flowering (90 days old) stages were performed for transcriptome analysis. Compared to wild type, a total of 938 differentially expressed genes (DEGs) were found to be up- or down-regulated in the two transgenic plants. GO and KEGG analysis revealed that tobacco of wild-type and transgenic lines were controlled by a complex gene network, which regulated multiple metabolic pathways. Phytohormone detection indicate GhSAMDC1 affect endogenous phytohormone content, ABA and JA content are remarkably increased in transgenic plants. Furthermore, transcript factor analysis indicated 18 transcript factor families, including stress response, development and flowering related transcript factor families, especially AP2-EREBP, WRKY, HSF and Tify are the most over-represented in those transcript factor families. In conclusion, transcriptome analysis provides insights into the molecular mechanism of GhSAMDC1 involving rapid vegetative growth and early flowering in tobacco.

Identification and Functional Analysis of lncRNA by CRISPR/Cas9 During the Cotton Response to Sap-Sucking Insect Infestation.

Sap-sucking insects cause severe damage to cotton production. Long non-coding RNAs (lncRNAs) play vital regulatory roles in various development processes and stress response, however, the function of lncRNAs during sap-sucking insect infection in cotton is largely unknown. In this study, the transcriptome profiles between resistant (HR) and susceptible (ZS) cotton cultivars under whitefly infestation at different time points (0, 4, 12, 24, and 48 h) were compared. A total of 6,651 lncRNAs transcript and 606 differentially expressed lncRNAs were identified from the RNA-seq data. A co-expression network indicated that lncA07 and lncD09 were potential hub genes that play a regulatory role in cotton defense against aphid infestation. Furthermore, CRISPR/Cas9 knock-out mutant of lncD09 and lncA07 showed a decrease of jasmonic acid (JA) content, which potentially lead to increased susceptibility toward insect infestation. Differentially expressed genes between wild type and lncRNA knock-out plants are enriched in modulating development and resistance to stimulus. Additionally, some candidate genes such as Ghir_A01G022270, Ghir_D04G014430, and Ghir_A01G022270 are involved in the regulation of the JA-mediated signaling pathway. This result provides a novel insight of the lncRNA role in the cotton defense system against pests.

Ectopic Expression of GhSAMDC3 Enhanced Salt Tolerance Due to Accumulated Spd Content and Activation of Salt Tolerance-Related Genes in Arabidopsis thaliana.

Polyamines (PAs), especially spermidine and spermine (which are involved in various types of abiotic stress tolerance), have been reported in many plant species. In this study, we identified 14 putative S-adenosylmethionine decarboxylase genes (GhSAMDC1-14) in upland cotton. Based on phylogenetic and expression analyses conducted under different abiotic stresses, we selected and transferred GhSAMDC3 into Arabidopsis thaliana. Compared to the wild type, transgenic plants displayed rapid growth and increases in average leaf area and leaf number of 52% and 36%, respectively. In transgenic plants, the germination vigor and rate were markedly enhanced under NaCl treatment, and the plant survival rate increased by 50% under 300 mM NaCl treatment. The spermidine content was significantly increased, possibly due to the synthesis of a series of PAs and oxidant and antioxidant genes, resulting in improved salinity tolerance in Arabidopsis. Various salinity resistance-related genes were upregulated in transgenic plants. Together, these results indicate that ectopic expression of GhSAMDC3 raised salinity tolerance by the accumulation of spermidine and activation of salinity tolerance-related genes in A. thaliana.

Ectopic expression of GhSAMDC1 improved plant vegetative growth and early flowering through conversion of spermidine to spermine in tobacco.

Polyamines play essential roles in plant development and various stress responses. In this study, one of the cotton S-adenosylmethionine decarboxylase (SAMDC) genes, GhSAMDC1, was constructed in the pGWB17 vector and overexpressed in tobacco. Leaf area and plant height increased 25.9-36.6% and 15.0-27.0%, respectively, compared to the wild type, and flowering time was advanced by 5 days in transgenic tobacco lines. Polyamine and gene expression analyses demonstrated that a decrease in spermidine and an increase in total polyamines and spermine might be regulated by NtSPDS4 and NtSPMS in transgenic plants. Furthermore, exogenous spermidine, spermine and spermidine synthesis inhibitor dicyclohexylamine were used for complementary tests, which resulted in small leaves and dwarf plants, big leaves and early flowering, and big leaves and dwarf plants, respectively. These results indicate that spermidine and spermine are mainly involved in the vegetative growth and early flowering stages, respectively. Expression analysis of flowering-related genes suggested that NtSOC1, NtAP1, NtNFL1 and NtFT4 were upregulated in transgenic plants. In conclusion, ectopic GhSAMDC1 is involved in the conversion of spermidine to spermine, resulting in rapid vegetative growth and early flowering in tobacco, which could be applied to genetically improve plants.

Simultaneous silencing of GhFAD2-1 and GhFATB enhances the quality of cottonseed oil with high oleic acid.

Cottonseed oil has become an important source of edible oil due to its significant cost advantage. However, there is a growing concern over its fatty acid composition and nutritional value. In Gossypium hirsutum, GhFAD2-1 and GhFATB encoding the microsomal oleate desaturase and palmitoyl-acyl carrier protein thioesterase, respectively, play critical roles in regulating the proportions of saturated and polyunsaturated fatty acids in cottonseed lipids. In this study, RNAi technology was used to simultaneously inhibit the expression levels of GhFAD2-1 and GhFATB to improve the quality of cottonseed oil by increasing oleic acid content. Transgenic cotton plants with reduced levels of both target genes were successfully generated. In mature seed kernels of transgenic plants, the content of oleic acid was 38.25%, accordingly increasing by 156.96%, while the content of palmitic acid and linoleic acid was 19.15% and 36.68%, decreasing by 21.28% and 33.92%, respectively, compared with that of the control. The total oil content in transgenic and control kernels was 22.48% and 29.83%, respectively. The reduced oil level in transgenic seeds was accompanied by a reduction in seed index, thereby causing disadvantageous effects on seed germination potentiality and seed vigor, particularly under cool stress conditions. Our results demonstrated the feasibility of simultaneous manipulation of multiple genes using RNAi technology and showed the important role of oil content in seed development and vigor. Our findings provide insight into the physiological significance of the fatty acid composition in cottonseeds.

Polyamine and Its Metabolite H2O2 Play a Key Role in the Conversion of Embryogenic Callus into Somatic Embryos in Upland Cotton (Gossypium hirsutum L.).

The objective of this study was to increase understanding about the mechanism by which polyamines (PAs) promote the conversion of embryogenic calli (EC) into somatic embryos in cotton (Gossypium hirsutum L.). We measured the levels of endogenous PAs and H2O2, quantified the expression levels of genes involved in the PAs pathway at various stages of cotton somatic embryogenesis (SE), and investigated the effects of exogenous PAs and H2O2 on differentiation and development of EC. Putrescine (Put), spermidine (Spd), and spermine (Spm) significantly increased from the EC stage to the early phase of embryo differentiation. The levels of Put then decreased until the somatic embryo stage whereas Spd and Spm remained nearly the same. The expression profiles of GhADC genes were consistent with changes in Put during cotton SE. The H2O2 concentrations began to increase significantly at the EC stage, during which time both GhPAO1 and GhPAO4 expressions were highest and PAO activity was significantly increased. Exogenous Put, Spd, Spm, and H2O2 not only enhanced embryogenic callus growth and embryo formation, but also alleviated the effects of D-arginine and 1, 8-diamino-octane, which are inhibitors of PA synthesis and PAO activity. Overall, the results suggest that both PAs and their metabolic product H2O2 are essential for the conversion of EC into somatic embryos in cotton.

Detection of somaclonal variation of cotton (Gossypium hirsutum) using cytogenetics, flow cytometry and molecular markers.

Two protocols of plant regeneration for cotton were adopted in this study, namely, 2, 4-D and kinetin hormone combination and IBA and kinetin hormone combination. Twenty-eight embryogenic cell lines via somatic embryogenesis and 67 regenerated plants from these embryogenic calli were selected and used for random amplified polymorphic DNA (RAPD), simple sequence repeat (SSR), chromosomal number counting, and flow cytometric analysis. The roles of RAPD and SSR markers in detecting somaclonal variation of cotton (Gossypium hirsutum L.) were evaluated. Two cluster analyses were performed to express, in the form of dendrograms, the relationships among the hormone combinations and the genetic variability. Both DNA-based techniques were able to amplify all of the cell clones and regenerated plantlets genomes and relative higher genetic variation could be detected in the culture type with 2, 4-D and kinetin hormone combination. The result suggested that 2, 4-D and kinetin hormone combination could induce relative high somaclonal variation and RAPD and SSR markers are useful in detecting somaclonal variation of regenerated cotton plants via somatic embryogenesis. Chromosome number counting and flow cytometry analysis revealed that the number of chromosomes and ploidy levels were nearly stable in all regenerated plants except two regenerated plantlets (lost 4 and 5 chromosomes, respectively) which meant that cytological changes were not correlated with the frequency of RAPD and SSR polymorphisms. This result also might mean that the cell lines with variation of chromosome numbers were difficult to regenerate plants.