A draft genome assembly of halophyte Suaeda aralocaspica, a plant that performs C4 photosynthesis within individual cells

Understanding the salinity resilience and productivity of halophytes in saline environments

The escalating salinity levels in cultivable soil pose a significant threat to agricultural productivity and, consequently, human sustenance. This problem is being exacerbated by natural processes and human activities, coinciding with a period of rapid population growth. Developing halophytic crops is needed to ensure food security is not impaired and land resources can be used sustainably. Evolution has created many close halophyte relatives of our major glycophytic crops, such as Puccinellia tenuiflora (relative of barley and wheat), Oryza coarctata (relative of rice) and Glycine soja (relative of soybean). There are also some halophytes have been subjected to semi-domestication and are considered as minor crops, such as Chenopodium quinoa. In this paper, we examine the prevailing comprehension of robust salinity resilience in halophytes. We summarize the existing strategies and technologies that equip researchers with the means to enhance the salt tolerance capabilities of primary crops and investigate the genetic makeup of halophytes.

Combined genome and transcriptome provides insight into the genetic evolution of an edible halophyte Suaeda salsa adaptation to high salinity

Suaeda salsa L. is a typical halophyte with high value as a vegetable. Here, we report a 447.98 Mb, chromosomal-level genome of S. salsa, assembled into nine pseudomolecules (contig N50 = 1.36 Mb) and annotated with 27,927 annotated protein-coding genes. Most of the assembled S. salsa genome, 58.03%, consists of transposable elements. Some gene families including HKT1, NHX, SOS and CASP related to salt resistance were significantly amplified. We also observed expansion of genes encoding protein that bind the trace elements Zn, Fe, Cu and Mn, and genes related to flavonoid and α-linolenic acid metabolism. Many expanded genes were significantly up-regulated under salinity, which might have contributed to the acquisition of salt tolerance in S. salsa. Transcriptomic data showed that high salinity markedly up-regulated salt-resistance related genes, compared to low salinity. Abundant metabolic pathways of secondary metabolites including flavonoid, unsaturated fatty acids and selenocompound were enriched, which indicates that the species is a nutrient-rich vegetable. Particularly worth mentioning is that there was no significant difference in the numbers of cis-elements in the promoters of salt-related and randomly selected genes in S. salsa when compared with Arabidopsis thaliana, which may affirm that plant salt tolerance is a quantitative rather than a qualitative trait in terms of promoter evolution. Our findings provide deep insight into the adaptation of halophytes to salinity from a genetic evolution perspective.

Identification and selection of reference genes for analysis of gene expression by quantitative real-time PCR in the euhalophyte Suaeda altissima (L.) Pall

Сoding sequences of seven housekeeping genes: actin SaACT7, ubiquitin-conjugating protein SaUBC10, glyceraldehyde-3-phosphate dehydrogenase SaGAPDH, protein of the large subunit of ribosomes SaL2, α-tubulin SaTUA, translation elongation factor SaeEF1α, and protein phosphatase SaPP2A were identified as candidate reference genes for expression analysis of target genes in the extremely salt tolerant plant Suaeda altissima (L.) Pall. The expression profiles of the genes differed. SaACT7 and SaeEF1α demonstrated the highest expression levels, while the lowest expression levels were found for SaPP2A and SaTUA. SaPP2A and SaeEF1α genes were the most stably expressed at different steady-state salinity levels and different nitrate concentrations in nutrient solutions (NSs). SaL2, SaPP2A, and SaeEF1α genes showed the greatest stability of expression when nitrate was added to nutrient solution of plants grown under conditions of nitrate deficiency. Less constant expression was demonstrated in this experiment by SaACT7 and SaTUA. SaL2, SaACT7, SaeEF1α, and SaUBC10 genes showed the smallest expression changes under salt shock. To validate the use of the most stably expressed genes for normalization of gene expression, we checked them as reference genes to study the expression of the nitrate transporter gene SaNPF6.3 in S. altissima roots under conditions of different salinity and different nitrate supply.

Genome-Wide Comparative Analysis of Five Amaranthaceae Species Reveals a Large Amount of Repeat Content

