The bracteatus pineapple genome and domestication of clonally propagated crops

Evolutionary signatures of the erosion of sexual reproduction genes in domesticated cassava (Manihot esculenta)

Centuries of clonal propagation in cassava (Manihot esculenta) have reduced sexual recombination, leading to the accumulation of deleterious mutations. This has resulted in both inbreeding depression affecting yield and a significant decrease in reproductive performance, creating hurdles for contemporary breeding programs. Cassava is a member of the Euphorbiaceae family, including notable species such as rubber tree (Hevea brasiliensis) and poinsettia (Euphorbia pulcherrima). Expanding upon preliminary draft genomes, we annotated 7 long-read genome assemblies and aligned a total of 52 genomes, to analyze selection across the genome and the phylogeny. Through this comparative genomic approach, we identified 48 genes under relaxed selection in cassava. Notably, we discovered an overrepresentation of floral expressed genes, especially focused at 6 pollen-related genes. Our results indicate that domestication and a transition to clonal propagation have reduced selection pressures on sexually reproductive functions in cassava leading to an accumulation of mutations in pollen-related genes. This relaxed selection and the genome-wide deleterious mutations responsible for inbreeding depression are potential targets for improving cassava breeding, where the generation of new varieties relies on recombining favorable alleles through sexual reproduction.

The Role of Aspergillus niger in Regulating Internal Browning Involves Flavonoid Biosynthesis and the Endophytic Fungal Community of Pineapple

Endophytic fungi are commonly used to control plant diseases, overcoming the drawbacks of chemical agents. The internal browning (IB) of postharvest pineapple fruit, a physiological disease, leads to quality losses and limits industrial development. This work investigated the relationship among the effects of Aspergillus niger (An) on IB controlling, flavonoid metabolism and the endophytic fungal community of pineapple through metabolomics, transcriptomics, microbiomics and microorganism mutagenesis technology. We obtained an endophyte An that can control the IB of pineapple and screened its mutant strain AnM, through chemical mutagenesis, that cannot control IB. The transcriptome of fungi showed that An and AnM were different in oxidative metabolism. Transcriptome and metabolome analyses of pineapple showed that An upregulated genes of flavonoid synthesis, including dihydroflavonol 4-reductase and flavonoid 3'-monooxygenase and increased the flavonoid content in pineapple fruit, i.e., Hispidulin, Hispidulin-7-O-Glucoside, and Diosmetin, while AnM could not. Microbiomics analysis identified an increase in the abundance of eight endophytic fungi in An-inoculated fruit, among which the abundance of six endophytic fungi (Filobasidium magnum, Naganishia albida, A. niger, Aureobasidium melanogenum, Kwoniella heveanensis and Lysurus cruciatus) was positively correlated with the content of three flavonoids mentioned above but not in AnM-inoculated fruit. Overall, this suggested, for the first time, that A. niger alleviated IB mainly by enhancing flavonoid synthesis and content and the abundance of endophytic fungi and by regulating the interaction between flavonoid content and endophytic fungi abundance in pineapple. This work adds to the understanding of the IB mechanism in postharvest pineapple and provides a new green approach for reducing postharvest losses and controlling physiological diseases.

Tracing the evolutionary and genetic footprints of atmospheric tillandsioids transition from land to air

Plant evolution is driven by key innovations of functional traits that enables their survivals in diverse ecological environments. However, plant adaptive evolution from land to atmospheric niches remains poorly understood. In this study, we use the epiphytic Tillandsioideae subfamily of Bromeliaceae as model plants to explore their origin, evolution and diversification. We provide a comprehensive phylogenetic tree based on nuclear transcriptomic sequences, indicating that core tillandsioids originated approximately 11.3 million years ago in the Andes. The geological uplift of the Andes drives the divergence of tillandsioids into tank-forming and atmospheric types. Our genomic and transcriptomic analyses reveal gene variations and losses associated with adaptive traits such as impounding tanks and absorptive trichomes. Furthermore, we uncover specific nitrogen-fixing bacterial communities in the phyllosphere of tillandsioids as potential source of nitrogen acquisition. Collectively, our study provides integrative multi-omics insights into the adaptive evolution of tillandsioids in response to elevated aerial habitats.

CAM evolution is associated with gene family expansion in an explosive bromeliad radiation

The subgenus Tillandsia (Bromeliaceae) belongs to one of the fastest radiating clades in the plant kingdom and is characterized by the repeated evolution of Crassulacean acid metabolism (CAM). Despite its complex genetic basis, this water-conserving trait has evolved independently across many plant families and is regarded as a key innovation trait and driver of ecological diversification in Bromeliaceae. By producing high-quality genome assemblies of a Tillandsia species pair displaying divergent photosynthetic phenotypes, and combining genome-wide investigations of synteny, transposable element (TE) dynamics, sequence evolution, gene family evolution, and temporal differential expression, we were able to pinpoint the genomic drivers of CAM evolution in Tillandsia. Several large-scale rearrangements associated with karyotype changes between the 2 genomes and a highly dynamic TE landscape shaped the genomes of Tillandsia. However, our analyses show that rewiring of photosynthetic metabolism is mainly obtained through regulatory evolution rather than coding sequence evolution, as CAM-related genes are differentially expressed across a 24-h cycle between the 2 species but are not candidates of positive selection. Gene orthology analyses reveal that CAM-related gene families manifesting differential expression underwent accelerated gene family expansion in the constitutive CAM species, further supporting the view of gene family evolution as a driver of CAM evolution.

The pineapple reference genome: Telomere-to-telomere assembly, manually curated annotation, and comparative analysis

Pineapple is the third most crucial tropical fruit worldwide and available in five varieties. Genomes of different pineapple varieties have been released to date; however, none of them are complete, with all exhibiting substantial gaps and representing only two of the five pineapple varieties. This significantly hinders the advancement of pineapple breeding efforts. In this study, we sequenced the genomes of three varieties: a wild pineapple variety, a fiber pineapple variety, and a globally cultivated edible pineapple variety. We constructed the first gap-free reference genome (Ref) for pineapple. By consolidating multiple sources of evidence and manually revising each gene structure annotation, we identified 26,656 protein-coding genes. The BUSCO evaluation indicated a completeness of 99.2%, demonstrating the high quality of the gene structure annotations in this genome. Utilizing these resources, we identified 7,209 structural variations across the three varieties. Approximately 30.8% of pineapple genes were located within ±5 kb of structural variations, including 30 genes associated with anthocyanin synthesis. Further analysis and functional experiments demonstrated that the high expression of AcMYB528 aligns with the accumulation of anthocyanins in the leaves, both of which may be affected by a 1.9-kb insertion fragment. In addition, we developed the Ananas Genome Database, which offers data browsing, retrieval, analysis, and download functions. The construction of this database addresses the lack of pineapple genome resource databases. In summary, we acquired a seamless pineapple reference genome with high-quality gene structure annotations, providing a solid foundation for pineapple genomics and a valuable reference for pineapple breeding.

Genomic variation across distribution of Micro-Tom, a model cultivar of tomato (Solanum lycopersicum)

