Wolfberry genome database: integrated genomic datasets for studying molecular biology.

Wolfberry, also known as goji berry or Lycium barbarum, is a highly valued fruit with significant health benefits and nutritional value. For more efficient and comprehensive usage of published L. barbarum genomic data, we established the Wolfberry database. The utility of the Wolfberry Genome Database (WGDB) is highlighted through the Genome browser, which enables the user to explore the L. barbarum genome, browse specific chromosomes, and access gene sequences. Gene annotation features provide comprehensive information about gene functions, locations, expression profiles, pathway involvement, protein domains, and regulatory transcription factors. The transcriptome feature allows the user to explore gene expression patterns using transcripts per kilobase million (TPM) and fragments per kilobase per million mapped reads (FPKM) metrics. The Metabolism pathway page provides insights into metabolic pathways and the involvement of the selected genes. In addition to the database content, we also introduce six analysis tools developed for the WGDB. These tools offer functionalities for gene function prediction, nucleotide and amino acid BLAST analysis, protein domain analysis, GO annotation, and gene expression pattern analysis. The WGDB is freely accessible at https://cosbi7.ee.ncku.edu.tw/Wolfberry/. Overall, WGDB serves as a valuable resource for researchers interested in the genomics and transcriptomics of L. barbarum. Its user-friendly web interface and comprehensive data facilitate the exploration of gene functions, regulatory mechanisms, and metabolic pathways, ultimately contributing to a deeper understanding of wolfberry and its potential applications in agronomy and nutrition.

Diploid and tetraploid genomes of Acorus and the evolution of monocots.

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

OrchidBase 5.0: updates of the orchid genome knowledgebase.

Containing the largest number of species, the orchid family provides not only materials for studying plant evolution and environmental adaptation, but economically and culturally important ornamental plants for human society. Previously, we collected genome and transcriptome information of Dendrobium catenatum, Phalaenopsis equestris, and Apostasia shenzhenica which belong to two different subfamilies of Orchidaceae, and developed user-friendly tools to explore the orchid genetic sequences in the OrchidBase 4.0. The OrchidBase 4.0 offers the opportunity for plant science community to compare orchid genomes and transcriptomes and retrieve orchid sequences for further study.In the year 2022, two whole-genome sequences of Orchidoideae species, Platanthera zijinensis and Platanthera guangdongensis, were de novo sequenced, assembled and analyzed. In addition, systemic transcriptomes from these two species were also established. Therefore, we included these datasets to develop the new version of OrchidBase 5.0. In addition, three new functions including synteny, gene order, and miRNA information were also developed for orchid genome comparisons and miRNA characterization.OrchidBase 5.0 extended the genetic information to three orchid subfamilies (including five orchid species) and provided new tools for orchid researchers to analyze orchid genomes and transcriptomes. The online resources can be accessed at https://cosbi.ee.ncku.edu.tw/orchidbase5/.

Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy.

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensis and holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

R2R3-MYB genes coordinate conical cell development and cuticular wax biosynthesis in Phalaenopsis aphrodite.

Petals of the monocot Phalaenopsis aphrodite (Orchidaceae) possess conical epidermal cells on their adaxial surfaces, and a large amount of cuticular wax is deposited on them to serve as a primary barrier against biotic and abiotic stresses. It has been widely reported that subgroup 9A members of the R2R3-MYB gene family, MIXTA and MIXTA-like in eudicots, act to regulate the differentiation of conical epidermal cells. However, the molecular pathways underlying conical epidermal cell development and cuticular wax biosynthesis in monocot petals remain unclear. Here, we characterized two subgroup 9A R2R3-MYB genes, PaMYB9A1 and PaMYB9A2 (PaMYB9A1/2), from P. aphrodite through the transient overexpression of their coding sequences and corresponding chimeric repressors in developing petals. We showed that PaMYB9A1/2 function to coordinate conical epidermal cell development and cuticular wax biosynthesis. In addition, we identified putative targets of PaMYB9A1/2 through comparative transcriptome analyses, revealing that PaMYB9A1/2 acts to regulate the expression of cell wall-associated and wax biosynthetic genes. Furthermore, a chemical composition analysis of cuticular wax showed that even-chain n-alkanes and odd-chain primary alcohols are the main chemical constituents of cuticular wax deposited on petals, which is inconsistent with the well-known biosynthetic pathways of cuticular wax, implying a distinct biosynthetic pathway occurring in P. aphrodite flowers. These results reveal that the function of subgroup 9A R2R3-MYB family genes in regulating the differentiation of epidermal cells is largely conserved in monocots and dicots. Furthermore, both PaMYB9A1/2 have evolved additional functions controlling the biosynthesis of cuticular wax.