Amaranthus is a genus of C4 dicotyledonous herbaceous plant species that are widely distributed in Asia, Africa, Australia, and Europe and are used as grain, vegetables, forages, and ornamental plants. Amaranth species have gained significant attention nowadays as potential sources of nutritious food and industrial products. In this study, we performed a comparative genome analysis of five amaranth species, namely, Amaranthus hypochondriacus, Amaranthus tuberculatus, Amaranthus hybridus, Amaranthus palmeri, and Amaranthus cruentus. The estimated repeat content ranged from 54.49% to 63.26% and was not correlated with the genome sizes. Out of the predicted repeat classes, the majority of repetitive sequences were Long Terminal Repeat (LTR) elements, which account for about 13.91% to 24.89% of all amaranth genomes. Phylogenetic analysis based on 406 single-copy orthologous genes revealed that A. hypochondriacus is most closely linked to A. hybridus and distantly related to A. cruentus. However, dioecious amaranth species, such as A. tuberculatus and A. palmeri, which belong to the subgenera Amaranthus Acnida, have formed their distinct clade. The comparative analysis of genomic data of amaranth species will be useful to identify and characterize agronomically important genes and their mechanisms of action. This will facilitate genomics-based, evolutionary studies, and breeding strategies to design faster, more precise, and predictable crop improvement programs.

Chromosome-scale genome sequence of Suaeda glauca sheds light on salt stress tolerance in halophytes

Soil salinity is a growing concern for global crop production and the sustainable development of humanity. Therefore, it is crucial to comprehend salt tolerance mechanisms and identify salt-tolerance genes to enhance crop tolerance to salt stress. Suaeda glauca, a halophyte species well adapted to the seawater environment, possesses a unique ability to absorb and retain high salt concentrations within its cells, particularly in its leaves, suggesting the presence of a distinct mechanism for salt tolerance. In this study, we performed de novo sequencing of the S. glauca genome. The genome has a size of 1.02 Gb (consisting of two sets of haplotypes) and contains 54 761 annotated genes, including alleles and repeats. Comparative genomic analysis revealed a strong synteny between the genomes of S. glauca and Beta vulgaris. Of the S. glauca genome, 70.56% comprises repeat sequences, with retroelements being the most abundant. Leveraging the allele-aware assembly of the S. glauca genome, we investigated genome-wide allele-specific expression in the analyzed samples. The results indicated that the diversity in promoter sequences might contribute to consistent allele-specific expression. Moreover, a systematic analysis of the ABCE gene families shed light on the formation of S. glauca's flower morphology, suggesting that dysfunction of A-class genes is responsible for the absence of petals in S. glauca. Gene family expansion analysis demonstrated significant enrichment of Gene Ontology (GO) terms associated with DNA repair, chromosome stability, DNA demethylation, cation binding, and red/far-red light signaling pathways in the co-expanded gene families of S. glauca and S. aralocaspica, in comparison with glycophytic species within the chenopodium family. Time-course transcriptome analysis under salt treatments revealed detailed responses of S. glauca to salt tolerance, and the enrichment of the transition-upregulated genes in the leaves associated with DNA repair and chromosome stability, lipid biosynthetic process, and isoprenoid metabolic process. Additionally, genome-wide analysis of transcription factors indicated a significant expansion of FAR1 gene family. However, further investigation is needed to determine the exact role of the FAR1 gene family in salt tolerance in S. glauca.

In-silico, evolutionary, and functional analysis of CHUP1 and its related proteins in Bienertia sinuspersici-a comparative study across C3, C4, CAM, and SCC4 model plants

