Assembled and annotated 26.5 Gbp coast redwood genome: a resource for estimating evolutionary adaptive potential and investigating hexaploid origin

Occurrence of aneuploidy across the range of coast redwood (Sequoia sempervirens)

Aneuploidy, a condition characterized by an abnormal number of chromosomes, can have significant consequences for fitness of an organism, often manifesting in reduced fertility and other developmental challenges. In plants, aneuploidy is particularly complex to study, especially in polyploid species such as coast redwood (Sequoia sempervirens (D. Don) Endl.), which is a hexaploid conifer (2n=6×=66). This study leverages a novel Markov Chain Monte Carlo method based on sequence depth to investigate the occurrence of aneuploidy across the range of coast redwood. We show that aneuploidy, defined here as a whole-chromosome gain or loss, is prevalent in second-growth redwoods, predominantly as additional chromosomes, while vegetatively propagated plants frequently experience chromosome loss. Although our study does not directly assess the fitness of aneuploids, the frequency of chromosomal instability observed in vegetatively propagated plants compared to second-growth and old-growth trees raises questions about their long-term developmental viability and potential to become established trees. These findings have significant implications for redwood conservation and restoration strategies, especially as methods such as tissue culture propagation becomes the primary mode of producing nursery stock plants used in reforestation.

Variant calling in polyploids for population and quantitative genetics

Advancements in genome assembly and sequencing technology have made whole genome sequence (WGS) data and reference genomes accessible to study polyploid species. Compared to popular reduced-representation sequencing approaches, the genome-wide coverage and greater marker density provided by WGS data can greatly improve our understanding of polyploid species and polyploid biology. However, biological features that make polyploid species interesting also pose challenges in read mapping, variant identification, and genotype estimation. Accounting for characteristics in variant calling like allelic dosage uncertainty, homology between subgenomes, and variance in chromosome inheritance mode can reduce errors. Here, I discuss the challenges of variant calling in polyploid WGS data and discuss where potential solutions can be integrated into a standard variant calling pipeline.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Piece and parcel of gymnosperm organellar genomes

MAIN CONCLUSION

Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity.

A genome sequence for the threatened whitebark pine

Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Plant secondary metabolites involved in the stress tolerance of long-lived trees

Ancient trees are natural wonders because of their longevity, having lived for hundreds or thousands of years, and their ability to withstand changing environments and a variety of stresses. These long-lived trees have sophisticated defense mechanisms, such as the production of specialized plant metabolites (SPMs). In this review, we provide an overview of the major biotic and abiotic stresses that long-lived trees often face, as well as an analysis of renowned ancient tree species and their unique protective SPMs against environmental stressors. We also discuss the synthesis and accumulation of defensive SPMs induced by environmental factors and endophytes in these trees. Furthermore, we conducted a comparative genomic analysis of 17 long-lived tree species and discovered significant expansions of SPM biosynthesis gene families in these species. Our comprehensive review reveals the crucial role of SPMs in high resistance in long-lived trees, providing a novel natural resource for plant defense, crop improvement and even the pharmaceutical industry.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Identification and high-throughput genotyping of single nucleotide polymorphism markers in a non-model conifer (Abies nordmanniana (Steven) Spach)

Single nucleotide polymorphism (SNP) markers are powerful tools for investigating population structures, linkage analysis, and genome-wide association studies, as well as for breeding and population management. The availability of SNP markers has been limited to the most commercially important timber species, primarily due to the cost of genome sequencing required for SNP discovery. In this study, a combination of reference-based and reference-free approaches were used to identify SNPs in Nordmann fir (Abies nordmanniana), a species previously lacking genomic sequence information. Using a combination of a genome assembly of the closely related Silver fir (Abies alba) species and a de novo assembly of low-copy regions of the Nordmann fir genome, we identified a high density of reliable SNPs. Reference-based approaches identified two million SNPs in common between the Silver fir genome and low-copy regions of Nordmann fir. A combination of one reference-free and two reference-based approaches identified 250 shared SNPs. A subset of 200 SNPs were used to genotype 342 individuals and thereby tested and validated in the context of identity analysis and/or clone identification. The tested SNPs successfully identified all ramets per clone and five mislabeled individuals via identity and genomic relatedness analysis. The identified SNPs will be used in ad hoc breeding of Nordmann fir in Denmark.

