wgd-simple command line tools for the analysis of ancient whole-genome duplications

Genome of Hippophae rhamnoides provides insights into a conserved molecular mechanism in actinorhizal and rhizobial symbioses

Sea buckthorn (Hippophae rhamnoides), a horticulturally multipurpose species in the family Elaeagnaceae, can build associations with Frankia actinomycetes to enable symbiotic nitrogen-fixing. Currently, no high-quality reference genome is available for an actinorhizal plant, which greatly hinders the study of actinorhizal symbiotic nodulation. Here, by combining short-read, long-read and Hi-C sequencing technologies, we generated a chromosome-level reference genome of H. rhamnoides (scaffold N50: 65 Mb, and genome size: 730 Mb) and predicted 30 812 protein-coding genes mainly on 12 pseudochromosomes. Hippophae rhamnoides was found to share a high proportion of symbiotic nodulation genes with Medicago truncatula, implying a shared molecular mechanism between actinorhizal and rhizobial symbioses. Phylogenetic analysis clustered the three paralogous NODULE INCEPTION (NIN) genes of H. rhamnoides with those of other nodulating species, forming the NIN group that most likely evolved from the ancestral NLP group. The genome of H. rhamnoides will help us to decipher the underlying genetic programming of actinorhizal symbiosis, and our high-quality genome and transcriptomic resources will make H. rhamnoides a new excellent model plant for actinorhizal symbiosis research.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

The genome of Corydalis reveals the evolution of benzylisoquinoline alkaloid biosynthesis in Ranunculales

Species belonging to the order Ranunculales have attracted much attention because of their phylogenetic position as a sister group to all other eudicot lineages and their ability to produce unique yet diverse benzylisoquinoline alkaloids (BIAs). The Papaveraceae family in Ranunculales is often used as a model system for studying BIA biosynthesis. Here, we report the chromosome-level genome assembly of Corydalis tomentella, a species of Fumarioideae, one of the two subfamilies of Papaveraceae. Based on comparisons of sequenced Ranunculalean species, we present clear evidence of a shared whole-genome duplication (WGD) event that has occurred before the divergence of Ranunculales but after its divergence from other eudicot lineages. The C. tomentella genome enabled us to integrate isotopic labeling and comparative genomics to reconstruct the BIA biosynthetic pathway for both sanguinarine biosynthesis shared by papaveraceous species and the cavidine biosynthesis that is specific to Corydalis. Also, our comparative analysis revealed that gene duplications, especially tandem gene duplications, underlie the diversification of BIA biosynthetic pathways in Ranunculales. In particular, tandemly duplicated berberine bridge enzyme-like genes appear to be involved in cavidine biosynthesis. In conclusion, our study of the C. tomentella genome provides important insights into the occurrence of WGDs during the early evolution of eudicots, as well as into the evolution of BIA biosynthesis in Ranunculales.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA

As the focus for CRISPR/Cas-edited plants moves from proof-of-concept to real-world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one guide RNA (gRNA) per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in hybrid poplar (Populus tremula × P. alba INRA 717-1B4). Unexpected duplications of MYB186 and MYB138 resulted in eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and polymerase chain reaction analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole-leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the nonglandular trichomes of poplar.

The flying spider-monkey tree fern genome provides insights into fern evolution and arborescence

To date, little is known about the evolution of fern genomes, with only two small genomes published from the heterosporous Salviniales. Here we assembled the genome of Alsophila spinulosa, known as the flying spider-monkey tree fern, onto 69 pseudochromosomes. The remarkable preservation of synteny, despite resulting from an ancient whole-genome duplication over 100 million years ago, is unprecedented in plants and probably speaks to the uniqueness of tree ferns. Our detailed investigations into stem anatomy and lignin biosynthesis shed new light on the evolution of stem formation in tree ferns. We identified a phenolic compound, alsophilin, that is abundant in xylem, and we provided the molecular basis for its biosynthesis. Finally, analysis of demographic history revealed two genetic bottlenecks, resulting in rapid demographic declines of A. spinulosa. The A. spinulosa genome fills a crucial gap in the plant genomic landscape and helps elucidate many unique aspects of tree fern biology.

The first genome for the Cape Primrose Streptocarpus rexii (Gesneriaceae), a model plant for studying meristem-driven shoot diversity

Cape Primroses (Streptocarpus, Gesneriaceae) are an ideal study system for investigating the genetics underlying species diversity in angiosperms. Streptocarpus rexii has served as a model species for plant developmental research for over five decades due to its unusual extended meristem activity present in the leaves. In this study, we sequenced and assembled the complete nuclear, chloroplast, and mitochondrial genomes of S. rexii using Oxford Nanopore Technologies long read sequencing. Two flow cells of PromethION sequencing resulted in 32 billion reads and were sufficient to generate a draft assembly including the chloroplast, mitochondrial and nuclear genomes, spanning 776 Mbp. The final nuclear genome assembly contained 5,855 contigs, spanning 766 Mbp of the 929-Mbp haploid genome with an N50 of 3.7 Mbp and an L50 of 57 contigs. Over 70% of the draft genome was identified as repeats. A genome repeat library of Gesneriaceae was generated and used for genome annotation, with a total of 45,045 genes annotated in the S. rexii genome. Ks plots of the paranomes suggested a recent whole genome duplication event, shared between S. rexii and Primulina huaijiensis. A new chloroplast and mitochondrial genome assembly method, based on contig coverage and identification, was developed, and successfully used to assemble both organellar genomes of S. rexii. This method was developed into a pipeline and proved widely applicable. The nuclear genome of S. rexii and other datasets generated and reported here will be invaluable resources for further research to aid in the identification of genes involved in morphological variation underpinning plant diversification.

The genome of hibiscus hamabo reveals its adaptation to saline and waterlogged habitat

Hibiscus hamabo is a semi-mangrove species with strong tolerance to salt and waterlogging stress. However, the molecular basis and mechanisms that underlie this strong adaptability to harsh environments remain poorly understood. Here, we assembled a high-quality, chromosome-level genome of this semi-mangrove plant and analyzed its transcriptome under different stress treatments to reveal regulatory responses and mechanisms. Our analyses suggested that H. hamabo has undergone two recent successive polyploidy events, a whole-genome duplication followed by a whole-genome triplication, resulting in an unusually large gene number (107 309 genes). Comparison of the H. hamabo genome with that of its close relative Hibiscus cannabinus, which has not experienced a recent WGT, indicated that genes associated with high stress resistance have been preferentially preserved in the H. hamabo genome, suggesting an underlying association between polyploidy and stronger stress resistance. Transcriptomic data indicated that genes in the roots and leaves responded differently to stress. In roots, genes that regulate ion channels involved in biosynthetic and metabolic processes responded quickly to adjust the ion concentration and provide metabolic products to protect root cells, whereas no such rapid response was observed from genes in leaves. Using co-expression networks, potential stress resistance genes were identified for use in future functional investigations. The genome sequence, along with several transcriptome datasets, provide insights into genome evolution and the mechanism of salt and waterlogging tolerance in H. hamabo, suggesting the importance of polyploidization for environmental adaptation.