Single-cell C₄ (SCC₄) plants with bienertioid anatomy carry out photosynthesis in a single cell. Chloroplast movement is the underlying phenomenon, where chloroplast unusual positioning 1 (CHUP1) plays a key role. This study aimed to characterize CHUP1 and CHUP1-like proteins in an SCC₄ photosynthetic plant, Bienertia sinuspersici. Also, a comparative analysis of SCC4 CHUP1 was made with C₃, C₄, and CAM model plants including an extant basal angiosperm, Amborella. The CHUP1 gene exists as a single copy from the basal angiosperms to SCC₄ plants. Our analysis identified that Chenopodium quinoa, a recently duplicated allotetraploid, has two copies of CHUP1. In addition, the numbers of CHUP1-like and its associated proteins such as CHUP1-like_a, CHUP1-like_b, HPR, TPR, and ABP varied between the species. Hidden Markov Model analysis showed that the gene size of CHUP1-like_a and CHUP1-like_b of SCC₄ species, Bienertia, and Suaeda were enlarged than other plants. Also, we identified that CHUP1-like_a and CHUP1-like_b are absent in Arabidopsis and Amborella, respectively. Motif analysis identified several conserved and variable motifs based on the orders (monocot and dicot) as well as photosynthetic pathways. For instance, CAM plants such as pineapple and cactus shared certain motifs of CHUP1-like_a irrespective of their distant phylogenetic relationship. The free ratio model showed that CHUP1 maintained purifying selection, whereas CHUP1-like_a and CHUP1-like_b have adaptive functions between SCC₄ plants and quinoa. Similarly, rice and maize branches displayed functional diversification on CHUP1-like_b. Relative gene expression data showed that during the subcellular compartmentalization process of Bienertia, CHUP1 and actin-binding proteins (ABP) genes showed a similar pattern of expression. Altogether, the results of this study provide insight into the evolutionary and functional details of CHUP1 and its associated proteins in the development of the SCC₄ system in comparison with other C₃, C₄, and CAM model plants.

Comparative transcriptomics reveals the role of altered energy metabolism in the establishment of single-cell C4 photosynthesis in Bienertia sinuspersici

Single-cell C₄ photosynthesis (SCC₄) in terrestrial plants without Kranz anatomy involves three steps: initial CO₂ fixation in the cytosol, CO₂ release in mitochondria, and a second CO₂ fixation in central chloroplasts. Here, we investigated how the large number of mechanisms underlying these processes, which occur in three different compartments, are orchestrated in a coordinated manner to establish the C₄ pathway in Bienertia sinuspersici, a SCC₄ plant. Leaves were subjected to transcriptome analysis at three different developmental stages. Functional enrichment analysis revealed that SCC₄ cycle genes are coexpressed with genes regulating cyclic electron flow and amino/organic acid metabolism, two key processes required for the production of energy molecules in C₃ plants. Comparative gene expression profiling of B. sinuspersici and three other species (Suaeda aralocaspica, Amaranthus hypochondriacus, and Arabidopsis thaliana) showed that the direction of metabolic flux was determined via an alteration in energy supply in peripheral chloroplasts and mitochondria via regulation of gene expression in the direction of the C₄ cycle. Based on these results, we propose that the redox homeostasis of energy molecules via energy metabolism regulation is key to the establishment of the SCC₄ pathway in B. sinuspersici.

Genome-Wide Development of Polymorphic Microsatellite Markers and Genetic Diversity Analysis for the Halophyte Suaeda aralocaspica (Amaranthaceae)

Suaeda aralocaspica, which is an annual halophyte, grows in saline deserts in Central Asia with potential use in saline soil reclamation and salt tolerance breeding. Studying its genetic diversity is critical for effective conservation and breeding programs. In this study, we aimed to develop a set of polymorphic microsatellite markers to analyze the genetic diversity of S. aralocaspica. We identified 177,805 SSRs from the S. aralocaspica genome, with an average length of 19.49 bp, which were present at a density of 393.37 SSR/Mb. Trinucleotide repeats dominated (75.74%) different types of motifs, and the main motif was CAA/TTG (44.25%). We successfully developed 38 SSR markers that exhibited substantial polymorphism, displaying an average of 6.18 alleles with accompanying average polymorphism information content (PIC) value of 0.516. The markers were used to evaluate the genetic diversity of 52 individuals collected from three populations of S. aralocaspica in Xinjiang, China. The results showed that the genetic diversity was moderate to high, with a mean expected heterozygosity (He) of 0.614, a mean Shannon's information index (I) of 1.23, and a mean genetic differentiation index (Fst) of 0.263. The SSR markers developed in this study provide a valuable resource for future genetic studies and breeding programs of S. aralocaspica, and even other species in Suaeda.