Three-dimensional chromatin architecture in plants - General features and novelties

Research on the three-dimensional (3D) structure of the genome and its distribution within the nuclear space has made a big leap in the last two decades. Work in the animal field has led to significant advances in our general understanding on eukaryotic genome organization. This did not only bring along insights into how the 3D genome interacts with the epigenetic landscape and the transcriptional machinery but also how 3D genome architecture is relevant for fundamental developmental processes, such as cell differentiation. In parallel, the 3D organization of plant genomes have been extensively studied, which resulted in both congruent and novel findings, contributing to a more complete view on how eukaryotic genomes are organized in multiple dimensions. Plant genomes are remarkably diverse in size, composition, and ploidy. Furthermore, as intrinsically sessile organisms without the possibility to relocate to more favorable environments, plants have evolved an elaborate epigenetic repertoire to rapidly respond to environmental challenges. The diversity in genome organization and the complex epigenetic programs make plants ideal study subjects to acquire a better understanding on universal features and inherent constraints of genome organization. Furthermore, considering a wide range of species allows us to study the evolutionary crosstalk between the various levels of genome architecture. In this article, we aim at summarizing important findings on 3D genome architecture obtained in various plant species. These findings cover many aspects of 3D genome organization on a wide range of levels, from gene loops to topologically associated domains and to global 3D chromosome configurations. We present an overview on plant 3D genome organizational features that resemble those in animals and highlight facets that have only been observed in plants to date.

A Genome Sequence for the Threatened Whitebark Pine

Whitebark pine (WBP, Pinus albicaulis ) is a white pine of subalpine regions in western contiguous US and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola ) and additional threats from mountain pine beetle ( Dendroctonus ponderosae ), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short-reads of haploid megametophyte tissue and Oxford Nanopore long-reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gbp of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gbp). Approximately 87.2% (24.0 Gbp) of total sequence was placed on the twelve WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich-repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the three subclasses of NLRs (TNL, CNL, RNL). Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo assembled transcriptomes.

The Metasequoia genome and evolutionary relationships among redwoods

Redwood trees (Sequoioideae), including Metasequoia glyptostroboides (dawn redwood), Sequoiadendron giganteum (giant sequoia), and Sequoia sempervirens (coast redwood), are threatened and widely recognized iconic tree species. Genomic resources for redwood trees could provide clues to their evolutionary relationships. Here, we report the 8-Gb reference genome of M. glyptostroboides and a comparative analysis with two related species. More than 62% of the M. glyptostroboides genome is composed of repetitive sequences. Clade-specific bursts of long terminal repeat retrotransposons may have contributed to genomic differentiation in the three species. The chromosomal synteny between M. glyptostroboides and S. giganteum is extremely high, whereas there has been significant chromosome reorganization in S. sempervirens. Phylogenetic analysis of marker genes indicates that S. sempervirens is an autopolyploid, and more than 48% of the gene trees are incongruent with the species tree. Results of multiple analyses suggest that incomplete lineage sorting (ILS) rather than hybridization explains the inconsistent phylogeny, indicating that genetic variation among redwoods may be due to random retention of polymorphisms in ancestral populations. Functional analysis of ortholog groups indicates that gene families of ion channels, tannin biosynthesis enzymes, and transcription factors for meristem maintenance have expanded in S. giganteum and S. sempervirens, which is consistent with their extreme height. As a wetland-tolerant species, M. glyptostroboides shows a transcriptional response to flooding stress that is conserved with that of analyzed angiosperm species. Our study offers insights into redwood evolution and adaptation and provides genomic resources to aid in their conservation and management.

Mabs, a suite of tools for gene-informed genome assembly

BACKGROUND

Despite constantly improving genome sequencing methods, error-free eukaryotic genome assembly has not yet been achieved. Among other kinds of problems of eukaryotic genome assembly are so-called "haplotypic duplications", which may manifest themselves as cases of alleles being mistakenly assembled as paralogues. Haplotypic duplications are dangerous because they create illusions of gene family expansions and, thus, may lead scientists to incorrect conclusions about genome evolution and functioning.

RESULTS