Micro-Tom is a cultivar of tomato (Solanum lycopersicum), which is known as a major crop and model plant in Solanaceae. Micro-Tom has phenotypic traits such as dwarfism, and substantial EMS-mutagenized lines have been reported. After Micro-Tom was generated in Florida, USA, it was distributed to research institutes worldwide and used as a genetic resource. In Japan, the Micro-Tom lines have been genetically fixed; currently, three lines have been re-distributed from three institutes, but many phenotypes among the lines have been observed. We have determined the genome sequence de novo of the Micro-Tom KDRI line, one of the Micro-Tom lines distributed from Kazusa DNA Research Institute (KDRI) in Japan, and have built chromosome-scale pseudomolecules. Genotypes among six Micro-Tom lines, including three in Japan, one in the United States, one in France, and one in Brazil showed phenotypic alternation. Here, we unveiled the swift emergence of genetic diversity in both phenotypes and genotypes within the Micro-Tom genome sequence during its propagation. These findings offer valuable insights crucial for the management of bioresources.

AcMYB266, a key regulator of the red coloration in pineapple peel: a case of subfunctionalization in tandem duplicated genes

Red fruit peel is an attractive target for pineapple breeding. Various pineapple accessions with distinct red coloration patterns exist; however, the precise molecular mechanism accounting for these differences remains unknown, which hinders the pineapple breeding process from combining high fruit quality with red peel. In this study, we characterized a transcription factor, AcMYB266, which is preferentially expressed in pineapple peel and positively regulates anthocyanin accumulation. Transgenic pineapple, Arabidopsis, and tobacco plants overexpressing AcMYB266 exhibited significant anthocyanin accumulation. Conversely, transient silencing of this gene led to decreased anthocyanin accumulation in pineapple red bracts. In-depth analysis indicated that variations of AcMYB266 sequences in the promoter instead of the protein-coding region seem to contribute to different red coloration patterns in peels of three representative pineapple varieties. In addition, we found that AcMYB266 was located in a cluster of four MYB genes exclusive to and conserved in Ananas species. Of this cluster, each was proved to regulate anthocyanin synthesis in different pineapple tissues, illustrating an interesting case of gene subfunctionalization after tandem duplication. In summary, we have characterized AcMYB266 as a key regulator of pineapple red fruit peel and identified an MYB cluster whose members were subfunctionalized to specifically regulate the red coloration of different pineapple tissues. The present study will assist in establishing a theoretical mechanism for pineapple breeding for red fruit peel and provide an interesting case for the investigation of gene subfunctionalization in plants.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

TCOD: an integrated resource for tropical crops