OrchidBase 4.0: a database for orchid genomics and molecular biology.

BACKGROUND: The Orchid family is the largest families of the monocotyledons and an economically important ornamental plant worldwide. Given the pivotal role of this plant to humans, botanical researchers and breeding communities should have access to valuable genomic and transcriptomic information of this plant. Previously, we established OrchidBase, which contains expressed sequence tags (ESTs) from different tissues and developmental stages of Phalaenopsis as well as biotic and abiotic stress-treated Phalaenopsis. The database includes floral transcriptomic sequences from 10 orchid species across all the five subfamilies of Orchidaceae. DESCRIPTION: Recently, the whole-genome sequences of Apostasia shenzhenica, Dendrobium catenatum, and Phalaenopsis equestris were de novo assembled and analyzed. These datasets were used to develop OrchidBase 4.0, including genomic and transcriptomic data for these three orchid species. OrchidBase 4.0 offers information for gene annotation, gene expression with fragments per kilobase of transcript per millions mapped reads (FPKM), KEGG pathways and BLAST search. In addition, assembled genome sequences and location of genes and miRNAs could be visualized by the genome browser. The online resources in OrchidBase 4.0 can be accessed by browsing or using BLAST. Users can also download the assembled scaffold sequences and the predicted gene and protein sequences of these three orchid species. CONCLUSIONS: OrchidBase 4.0 is the first database that contain the whole-genome sequences and annotations of multiple orchid species. OrchidBase 4.0 is available at http://orchidbase.itps.ncku.edu.tw/.

Corrigendum to "Comparative analysis of Phytophthora genomes reveals oomycete pathogenesis in crops" [Heliyon 7 (2) (February 2021) e06317].

[This corrects the article DOI: 10.1016/j.heliyon.2021.e06317.].

Wolfberry genomes and the evolution of Lycium (Solanaceae).

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

The ancestral duplicated DL/CRC orthologs, PeDL1 and PeDL2, function in orchid reproductive organ innovation.

Orchid gynostemium, the fused organ of the androecium and gynoecium, and ovule development are unique developmental processes. Two DROOPING LEAF/CRABS CLAW (DL/CRC) genes, PeDL1 and PeDL2, were identified from the Phalaenopsis orchid genome and functionally characterized. Phylogenetic analysis indicated that the most recent common ancestor of orchids contained the duplicated DL/CRC-like genes. Temporal and spatial expression analysis indicated that PeDL genes are specifically expressed in the gynostemium and at the early stages of ovule development. Both PeDLs could partially complement an Arabidopsis crc-1 mutant. Virus-induced gene silencing (VIGS) of PeDL1 and PeDL2 affected the number of protuberant ovule initials differentiated from the placenta. Transient overexpression of PeDL1 in Phalaenopsis orchids caused abnormal development of ovule and stigmatic cavity of gynostemium. PeDL1, but not PeDL2, could form a heterodimer with Phalaenopsis equestris CINCINNATA 8 (PeCIN8). Paralogous retention and subsequent divergence of the gene sequences of PeDL1 and PeDL2 in P. equestris might result in the differentiation of function and protein behaviors. These results reveal that the ancestral duplicated DL/CRC-like genes play important roles in orchid reproductive organ innovation.

Comparative analysis of Phytophthora genomes reveals oomycete pathogenesis in crops.