Here, I present Mabs, a suite of tools that serve as parameter optimizers of the popular genome assemblers Hifiasm and Flye. By optimizing the parameters of Hifiasm and Flye, Mabs tries to create genome assemblies with the genes assembled as accurately as possible. Tests on 6 eukaryotic genomes showed that in 6 out of 6 cases, Mabs created assemblies with more accurately assembled genes than those generated by Hifiasm and Flye when they were run with default parameters. When assemblies of Mabs, Hifiasm and Flye were postprocessed by a popular tool for haplotypic duplication removal, Purge_dups, genes were better assembled by Mabs in 5 out of 6 cases.

CONCLUSIONS

Mabs is useful for making high-quality genome assemblies. It is available at https://github.com/shelkmike/Mabs.

A gap-free genome assembly of Chlamydomonas reinhardtii and detection of translocations induced by CRISPR-mediated mutagenesis

Genomic assemblies of the unicellular green alga Chlamydomonas reinhardtii have provided important resources for researchers. However, assembly errors, large gaps, and unplaced scaffolds as well as strain-specific variants currently impede many types of analysis. By combining PacBio HiFi and Oxford Nanopore long-read technologies, we generated a de novo genome assembly for strain CC-5816, derived from crosses of strains CC-125 and CC-124. Multiple methods of evaluating genome completeness and base-pair error rate suggest that the final telomere-to-telomere assembly is highly accurate. The CC-5816 assembly enabled previously difficult analyses that include characterization of the 17 centromeres, rDNA arrays on three chromosomes, and 56 insertions of organellar DNA into the nuclear genome. Using Nanopore sequencing, we identified sites of cytosine (CpG) methylation, which are enriched at centromeres. We analyzed CRISPR-Cas9 insertional mutants in the PF23 gene. Two of the three alleles produced progeny that displayed patterns of meiotic inviability that suggested the presence of a chromosomal aberration. Mapping Nanopore reads from pf23-2 and pf23-3 onto the CC-5816 genome showed that these two strains each carry a translocation that was initiated at the PF23 gene locus on chromosome 11 and joined with chromosomes 5 or 3, respectively. The translocations were verified by demonstrating linkage between loci on the two translocated chromosomes in meiotic progeny. The three pf23 alleles display the expected short-cilia phenotype, and immunoblotting showed that pf23-2 lacks the PF23 protein. Our CC-5816 genome assembly will undoubtedly provide an important tool for the Chlamydomonas research community.

Comprehensive collection of genes and comparative analysis of full-length transcriptome sequences from Japanese larch (Larix kaempferi) and Kuril larch (Larix gmelinii var. japonica)

BACKGROUND

Japanese larch (Larix kaempferi) is an economically important deciduous conifer species that grows in cool-temperate forests and is endemic to Japan. Kuril larch (L. gmelinii var. japonica) is a variety of Dahurian larch that is naturally distributed in the Kuril Islands and Sakhalin. The hybrid larch (L. gmelinii var. japonica × L. kaempferi) exhibits heterosis, which manifests as rapid juvenile growth and high resistance to vole grazing. Since these superior characteristics have been valued by forestry managers, the hybrid larch is one of the most important plantation species in Hokkaido. To accelerate molecular breeding in these species, we collected and compared full-length cDNA isoforms (Iso-Seq) and RNA-Seq short-read, and merged them to construct candidate gene as reference for both Larix species. To validate the results, candidate protein-coding genes (ORFs) related to some flowering signal-related genes ​were screened from the reference sequences, and the phylogenetic relationship with closely related species was elucidated.

RESULTS

Using the isoform sequencing of PacBio RS ll and the de novo assembly of RNA-Seq short-read sequences, we identified 50,690 and 38,684 ORFs in Japanese larch and Kuril larch, respectively. BUSCO completeness values were 90.5% and 92.1% in the Japanese and Kuril larches, respectively. After comparing the collected ORFs from the two larch species, a total of 19,813 clusters, comprising 22,571 Japanese larch ORFs and 22,667 Kuril larch ORFs, were contained in the intersection of the Venn diagram. In addition, we screened several ORFs related to flowering signals (SUPPRESSER OF OVEREXPRESSION OF CO1: SOC1, LEAFY: LFY, FLOWERING Locus T: FT, CONSTANCE: CO) from both reference sequences, and very similar found in other species.

CONCLUSIONS