Tropical crops are vital for tropical agriculture, with resource scarcity, functional diversity and extensive market demand, providing considerable economic benefits for the world's tropical agriculture-producing countries. The rapid development of sequencing technology has promoted a milestone in tropical crop research, resulting in the generation of massive amount of data, which urgently needs an effective platform for data integration and sharing. However, the existing databases cannot fully satisfy researchers' requirements due to the relatively limited integration level and untimely update. Here, we present the Tropical Crop Omics Database (TCOD, https://ngdc.cncb.ac.cn/tcod), a comprehensive multi-omics data platform for tropical crops. TCOD integrates diverse omics data from 15 species, encompassing 34 chromosome-level de novo assemblies, 1 255 004 genes with functional annotations, 282 436 992 unique variants from 2048 WGS samples, 88 transcriptomic profiles from 1997 RNA-Seq samples and 13 381 germplasm items. Additionally, TCOD not only employs genes as a bridge to interconnect multi-omics data, enabling cross-species comparisons based on homology relationships, but also offers user-friendly online tools for efficient data mining and visualization. In short, TCOD integrates multi-species, multi-omics data and online tools, which will facilitate the research on genomic selective breeding and trait biology of tropical crops.

Identification of bromelain subfamily proteases encoded in the pineapple genome

Papain (aka C1A) family proteases, including bromelain enzymes, are widespread across the plant kingdom and play critical regulatory functions in protein turnover during development. The proteolytic activity exhibited by papain family proteases has led to their increased usage for a wide range of cosmetic, therapeutic, and medicinal purposes. Bromelain enzymes, or bromelains in short, are members of the papain family that are specific to the bromeliad plant family. The only major commercial extraction source of bromelain is pineapple. The importance of C1A family and bromelain subfamily proteases in pineapple development and their increasing economic importance led several researchers to utilize available genomic resources to identify protease-encoding genes in the pineapple genome. To date, studies are lacking in screening bromelain genes for targeted use in applied science studies. In addition, the bromelain genes coding for the enzymes present in commercially available bromelain products have not been identified and their evolutionary origin has remained unclear. Here, using the newly developed MD2 v2 pineapple genome, we aimed to identify bromelain-encoding genes and elucidate their evolutionary origin. Orthologous and phylogenetic analyses of all papain-family proteases encoded in the pineapple genome revealed a single orthogroup (189) and phylogenetic clade (XIII) containing the bromelain subfamily. Duplication mode and synteny analyses provided insight into the origin and expansion of the bromelain subfamily in pineapple. Proteomic analysis identified four bromelain enzymes present in two commercially available bromelain products derived from pineapple stem, corresponding to products of four putative bromelain genes. Gene expression analysis using publicly available transcriptome data showed that 31 papain-family genes identified in this study were up-regulated in specific tissues, including stem, fruit, and floral tissues. Some of these genes had higher expression in earlier developmental stages of different tissues. Similar expression patterns were identified by RT-qPCR analysis with leaf, stem, and fruit. Our results provide a strong foundation for future applicable studies on bromelain, such as transgenic approaches to increase bromelain content in pineapple, development of bromelain-producing bioreactors, and studies that aim to determine the medicinal and/or therapeutic viability of individual bromelain enzymes.

Loss of Heterozygosity and Its Importance in Evolution

Loss of heterozygosity (LOH) is a mitotic recombination event that converts heterozygous loci to homozygous loci. This mutation event is widespread in organisms that have asexual reproduction like budding yeasts, and is also an important and frequent mutation event in tumorigenesis. Mutation accumulation studies have demonstrated that LOH occurs at a rate higher than the point mutation rate, and can impact large portions of the genome. Laboratory evolution experiments of heterozygous yeasts have revealed that LOH often unmasks beneficial recessive alleles that can confer large fitness advantages. Here, I highlight advances in understanding dominance, fitness, and phenotypes in laboratory evolved heterozygous yeast strains. I discuss best practices for detecting LOH in intraspecific and interspecific evolved clones and populations. Utilizing heterozygous strain backgrounds in laboratory evolution experiments offers an opportunity to advance our understanding of this important mutation type in shaping adaptation and genome evolution in wild, domesticated, and clinical populations.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

Whole-genome Duplication Reshaped Adaptive Evolution in A Relict Plant Species, Cyclocarya paliurus

Cyclocarya paliurus is a relict plant species that survived the last glacial period and shows a population expansion recently. Its leaves have been traditionally used to treat obesity and diabetes with the well-known active ingredient cyclocaric acid B. Here, we presented three C. paliurus genomes from two diploids with different flower morphs and one haplotype-resolved tetraploid assembly. Comparative genomic analysis revealed two rounds of recent whole-genome duplication events and identified 691 genes with dosage effects that likely contribute to adaptive evolution through enhanced photosynthesis and increased accumulation of triterpenoids. Resequencing analysis of 45 C. paliurus individuals uncovered two bottlenecks, consistent with the known events of environmental changes, and many selectively swept genes involved in critical biological functions, including plant defense and secondary metabolite biosynthesis. We also proposed the biosynthesis pathway of cyclocaric acid B based on multi-omics data and identified key genes, in particular gibberellin-related genes, associated with the heterodichogamy in C. paliurus species. Our study sheds light on evolutionary history of C. paliurus and provides genomic resources to study the medicinal herbs.

Molecular regulation of tomato male reproductive development

The reproductive success of flowering plants, which directly affects crop yield, is sensitive to environmental changes. A thorough understanding of how crop reproductive development adapts to climate changes is vital for ensuring global food security. In addition to being a high-value vegetable crop, tomato is also a model plant used for research on plant reproductive development. Tomato crops are cultivated under highly diverse climatic conditions worldwide. Targeted crosses of hybrid varieties have resulted in increased yields and abiotic stress resistance; however, tomato reproduction, especially male reproductive development, is sensitive to temperature fluctuations, which can lead to aborted male gametophytes, with detrimental effects on fruit set. We herein review the cytological features as well as genetic and molecular pathways influencing tomato male reproductive organ development and responses to abiotic stress. We also compare the shared features among the associated regulatory mechanisms of tomato and other plants. Collectively, this review highlights the opportunities and challenges related to characterizing and exploiting genic male sterility in tomato hybrid breeding programs.

Genome-Wide Identification and Characterization of R2R3-MYB Provide Insight into Anthocyanin Biosynthesis Regulation Mechanism of Ananas comosus var. bracteatus

The R2R3-MYB proteins comprise the largest class of MYB transcription factors, which play an essential role in regulating anthocyanin synthesis in various plant species. Ananas comosus var. bracteatus is an important colorful anthocyanins-rich garden plant. The spatio-temporal accumulation of anthocyanins in chimeric leaves, bracts, flowers, and peels makes it an important plant with a long ornamental period and highly improves its commercial value. We conducted a comprehensive bioinformatic analysis of the R2R3-MYB gene family based on genome data from A. comosus var. bracteatus. Phylogenetic analysis, gene structure and motif analysis, gene duplication, collinearity, and promoter analysis were used to analyze the characteristics of this gene family. In this work, a total of 99 R2R3-MYB genes were identified and classified into 33 subfamilies according to phylogenetic analysis, and most of them were localized in the nucleus. We found these genes were mapped to 25 chromosomes. Gene structure and protein motifs were conserved among AbR2R3-MYB genes, especially within the same subfamily. Collinearity analysis revealed four pairs of tandem duplicated genes and 32 segmental duplicates in AbR2R3-MYB genes, indicating that segmental duplication contributed to the amplification of the AbR2R3-MYB gene family. A total of 273 ABRE responsiveness, 66 TCA elements, 97 CGTCA motifs, and TGACG motifs were the main cis elements in the promoter region under response to ABA, SA, and MEJA. These results revealed the potential function of AbR2R3-MYB genes in response to hormone stress. Ten R2R3-MYBs were found to have high homology to MYB proteins reported to be involved in anthocyanin biosynthesis from other plants. RT-qPCR results revealed the 10 AbR2R3-MYB genes showed tissue-specific expression patterns, six of them expressed the highest in the flower, two genes in the bract, and two genes in the leaf. These results suggested that these genes may be the candidates that regulate anthocyanin biosynthesis of A. comosus var. bracteatus in the flower, leaf, and bract, respectively. In addition, the expressions of these 10 AbR2R3-MYB genes were differentially induced by ABA, MEJA, and SA, implying that these genes may play crucial roles in hormone-induced anthocyanin biosynthesis. Our study provided a comprehensive and systematic analysis of AbR2R3-MYB genes and identified the AbR2R3-MYB genes regulating the spatial-temporal anthocyanin biosynthesis in A. comosus var. bracteatus, which would be valuable for further study on the anthocyanin regulation mechanism of A. comosus var. bracteatus.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

Genomic and transcriptomic-based analysis of agronomic traits in sugar beet (Beta vulgaris L.) pure line IMA1

Sugar beet (Beta vulgaris L.) is an important sugar-producing and energy crop worldwide. The sugar beet pure line IMA1 independently bred by Chinese scientists is a standard diploid parent material that is widely used in hybrid-breeding programs. In this study, a high-quality, chromosome-level genome assembly for IMA1was conducted, and 99.1% of genome sequences were assigned to nine chromosomes. A total of 35,003 protein-coding genes were annotated, with 91.56% functionally annotated by public databases. Compared with previously released sugar beet assemblies, the new genome was larger with at least 1.6 times larger N50 size, thereby substantially improving the completeness and continuity of the sugar beet genome. A Genome-Wide Association Studies analysis identified 10 disease-resistance genes associated with three important beet diseases and five genes associated with sugar yield per hectare, which could be key targets to improve sugar productivity. Nine highly expressed genes associated with pollen fertility of sugar beet were also identified. The results of this study provide valuable information to identify and dissect functional genes affecting sugar beet agronomic traits, which can increase sugar beet production and help screen for excellent sugar beet breeding materials. In addition, information is provided that can precisely incorporate biotechnology tools into breeding efforts.

Chromosome-scale genome assemblies and annotations for Poales species Carex cristatella, Carex scoparia, Juncus effusus, and Juncus inflexus

The majority of sequenced genomes in the monocots are from species belonging to Poaceae, which include many commercially important crops. Here, we expand the number of sequenced genomes from the monocots to include the genomes of 4 related cyperids: Carex cristatella and Carex scoparia from Cyperaceae and Juncus effusus and Juncus inflexus from Juncaceae. The high-quality, chromosome-scale genome sequences from these 4 cyperids were assembled by combining whole-genome shotgun sequencing of Nanopore long reads, Illumina short reads, and Hi-C sequencing data. Some members of the Cyperaceae and Juncaceae are known to possess holocentric chromosomes. We examined the repeat landscapes in our sequenced genomes to search for potential repeats associated with centromeres. Several large satellite repeat families, comprising 3.2–9.5% of our sequenced genomes, showed dispersed distribution of large satellite repeat clusters across all Carex chromosomes, with few instances of these repeats clustering in the same chromosomal regions. In contrast, most large Juncus satellite repeats were clustered in a single location on each chromosome, with sporadic instances of large satellite repeats throughout the Juncus genomes. Recognizable transposable elements account for about 20% of each of the 4 genome assemblies, with the Carex genomes containing more DNA transposons than retrotransposons while the converse is true for the Juncus genomes. These genome sequences and annotations will facilitate better comparative analysis within monocots.

Evolution of the DEHYDRATION-RESPONSIVE ELEMENT-BINDING PROTEIN subfamily in green plants

Adapting to unfavorable environments is a necessary step in plant terrestrialization and radiation. The dehydration-responsive element-binding (DREB) protein subfamily plays a pivotal role in plant abiotic stress regulation. However, relationships between the origin and expansion of the DREB subfamily and adaptive evolution of land plants are still being elucidated. Here, we constructed the evolutionary history of the DREB subfamily by compiling APETALA2/ethylene-responsive element-binding protein superfamily genes from 169 representative species of green plants. Through extensive phylogenetic analyses and comparative genomic analysis, our results revealed that the DREB subfamily diverged from the ethylene-responsive factor (ERF) subfamily in the common ancestor of Zygnemophyceae and Embryophyta during the colonization of land by plants, followed by expansions to form three different ancient archetypal genes in Zygnemophyceae species, designated as groups archetype-I, archetype-II/III, and archetype-IV. Four large-scale expansions paralleling the evolution of land plants led to the nine-subgroup divergence of group archetype-II/III in angiosperms, and five whole-genome duplications during Brassicaceae and Poaceae radiation shaped the diversity of subgroup IIb-1. We identified a Poaceae-specific gene in subgroup IIb-1, ERF014, remaining in a Poaceae-specific microsynteny block and co-evolving with a small heat shock protein cluster. Expression analyses demonstrated that heat acclimation may have driven the neofunctionalization of ERF014s in Pooideae by engaging in the conserved heat-responsive module in Poaceae. This study provides insights into lineage-specific expansion and neofunctionalization in the DREB subfamily, together with evolutionary information valuable for future functional studies of plant stress biology.

SunUp and Sunset genomes revealed impact of particle bombardment mediated transformation and domestication history in papaya

Transgenic papaya is widely publicized for controlling papaya ringspot virus. However, the impact of particle bombardment on the genome remains unknown. The transgenic SunUp and its progenitor Sunset genomes were assembled into 351.5 and 350.3 Mb in nine chromosomes, respectively. We identified a 1.64 Mb insertion containing three transgenic insertions in SunUp chromosome 5, consisting of 52 nuclear-plastid, 21 nuclear-mitochondrial and 1 nuclear genomic fragments. A 591.9 kb fragment in chromosome 5 was translocated into the 1.64 Mb insertion. We assembled a gapless 9.8 Mb hermaphrodite-specific region of the Y^h chromosome and its 6.0 Mb X counterpart. Resequencing 86 genomes revealed three distinct groups, validating their geographic origin and breeding history. We identified 147 selective sweeps and defined the essential role of zeta-carotene desaturase in carotenoid accumulation during domestication. Our findings elucidated the impact of particle bombardment and improved our understanding of sex chromosomes and domestication to expedite papaya improvement.

Gene prediction for leaf margin phenotype and fruit flesh color in pineapple (Ananas comosus) using haplotype-resolved genome sequencing

Pineapple (Ananas comosus (L.) Merr.) is one of the most economically important tropical fruit species. The major aim of the breeding programs in several countries, including Japan, is quality improvement, mainly for the fresh market. ‘Yugafu’, a Japanese cultivar with distinctive pipe-type leaf margin phenotype and white flesh color, is popular for fresh consumption. Therefore, genome sequencing of ‘Yugafu’ is expected to assist pineapple breeding. Here, we developed a haplotype-resolved assembly for the heterozygous genome of ‘Yugafu’ using long-read sequencing technology and obtained a pair of 25 pseudomolecule sequences inherited from the parental accessions ‘Cream pineapple’ and ‘HI101’. The causative genes for leaf margin and fruit flesh color were identified. Fine mapping revealed a 162-kb region on CLG23 for the leaf margin phenotype. In this region, 20 kb of inverted repeat was specifically observed in the ‘Cream pineapple’ derived allele, and the WUSCHEL-related homeobox 3 (AcWOX3) gene was predicted as the key gene for leaf margin morphogenesis. Dominantly repressed AcWOX3 via RNAi was suggested to be the cause of the pipe-type leaf margin phenotype. Quantitative trait locus (QTL) analysis revealed that the terminal region of CLG08 contributed to white flesh and low carotenoid content. Carotenoid cleaved dioxygenase 4 (AcCCD4), a key gene for carotenoid degradation underlying this QTL, was predicted as the key gene for white flesh color through expression analysis. These findings could assist in modern pineapple breeding and facilitate marker-assisted selection for important traits.

The Gastrodia menghaiensis (Orchidaceae) genome provides new insights of orchid mycorrhizal interactions

BACKGROUND

To illustrate the molecular mechanism of mycoheterotrophic interactions between orchids and fungi, we assembled chromosome-level reference genome of Gastrodia menghaiensis (Orchidaceae) and analyzed the genomes of two species of Gastrodia.

RESULTS

Our analyses indicated that the genomes of Gastrodia are globally diminished in comparison to autotrophic orchids, even compared to Cuscuta (a plant parasite). Genes involved in arbuscular mycorrhizae colonization were found in genomes of Gastrodia, and many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids. The highly expressed genes for fatty acid and ammonium root transporters suggest that fungi receive material from orchids, although most raw materials flow from the fungi. Many nuclear genes (e.g. biosynthesis of aromatic amino acid L-tryptophan) supporting plastid functions are expanded compared to photosynthetic orchids, an indication of the importance of plastids even in totally mycoheterotrophic species.

CONCLUSION

Gastrodia menghaiensis has the smallest proteome thus far among angiosperms. Many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids.

Reproductive biology of wild and domesticated Ensete ventricosum: Further evidence for maintenance of sexual reproductive capacity in a vegetatively propagated perennial crop

Loss of sexual reproductive capacity has been proposed as a syndrome of domestication in vegetatively propagated crops, but there are relatively few examples from agricultural systems. In this study, we compare sexual reproductive capacity in wild (sexual) and domesticated (vegetative) populations of enset (Ensete ventricosum (Welw.) Cheesman), a tropical banana relative and Ethiopian food security crop. We examined floral and seed morphology and germination ecology across 35 wild and domesticated enset. We surveyed variation in floral and seed traits, including seed weight, viability and internal morphology, and germinated seeds across a range of constant and alternating temperature regimes to characterize optimum germination requirements. We report highly consistent floral allometry, seed viability, internal morphology and days to germination in wild and domesticated enset. However, seeds from domesticated plants responded to cooler temperatures with greater diurnal range. Shifts in germination behaviour appear concordant with a climatic envelope shift in the domesticated distribution. Our findings provide evidence that sexual reproductive capacity has been maintained despite long-term near-exclusive vegetative propagation in domesticated enset. Furthermore, certain traits such as germination behaviour and floral morphology may be under continued selection, presumably through rare sexually reproductive events. Compared to sexually propagated crops banked as seeds, vegetative crop diversity is typically conserved in living collections that are more costly and insecure. Improved understanding of sexual propagation in vegetative crops may have applications in germplasm conservation and plant breeding.

The highly continuous reference genome of a leaf-chimeric red pineapple (Ananas comosus var. bracteatus f. tricolor) provides insights into elaboration of leaf color

Ananas comosus var. bracteatus f. tricolor (GL1) is a red pineapple accession whose mostly green leaves with chimeric white leaf margins turn red in spring and autumn and during flowering. It is an important ornamental plant and ideal plant research model for anthocyanin metabolism, chimeric leaf development, and photosynthesis. Here, we generated a highly contiguous chromosome-scale genome assembly for GL1 and compared it with other 3 published pineapple assemblies (var. comosus accessions MD2 and F153, and var. bracteatus accession CB5). The GL1 assembly has a total size of ∼461 Mb, with a contig N50 of ∼2.97 Mb and Benchmarking Universal Single-Copy Ortholog score of 97.3%. More than 99% of the contigs are anchored to 25 pseudochromosomes. Compared with the other 3 published pineapple assemblies, the GL1 assembly was confirmed to be more continuous. Our evolutionary analysis showed that the Bromeliaceae and Poaceae diverged from their nearest common ancestor ∼82.36 million years ago (MYA). Population structure analysis showed that while GL1 has not undergone admixture, bracteatus accession CB5 has resulted from admixture of 3 species of Ananas. Through classification of orthogroups, analysis of genes under positive selection, and analysis of presence/absence variants, we identified a series of genes related to anthocyanin metabolism and development of chimeric leaves. The structure and evolution of these genes were compared among the published pineapple assemblies with reveal candidate genes for these traits. The GL1 genome assembly and its comparisons with other 3 pineapple genome assemblies provide a valuable resource for the genetic improvement of pineapple and serve as a model for understanding the genomic basis of important traits in different pineapple varieties and other pan-cereal crops.

Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales

The Poales is one of the largest orders of flowering plants with significant economic and ecological values. Reconstructing the phylogeny of the Poales is important for understanding its evolutionary history that forms the basis for biological studies. However, due to sparse taxon sampling and limited molecular data, previous studies have resulted in a variety of contradictory topologies. In particular, there are three nodes surrounded by incongruence: the phylogenetic ambiguity near the root of the Poales tree, the sister family of Poaceae, and the delimitation of the xyrid clade. We conducted a comprehensive sampling and reconstructed the phylogenetic tree using plastid and mitochondrial genomic data from 91 to 66 taxa, respectively, representing all the 16 families of Poales. Our analyses support the finding of Bromeliaceae and Typhaceae as the earliest diverging groups within the Poales while having phylogenetic relationships with the polytomy. The clade of Ecdeiocoleaceae and Joinvilleaceae is recovered as the sister group of Poaceae. The three families, Mayacaceae, Eriocaulaceae, and Xyridaceae, of the xyrid assembly diverged successively along the backbone of the Poales phylogeny, and thus this assembly is paraphyletic. Surprisingly, we find substantial phylogenetic conflicts within the plastid genomes of the Poales, as well as among the plastid, mitochondrial, and nuclear data. These conflicts suggest that the Poales could have a complicated evolutionary history, such as rapid radiation and polyploidy, particularly allopolyploidy through hybridization. In sum, our study presents a new perspicacity into the complex phylogenetic relationships and the underlying phylogenetic conflicts within the Poales.

Origin, loss, and regain of self-incompatibility in angiosperms

The self-incompatibility (SI) system with the broadest taxonomic distribution in angiosperms is based on multiple S-locus F-box genes (SLFs) tightly linked to an S-RNase termed type-1. Multiple SLFs collaborate to detoxify nonself S-RNases while being unable to detoxify self S-RNases. However, it is unclear how such a system evolved, because in an ancestral system with a single SLF, many nonself S-RNases would not be detoxified, giving low cross-fertilization rates. In addition, how the system has been maintained in the face of whole-genome duplications (WGDs) or lost in other lineages remains unclear. Here we show that SLFs from a broad range of species can detoxify S-RNases from Petunia with a high detoxification probability, suggestive of an ancestral feature enabling cross-fertilization and subsequently modified as additional SLFs evolved. We further show, based on its genomic signatures, that type-1 was likely maintained in many lineages, despite WGD, through deletion of duplicate S-loci. In other lineages, SI was lost either through S-locus deletions or by retaining duplications. Two deletion lineages regained SI through type-2 (Brassicaceae) or type-4 (Primulaceae), and one duplication lineage through type-3 (Papaveraceae) mechanisms. Thus, our results reveal a highly dynamic process behind the origin, maintenance, loss, and regain of SI.

Signatures of selection in recently domesticated macadamia

Macadamia is a high value nut crop that is recently domesticated, ideal for testing the effect of artificial selection. Here, we sequence the genome of Hawaiian cultivar ‘Kau’ and assemble into 794 Mb in 14 pseudo-chromosomes with 37,728 genes. Genome analysis reveals a whole-genome duplication event, occurred 46.8 million years ago. Gene expansions occurred in gene families involves in fatty acid biosynthesis. Gene duplication of MADS-Box transcription factors in proanthocyanidin biosynthesis are relevant for seed coat development. Genome re-sequencing of 112 accessions reveals the origin of Hawaiian cultivars from Mount Bauple in southeast Queensland in Australia. Selective sweeps are detected in macadamia cultivars, including genes involved in fatty acid biosynthesis, seed coat development, and heat stress response. Such strong effects of artificial selection in few generations reveals the genomic basis for ‘one-step operation’ for clonal crop domestication. The knowledge gained could accelerate domestication of new crops from wild species.

Macadamia is a recently domesticated nut crop. Here, the authors report the genome assembly of Hawaiian cultivar ‘Kau’ and conduct population genomic analyses to reveal the origin of Hawaiian cultivars and the genomic basis for one-step operation for the clonal crop domestication.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated bioinformatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for over 90 plant genomes including agronomically important cereals (e.g., maize, sorghum, wheat, teff), fruits and vegetables (e.g., apple, watermelon, clementine, tomato, cassava), specialty crops (e.g., coffee, olive tree, pistachio, almond), and plants of special or emerging interest (e.g., cotton, tobacco, cannabis, or hemp). For some species, the resource includes multiple varieties of the same species, which has paved the road for the creation of species-specific pan-genome browsers. The resource also features plant research models, including Arabidopsis and C4 warm-season grasses and brassicas, as well as other species that fill phylogenetic gaps for plant evolution studies. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end-users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) explore Gramene's search results, (2) explore gene-centric comparative genomics data visualizations in Gramene, and (3) explore genetic variation associated with a gene locus. This is the first publication describing in detail Gramene's integrated search interface-intended to provide a simplified entry portal for the resource's main data categories (genomic location, phylogeny, gene expression, pathways, and external references) to the most complete and up-to-date set of plant genome and pathway annotations.