The oomycete genus Phytophthora includes devastating plant pathogens that are found in almost all ecosystems. We sequenced the genomes of two quarantined Phytophthora species-P. fragariae and P. rubi. Comparing these Phytophthora species and related genera allowed reconstruction of the phylogenetic relationships within the genus Phytophthora and revealed Phytophthora genomic features associated with infection and pathogenicity. We found that several hundred Phytophthora genes are putatively inherited from red algae, but Phytophthora does not have vestigial plastids originating from phototrophs. The horizontally-transferred Phytophthora genes are abundant transposons that "transmit" exogenous gene to Phytophthora species thus bring about the gene recombination possibility. Several expansion events of Phytophthora gene families associated with cell wall biogenesis can be used as mutational targets to elucidate gene function in pathogenic interactions with host plants. This work enhanced the understanding of Phytophthora evolution and will also be helpful for the design of phytopathological control strategies.

Orchid Bsister gene PeMADS28 displays conserved function in ovule integument development.

The ovules and egg cells are well developed to be fertilized at anthesis in many flowering plants. However, ovule development is triggered by pollination in most orchids. In this study, we characterized the function of a Bsister gene, named PeMADS28, isolated from Phalaenopsis equestris, the genome-sequenced orchid. Spatial and temporal expression analysis showed PeMADS28 predominantly expressed in ovules between 32 and 48 days after pollination, which synchronizes with integument development. Subcellular localization and protein-protein interaction analyses revealed that PeMADS28 could form a homodimer as well as heterodimers with D-class and E-class MADS-box proteins. In addition, ectopic expression of PeMADS28 in Arabidopsis thaliana induced small curled rosette leaves, short silique length and few seeds, similar to that with overexpression of other species' Bsister genes in Arabidopsis. Furthermore, complementation test revealed that PeMADS28 could rescue the phenotype of the ABS/TT16 mutant. Together, these results indicate the conserved function of Bsister PeMADS28 associated with ovule integument development in orchid.

Identification of high-copy number long terminal repeat retrotransposons and their expansion in Phalaenopsis orchids.

BACKGROUND: Transposable elements (TEs) are fragments of DNA that can insert into new chromosomal locations. They represent a great proportion of eukaryotic genomes. The identification and characterization of TEs facilitates understanding the transpositional activity of TEs with their effects on the orchid genome structure. RESULTS: We combined the draft whole-genome sequences of Phalaenopsis equestris with BAC end sequences, Roche 454, and Illumina/Solexa, and identified long terminal repeat (LTR) retrotransposons in these genome sequences by using LTRfinder and classified by using Gepard software. Among the 10 families Gypsy-like retrotransposons, three families Gypsy1, Gypsy2, and Gypsy3, contained the most copies among these predicted elements. In addition, six high-copy retrotransposons were identified according to their reads in the sequenced raw data. The 12-kb Orchid-rt1 contains 18,000 copies representing 220 Mbp of the P. equestris genome. Southern blot and slot blot assays showed that these four retrotransposons Gypsy1, Gypsy2, Gypsy3, and Orchid-rt1 contained high copies in the large-genome-size/large-chromosome species P. violacea and P. bellina. Both Orchid-rt1 and Gypsy1 displayed various ratios of copy number for the LTR sequences versus coding sequences among four Phalaenopsis species, including P. violacea and P. bellina and small-genome-size/small-chromosome P. equestris and P. ahprodite subsp. formosana, which suggests that Orchid-rt1 and Gypsy1 have been through various mutations and homologous recombination events. FISH results showed amplification of Orchid-rt1 in the euchromatin regions among the four Phalaenopsis species. The expression levels of Peq018599 encoding copper transporter 1 is highly upregulated with the insertion of Orchid-rt1, while it is down regulated for Peq009948 and Peq014239 encoding for a 26S proteasome non-ATP regulatory subunit 4 homolog and auxin-responsive factor AUX/IAA-related. In addition, insertion of Orchid-rt1 in these three genes are all in their intron regions. CONCLUSION: Orchid-rt1 and Gypsy1-3 have amplified within Phalaenopsis orchids concomitant with the expanded genome sizes, and Orchid-rt1 and Gypsy1 may have gone through various mutations and homologous recombination events. Insertion of Orchid-rt1 is in the introns and affects gene expression levels.

Genome-Wide Identification of YABBY Genes in Orchidaceae and Their Expression Patterns in Phalaenopsis Orchid.