Phylotranscriptomics of Theaceae: generic-level relationships, reticulation and whole-genome duplication

BACKGROUND AND AIMS

Theaceae, with three tribes, nine genera and more than 200 species, are of great economic and ecological importance. Recent phylogenetic analyses based on plastomic data resolved the relationships among the three tribes and the intergeneric relationships within two of those tribes. However, generic-level relationships within the largest tribe, Theeae, were not fully resolved. The role of putative whole-genome duplication (WGD) events in the family and possible hybridization events among genera within Theeae also remain to be tested further.

METHODS

Transcriptomes or low-depth whole-genome sequencing of 57 species of Theaceae, as well as additional plastome sequence data, were generated. Using a dataset of low-copy nuclear genes, we reconstructed phylogenetic relationships using concatenated, species tree and phylogenetic network approaches. We further conducted molecular dating analyses and inferred possible WGD events by examining the distribution of the number of synonymous substitutions per synonymous site (Ks) for paralogues in each species. For plastid protein-coding sequences , phylogenies were reconstructed for comparison with the results obtained from analysis of the nuclear dataset.

RESULTS

Based on the 610 low-copy nuclear genes (858 606 bp in length) investigated, Stewartieae was resolved as sister to the other two tribes. Within Theeae, the Apterosperma-Laplacea clade grouped with Pyrenaria, leaving Camellia and Polyspora as sister. The estimated ages within Theaceae were largely consistent with previous studies based mainly on plastome data. Two reticulation events within Camellia and one between the common ancestor of Gordonia and Schima were found. All members of the tea family shared two WGD events, an older At-γ and a recent Ad-β; both events were also shared with the outgroups (Diapensiaceae, Pentaphylacaceae, Styracaceae and Symplocaceae).

CONCLUSIONS

Our analyses using low-copy nuclear genes improved understanding of phylogenetic relationships at the tribal and generic levels previously proposed based on plastome data, but the phylogenetic position of the Apterosperma-Laplacea clade needs more attention. There is no evidence for extensive intergeneric hybridization within Theeae or for a Theaceae-specific WGD event. Land bridges (e.g. the Bering land bridge) during the Late Oligocene may have permitted the intercontinental plant movements that facilitated the putative ancient introgression between the common ancestor of Gordonia and Schima.

Comparative Genome Analysis Across 128 Phytophthora Isolates Reveal Species-Specific Microsatellite Distribution and Localized Evolution of Compartmentalized Genomes

Phytophthora sp. are invasive groups of pathogens belonging to class Oomycetes. In order to contain and control them, a deep knowledge of their biology and infection strategy is imperative. With the availability of large-scale sequencing data, it has been possible to look directly into their genetic material and understand the strategies adopted by them for becoming successful pathogens. Here, we have studied the genomes of 128 Phytophthora species available publicly with reasonable quality. Our analysis reveals that the simple sequence repeats (SSRs) of all Phytophthora sp. follow distinct isolate specific patterns. We further show that TG/CA dinucleotide repeats are far more abundant in Phytophthora sp. than other classes of repeats. In case of tri- and tetranucleotide SSRs also, TG/CA-containing motifs always dominate over others. The GC content of the SSRs are stable without much variation across the isolates of Phytophthora. Telomeric repeats of Phytophthora follow a pattern of (TTTAGGG)n or (TTAGGGT)n rather than the canonical (TTAGGG)n. RxLR (arginine-any amino acid-leucine-arginine) motifs containing effectors diverge rapidly in Phytophthora and do not show any core common group. The RxLR effectors of some Phytophthora isolates have a tendency to form clusters with RxLRs from other species than within the same species. An analysis of the flanking intergenic distance clearly indicates a two-speed genome organization for all the Phytophthora isolates. Apart from effectors and the transposons, a large number of other virulence genes such as carbohydrate-active enzymes (CAZymes), transcriptional regulators, signal transduction genes, ATP-binding cassette transporters (ABC), and ubiquitins are also present in the repeat-rich compartments. This indicates a rapid co-evolution of this powerful arsenal for successful pathogenicity. Whole genome duplication studies indicate that the pattern followed is more specific to a geographic location. To conclude, the large-scale genomic studies of Phytophthora have thrown light on their adaptive evolution, which is largely guided by the localized host-mediated selection pressure.

The chromosome-scale genomes of Dipterocarpus turbinatus and Hopea hainanensis (Dipterocarpaceae) provide insights into fragrant oleoresin biosynthesis and hardwood formation

Dipterocarpaceae are typical tropical plants (dipterocarp forests) that are famous for their high economic value because of their production of fragrant oleoresins, top-quality timber and usage in traditional Chinese medicine. Currently, the lack of Dipterocarpaceae genomes has been a limiting factor to decipher the fragrant oleoresin biosynthesis and gain evolutionary insights into high-quality wood formation in Dipterocarpaceae. We generated chromosome-level genome assemblies for two representative Dipterocarpaceae species viz. Dipterocarpus turbinatus Gaertn. f. and Hopea hainanensis Merr. et Chun. Our whole-genome duplication (WGD) analysis revealed that Dipterocarpaceae underwent a shared WGD event, which showed significant impacts on increased copy numbers of genes related to the biosynthesis of terpene, BAHD acyltransferases, fatty acid and benzenoid/phenylpropanoid, which probably confer to the formation of their characteristic fragrant oleoresin. Additionally, compared with common soft wood plants, the expansion of gene families was also found to be associated with wood formation, such as in CESA (cellulose synthase), CSLE (cellulose synthase-like protein E), laccase and peroxidase in Dipterocarpaceae genomes, which might also contribute to the formation of harder, stronger and high-density timbers. Finally, an integrative analysis on a combination of genomic, transcriptomic and metabolic data from different tissues provided further insights into the molecular basis of fragrant oleoresins biosynthesis and high-quality wood formation of Dipterocarpaceae. Our study contributes the first two representative genomes for Dipterocarpaceae, which are valuable genetic resources for further researches on the fragrant oleoresins and superior-quality timber, genome-assisted breeding and improvement, and conservation biology of this family.

Chromosome-level genome assembly of Aristolochia contorta provides insights into the biosynthesis of benzylisoquinoline alkaloids and aristolochic acids

Aristolochic acids (AAs) and their derivatives exist in multiple Aristolochiaceae species which had been or are being used as medicinal materials. During the past decades, AAs have received increasing attention due to their nephrotoxicity and carcinogenecity. Elimination of AAs in medicinal materials using biotechnological approaches is important to improve medication safety. However, it has not been achieved because of the limited information of AA biosynthesis available. Here, we report a high-quality reference-grade genome assembly of the AA-containing vine, Aristolochia contorta. Total size of the assembly is 209.27 Mb, which is assembled into 7 pseudochromosomes. Synteny analysis, Ks distribution and 4DTv suggest absences of whole-genome duplication events in A. contorta after the angiosperm-wide WGD. Based on genomic, transcriptomic and metabolic data, pathways and candidate genes of benzylisoquinoline alkaloid (BIA) and AA biosynthesis in A. contorta were proposed. Five O-methyltransferase genes, including AcOMT1-3, AcOMT5 and AcOMT7, were cloned and functionally characterized. The results provide a high-quality reference genome for AA-containing species of Aristolochiaceae. It lays a solid foundation for further elucidation of AA biosynthesis and regulation and molecular breeding of Aristolochiaceae medicinal materials.