Detecting macroevolutionary genotype-phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype-phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω_C, which measures the error-corrected convergence rate of protein evolution. While ω_C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype-phenotype associations, even in lineages that diverged for hundreds of millions of years.

Genome-Wide Characterization and Analysis of the bHLH Transcription Factor Family in Suaeda aralocaspica, an Annual Halophyte With Single-Cell C4 Anatomy

Basic helix-loop-helix (bHLH) transcription factors play important roles in plant growth, development, metabolism, hormone signaling pathways, and responses to abiotic stresses. However, comprehensive genomic and functional analyses of bHLH genes have not yet been reported in desert euhalophytes. Suaeda aralocaspica, an annual C4 halophyte without Kranz anatomy, presents high photosynthetic efficiency in harsh natural habitats and is an ideal plant for identifying transcription factors involved in stress resistance. In this study, 83 bHLH genes in S. aralocaspica were identified and categorized into 21 subfamilies based on conserved motifs, gene structures, and phylogenetic analysis. Functional annotation enrichment revealed that the majority of SabHLHs were enriched in Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways involved in the response to stress conditions, as transcription factors. A number of cis-acting elements related to plant hormones and stress responses were also predicted in the promoter regions of SabHLHs, which were confirmed by expression analysis under various abiotic stress conditions (NaCl, mannitol, low temperature, ABA, GA3, MeJA, and SA); most were involved in tolerance to drought and salinity. SabHLH169 (076) protein localized in the nucleus was involved in transcriptional activity, and gene expression could be affected by different light qualities. This study is the first comprehensive analysis of the bHLH gene family in S. aralocaspica. These data will facilitate further characterization of their molecular functions in the adaptation of desert plants to abiotic stress.

The genome of the recretohalophyte Limonium bicolor provides insights into salt gland development and salinity adaptation during terrestrial evolution

Halophytes have evolved specialized strategies to cope with high salinity. The extreme halophyte sea lavender (Limonium bicolor) lacks trichomes but possesses salt glands on its epidermis that can excrete harmful ions, such as sodium, to avoid salt damage. Here, we report a high-quality, 2.92-Gb, chromosome-scale L. bicolor genome assembly based on a combination of Illumina short reads, single-molecule, real-time long reads, chromosome conformation capture (Hi-C) data, and Bionano genome maps, greatly enriching the genomic information on recretohalophytes with multicellular salt glands. Although the L. bicolor genome contains genes that show similarity to trichome fate genes from Arabidopsis thaliana, it lacks homologs of the decision fate genes GLABRA3, ENHANCER OF GLABRA3, GLABRA2, TRANSPARENT TESTA GLABRA2, and SIAMESE, providing a molecular explanation for the absence of trichomes in this species. We identified key genes (LbHLH and LbTTG1) controlling salt gland development among classical trichome homologous genes and confirmed their roles by showing that their mutations markedly disrupted salt gland initiation, salt secretion, and salt tolerance, thus offering genetic support for the long-standing hypothesis that salt glands and trichomes may share a common origin. In addition, a whole-genome duplication event occurred in the L. bicolor genome after its divergence from Tartary buckwheat and may have contributed to its adaptation to high salinity. The L. bicolor genome resource and genetic evidence reported in this study provide profound insights into plant salt tolerance mechanisms that may facilitate the engineering of salt-tolerant crops.

The Developmental Delay of Seedlings With Cotyledons Only Confers Stress Tolerance to Suaeda aralocaspica (Chenopodiaceae) by Unique Performance on Morphology, Physiology, and Gene Expression