Improved High-Quality Genome Assembly and Annotation of Pineapple (Ananas comosus) Cultivar MD2 Revealed Extensive Haplotype Diversity and Diversified FRS/FRF Gene Family

Pineapple (Ananas comosus (L.) Merr.) is the second most important tropical fruit crop globally, and ‘MD2’ is the most important cultivated variety. A high-quality genome is important for molecular-based breeding, but available pineapple genomes still have some quality limitations. Here, PacBio and Hi-C data were used to develop a new high-quality MD2 assembly and gene prediction. Compared to the previous MD2 assembly, major improvements included a 26.6-fold increase in contig N50 length, phased chromosomes, and >6000 new genes. The new MD2 assembly also included 161.6 Mb additional sequences and >3000 extra genes compared to the F153 genome. Over 48% of the predicted genes harbored potential deleterious mutations, indicating that the high level of heterozygosity in this species contributes to maintaining functional alleles. The genome was used to characterize the FAR1-RELATED SEQUENCE (FRS) genes that were expanded in pineapple and rice. Transposed and dispersed duplications contributed to expanding the numbers of these genes in the pineapple lineage. Several AcFRS genes were differentially expressed among tissue-types and stages of flower development, suggesting that their expansion contributed to evolving specialized functions in reproductive tissues. The new MD2 assembly will serve as a new reference for genetic and genomic studies in pineapple.

Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity

Background

Heterozygous genomes are widespread in outcrossing and clonally propagated crops. However, the variation in heterozygosity underlying key agronomic traits and crop domestication remains largely unknown. Cassava is a staple crop in Africa and other tropical regions and has a highly heterozygous genome.

Results

We describe a genomic variation map from 388 resequenced genomes of cassava cultivars and wild accessions. We identify 52 loci for 23 agronomic traits through a genome-wide association study. Eighteen allelic variations in heterozygosity for nine candidate genes are significantly associated with seven key agronomic traits. We detect 81 selective sweeps with decreasing heterozygosity and nucleotide diversity, harboring 548 genes, which are enriched in multiple biological processes including growth, development, hormone metabolisms and responses, and immune-related processes. Artificial selection for decreased heterozygosity has contributed to the domestication of the large starchy storage root of cassava. Selection for homozygous GG allele in MeTIR1 during domestication contributes to increased starch content. Selection of homozygous AA allele in MeAHL17 is associated with increased storage root weight and cassava bacterial blight (CBB) susceptibility. We have verified the positive roles of MeTIR1 in increasing starch content and MeAHL17 in resistance to CBB by transient overexpression and silencing analysis. The allelic combinations in MeTIR1 and MeAHL17 may result in high starch content and resistance to CBB.