Chromosome-level genome assembly of Zizania latifolia provides insights into its seed shattering and phytocassane biosynthesis

Chinese wild rice (Zizania latifolia; family: Gramineae) is a valuable medicinal homologous grain in East and Southeast Asia. Here, using Nanopore sequencing and Hi-C scaffolding, we generated a 547.38 Mb chromosome-level genome assembly comprising 332 contigs and 164 scaffolds (contig N50 = 4.48 Mb; scaffold N50 = 32.79 Mb). The genome harbors 38,852 genes, with 52.89% of the genome comprising repetitive sequences. Phylogenetic analyses revealed close relation of Z. latifolia to Leersia perrieri and Oryza species, with a divergence time of 19.7-31.0 million years. Collinearity and transcriptome analyses revealed candidate genes related to seed shattering, providing basic information on abscission layer formation and degradation in Z. latifolia. Moreover, two genomic blocks in the Z. latifolia genome showed good synteny with the rice phytocassane biosynthetic gene cluster. The updated genome will support future studies on the genetic improvement of Chinese wild rice and comparative analyses between Z. latifolia and other plants.

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.

Genome sequence and transcriptome of Sorbus pohuashanensis provides insights into population evolution and leaf sunburn response

Sorbus pohuashanensis is a potential horticulture and medicinal plant, but its genomic and genetic backgrounds remain unknown. Here, we sequenced and assembled the S. pohuashanensis (Hance) Hedl. reference genome using PacBio long reads. Based on the new reference genome, we resequenced a core collection of 22 Sorbus spp. samples, which were divided into 2 groups (G1 and G2) based on phylogenetic and PCA analyses. These phylogenetic clusters were highly consistent with their classification based on leaf shape. Natural hybridization between the G1 and G2 groups was evidenced by a sample (R21) with a highly heterozygous genotype. Nucleotide diversity (π) analysis showed that G1 had a higher diversity than G2 and that G2 originated from G1. During the evolution process, the gene families involved in photosynthesis pathways expanded and the gene families involved in energy consumption contracted. RNA-seq data suggested that flavonoid biosynthesis and heat-shock protein (HSP)-heat-shock factor (HSF) pathways play important roles in protection against sunburn. This study provides new insights into the evolution of Sorbus spp. genomes. In addition, the genomic resources, and the identified genetic variations, especially those related to stress resistance, will help future efforts to produce and breed Sorbus spp.

ksrates: positioning whole-genome duplications relative to speciation events in KS distributions

SUMMARY

We present ksrates, a user-friendly command-line tool to position ancient whole-genome duplication events with respect to speciation events in a phylogeny by comparing paralog and ortholog KS distributions derived from genomic or transcriptomic sequences, while adjusting for substitution rate differences among the lineages involved.

AVAILABILITY AND IMPLEMENTATION

ksrates is implemented in Python 3 and as a Nextflow pipeline. The source code, Singularity and Docker containers, documentation and tutorial are available via https://github.com/VIB-PSB/ksrates.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Antarctic Moss Pohlia nutans Genome Provides Insights Into the Evolution of Bryophytes and the Adaptation to Extreme Terrestrial Habitats

The Antarctic continent has extreme natural environment and fragile ecosystem. Mosses are one of the dominant floras in the Antarctic continent. However, their genomic features and adaptation processes to extreme environments remain poorly understood. Here, we assembled the high-quality genome sequence of the Antarctic moss (Pohlia nutans) with 698.20 Mb and 22 chromosomes. We found that the high proportion of repeat sequences and a recent whole-genome duplication (WGD) contribute to the large size genome of P. nutans when compared to other bryophytes. The genome of P. nutans harbors the signatures of massive segmental gene duplications and large expansions of gene families, likely facilitating neofunctionalization. Genomic characteristics that may support the Antarctic lifestyle of this moss comprise expanded gene families involved in phenylpropanoid biosynthesis, unsaturated fatty acid biosynthesis, and plant hormone signal transduction. Additional contributions include the significant expansion and upregulation of several genes encoding DNA photolyase, antioxidant enzymes, flavonoid biosynthesis enzymes, possibly reflecting diverse adaptive strategies. Notably, integrated multi-omic analyses elucidate flavonoid biosynthesis may function as the reactive oxygen species detoxification under UV-B radiation. Our studies provide insight into the unique features of the Antarctic moss genome and their molecular responses to extreme terrestrial environments.

A chromosome-level genome assembly of Amorphophallus konjac provides insights into konjac glucomannan biosynthesis

Amorphophallus konjac, a perennial herb in the Araceae family, is a cash crop that can produce a large amount of konjac glucomannan. To explore mechanisms underlying such large genomes in the genus Amorphophallus as well as the gene regulation of glucomannan biosynthesis, we present a chromosome-level genome assembly of A. konjac with a total genome size of 5.60 Gb and a contig N50 of 1.20 Mb. Comparative genomic analysis reveals that A. konjac has undergone two whole-genome duplication (WGD) events in quick succession. Two recent bursts of transposable elements are identified in the A. konjac genome, which contribute greatly to the large genome size. Our transcriptomic analysis of the developmental corms characterizes key genes involved in the biosynthesis of glucomannan and related starches. High expression of cellulose synthase-like A, Cellulose synthase-like D, mannan-synthesis related 1, GDP-mannose pyrophosphorylase and phosphomannomutase fructokinase contributes to glucomannan synthesis during the corm expansion period while high expression of starch synthase, starch branching enzyme and phosphoglucomutase is responsible for starch synthesis in the late corm development stage. In conclusion, we generate a high-quality genome of A. konjac with different sequencing technologies. The expansion of transposable elements has caused the large genome of this species. And the identified key genes in the glucomannan biosynthesis provide valuable candidates for molecular breeding of this crop in the future.

The genome of the endangered Macadamia jansenii displays little diversity but represents an important genetic resource for plant breeding

Macadamia, a recently domesticated expanding nut crop in the tropical and subtropical regions of the world, is one of the most economically important genera in the diverse and widely adapted Proteaceae family. All four species of Macadamia are rare in the wild with the most recently discovered, M. jansenii, being endangered. The M. jansenii genome has been used as a model for testing sequencing methods using a wide range of long read sequencing techniques. Here, we report a chromosome level genome assembly, generated using a combination of Pacific Biosciences sequencing and Hi-C, comprising 14 pseudo-molecules, with a N50 of 52 Mb and a total genome assembly size of 758 Mb of which 56% is repetitive. Completeness assessment revealed that the assembly covered -97.1% of the conserved single copy genes. Annotation predicted 31,591 protein coding genes and allowed the characterization of genes encoding biosynthesis of cyanogenic glycosides, fatty acid metabolism, and anti-microbial proteins. Re-sequencing of seven other genotypes confirmed low diversity and low heterozygosity within this endangered species. Important morphological characteristics of this species such as small tree size and high kernel recovery suggest that M. jansenii is an important source of these commercial traits for breeding. As a member of a small group of families that are sister to the core eudicots, this high-quality genome also provides a key resource for evolutionary and comparative genomics studies.