Cotyledons play an important role in seedling establishment, although they may just exist for a short time and become senescent upon the emergence of euphylla. So far, the detailed function of cotyledons has not been well understood. Suaeda aralocaspica is an annual halophyte distributed in cold deserts; its cotyledons could exist for a longer time, even last until maturity, and they must exert a unique function in seedling development. Therefore, in this study, we conducted a series of experiments to investigate the morphological and physiological performances of cotyledons under salt stress at different developmental stages. The results showed that the cotyledons kept growing slowly to maintain the normal physiological activities of seedlings by balancing phytohormone levels, accumulating osmoprotectants and antioxidants, and scavenging reactive oxygen species (ROS). Salt stress activated the expression of osmoprotectant-related genes and enhanced the accumulation of related primary metabolites. Furthermore, differentially expressed transcriptional profiles of the cotyledons were also analyzed by cDNA-AFLP to gain an understanding of cotyledons in response to development and salt stress, and the results revealed a progressive increase in the expression level of development-related genes, which accounted for a majority of the total tested TDFs. Meanwhile, key photosynthetic and important salt stress-related genes also actively responded. All these performances suggest that "big cotyledons" are experiencing a delayed but active developmental process, by which S. aralocaspica may survive the harsh condition of the seedling stage.

Optimum Sterilization Method for In Vitro Cultivation of Dimorphic Seeds of the Succulent Halophyte Suaeda aralocaspica

Suaeda aralocaspica is an annual halophyte in the Amaranthaceae in the saline deserts of central Asia. This plant has succulent leaves and grape-like fruits and is a potential horticultural plant. To obtain the efficient sterilization method and optimal culture conditions, two types of seeds produced from a single plant of S. aralocaspica were treated with 75% ethanol for different time durations first, and then sodium hypochlorite (NaClO) or mercury chloride (HgCl2), with five different timing treatments were used for second seed surface sterilization. Sterilized seeds were germinated on a Murashige and Skoog (MS) medium at different potential hydrogenation (pH) levels, to examine germination and seedling performance. The results showed that the highest germination percentage of brown seeds was 100% and that of black seeds was 17%. Thus, brown seeds were more suitable for further culture experiments than black seeds. For brown seeds, the sterilization effect of NaClO was better than that of HgCl2, based on the results of seed germination, contamination, and seedling survival. Rinsing with 75% ethanol for 60 s, sterilizing with NaClO for 8 min, and cultivating at pH 8.0 MS for 7 days was the best of all sterilization procedures and cultivation methods tested, which has been successfully applied to S. aralocaspica in vitro culture. The optimized protocol described here can be used as the reference for the Suaeda genus.

Chromosome-level genome assembly of a xerophytic plant, Haloxylon ammodendron

Haloxylon ammodendron is a xerophytic perennial shrub or small tree that has a high ecological value in anti-desertification due to its high tolerance to drought and salt stress. Here, we report a high-quality, chromosome-level genome assembly of H. ammodendron by integrating PacBio's high-fidelity sequencing and Hi-C technology. The assembled genome size was 685.4 Mb, of which 99.6% was assigned to nine pseudochromosomes with a contig N50 value of 23.6 Mb. Evolutionary analysis showed that both the recent substantial amplification of long terminal repeat retrotransposons and tandem gene duplication may have contributed to its genome size expansion and arid adaptation. An ample amount of low-GC genes was closely related to functions that may contribute to the desert adaptation of H. ammodendron. Gene family clustering together with gene expression analysis identified differentially expressed genes that may play important roles in the direct response of H. ammodendron to water-deficit stress. We also identified several genes possibly related to the degraded scaly leaves and well-developed root system of H. ammodendron. The reference-level genome assembly presented here will provide a valuable genomic resource for studying the genome evolution of xerophytic plants, as well as for further genetic breeding studies of H. ammodendron.

Rewilding staple crops for the lost halophytism: Toward sustainability and profitability of agricultural production systems