Conclusions

This study provides insights into allelic variation in heterozygosity associated with key agronomic traits and cassava domestication. It also offers valuable resources for the improvement of cassava and other highly heterozygous crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02524-7.

Characterization of germline development and identification of genes associated with germline specification in pineapple

Understanding germline specification in plants could be advantageous for agricultural applications. In recent decades, substantial efforts have been made to understand germline specification in several plant species, including Arabidopsis, rice, and maize. However, our knowledge of germline specification in many agronomically important plant species remains obscure. Here, we characterized the female germline specification and subsequent female gametophyte development in pineapple using callose staining, cytological, and whole-mount immunolocalization analyses. We also determined the male germline specification and gametophyte developmental timeline and observed male meiotic behavior using chromosome spreading assays. Furthermore, we identified 229 genes that are preferentially expressed at the megaspore mother cell (MMC) stage during ovule development and 478 genes that are preferentially expressed at the pollen mother cell (PMC) stage of anther development using comparative transcriptomic analysis. The biological functions, associated regulatory pathways and expression patterns of these genes were also analyzed. Our study provides a convenient cytological reference for exploring pineapple germline development and a molecular basis for the future functional analysis of germline specification in related plant species.

Recent amplification of microsatellite-associated miniature inverted-repeat transposable elements in the pineapple genome

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous DNA transposable elements that play important roles in genome organization and evolution. Genome-wide identification and characterization of MITEs provide essential information for understanding genome structure and evolution.

RESULTS

We performed genome-wide identification and characterization of MITEs in the pineapple genome. The top two MITE families, accounting for 29.39% of the total MITEs and 3.86% of the pineapple genome, have insertion preference in (TA) n dinucleotide microsatellite regions. We therefore named these MITEs A. comosus microsatellite-associated MITEs (Ac-mMITEs). The two Ac-mMITE families, Ac-mMITE-1 and Ac-mMITE-2, shared sequence similarity in the terminal inverted repeat (TIR) regions, suggesting that these two Ac-mMITE families might be derived from a common or closely related autonomous elements. The Ac-mMITEs are frequently clustered via adjacent insertions. Among the 21,994 full-length Ac-mMITEs, 46.1% of them were present in clusters. By analyzing the Ac-mMITEs without (TA) n microsatellite flanking sequences, we found that Ac-mMITEs were likely derived from Mutator-like DNA transposon. Ac-MITEs showed highly polymorphic insertion sites between cultivated pineapples and their wild relatives. To better understand the evolutionary history of Ac-mMITEs, we filtered and performed comparative analysis on the two distinct groups of Ac-mMITEs, microsatellite-targeting MITEs (mt-MITEs) that are flanked by dinucleotide microsatellites on both sides and mutator-like MITEs (ml-MITEs) that contain 9/10 bp TSDs. Epigenetic analysis revealed a lower level of host-induced silencing on the mt-MITEs in comparison to the ml-MITEs, which partially explained the significantly higher abundance of mt-MITEs in pineapple genome. The mt-MITEs and ml-MITEs exhibited differential insertion preference to gene-related regions and RNA-seq analysis revealed their differential influences on expression regulation of nearby genes.

CONCLUSIONS

Ac-mMITEs are the most abundant MITEs in the pineapple genome and they were likely derived from Mutator-like DNA transposon. Preferential insertion in (TA) n microsatellite regions of Ac-mMITEs occurred recently and is likely the result of damage-limiting strategy adapted by Ac-mMITEs during co-evolution with their host. Insertion in (TA) n microsatellite regions might also have promoted the amplification of mt-MITEs. In addition, mt-MITEs showed no or negligible impact on nearby gene expression, which may help them escape genome control and lead to their amplification.

Construction of a High-Density Genetic Map and Identification of Leaf Trait-Related QTLs in Chinese Bayberry (Myrica rubra)

Chinese bayberry (Myrica rubra) is an economically important fruit tree that is grown in southern China. Owing to its over 10-year seedling period, the crossbreeding of bayberry is challenging. The characteristics of plant leaves are among the primary factors that control plant architecture and potential yields, making the analysis of leaf trait-related genetic factors crucial to the hybrid breeding of any plant. In the present study, molecular markers associated with leaf traits were identified via a whole-genome re-sequencing approach, and a genetic map was thereby constructed. In total, this effort yielded 902.11 Gb of raw data that led to the identification of 2,242,353 single nucleotide polymorphisms (SNPs) in 140 F1 individuals and parents (Myrica rubra cv. Biqizhong × Myrica rubra cv. 2012LXRM). The final genetic map ultimately incorporated 31,431 SNPs in eight linkage groups, spanning 1,351.85 cM. This map was then used to assemble and update previous scaffold genomic data at the chromosomal level. The genome size of M. rubra was thereby established to be 275.37 Mb, with 94.98% of sequences being assembled into eight pseudo-chromosomes. Additionally, 18 quantitative trait loci (QTLs) associated with nine leaf and growth-related traits were identified. Two QTL clusters were detected (the LG3 and LG5 clusters). Functional annotations further suggested two chlorophyll content-related candidate genes being identified in the LG5 cluster. Overall, this is the first study on the QTL mapping and identification of loci responsible for the regulation of leaf traits in M. rubra, offering an invaluable scientific for future marker-assisted selection breeding and candidate gene analyses.