The genome of the endangered Macadamia jansenii displays little diversity but represents an important genetic resource for plant breeding

Macadamia, a recently domesticated expanding nut crop in the tropical and subtropical regions of the world, is one of the most economically important genera in the diverse and widely adapted Proteaceae family. All four species of Macadamia are rare in the wild with the most recently discovered, M. jansenii, being endangered. The M. jansenii genome has been used as a model for testing sequencing methods using a wide range of long read sequencing techniques. Here, we report a chromosome level genome assembly, generated using a combination of Pacific Biosciences sequencing and Hi-C, comprising 14 pseudo-molecules, with a N50 of 52 Mb and a total genome assembly size of 758 Mb of which 56% is repetitive. Completeness assessment revealed that the assembly covered -97.1% of the conserved single copy genes. Annotation predicted 31,591 protein coding genes and allowed the characterization of genes encoding biosynthesis of cyanogenic glycosides, fatty acid metabolism, and anti-microbial proteins. Re-sequencing of seven other genotypes confirmed low diversity and low heterozygosity within this endangered species. Important morphological characteristics of this species such as small tree size and high kernel recovery suggest that M. jansenii is an important source of these commercial traits for breeding. As a member of a small group of families that are sister to the core eudicots, this high-quality genome also provides a key resource for evolutionary and comparative genomics studies.

Chromosome-scale assembly reveals asymmetric paleo-subgenome evolution and targets for the acceleration of fungal resistance breeding in the nut crop, pecan

Pecan (Carya illinoinensis) is a tree nut crop of worldwide economic importance that is rich in health-promoting factors. However, pecan production and nut quality are greatly challenged by environmental stresses such as the outbreak of severe fungal diseases. Here, we report a high-quality, chromosome-scale genome assembly of the controlled-cross pecan cultivar 'Pawnee' constructed by integrating Nanopore sequencing and Hi-C technologies. Phylogenetic and evolutionary analyses reveal two whole-genome duplication (WGD) events and two paleo-subgenomes in pecan and walnut. Time estimates suggest that the recent WGD event and considerable genome rearrangements in pecan and walnut account for expansions in genome size and chromosome number after the divergence from bayberry. The two paleo-subgenomes differ in size and protein-coding gene sets. They exhibit uneven ancient gene loss, asymmetrical distribution of transposable elements (especially LTR/Copia and LTR/Gypsy), and expansions in transcription factor families (such as the extreme pecan-specific expansion in the far-red impaired response 1 family), which are likely to reflect the long evolutionary history of species in the Juglandaceae. A whole-genome scan of resequencing data from 86 pecan scab-associated core accessions identified 47 chromosome regions containing 185 putative candidate genes. Significant changes were detected in the expression of candidate genes associated with the chitin response pathway under chitin treatment in the scab-resistant and scab-susceptible cultivars 'Excell' and 'Pawnee'. These findings enable us to identify key genes that may be important susceptibility factors for fungal diseases in pecan. The high-quality sequences are valuable resources for pecan breeders and will provide a foundation for the production and quality improvement of tree nut crops.

Underwater CAM photosynthesis elucidated by Isoetes genome

To conserve water in arid environments, numerous plant lineages have independently evolved Crassulacean Acid Metabolism (CAM). Interestingly, Isoetes, an aquatic lycophyte, can also perform CAM as an adaptation to low CO2 availability underwater. However, little is known about the evolution of CAM in aquatic plants and the lack of genomic data has hindered comparison between aquatic and terrestrial CAM. Here, we investigate underwater CAM in Isoetes taiwanensis by generating a high-quality genome assembly and RNA-seq time course. Despite broad similarities between CAM in Isoetes and terrestrial angiosperms, we identify several key differences. Notably, Isoetes may have recruited the lesser-known 'bacterial-type' PEPC, along with the 'plant-type' exclusively used in other CAM and C4 plants for carboxylation of PEP. Furthermore, we find that circadian control of key CAM pathway genes has diverged considerably in Isoetes relative to flowering plants. This suggests the existence of more evolutionary paths to CAM than previously recognized.

Rambutan genome revealed gene networks for spine formation and aril development

Rambutan is a popular tropical fruit known for its exotic appearance, has long flexible spines on shells, extraordinary aril growth, desirable nutrition, and a favorable taste. The genome of an elite rambutan cultivar Baoyan 7 was assembled into 328 Mb in 16 pseudo-chromosomes. Comparative genomics analysis between rambutan and lychee revealed that rambutan chromosomes 8 and 12 are collinear with lychee chromosome 1, which resulted in a chromosome fission event in rambutan (n = 16) or a fusion event in lychee (n = 15) after their divergence from a common ancestor 15.7 million years ago. Root development genes played a crucial role in spine development, such as endoplasmic reticulum pathway genes, jasmonic acid response genes, vascular bundle development genes, and K⁺ transport genes. Aril development was regulated by D-class genes (STK and SHP1), plant hormone and phenylpropanoid biosynthesis genes, and sugar metabolism genes. The lower rate of male sterility of hermaphroditic flowers appears to be regulated by MYB24. Population genomic analyses revealed genes in selective sweeps during domestication that are related to fruit morphology and environment stress response. These findings enhance our understanding of spine and aril development and provide genomic resources for rambutan improvement.

Improved chromosome-level genome assembly and annotation of the seagrass, Zostera marina (eelgrass)