The plant YABBY transcription factors are key regulators in the lamina development of lateral organs. Orchid is one of the largest families in angiosperm and known for their unique floral morphology, reproductive biology, and diversified lifestyles. However, nothing is known about the role of YABBY genes in orchids, although biologists have never lost their fascination with orchids. In this study, a total of 54 YABBY genes, including 15 genes in CRC/DL, eight in INO, 17 in YAB2, and 14 in FIL clade, were identified from the eight orchid species. A sequence analysis showed that all protein sequences encoded by these YABBY genes share the highly conserved C2C2 zinc-finger domain and YABBY domain (a helix-loop-helix motif). A gene structure analysis showed that the number of exons is highly conserved in the same clades. The genes in YAB2 clade have six exons, and genes in CRC/DL, INO, and FIL have six or seven exons. A phylogenetic analysis showed all 54 orchid YABBY genes could be classified into four major clades, including CRC/DL, INO, FIL, and YAB2. Many of orchid species maintain more than one member in CRC/DL, FIL, and YAB2 clades, implying functional differentiation among these genes, which is supported by sequence diversification and differential expression. An expression analysis of PhalaenopsisYABBY genes revealed that members in the CRC/DL clade have concentrated expressions in the early floral development stage and gynostemium, the fused male and female reproductive organs. The expression of PeINO is consistent with the biological role it played in ovule integument morphogenesis. Transcripts of members in the FIL clade could be obviously detected at the early developmental stage of the flowers. The expression of three genes, PeYAB2,PeYAB3, and PeYAB4, in the YAB2 clade could be revealed both in vegetative and reproductive tissues, and PeYAB4 was transcribed at a relatively higher level than that of PeYAB2 and PeYAB3. Together, this comprehensive analysis provides the basic information for understanding the function of the YABBY gene in Orchidaceae.

Author Correction: The Apostasia genome and the evolution of orchids.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Expression regulation of MALATE SYNTHASE involved in glyoxylate cycle during protocorm development in Phalaenopsis aphrodite (Orchidaceae).

Orchid (Orchidaceae) is one of the largest families in angiosperms and presents exceptional diversity in lifestyle. Their unique reproductive characteristics of orchid are attracted by scientist for centuries. One of the synapomorphies of orchid plants is that their seeds do not contain endosperm. Lipids are used as major energy storage in orchid seeds. However, regulation and mobilization of lipid usage during early seedling (protocorm) stage of orchid is not understood. In this study, we compared transcriptomes from developing Phalaenopsis aphrodite protocorms grown on 1/2-strength MS medium with sucrose. The expression of P. aphrodite MALATE SYNTHASE (PaMLS), involved in the glyoxylate cycle, was significantly decreased from 4 days after incubation (DAI) to 7 DAI. On real-time RT-PCR, both P. aphrodite ISOCITRATE LYASE (PaICL) and PaMLS were down-regulated during protocorm development and suppressed by sucrose treatment. In addition, several genes encoding transcription factors regulating PaMLS expression were identified. A gene encoding homeobox transcription factor (named PaHB5) was involved in positive regulation of PaMLS. This study showed that sucrose regulates the glyoxylate cycle during orchid protocorm development in asymbiotic germination and provides new insights into the transcription factors involved in the regulation of malate synthase expression.

The genome sequence of star fruit (Averrhoa carambola).

Oxalidaceae is one of the most important plant families in horticulture, and its key commercially relevant genus, Averrhoa, has diverse growth habits and fruit types. Here, we describe the assembly of a high-quality chromosome-scale genome sequence for Averrhoa carambola (star fruit). Ks distribution analysis showed that A. carambola underwent a whole-genome triplication event, i.e., the gamma event shared by most eudicots. Comparisons between A. carambola and other angiosperms also permitted the generation of Oxalidaceae gene annotations. We identified unique gene families and analyzed gene family expansion and contraction. This analysis revealed significant changes in MADS-box gene family content, which might be related to the cauliflory of A. carambola. In addition, we identified and analyzed a total of 204 nucleotide-binding site, leucine-rich repeat receptor (NLR) genes and 58 WRKY genes in the genome, which may be related to the defense response. Our results provide insights into the origin, evolution and diversification of star fruit.