Abiotic stress tolerance has been weakened during the domestication of all major staple crops. Soil salinity is a major environmental constraint that impacts over half of the world population; however, given the increasing reliance on irrigation and the lack of available freshwater, agriculture in the 21st century will increasingly become saline. Therefore, global food security is critically dependent on the ability of plant breeders to create high-yielding staple crop varieties that will incorporate salinity tolerance traits and account for future climate scenarios. Previously, we have argued that the current agricultural practices and reliance on crops that exclude salt from uptake is counterproductive and environmentally unsustainable, and thus called for a need for a major shift in a breeding paradigm to incorporate some halophytic traits that were present in wild relatives but were lost in modern crops during domestication. In this review, we provide a comprehensive physiological and molecular analysis of the key traits conferring crop halophytism, such as vacuolar Na⁺ sequestration, ROS desensitization, succulence, metabolic photosynthetic switch, and salt deposition in trichomes, and discuss the strategies for incorporating them into elite germplasm, to address a pressing issue of boosting plant salinity tolerance.

Genome-Wide Identification and Analysis of the Phosphoenolpyruvate Carboxylase Gene Family in Suaeda aralocaspica, an Annual Halophyte With Single-Cellular C4 Anatomy

Phosphoenolpyruvate carboxylase (PEPC) plays pivotal roles in the carbon fixation of photosynthesis and a variety of metabolic and stress pathways. Suaeda aralocaspica belongs to a single-cellular C₄ species and carries out a photosynthetic pathway in an unusually elongated chlorenchyma cell, which is expected to have PEPCs with different characteristics. To identify the different isoforms of PEPC genes in S. aralocaspica and comparatively analyze their expression and regulation patterns as well as the biochemical and enzymatic properties in this study, we characterized a bacterial-type PEPC (BTPC; SaPEPC-4) in addition to the two plant-type PEPCs (PTPCs; SaPEPC-1 and SaPEPC-2) using a genome-wide identification. SaPEPC-4 presented a lower expression level in all test combinations with an unknown function; two SaPTPCs showed distinct subcellular localizations and different spatiotemporal expression patterns but positively responded to abiotic stresses. Compared to SaPEPC-2, the expression of SaPEPC-1 specifically in chlorenchyma cell tissues was much more active with the progression of development and under various stresses, particularly sensitive to light, implying the involvement of SaPEPC-1 in a C₄ photosynthetic pathway. In contrast, SaPEPC-2 was more like a non-photosynthetic PEPC. The expression trends of two SaPTPCs in response to light, development, and abiotic stresses were also matched with the changes in PEPC activity in vivo (native) or in vitro (recombinant), and the biochemical properties of the two recombinant SaPTPCs were similar in response to various effectors while the catalytic efficiency, substrate affinity, and enzyme activity of SaPEPC-2 were higher than that of SaPEPC-1 in vitro. All the different properties between these two SaPTPCs might be involved in transcriptional (e.g., specific cis-elements), posttranscriptional [e.g., 5'-untranslated region (5'-UTR) secondary structure], or translational (e.g., PEPC phosphorylation/dephosphorylation) regulatory events. The comparative studies on the different isoforms of the PEPC gene family in S. aralocaspica may help to decipher their exact role in C₄ photosynthesis, plant growth/development, and stress resistance.

The first draft genome of Picrorhiza kurrooa, an endangered medicinal herb from Himalayas

Picrorhiza kurrooa is an endangered medicinal herb which is distributed across the Himalayan region at an altitude between 3000–5000 m above mean sea level. The medicinal properties of P. kurrooa are attributed to monoterpenoid picrosides present in leaf, rhizome and root of the plant. However, no genomic information is currently available for P. kurrooa, which limits our understanding about its molecular systems and associated responses. The present study brings the first assembled draft genome of P. kurrooa by using 227 Gb of raw data generated by Illumina and PacBio RS II sequencing platforms. The assembled genome has a size of n = ~ 1.7 Gb with 12,924 scaffolds. Four pronged assembly quality validations studies, including experimentally reported ESTs mapping and directed sequencing of the assembled contigs, confirmed high reliability of the assembly. About 76% of the genome is covered by complex repeats alone. Annotation revealed 24,798 protein coding and 9789 non-coding genes. Using the assembled genome, a total of 710 miRNAs were discovered, many of which were found responsible for molecular response against temperature changes. The miRNAs and targets were validated experimentally. The availability of draft genome sequence will aid in genetic improvement and conservation of P. kurrooa. Also, this study provided an efficient approach for assembling complex genomes while dealing with repeats when regular assemblers failed to progress due to repeats.