Draft genome of Puya raimondii (Bromeliaceae), the Queen of the Andes

Puya raimondii, the Queen of the Andes, is an endangered high Andean species in the Bromeliaceae family. Here, we report its first genome to promote its conservation and evolutionary study. Comparative genomics showed P. raimondii diverged from Ananas comosus about 14.8 million years ago, and the long terminal repeats were likely to contribute to the genus diversification in last 3.5 million years. The gene families related to plant reproductive development and stress responses significantly expanded in the genome. At the same time, gene families involved in disease defense, photosynthesis and carbohydrate metabolism significantly contracted, which may be an evolutionary strategy to adapt to the harsh conditions in high Andes. The demographic history analysis revealed the P. raimondii population size sharply declined in the Pleistocene and then increased in the Holocene. We also designed and tested 46 pairs of universal primers for amplifying orthologous single-copy nuclear genes in Puya species.

Metabolome and transcriptome profiling reveals anthocyanin contents and anthocyanin-related genes of chimeric leaves in Ananas comosus var. bracteatus

BACKGROUND

Ananas comosus var. bracteatus is a colorful plant used as a cut flower or landscape ornamental. The unique foliage color of this plant includes both green and red leaves and, as a trait of interest, deserves investigation. In order to explore the pigments behind the red section of the chimeric leaves, the green and red parts of chimeric leaves of Ananas comosus var. bracteatus were sampled and analyzed at phenotypic, cellular and molecular levels in this study.

RESULTS

The CIELAB results indicated that the a* values and L* values samples had significant differences between two parts. Freehand sections showed that anthocyanin presented limited accumulation in the green leaf tissues but obviously accumulation in the epidermal cells of red tissues. Transcriptomic and metabolomic analyses were performed by RNA-seq and LC-ESI-MS/MS. Among the 508 identified metabolites, 10 kinds of anthocyanins were detected, with 6 significantly different between the two samples. The cyanidin-3,5-O-diglucoside content that accounts for nearly 95.6% in red samples was significantly higher than green samples. RNA-Seq analyses showed that 11 out of 40 anthocyanin-related genes were differentially expressed between the green and red samples. Transcriptome and metabolome correlations were determined by nine quadrant analyses, and 9 anthocyanin-related genes, including MYB5 and MYB82, were correlated with 7 anthocyanin-related metabolites in the third quadrant in which genes and metabolites showing consistent change. Particularly, the PCCs between these two MYB genes and cyanidin-3,5-O-diglucoside were above 0.95.

CONCLUSION

Phenotypic colors are closely related to the tissue structures of different leaf parts of Ananas comosus var. bracteatus, and two MYB transcription factors might contribute to differences of anthocyanin accumulation in two parts of Ananas comosus var. bracteatus chimeric leaves. This study lay a foundation for further researches on functions of MYBs in Ananas comosus var. bracteatus and provides new insights to anthocyanin accumulation in different parts of chimeric leaves.

The population genomics of adaptive loss of function

Discoveries of adaptive gene knockouts and widespread losses of complete genes have in recent years led to a major rethink of the early view that loss-of-function alleles are almost always deleterious. Today, surveys of population genomic diversity are revealing extensive loss-of-function and gene content variation, yet the adaptive significance of much of this variation remains unknown. Here we examine the evolutionary dynamics of adaptive loss of function through the lens of population genomics and consider the challenges and opportunities of studying adaptive loss-of-function alleles using population genetics models. We discuss how the theoretically expected existence of allelic heterogeneity, defined as multiple functionally analogous mutations at the same locus, has proven consistent with empirical evidence and why this impedes both the detection of selection and causal relationships with phenotypes. We then review technical progress towards new functionally explicit population genomic tools and genotype-phenotype methods to overcome these limitations. More broadly, we discuss how the challenges of studying adaptive loss of function highlight the value of classifying genomic variation in a way consistent with the functional concept of an allele from classical population genetics.

Breeding Potentials of Bambara Groundnut for Food and Nutrition Security in the Face of Climate Change

Constant production of quality food should be a norm in any community, but climate change, increasing population, and unavailability of land for farming affect food production. As a result, food scarcity is affecting some communities, especially in the developing world. Finding a stable solution to this problem is a major cause of concern for researchers. Synergistic application of molecular marker techniques with next generation sequencing (NGS) technologies can unlock the potentials hidden in most crop genomes for improving yield and food availability. Most crops such as Bambara groundnut (BGN), Winged bean, and African yam bean are underutilized. These underutilized crops can compete with the major crops such as cowpea, soybean, maize, and rice, in areas of nutrition, ability to withstand drought stress, economic importance, and food production. One of these underutilized crops, BGN [Vigna subterranea (L.), Verdc.], is an indigenous African legume and can survive in tropical climates and marginal soils. In this review, we focus on the roles of BGN and the opportunities it possesses in tackling food insecurity and its benefits to local farmers. We will discuss BGN's potential impact on global food production and how the advances in NGS technologies can enhance its production.

Floral transcriptomes reveal gene networks in pineapple floral growth and fruit development

Proper flower development is essential for sexual reproductive success and the setting of fruits and seeds. The availability of a high quality genome sequence for pineapple makes it an excellent model for studying fruit and floral organ development. In this study, we sequenced 27 different pineapple floral samples and integrated nine published RNA-seq datasets to generate tissue- and stage-specific transcriptomic profiles. Pairwise comparisons and weighted gene co-expression network analysis successfully identified ovule-, stamen-, petal- and fruit-specific modules as well as hub genes involved in ovule, fruit and petal development. In situ hybridization confirmed the enriched expression of six genes in developing ovules and stamens. Mutant characterization and complementation analysis revealed the important role of the subtilase gene AcSBT1.8 in petal development. This work provides an important genomic resource for functional analysis of pineapple floral organ growth and fruit development and sheds light on molecular networks underlying pineapple reproductive organ growth.

On the necessity of combining ethnobotany and genetics to assess agrobiodiversity and its evolution in crops: A case study on date palms (Phoenix dactylifera L.) in Siwa Oasis, Egypt

