The first homosporous lycophyte genome revealed the association between the recent dynamic accumulation of LTR-RTs and genome size variation.

The contrasting genome size between homosporous and heterosporous plants is fascinating. Different from the heterosporous seed plants and mainly homosporous ferns, the lycophytes are either heterosporous (Isoetales and Selaginellales) or homosporous (Lycopodiales). Many lycophytes are the resource plants of Huperzine A (HupA) which is invaluable for treating Alzheimer's disease. For the seed-free vascular plants, several high-quality genomes of heterosporous Selaginella, homosporous ferns (maidenhair fern, monkey spider tree fern), and heterosporous ferns (Azolla) have been published and provided important insights into the origin and evolution of early land plants. However, the homosporous lycophyte genome has not been decoded. Here, we assembled the first homosporous lycophyte genome and conducted comparative genomic analyses by applying a reformed pipeline for filtering out non-plant sequences. The obtained genome size of Lycopodium clavatum is 2.30 Gb, distinguished in more than 85% repetitive elements of which 62% is long terminal repeat (LTR). This study disclosed a high birth rate and a low death rate of the LTR-RTs in homosporous lycophytes, but the opposite occurs in heterosporous lycophytes. we propose that the recent activity of LTR-RT is responsible for the immense genome size variation between homosporous and heterosporous lycophytes. By combing Ks analysis with a phylogenetic approach, we discovered two whole genome duplications (WGD). Morover, we identified all the five recognized key enzymes for the HupA biosynthetic pathway in the L. clavatum genome, but found this pathway incomplete in other major lineages of land plants. Overall, this study is of great importance for the medicinal utilization of lycophytes and the decoded genome data will be a key cornerstone to elucidate the evolution and biology of early vascular land plants.

Genome structure and evolutionary history of frankincense producing Boswellia sacra.

Boswellia sacra Flueck (family Burseraceae) tree is wounded to produce frankincense. We report its de novo assembled genome (667.8 Mb) comprising 18,564 high-confidence protein-encoding genes. Comparing conserved single-copy genes across eudicots suggest >97% gene space assembly of B. sacra genome. Evolutionary history shows B. sacra gene-duplications derived from recent paralogous events and retained from ancient hexaploidy shared with other eudicots. The genome indicated a major expansion of Gypsy retroelements in last 2 million years. The B. sacra genetic diversity showed four clades intermixed with a primary genotype-dominating most resin-productive trees. Further, the stem transcriptome revealed that wounding concurrently activates phytohormones signaling, cell wall fortification, and resin terpenoid biosynthesis pathways leading to the synthesis of boswellic acid-a key chemotaxonomic marker of Boswellia. The sequence datasets reported here will serve as a foundation to investigate the genetic determinants of frankincense and other resin-producing species in Burseraceae.

The evolution of extremely diverged plastomes in Selaginellaceae (lycophyte) is driven by repeat patterns and the underlying DNA maintenance machinery.

Two factors are proposed to account for the unusual features of organellar genomes: the disruptions of organelle-targeted DNA replication, repair, and recombination (DNA-RRR) systems in the nuclear genome and repetitive elements in organellar genomes. Little is known about how these factors affect organellar genome evolution. The deep-branching vascular plant family Selaginellaceae is known to have a deficient DNA-RRR system and convergently evolved organellar genomes. However, we found that the plastid genome (plastome) of Selaginella sinensis has extremely accelerated substitution rates, a low GC content, pervasive repeat elements, a dynamic network structure, and it lacks direct or inverted repeats. Unexpectedly, its organelle DNA-RRR system is short of a plastid-targeted Recombinase A1 (RecA1) and a mitochondrion-targeted RecA3, in line with other explored Selaginella species. The plastome contains a large collection of short- and medium-sized repeats. Given the absence of RecA1 surveillance, we propose that these repeats trigger illegitimate recombination, accelerated mutation rates, and structural instability. The correlations between repeat quantity and architectural complexity in the Selaginella plastomes support these conclusions. We, therefore, hypothesize that the interplay of the deficient DNA-RRR system and the high repeat content has led to the extraordinary divergence of the S. sinensis plastome. Our study not only sheds new light on the mechanism of plastome divergence by emphasizing the power of cytonuclear integration, but it also reconciles the longstanding contradiction on the effects of DNA-RRR system disruption on genome structure evolution.

Preferential gene retention increases the robustness of cold regulation in Brassicaceae and other plants after polyploidization.

Cold stress profoundly affects plant growth and development and is a key factor affecting the geographic distribution and evolution of plants. Plants have evolved adaptive mechanisms to cope with cold stress. Here, through the genomic analysis of Arabidopsis, three Brassica species and 17 other representative species, we found that both cold-related genes (CRGs) and their collinearity were preferentially retained after polyploidization followed by genome instability, while genome-wide gene sets exhibited a variety of other expansion mechanisms. The cold-related regulatory network was increased in Brassicaceae genomes, which were recursively affected by polyploidization. By combining our findings regarding the selective retention of CRGs from this ecological genomics study with the available knowledge of cold-induced chromosome doubling, we hypothesize that cold stress may have contributed to the success of polyploid plants through both increasing polyploidization and selectively maintaining CRGs during evolution. This hypothesis requires further biological and ecological exploration to obtain solid supporting evidence, which will potentially contribute to understanding the generation of polyploids and to the field of ecological genomics.

Two Likely Auto-Tetraploidization Events Shaped Kiwifruit Genome and Contributed to Establishment of the Actinidiaceae Family.

The genome of kiwifruit (Actinidia chinensis) was sequenced previously, the first in the Actinidiaceae family. It was shown to have been affected by polyploidization events, the nature of which has been elusive. Here, we performed a reanalysis of the genome and found clear evidence of 2 tetraploidization events, with one occurring ∼50-57 million years ago (Mya) and the other ∼18-20 Mya. Two subgenomes produced by each event have been under balanced fractionation. Moreover, genes were revealed to express in a balanced way between duplicated copies of chromosomes. Besides, lowered evolutionary rates of kiwifruit genes were observed. These findings could be explained by the likely auto-tetraploidization nature of the polyploidization events. Besides, we found that polyploidy contributed to the expansion of key functional genes, e.g., vitamin C biosynthesis genes. The present work also provided an important comparative genomics resource in the Actinidiaceae and related families.