Plant Biosystems Design Research Roadmap 1.0

Human life intimately depends on plants for food, biomaterials, health, energy, and a sustainable environment. Various plants have been genetically improved mostly through breeding, along with limited modification via genetic engineering, yet they are still not able to meet the ever-increasing needs, in terms of both quantity and quality, resulting from the rapid increase in world population and expected standards of living. A step change that may address these challenges would be to expand the potential of plants using biosystems design approaches. This represents a shift in plant science research from relatively simple trial-and-error approaches to innovative strategies based on predictive models of biological systems. Plant biosystems design seeks to accelerate plant genetic improvement using genome editing and genetic circuit engineering or create novel plant systems through de novo synthesis of plant genomes. From this perspective, we present a comprehensive roadmap of plant biosystems design covering theories, principles, and technical methods, along with potential applications in basic and applied plant biology research. We highlight current challenges, future opportunities, and research priorities, along with a framework for international collaboration, towards rapid advancement of this emerging interdisciplinary area of research. Finally, we discuss the importance of social responsibility in utilizing plant biosystems design and suggest strategies for improving public perception, trust, and acceptance.

Genome of extreme halophyte Puccinellia tenuiflora

Background

Puccinellia tenuiflora, a forage grass, is considered a model halophyte given its strong tolerance for multiple stress conditions and its close genetic relationship with cereals. This halophyte has enormous values for improving our understanding of salinity tolerance mechanisms. The genetic information of P. tenuiflora also is a potential resource that can be used for improving the salinity tolerance of cereals.

Results

Here, we sequenced and assembled the P. tenuiflora genome (2n = 14) through the combined strategy of Illumina, PacBio, and 10× genomic technique. We generated 43.2× PacBio long reads, 123.87× 10× genomic reads, and 312.6× Illumina reads. Finally, we assembled 2638 scaffolds with a total size of 1.107 Gb, contig N50 of 117 kb, and scaffold N50 of 950 kb. We predicted 39,725 protein-coding genes, and identified 692 tRNAs, 68 rRNAs, 702 snRNAs, 1376 microRNAs, and 691 Mb transposable elements.

Conclusions

We deposited the genome sequence in NCBI and the Genome Warehouse in National Genomics Data Center. Our work may improve current understanding of plant salinity tolerance, and provides extensive genetic resources necessary for improving the salinity and drought tolerance of cereals.

A draft genome assembly of halophyte Suaeda aralocaspica, a plant that performs C4 photosynthesis within individual cells

BACKGROUND

The halophyte Suaeda aralocaspica performs complete C4 photosynthesis within individual cells (SCC4), which is distinct from typical C4 plants, which require the collaboration of 2 types of photosynthetic cells. However, despite SCC4 plants having features that are valuable in engineering higher photosynthetic efficiencies in agriculturally important C3 species such as rice, there are no reported sequenced SCC4 plant genomes, limiting our understanding of the mechanisms involved in, and evolution of, SCC4 photosynthesis.

FINDINGS

Using Illumina and Pacific Biosciences sequencing platforms, we generated ∼202 Gb of clean genomic DNA sequences having a 433-fold coverage based on the 467 Mb estimated genome size of S. aralocaspica. The final genome assembly was 452 Mb, consisting of 4,033 scaffolds, with a scaffold N50 length of 1.83 Mb. We annotated 29,604 protein-coding genes using Evidence Modeler based on the gene information from ab initio predictions, homology levels with known genes, and RNA sequencing-based transcriptome evidence. We also annotated noncoding genes, including 1,651 long noncoding RNAs, 21 microRNAs, 382 transfer RNAs, 88 small nuclear RNAs, and 325 ribosomal RNAs. A complete (circular with no gaps) chloroplast genome of S. aralocaspica 146,654 bp in length was also assembled.

CONCLUSIONS

We have presented the genome sequence of the SCC4 plant S. aralocaspica. Knowledge of the genome of S. aralocaspica should increase our understanding of the evolution of SCC4 photosynthesis and contribute to the engineering of C4 photosynthesis into economically important C3 crops.