Background: Seagrasses (Alismatales) are the only fully marine angiosperms.  Zostera marina (eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system for exploring adaptation under rapid climate change. The original draft genome (v.1.0) of the seagrass  Z. marina (L.) was based on a combination of Illumina mate-pair libraries and fosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7× genomic coverage. The assembly resulted in ~2000 unordered scaffolds (L50 of 486 Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here, we present an upgraded chromosome-scale genome assembly and compare v.1.0 and the new v.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings.  

Methods: The same high molecular weight DNA used in the original sequencing of the Finnish clone was used. A high-quality reference genome was assembled with the MECAT assembly pipeline combining PacBio long-read sequencing and Hi-C scaffolding. 

Results: In total, 75.97 Gb PacBio data was produced. The final assembly comprises six pseudo-chromosomes and 304 unanchored scaffolds with a total length of 260.5Mb and an N50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%) obtained at least one functional assignment based on similarity to known proteins. 

Conclusions: As an important marine angiosperm, the improved  Z. marina genome assembly will further assist evolutionary, ecological, and comparative genomics at the chromosome level. The new genome assembly will further our understanding into the structural and physiological adaptations from land to marine life.

The chromosome-scale genome of Magnolia officinalis provides insight into the evolutionary position of magnoliids

Magnolia officinalis, a representative tall aromatic tree of the Magnoliaceae family, is a medicinal plant that is widely used in diverse industries from medicine to cosmetics. We report a chromosome-scale draft genome of M. officinalis, in which ∼99.66% of the sequences were anchored onto 19 chromosomes with the scaffold N50 of 76.62 Mb. We found that a high proportion of repetitive sequences was a common feature of three Magnoliaceae with known genomic data. Magnoliids were a sister clade to eudicots-monocots, which provided more support for understanding the phylogenetic position among angiosperms. An ancient duplication event occurred in the genome of M. officinalis and was shared with Lauraceae. Based on RNA-seq analysis, we identified several key enzyme-coding gene families associated with the biosynthesis of lignans in the genome. The construction of the M. officinalis genome sequence will serve as a reference for further studies of Magnolia, as well as other Magnoliaceae.

The Clausena lansium (Wampee) genome reveal new insights into the carbazole alkaloids biosynthesis pathway

Clausena lansium (Lour.) Skeels (Rutaceae), recognized as wampee, is a widely distributed fruit tree which is utilized as a folk-medicine for treatment of several common diseases. However, the genomic information about this medicinally important species is still lacking. Therefore, we assembled the first genome of Clausena genus with a total length of 310.51 Mb and scaffold N50 of 2.24 Mb by using 10× Genomics technology. Further annotation revealed a total of 34,419 protein-coding genes, while repetitive elements covered 39.08% (121.36 Mb) of the genome. The Clausena and Citrus genus were found to diverge around 22 MYA, and also shared an ancient whole-genome triplication event with Vitis. Furthermore, multi-tissue transcriptomic analysis enabled the identification of genes involved in the synthesis of carbazole alkaloids. Altogether, these findings provided new insights into the genome evolution of Wampee species and highlighted the possible role of key genes involved in the carbazole alkaloids biosynthetic pathway.

Whole-genome sequencing of brown-marbled grouper (Epinephelus fuscoguttatus) provides insights into adaptive evolution and growth differences

The brown-marbled grouper (Epinephelus fuscoguttatus) is an important species of fish in the coral reef ecosystem and marine aquaculture industry. In this study, a high-quality chromosome-level genome of brown-marbled grouper was assembled using Oxford Nanopore technology and Hi-C technology. The GC content and heterozygosity were approximately 42% and 0.35%, respectively. A total of 230 contigs with a total length of 1047 Mb and contig N50 of 13.8 Mb were assembled, and 228 contigs (99.13%) were anchored into 24 chromosomes. A total of 24,005 protein-coding genes were predicted, among which 23,862 (99.4%) predicted genes were annotated. Phylogenetic analysis showed that brown-marbled grouper and humpback grouper were clustered into one clade that separated approximately 11-23 million years ago. Collinearity analyses showed that there was no obvious duplication of large fragments between chromosomes in the brown-marbled grouper. Genomes of the humpback grouper and giant grouper showed a high collinearity with that of the brown-marbled grouper. A total of 305 expanded gene families were detected in the brown-marbled grouper genome, which is mainly involved in disease resistance. In addition, a genetic linkage map with 3061.88 cM was constructed. Based on the physical and genetic map, one growth-related quantitative trait loci was detected in 32,332,447 bp of chromosome 20, and meox1 and etv4 were considered candidate growth-related genes. This study provides pivotal genetic resources for further evolutionary analyses and artificial breeding of groupers.

Chromosome-level genome assembly of Welwitschia mirabilis, a unique Namib Desert species

Welwitschia mirabilis, which is endemic to the Namib Desert, is the only living species within the family Welwitschiaceae. This species has an extremely long lifespan of up to 2,000 years and bears a single pair of opposite leaves that persist whilst alive. However, the underlying genetic mechanisms and evolution of the species remain poorly elucidated. Here, we report on a chromosome-level genome assembly for W. mirabilis, with a 6.30-Gb genome sequence and contig N50 of 27.50 Mb. In total, 39,019 protein-coding genes were predicted from the genome. Two brassinosteroid-related genes (BRI1 and CYCD3), key regulators of cell division and elongation, were strongly selected in W. mirabilis and may contribute to their long ever-growing leaves. Furthermore, 29 gene families in the mitogen-activated protein kinase signalling pathway showed significant expansion, which may contribute to the desert adaptations of the plant. Three positively selected genes (EHMT1, EIF4E, SOD2) may be involved in the mechanisms leading to long lifespan. Based on molecular clock dating and fossil calibrations, the divergence time of W. mirabilis and Gnetum montanum was estimated at ~123.5 million years ago. Reconstruction of population dynamics from genome data coincided well with the aridification of the Namib Desert. The genome sequence detailed in the current study provides insight into the evolution of W. mirabilis and should be an important resource for further study on gnetophyte and gymnosperm evolution.

Chromosomal-level assembly of the Leptodermis oblonga (Rubiaceae) genome and its phylogenetic implications

Rubiaceae is the fourth largest and a taxonomically complex family of angiosperms. Many species in this family harbor low reproductive isolation and frequently exhibit inconsistent phenotypic characteristics. Therefore, taxonomic classification and their phylogenetic relationships in the Rubiaceae family is challenging, especially in the genus Leptodermis. Considering the low taxonomic confusion and wide distribution, Leptodermis oblonga is selected as a representative Leptodermis for genome sequencing. The assemblies resulted in 497 Mbp nuclear and 155,100 bp chloroplast genomes, respectively. Using the nuclear genome as a reference, SNPs were called from 37 Leptodermis species or varieties. The phylogenetic tree based on SNPs exhibited high resolution for species delimitation of the complex and well-resolved phylogenetic relationships in the genus. Moreover, 28,987 genes were predicted in the nuclear genome and used for comparative genomics study. As the first chromosomal-level genome of the subfamily Rubioideae in Rubiaceae, it will provide fruitfully evolutionary understanding in the family.

A chromosome-level Amaranthus cruentus genome assembly highlights gene family evolution and biosynthetic gene clusters that may underpin the nutritional value of this traditional crop

Traditional crops have historically provided accessible and affordable nutrition to millions of rural dwellers but have been neglected, with most modern agricultural systems over-reliant on a small number of internationally traded crops. Traditional crops are typically well-adapted to local agro-ecological conditions and many are nutrient-dense. They can play a vital role in local food systems through enhanced nutrition (particularly where diets are dominated by starch crops), food security and livelihoods for smallholder farmers, and a climate-resilient and biodiverse agriculture. Using short-read, long-read and phased sequencing technologies, we generated a high-quality chromosome-level genome assembly for Amaranthus cruentus, an under-researched crop with micronutrient- and protein-rich leaves and gluten-free seed, but lacking improved varieties, with respect to productivity and quality traits. The 370.9 Mb genome demonstrates a shared whole genome duplication with a related species, Amaranthus hypochondriacus. Comparative genome analysis indicates chromosomal loss and fusion events following genome duplication that are common to both species, as well as fission of chromosome 2 in A. cruentus alone, giving rise to a haploid chromosome number of 17 (versus 16 in A. hypochondriacus). Genomic features potentially underlying the nutritional value of this crop include two A. cruentus-specific genes with a likely role in phytic acid synthesis (an anti-nutrient), expansion of ion transporter gene families, and identification of biosynthetic gene clusters conserved within the amaranth lineage. The A. cruentus genome assembly will underpin much-needed research and global breeding efforts to develop improved varieties for economically viable cultivation and realization of the benefits to global nutrition security and agrobiodiversity.

A chromosome-level genome assembly of rugged rose (Rosa rugosa) provides insights into its evolution, ecology, and floral characteristics

Rosa rugosa, commonly known as rugged rose, is a perennial ornamental shrub. It produces beautiful flowers with a mild fragrance and colorful seed pods. Unlike many other cultivated roses, R. rugosa adapts to a wide range of habitat types and harsh environmental conditions such as salinity, alkaline, shade, drought, high humidity, and frigid temperatures. Here, we produced and analyzed a high-quality genome sequence for R. rugosa to understand its ecology, floral characteristics and evolution. PacBio HiFi reads were initially used to construct the draft genome of R. rugosa, and then Hi-C sequencing was applied to assemble the contigs into 7 chromosomes. We obtained a 382.6 Mb genome encoding 39,704 protein-coding genes. The genome of R. rugosa appears to be conserved with no additional whole-genome duplication after the gamma whole-genome triplication (WGT), which occurred ~100 million years ago in the ancestor of core eudicots. Based on a comparative analysis of the high-quality genome assembly of R. rugosa and other high-quality Rosaceae genomes, we found a unique large inverted segment in the Chinese rose R. chinensis and a retroposition in strawberry caused by post-WGT events. We also found that floral development- and stress response signaling-related gene modules were retained after the WGT. Two MADS-box genes involved in floral development and the stress-related transcription factors DREB2A-INTERACTING PROTEIN 2 (DRIP2) and PEPTIDE TRANSPORTER 3 (PTR3) were found to be positively selected in evolution, which may have contributed to the unique ability of this plant to adapt to harsh environments. In summary, the high-quality genome sequence of R. rugosa provides a map for genetic studies and molecular breeding of this plant and enables comparative genomic studies of Rosa in the near future.

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.

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.

Gene-rich UV sex chromosomes harbor conserved regulators of sexual development

Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants.

The Pharus latifolius genome bridges the gap of early grass evolution

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

Linked by Ancestral Bonds: Multiple Whole-Genome Duplications and Reticulate Evolution in a Brassicaceae Tribe

Pervasive hybridization and whole-genome duplications (WGDs) influenced genome evolution in several eukaryotic lineages. Although frequent and recurrent hybridizations may result in reticulate phylogenies, the evolutionary events underlying these reticulations, including detailed structure of the ancestral diploid and polyploid genomes, were only rarely reconstructed. Here, we elucidate the complex genomic history of a monophyletic clade from the mustard family (Brassicaceae), showing contentious relationships to the early-diverging clades of this model plant family. Genome evolution in the crucifer tribe Biscutelleae (∼60 species, 5 genera) was dominated by pervasive hybridizations and subsequent genome duplications. Diversification of an ancestral diploid genome into several divergent but crossable genomes was followed by hybridizations between these genomes. Whereas a single genus (Megadenia) remained diploid, the four remaining genera originated by allopolyploidy (Biscutella, Lunaria, Ricotia) or autopolyploidy (Heldreichia). The contentious relationships among the Biscutelleae genera, and between the tribe and other early diverged crucifer lineages, are best explained by close genomic relatedness among the recurrently hybridizing ancestral genomes. By using complementary cytogenomics and phylogenomics approaches, we demonstrate that the origin of a monophyletic plant clade can be more complex than a parsimonious assumption of a single WGD spurring postpolyploid cladogenesis. Instead, recurrent hybridization among the same and/or closely related parental genomes may phylogenetically interlink diploid and polyploid genomes despite the incidence of multiple independent WGDs. Our results provide new insights into evolution of early-diverging Brassicaceae lineages and elucidate challenges in resolving the contentious relationships within and between land plant lineages with pervasive hybridization and WGDs.

A chromosome-level genome assembly provides insights into ascorbic acid accumulation and fruit softening in guava (Psidium guajava)

Guava (Psidium guajava) is an important fleshy-fruited tree of the Myrtaceae family that is widely cultivated in tropical and subtropical areas of the world and has attracted considerable attention for the richness of ascorbic acid in its fruits. However, studies on the evolution and genetic breeding potential of guava are hindered by the lack of a reference genome. Here, we present a chromosome-level genomic assembly of guava using PacBio sequencing and Hi-C technology. We found that the genome assembly size was 443.8 Mb with a contig N50 of ~15.8 Mb. We annotated a total of 25 601 genes and 193.2 Mb of repetitive sequences for this genome. Comparative genomic analysis revealed that guava has undergone a recent whole-genome duplication (WGD) event shared by all species in Myrtaceae. In addition, through metabolic analysis, we determined that the L-galactose pathway plays a major role in ascorbic acid biosynthesis in guava fruits. Moreover, the softening of fruits of guava may result from both starch and cell wall degradation according to analyses of gene expression profiles and positively selected genes. Our data provide a foundational resource to support molecular breeding of guava and represent new insights into the evolution of soft, fleshy fruits in Myrtaceae.

SWEET Transporters and the Potential Functions of These Sequences in Tea (Camellia sinensis)

Tea (Camellia sinensis) is an important economic beverage crop. Its flowers and leaves could be used as healthcare tea for its medicinal value. SWEET proteins were recently identified in plants as sugar transporters, which participate in diverse physiological processes, including pathogen nutrition, seed filling, nectar secretion, and phloem loading. Although SWEET genes have been characterized and identified in model plants, such as Arabidopsis thaliana and Oryza sativa, there is very little knowledge of these genes in C. sinensis. In this study, 28 CsSWEETs were identified in C. sinensis and further phylogenetically divided into four subfamilies with A. thaliana. These identified CsSWEETs contained seven transmembrane helixes (TMHs) which were generated by an ancestral three-TMH unit with an internal duplication experience. Microsynteny analysis revealed that the large-scale duplication events were the main driving forces for members from CsSWEET family expansion in C. sinensis. The expression profiles of the 28 CsSWEETs revealed that some genes were highly expressed in reproductive tissues. Among them, CsSWEET1a might play crucial roles in the efflux of sucrose, and CsSWEET17b could control fructose content as a hexose transporter in C. sinensis. Remarkably, CsSWEET12 and CsSWEET17c were specifically expressed in flowers, indicating that these two genes might be involved in sugar transport during flower development. The expression patterns of all CsSWEETs were differentially regulated under cold and drought treatments. This work provided a systematic understanding of the members from the CsSWEET gene family, which would be helpful for further functional studies of CsSWEETs in C. sinensis.

Datura genome reveals duplications of psychoactive alkaloid biosynthetic genes and high mutation rate following tissue culture

BACKGROUND

Datura stramonium (Jimsonweed) is a medicinally and pharmaceutically important plant in the nightshade family (Solanaceae) known for its production of various toxic, hallucinogenic, and therapeutic tropane alkaloids. Recently, we published a tissue-culture based transformation protocol for D. stramonium that enables more thorough functional genomics studies of this plant. However, the tissue culture process can lead to undesirable phenotypic and genomic consequences independent of the transgene used. Here, we have assembled and annotated a draft genome of D. stramonium with a focus on tropane alkaloid biosynthetic genes. We then use mRNA sequencing and genome resequencing of transformants to characterize changes following tissue culture.

RESULTS

Our draft assembly conforms to the expected 2 gigabasepair haploid genome size of this plant and achieved a BUSCO score of 94.7% complete, single-copy genes. The repetitive content of the genome is 61%, with Gypsy-type retrotransposons accounting for half of this. Our gene annotation estimates the number of protein-coding genes at 52,149 and shows evidence of duplications in two key alkaloid biosynthetic genes, tropinone reductase I and hyoscyamine 6 β-hydroxylase. Following tissue culture, we detected only 186 differentially expressed genes, but were unable to correlate these changes in expression with either polymorphisms from resequencing or positional effects of transposons.

CONCLUSIONS

We have assembled, annotated, and characterized the first draft genome for this important model plant species. Using this resource, we show duplications of genes leading to the synthesis of the medicinally important alkaloid, scopolamine. Our results also demonstrate that following tissue culture, mutation rates of transformed plants are quite high (1.16 × 10- 3 mutations per site), but do not have a drastic impact on gene expression.

The chromosome-level genome of dragon fruit reveals whole-genome duplication and chromosomal co-localization of betacyanin biosynthetic genes

Dragon fruits are tropical fruits economically important for agricultural industries. As members of the family of Cactaceae, they have evolved to adapt to the arid environment. Here we report the draft genome of Hylocereus undatus, commercially known as the white-fleshed dragon fruit. The chromosomal level genome assembly contains 11 longest scaffolds corresponding to the 11 chromosomes of H. undatus. Genome annotation of H. undatus found ~29,000 protein-coding genes, similar to Carnegiea gigantea (saguaro). Whole-genome duplication (WGD) analysis revealed a WGD event in the last common ancestor of Cactaceae followed by extensive genome rearrangements. The divergence time between H. undatus and C. gigantea was estimated to be 9.18 MYA. Functional enrichment analysis of orthologous gene clusters (OGCs) in six Cactaceae plants found significantly enriched OGCs in drought resistance. Fruit flavor-related functions were overrepresented in OGCs that are significantly expanded in H. undatus. The H. undatus draft genome also enabled the discovery of carbohydrate and plant cell wall-related functional enrichment in dragon fruits treated with trypsin for a longer storage time. Lastly, genes of the betacyanin (a red-violet pigment and antioxidant with a very high concentration in dragon fruits) biosynthetic pathway were found to be co-localized on a 12 Mb region of one chromosome. The consequence may be a higher efficiency of betacyanin biosynthesis, which will need experimental validation in the future. The H. undatus draft genome will be a great resource to study various cactus plants.

The 'Tommy Atkins' mango genome reveals candidate genes for fruit quality

BACKGROUND

Mango, Mangifera indica L., an important tropical fruit crop, is grown for its sweet and aromatic fruits. Past improvement of this species has predominantly relied on chance seedlings derived from over 1000 cultivars in the Indian sub-continent with a large variation for fruit size, yield, biotic and abiotic stress resistance, and fruit quality among other traits. Historically, mango has been an orphan crop with very limited molecular information. Only recently have molecular and genomics-based analyses enabled the creation of linkage maps, transcriptomes, and diversity analysis of large collections. Additionally, the combined analysis of genomic and phenotypic information is poised to improve mango breeding efficiency.

RESULTS

This study sequenced, de novo assembled, analyzed, and annotated the genome of the monoembryonic mango cultivar 'Tommy Atkins'. The draft genome sequence was generated using NRGene de-novo Magic on high molecular weight DNA of 'Tommy Atkins', supplemented by 10X Genomics long read sequencing to improve the initial assembly. A hybrid population between 'Tommy Atkins' x 'Kensington Pride' was used to generate phased haplotype chromosomes and a highly resolved phased SNP map. The final 'Tommy Atkins' genome assembly was a consensus sequence that included 20 pseudomolecules representing the 20 chromosomes of mango and included ~ 86% of the ~ 439 Mb haploid mango genome. Skim sequencing identified ~ 3.3 M SNPs using the 'Tommy Atkins' x 'Kensington Pride' mapping population. Repeat masking identified 26,616 genes with a median length of 3348 bp. A whole genome duplication analysis revealed an ancestral 65 MYA polyploidization event shared with Anacardium occidentale. Two regions, one on LG4 and one on LG7 containing 28 candidate genes, were associated with the commercially important fruit size characteristic in the mapping population.

CONCLUSIONS

The availability of the complete 'Tommy Atkins' mango genome will aid global initiatives to study mango genetics.

Combined genomic, transcriptomic, and metabolomic analyses provide insights into chayote (Sechium edule) evolution and fruit development

Chayote (Sechium edule) is an agricultural crop in the Cucurbitaceae family that is rich in bioactive components. To enhance genetic research on chayote, we used Nanopore third-generation sequencing combined with Hi-C data to assemble a draft chayote genome. A chromosome-level assembly anchored on 14 chromosomes (N50 contig and scaffold sizes of 8.40 and 46.56 Mb, respectively) estimated the genome size as 606.42 Mb, which is large for the Cucurbitaceae, with 65.94% (401.08 Mb) of the genome comprising repetitive sequences; 28,237 protein-coding genes were predicted. Comparative genome analysis indicated that chayote and snake gourd diverged from sponge gourd and that a whole-genome duplication (WGD) event occurred in chayote at 25 ± 4 Mya. Transcriptional and metabolic analysis revealed genes involved in fruit texture, pigment, flavor, flavonoids, antioxidants, and plant hormones during chayote fruit development. The analysis of the genome, transcriptome, and metabolome provides insights into chayote evolution and lays the groundwork for future research on fruit and tuber development and genetic improvements in chayote.

Long read sequencing of Toona sinensis (A. Juss) Roem: A chromosome-level reference genome for the family Meliaceae

Chinese mahogany (Toona sinensis) is a woody plant that is widely cultivated in China and Malaysia. Toona sinensis is important economically, including as a nutritious food source, as material for traditional Chinese medicine and as a high-quality hardwood. However, the absence of a reference genome has hindered in-depth molecular and evolutionary studies of this plant. In this study, we report a high-quality T. sinensis genome assembly, with scaffolds anchored to 28 chromosomes and a total assembled length of 596 Mb (contig N50 = 1.5 Mb and scaffold N50 = 21.5 Mb). A total of 34,345 genes were predicted in the genome after homology-based and de novo annotation analyses. Evolutionary analysis showed that the genomes of T. sinensis and Populus trichocarpa diverged ~99.1-103.1 million years ago, and the T. sinensis genome underwent a recent genome-wide duplication event at ~7.8 million years and one more ancient whole genome duplication event at ~71.5 million years. These results provide a high-quality chromosome-level reference genome for T. sinensis and confirm its evolutionary position at the genomic level. Such information will offer genomic resources to study the molecular mechanism of terpenoid biosynthesis and the formation of flavour compounds, which will further facilitate its molecular breeding. As the first chromosome-level genome assembled in the family Meliaceae, it will provide unique insights into the evolution of members of the Meliaceae.

The genome and transcriptome analysis of snake gourd provide insights into its evolution and fruit development and ripening

Snake gourd (Trichosanthes anguina L.), which belongs to the Cucurbitaceae family, is a popular ornamental and food crop species with medicinal value and is grown in many parts of the world. Although progress has been made in its genetic improvement, the organization, composition, and evolution of the snake gourd genome remain largely unknown. Here, we report a high-quality genome assembly for snake gourd, comprising 202 contigs, with a total size of 919.8 Mb and an N50 size of 20.1 Mb. These findings indicate that snake gourd has one of the largest genomes of Cucurbitaceae species sequenced to date. The snake gourd genome assembly harbors 22,874 protein-coding genes and 80.0% of the genome consists of repetitive sequences. Phylogenetic analysis reveals that snake gourd is closely related to sponge gourd but diverged from their common ancestor ~33-47 million years ago. The genome sequence reported here serves as a valuable resource for snake gourd genetic research and comparative genomic studies in Cucurbitaceae and other plant species. In addition, fruit transcriptome analysis reveals the candidate genes related to quality traits during snake gourd fruit development and provides a basis for future research on snake gourd fruit development and ripening at the transcript level.

The genome of a cave plant, Primulina huaijiensis, provides insights into adaptation to limestone karst habitats

Although whole genome duplication (WGD) has been suggested to facilitate adaptive evolution and diversification, the role of specific WGD events in promoting diversification and adaptation in angiosperms remains poorly understood. Primulina, a species-rich genus with > 180 species associated with limestone karst habitat, constitutes an ideal system for studying the impact of WGD events on speciation and evolutionary adaptation. We sequenced and assembled a chromosome-level genome of the cave-dwelling species P. huaijiensis to study gene family expansion and gene retention following WGDs. We provide evidence that P. huaijiensis has undergone two WGDs since the γ triplication event shared by all eudicots. In addition to a WGD shared by almost all Lamiales (L event), we identified a lineage-specific WGD (D event) that occurred in the early Miocene around 20.6-24.2 Myr ago and that is shared by almost the entire subtribe Didymocarpinae. We found that gene retentions following the D event led to gene family proliferation (e.g. WRKYs) that probably facilitated adaptation to the high salinity and drought stress in limestone karst. Our study highlights the role of lineage-specific WGD in species diversification and adaptation of plants from special habitats.

The Ceratodon purpureus genome uncovers structurally complex, gene rich sex chromosomes.. [PMID: 0 DOI: 10.1101/2020.07.03.163634]

Non-recombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration1,2. The correlation between suppressed recombination and degeneration is clear in animal XY systems1,2, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions3,4. Here we test for degeneration in the bryophyte UV sex chromosome system through genomic comparisons with new female and male chromosome-scale reference genomes of the moss Ceratodon purpureus. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes show signs of weaker purifying selection than autosomes, we find suppressed recombination alone is insufficient to drive gene loss on sex-specific chromosomes. Instead, the U and V sex chromosomes harbor thousands of broadly-expressed genes, including numerous key regulators of sexual development across land plants.

Dissecting the genome of star fruit (Averrhoa carambola L.)

Averrhoa carambola is commonly known as star fruit because of its peculiar shape, and its fruit is a rich source of minerals and vitamins. It is also used in traditional medicines in countries such as India, China, the Philippines, and Brazil for treating various ailments, including fever, diarrhea, vomiting, and skin disease. Here, we present the first draft genome of the Oxalidaceae family, with an assembled genome size of 470.51 Mb. In total, 24,726 protein-coding genes were identified, and 16,490 genes were annotated using various well-known databases. The phylogenomic analysis confirmed the evolutionary position of the Oxalidaceae family. Based on the gene functional annotations, we also identified enzymes that may be involved in important nutritional pathways in the star fruit genome. Overall, the data from this first sequenced genome in the Oxalidaceae family provide an essential resource for nutritional, medicinal, and cultivational studies of the economically important star-fruit plant.

Multi-faceted analysis provides little evidence for recurrent whole-genome duplications during hexapod evolution

BACKGROUND

Gene duplication events play an important role in the evolution and adaptation of organisms. Duplicated genes can arise through different mechanisms, including whole-genome duplications (WGDs). Recently, WGD was suggested to be an important driver of evolution, also in hexapod animals.

RESULTS

Here, we analyzed 20 high-quality hexapod genomes using whole-paranome distributions of estimated synonymous distances (KS), patterns of within-genome co-linearity, and phylogenomic gene tree-species tree reconciliation methods. We observe an abundance of gene duplicates in the majority of these hexapod genomes, yet we find little evidence for WGD. The majority of gene duplicates seem to have originated through small-scale gene duplication processes. We did detect segmental duplications in six genomes, but these lacked the within-genome co-linearity signature typically associated with WGD, and the age of these duplications did not coincide with particular peaks in KS distributions. Furthermore, statistical gene tree-species tree reconciliation failed to support all but one of the previously hypothesized WGDs.

CONCLUSIONS

Our analyses therefore provide very limited evidence for WGD having played a significant role in the evolution of hexapods and suggest that alternative mechanisms drive gene duplication events in this group of animals. For instance, we propose that, along with small-scale gene duplication events, episodes of increased transposable element activity could have been an important source for gene duplicates in hexapods.

Genome and single-cell RNA-sequencing of the earthworm Eisenia andrei identifies cellular mechanisms underlying regeneration

The earthworm is particularly fascinating to biologists because of its strong regenerative capacity. However, many aspects of its regeneration in nature remain elusive. Here we report chromosome-level genome, large-scale transcriptome and single-cell RNA-sequencing data during earthworm (Eisenia andrei) regeneration. We observe expansion of LINE2 transposable elements and gene families functionally related to regeneration (for example, EGFR, epidermal growth factor receptor) particularly for genes exhibiting differential expression during earthworm regeneration. Temporal gene expression trajectories identify transcriptional regulatory factors that are potentially crucial for initiating cell proliferation and differentiation during regeneration. Furthermore, early growth response genes related to regeneration are transcriptionally activated in both the earthworm and planarian. Meanwhile, single-cell RNA-sequencing provides insight into the regenerative process at a cellular level and finds that the largest proportion of cells present during regeneration are stem cells.

The mechanisms regulating regeneration of the earthworm are unclear. Here, the authors use genomic and transcriptomic analysis of the earthworm Eisenia andrei together with Hi-C analysis to identify genes involved and show activation of LINE2 transposable elements on regeneration.