The Litsea genome and the evolution of the laurel family.

The laurel family within the Magnoliids has attracted attentions owing to its scents, variable inflorescences, and controversial phylogenetic position. Here, we present a chromosome-level assembly of the Litsea cubeba genome, together with low-coverage genomic and transcriptomic data for many other Lauraceae. Phylogenomic analyses show phylogenetic discordance at the position of Magnoliids, suggesting incomplete lineage sorting during the divergence of monocots, eudicots, and Magnoliids. An ancient whole-genome duplication (WGD) event occurred just before the divergence of Laurales and Magnoliales; subsequently, independent WGDs occurred almost simultaneously in the three Lauralean lineages. The phylogenetic relationships within Lauraceae correspond to the divergence of inflorescences, as evidenced by the phylogeny of FUWA, a conserved gene involved in determining panicle architecture in Lauraceae. Monoterpene synthases responsible for production of specific volatile compounds in Lauraceae are functionally verified. Our work sheds light on the evolution of the Lauraceae, the genetic basis for floral evolution and specific scents.

PePIF1, a P-lineage of PIF-like transposable element identified in protocorm-like bodies of Phalaenopsis orchids.

BACKGROUND: Orchids produce a colorless protocorm by symbiosis with fungi upon seed germination. For mass production of orchids, the prevailing approaches are both generation of protocorm-like bodies (PLBs) from callus and multiplication of adventitious buds on inflorescence. However, somaclonal variations occur during micropropagation. RESULTS: We isolated the two most expressed transposable elements belonging to P Instability Factor (PIF)-like transposons. Among them, a potential autonomous element was identified by similarity analysis against the whole-genome sequence of Phalaenopsis equestris and named PePIF1. It contains a 19-bp terminal inverted repeat flanked by a 3-bp target site duplication and two coding regions encoding ORF1- and transposase-like proteins. Phylogenetic analysis revealed that PePIF1 belongs to a new P-lineage of PIF. Furthermore, two distinct families, PePIF1a and PePIF1b, with 29 and 37 putative autonomous elements, respectively, were isolated, along with more than 3000 non-autonomous and miniature inverted-repeat transposable element (MITE)-like elements. Among them, 828 PePIF1-related elements were inserted in 771 predicted genes. Intriguingly, PePIF1 was transposed in the somaclonal variants of Phalaenopsis cultivars, as revealed by transposon display, and the newly inserted genes were identified and sequenced. CONCLUSION: A PIF-like element, PePIF1, was identified in the Phalaenopsis genome and actively transposed during micropropagation. With the identification of PePIF1, we have more understanding of the Phalaenopsis genome structure and somaclonal variations during micropropagation for use in orchid breeding and production.

Post genomics era for orchid research.

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.

Critical Intermediate Structure That Directs the Crystalline Texture and Surface Morphology of Organo-Lead Trihalide Perovskite.

We have identified an often observed yet unresolved intermediate structure in a popular processing with dimethylformamide solutions of lead chloride and methylammonium iodide for perovskite solar cells. With subsecond time-resolved grazing-incidence X-ray scattering and X-ray photoemission spectroscopy, supplemental with ab initio calculation, the resolved intermediate structure (CH3NH3)2PbI2Cl2·CH3NH3I features two-dimensional (2D) perovskite bilayers of zigzagged lead-halide octahedra and sandwiched CH3NH3I layers. Such intermediate structure reveals a hidden correlation between the intermediate phase and the composition of the processing solution. Most importantly, the 2D perovskite lattice of the intermediate phase is largely crystallographically aligned with the [110] planes of the three-dimensional perovskite cubic phase; consequently, with sublimation of Cl ions from the organo-lead octahedral terminal corners in prolonged annealing, the zigzagged octahedral layers of the intermediate phase can merge with the intercalated methylammonium iodide layers for templated growth of perovskite crystals. Regulated by annealing temperature and the activation energies of the intermediate and perovskite, deduced from analysis of temperature-dependent structural kinetics, the intermediate phase is found to selectively mature first and then melt along the layering direction for epitaxial conversion into perovskite crystals. The unveiled epitaxial conversion under growth kinetics controls might be general for solution-processed and intermediate-templated perovskite formation.