Crop diversity is shaped by biological and social processes interacting at different spatiotemporal scales. Here, we combined population genetics and ethnobotany to investigate date palm (Phoenix dactylifera L.) diversity in Siwa Oasis, Egypt. Based on interviews with farmers and observation of practices in the field, we collected 149 date palms from Siwa Oasis and 27 uncultivated date palms from abandoned oases in the surrounding desert. Using genotyping data from 18 nuclear and plastid microsatellite loci, we confirmed that some named types each constitute a clonal line, that is, a true-to-type cultivar. We also found that others are collections of clonal lines, that is, ethnovarieties, or even unrelated samples, that is, local categories. This alters current assessments of agrobiodiversity, which are visibly underestimated, and uncovers the impact of low-intensity, but highly effective, farming practices on biodiversity. These hardly observable practices, hypothesized by ethnographic survey and confirmed by genetic analysis, are enabled by the way Isiwans conceive and classify living beings in their oasis, which do not quite match the way biologists do: a classic disparity of etic versus. emic categorizations. In addition, we established that Siwa date palms represent a unique and highly diverse genetic cluster, rather than a subset of North African and Middle Eastern palm diversity. As previously shown, North African date palms display evidence of introgression by the wild relative Phoenix theophrasti, and we found that the uncultivated date palms from the abandoned oases share even more alleles with this species than cultivated palms in this region. The study of Siwa date palms could hence be a key to the understanding of date palm diversification in North Africa. Integration of ethnography and population genetics promoted the understanding of the interplay between diversity management in the oasis (short-time scale), and the origins and dynamic of diversity through domestication and diversification (long-time scale).

Intertwined signatures of desiccation and drought tolerance in grasses

Grasses are among the most resilient plants, and some can survive prolonged desiccation in semiarid regions with seasonal rainfall. However, the genetic elements that distinguish grasses that are sensitive versus tolerant to extreme drying are largely unknown. Here, we leveraged comparative genomic approaches with the desiccation-tolerant grass Eragrostis nindensis and the related desiccation-sensitive cereal Eragrostis tef to identify changes underlying desiccation tolerance. These analyses were extended across C4 grasses and cereals to identify broader evolutionary conservation and divergence. Across diverse genomic datasets, we identified changes in chromatin architecture, methylation, gene duplications, and expression dynamics related to desiccation in E. nindensis It was previously hypothesized that transcriptional rewiring of seed desiccation pathways confers vegetative desiccation tolerance. Here, we demonstrate that the majority of seed-dehydration-related genes showed similar expression patterns in leaves of both desiccation-tolerant and -sensitive species. However, we identified a small set of seed-related orthologs with expression specific to desiccation-tolerant species. This supports a broad role for seed-related genes, where many are involved in typical drought responses, with only a small subset of crucial genes specifically induced in desiccation-tolerant plants.

Simulation-Based Evaluation of Three Methods for Local Ancestry Deconvolution of Non-model Crop Species Genomes

Hybridizations between species and subspecies represented major steps in the history of many crop species. Such events generally lead to genomes with mosaic patterns of chromosomal segments of various origins that may be assessed by local ancestry inference methods. However, these methods have mainly been developed in the context of human population genetics with implicit assumptions that may not always fit plant models. The purpose of this study was to evaluate the suitability of three state-of-the-art inference methods (SABER, ELAI and WINPOP) for local ancestry inference under scenarios that can be encountered in plant species. For this, we developed an R package to simulate genotyping data under such scenarios. The tested inference methods performed similarly well as far as representatives of source populations were available. As expected, the higher the level of differentiation between ancestral source populations and the lower the number of generations since admixture, the more accurate were the results. Interestingly, the accuracy of the methods was only marginally affected by i) the number of ancestries (up to six tested); ii) the sample design (i.e., unbalanced representation of source populations); and iii) the reproduction mode (e.g., selfing, vegetative propagation). If a source population was not represented in the data set, no bias was observed in inference accuracy for regions originating from represented sources and regions from the missing source were assigned differently depending on the methods. Overall, the selected ancestry inference methods may be used for crop plant analysis if all ancestral sources are known.

A role for MIR828 in pineapple fruit development

Chen et al. ( Nature Genet. 51: 1549-1558; Oct. 2019) sequenced Ananas comosus var. bracteatus accession CB5, cultivated for its bright pink-to-red colored fruit, and yellow-fleshed A. comosus accession F153, reporting an improved F153 reference assembly while annotating MICRORNA (MIRNA) loci and gene family expressions relevant to lignin and anthocyanin biosynthesis. An independent article (Xiong et al.Sci. Rep. 8: 1947; 2018) reported var. bracteatus MIRNAs but not MIR828, a negative regulator of anthocyanin and polyphenolics biosynthesis by targeting MYB transcription factors associated with UV light- and sugar-signaling in dicots. MIR828 has been reported in gymnosperms, Amborella (sister to flowering plants), and basal monocot orders Liliales, Asparagales, Zingiberales, Arecales, but not in the Poales, a sister order comprising grasses and ~3,000 species of bromeliads including pineapple. Here I show MIR828 exists in pineapple and directs post-transcriptional gene silencing of mRNAs encoding MYB family members with inferred function to regulate the conspicuous red fruit trait in var. bracteatus. MIR828 plesiomorphy (an ancient basal trait) may shed light on monocot apomorphic fruit development, postulated for 21 monocot families with fleshy fruits as due to homoplasy/convergence driven by tropical climate and/or enticements to vertebrate endozoic seed dispersers.

A role for MIR828 in pineapple fruit development

Chen et al. ( Nature Genet. 51: 1549-1558; Oct. 2019) sequenced Ananas comosus var. bracteatus accession CB5, cultivated for its bright pink-to-red colored fruit, and yellow-fleshed A. comosus accession F153, reporting an improved F153 reference assembly while annotating MICRORNA (MIRNA) loci and gene family expressions relevant to lignin and anthocyanin biosynthesis. An independent article (Xiong et al.Sci. Rep. 8: 1947; 2018) reported var. bracteatus MIRNAs but not MIR828, a negative regulator of anthocyanin and polyphenolics biosynthesis by targeting MYB transcription factors associated with UV light- and sugar-signaling in dicots. MIR828 has been reported in gymnosperms, Amborella (sister to flowering plants), and basal monocot orders Liliales, Asparagales, Zingiberales, Arecales, but not in the Poales, a sister order comprising grasses and ~3,000 species of bromeliads including pineapple. Here I show MIR828 exists in pineapple and directs post-transcriptional gene silencing of mRNAs encoding MYB family members with inferred function to regulate the conspicuous red fruit trait in var. bracteatus. MIR828 plesiomorphy (an ancient basal trait) may shed light on monocot apomorphic fruit development, postulated for 21 monocot families with fleshy fruits as due to homoplasy/convergence driven by tropical climate and/or enticements to vertebrate endozoic seed dispersers.

The role of cis-elements in the evolution of crassulacean acid metabolism photosynthesis

Crassulacean acid metabolism (CAM) photosynthesis is an innovation of carbon concentrating mechanism that is characterized by nocturnal CO2 fixation. Recent progresses in genomics, transcriptomics, proteomics, and metabolomics of CAM species yielded new knowledge and abundant genomic resources. In this review, we will discuss the pattern of cis-elements in stomata movement-related genes and CAM CO2 fixation genes, and analyze the expression dynamic of CAM related genes in green leaf tissues. We propose that CAM photosynthesis evolved through the re-organization of existing enzymes and associated membrane transporters in central metabolism and stomatal movement-related genes, at least in part by selection of existing circadian clock cis-regulatory elements in their promoter regions. Better understanding of CAM evolution will help us to design crops that can thrive in arid or semi-arid regions, which are likely to expand due to global climate change.