The collected ORFs will be useful as reference sequences for molecular breeding of Japanese and Kuril larches, and also for clarifying the evolution of the conifer genome and investigating functional genomics.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

Current status and trends in forest genomics

Forests are not only the most predominant of the Earth's terrestrial ecosystems, but are also the core supply for essential products for human use. However, global climate change and ongoing population explosion severely threatens the health of the forest ecosystem and aggravtes the deforestation and forest degradation. Forest genomics has great potential of increasing forest productivity and adaptation to the changing climate. In the last two decades, the field of forest genomics has advanced quickly owing to the advent of multiple high-throughput sequencing technologies, single cell RNA-seq, clustered regularly interspaced short palindromic repeats (CRISPR)-mediated genome editing, and spatial transcriptomes, as well as bioinformatics analysis technologies, which have led to the generation of multidimensional, multilayered, and spatiotemporal gene expression data. These technologies, together with basic technologies routinely used in plant biotechnology, enable us to tackle many important or unique issues in forest biology, and provide a panoramic view and an integrative elucidation of molecular regulatory mechanisms underlying phenotypic changes and variations. In this review, we recapitulated the advancement and current status of 12 research branches of forest genomics, and then provided future research directions and focuses for each area. Evidently, a shift from simple biotechnology-based research to advanced and integrative genomics research, and a setup for investigation and interpretation of many spatiotemporal development and differentiation issues in forest genomics have just begun to emerge.

Evolution of complex genome architecture in gymnosperms

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5–15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants.

Repetitive Elements, Sequence Turnover and Cyto-Nuclear Gene Transfer in Gymnosperm Mitogenomes

Among the three genomes in plant cells, the mitochondrial genome (mitogenome) is the least studied due to complex recombination and intergenomic transfer. In gymnosperms only ∼20 mitogenomes have been released thus far, which hinders a systematic investigation into the tempo and mode of mitochondrial DNA evolution in seed plants. Here, we report the complete mitogenome sequence of Platycladus orientalis (Cupressaceae). This mitogenome is assembled as two circular-mapping chromosomes with a size of ∼2.6 Mb and which contains 32 protein-coding genes, three rRNA and seven tRNA genes, and 1,068 RNA editing sites. Repetitive sequences, including dispersed repeats, transposable elements (TEs), and tandem repeats, made up 23% of the genome. Comparative analyses with 17 other mitogenomes representing the five gymnosperm lineages revealed a 30-fold difference in genome size, 80-fold in repetitive content, and 230-fold in substitution rate. We found dispersed repeats are highly associated with mitogenome expansion (r = 0.99), and most of them were accumulated during recent duplication events. Syntenic blocks and shared sequences between mitogenomes decay rapidly with divergence time (r = 0.53), with the exceptions of Ginkgo and Cycads which retained conserved genome structure over long evolutionary time. Our phylogenetic analysis supports a sister group relationship of Cupressophytes and Gnetophytes; both groups are unique in that they lost 8–12 protein-coding genes, of which 4–7 intact genes are likely transferred to nucleus. These two clades also show accelerated and highly variable substitution rates relative to other gymnosperms. Our study highlights the dynamic and enigmatic evolution of gymnosperm mitogenomes.

Wood Formation under Changing Environment: Omics Approaches to Elucidate the Mechanisms Driving the Early-to-Latewood Transition in Conifers

The global change scenarios highlight the urgency of clarifying the mechanisms driving the determination of wood traits in forest trees. Coniferous xylem is characterized by the alternation between earlywood (EW) and latewood (LW), on which proportions the wood density depend, one of the most important mechanical xylem qualities. However, the molecular mechanisms triggering the transition between the production of cells with the typical features of EW to the LW are still far from being completely elucidated. The increasing availability of omics resources for conifers, e.g., genomes and transcriptomes, would lay the basis for the comprehension of wood formation dynamics, boosting both breeding and gene-editing approaches. This review is intended to introduce the importance of wood formation dynamics and xylem traits of conifers in a changing environment. Then, an up-to-date overview of the omics resources available for conifers was reported, focusing on both genomes and transcriptomes. Later, an analysis of wood formation studies using omics approaches was conducted, with the aim of elucidating the main metabolic pathways involved in EW and LW determination. Finally, the future perspectives and the urgent needs on this research topic were highlighted.