Genome Size Expansion Research Articles

The genome of Vicia sativa ssp. amphicarpa provides insights into the role of terpenoids in antimicrobial resistance within subterranean fruits.

Vicia sativa ssp. amphicarpa is a unique forage crop capable of simultaneously producing fruits above and below ground, representing a typical amphicarpic plant. In this study, we sequenced and assembled seven pseudo-chromosomes of the genome of V. sativa ssp. amphicarpa (n = 7) yielding a genome size of 1.59 Gb, with a total annotation of 48 932 protein-coding genes. Long terminal repeat (LTR) elements constituted 62.28% of the genome, significantly contributing to the expansion of genome size. Phylogenetic analysis revealed that the divergence between V. sativa ssp. amphicarpa and V. sativa was around 0.88 million years ago (MYA). Comparative transcriptomic and metabolomic analysis of aerial and subterranean pod shells showed biosynthesis of terpenoids in the subterranean pod shells indicating a correlation between the antimicrobial activity of subterranean pod shells and the biosynthesis of terpenoids. Furthermore, functional validation indicates that overexpression of VsTPS5 and VsTPS16 enhances terpenoid biosynthesis for antibacterial activity. Metabolomic analysis suggests the involvement of terpenoids in the antimicrobial properties of subterranean pod shells. Deciphering the genome of V. sativa ssp. amphicarpa elucidated the molecular mechanisms behind the antimicrobial properties of subterranean fruits in amphicarpic plants, providing valuable insights for the study of amphicarpic plant biology.

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana.

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

Chromosome-level genome assembly of Pedicularis kansuensis illuminates genome evolution of facultative parasitic plant.

Parasitic plants have a heterotrophic lifestyle, in which they withdraw all or part of their nutrients from their host through the haustorium. Despite the release of many draft genomes of parasitic plants, the genome evolution related to the parasitism feature of facultative parasites remains largely unknown. In this study, we present a high-quality chromosomal-level genome assembly for the facultative parasite Pedicularis kansuensis (Orobanchaceae), which invades both legume and grass host species in degraded grasslands on the Qinghai-Tibet Plateau. This species has the largest genome size compared with other parasitic species, and expansions of long terminal repeat retrotransposons accounting for 62.37% of the assembly greatly contributed to the genome size expansion of this species. A total of 42,782 genes were annotated, and the patterns of gene loss in P. kansuensis differed from other parasitic species. We also found many mobile mRNAs between P. kansuensis and one of its host species, but these mobile mRNAs could not compensate for the functional losses of missing genes in P. kansuensis. In addition, we identified nine horizontal gene transfer (HGT) events from rosids and monocots, as well as one single-gene duplication events from HGT genes, which differ distinctly from that of other parasitic species. Furthermore, we found evidence for HGT through transferring genomic fragments from phylogenetically remote host species. Taken together, these findings provide genomic insights into the evolution of facultative parasites and broaden our understanding of the diversified genome evolution in parasitic plants and the molecular mechanisms of plant parasitism.

Apomixis and the paradox of sex in plants.

The predominance of sex in eukaryotes, despite the high costs of meiosis and mating, is still an evolutionary enigma. Many theories have been proposed, none of them being conclusive on its own, and they are partly not well applicable to land plants. Sexual reproduction is obligate in embryophytes for the great majority of species. This review will compare the main forms of sexual and asexual reproduction in ferns and angiosperms, based on the generation cycling of sporophyte and gametophyte (leaving vegetative propagation aside). The benefits of sexual reproduction for maintenance of genomic integrity compared to asexuality will be discussed in the light of developmental, evolutionary, genetic and phylogenetic studies. Asexual reproduction represents modifications of the sexual pathway, with various forms of facultative sexuality. For sexual land plants, meiosis provides direct DNA repair mechanisms of oxidative damage in reproductive tissues. The ploidy alternations of meiosis-syngamy cycles, and prolonged, multicellular stages in the haploid phase in the gametophytes provide a high efficiency of purifying selection against recessive deleterious mutations. Asexual lineages might buffer effects of such mutations via polyploidy, and can purge the mutational load via facultative sexuality. The role of organelle-nuclear genome compatibility for maintenance of genome integrity is still not well understood. The costs of mating are in plants in general low because of predominant hermaphroditism. Phylogenetic patterns in the Archaeplastid clade suggest that high frequencies of sexuality in land plants are concomitant with a stepwise increase of intrinsic and extrinsic stress factors. Furthermore, expansion of genome size in land plants would increase the potential mutational load. Sexual reproduction appears to be essential for keeping long term genomic integrity, and only rare combinations of extrinsic and intrinsic factors allow for shifts to asexuality.

Daidara: A gigantic Gypsy LTR retrotransposon lineage in the springtail Allacma fusca genome.

Long terminal repeat (LTR) retrotransposons are the major contributor to genome size expansion, as in the cases of the maize genome or the axolotl genome. Despite their impact on the genome size, the length of each retrotransposon is limited, compared to DNA transposons, which sometimes exceed over 100 kb. The longest LTR retrotransposon known to date is Burro-1 from the planarian Schmidtea medierranea, which is around 35.7 kb long. Here through bioinformatics analysis, a new lineage of gigantic LTR retrotransposons, designated Daidara, is reported from the springtail Allacma fusca genome. Their entire length (25-33 kb) rivals Burro families, while their LTRs are shorter than 1.5 kb, in contrast to other gigantic LTR retrotransposon lineages Burro and Ogre, whose LTRs are around 5 kb long. Daidara encodes three core proteins corresponding to gag, pol, and an additional protein of unknown function. The phylogenetic analysis supports the independent gigantification of Daidara from Burro or Ogre.

Integrated global analysis in spider flowers illuminates features underlying the evolution and maintenance of C4 photosynthesis.

The carbon concentrating mechanism-C4 photosynthesis-represents a classic example of convergent evolution, but how this important trait originated and evolved remains largely enigmatic. The spider flower Gynandropsis gynandra is a valuable leafy vegetable crop and medicinal plant that has also been recognized as a C4 model species. Here we present a high-quality chromosome-scale annotated genome assembly of G. gynandra through a combination of Oxford Nanopore Technology (ONT), HiFi and Hi-C technology. The 17 super-scaffolds cover 98.66% of the estimated genome (997.61Mb), with a contig N50 of 11.43Mb and a scaffold N50 of 51.02Mb. Repetitive elements occupy up to 71.91% of its genome, and over half are long terminal repeat retrotransposons (LTR-RTs) derived from recent bursts, contributing to genome size expansion. Strikingly, LTR-RT explosion also played a critical role in C4 evolution by altering expression features of photosynthesis-associated genes via preferential insertion in promoters. Integrated multiomics analyses of G. gynandra and the ornamental horticulture C3 relative Tarenaya hassleriana reveal that species-specific whole-genome duplication, gene family expansion, recent LTR-RT amplification, and more recent tandem duplication events have all facilitated the evolution of C4 photosynthesis, revealing uniqueness of C4 evolution in the Cleome genus. Moreover, high leaf vein density and heat stress resilience are associated with shifted gene expression patterns. The mode of C3-to-C4 transition found here yields new insights into evolutionary convergence of a complex plant trait. The availability of this reference-grade genomic resource makes G. gynandra an ideal model system facilitating efforts toward C4-aimed crop engineering.

MHC architecture in amphibians - ancestral reconstruction, gene rearrangements and duplication patterns.

The hypervariable major histocompatibility complex (MHC) is a crucial component of vertebrate adaptive immunity, but large-scale studies on MHC macroevolution in non-model vertebrates have long been constrained by methodological limitations. Here, we used rapidly accumulating genomic data to reconstruct macroevolution of the MHC region in amphibians. We retrieved contigs containing the MHC region from genome assemblies of 32 amphibian species and examined major structural rearrangements, duplication patterns and gene structure across the amphibian phylogeny. Based on the few available caecilian and urodele genomes we showed that the structure of ancestral MHC region in amphibians was probably relatively simple and compact, with a close physical linkage between MHC-I and MHC-II regions. This ancestral MHC architecture was generally conserved in anurans, although the evolution of class I subregion proceeded towards more extensive duplication and rapid expansion of gene copy number, providing evidence for dynamic evolutionary trajectories. Although in anurans we recorded tandems of duplicated MHC-I genes outside the core subregion, our phylogenetic analyses of MHC-I sequences provided little support for an expansion of nonclassical MHC-Ib genes across amphibian families. Finally, we found that intronic regions of amphibian classical MHC genes were much longer when compared to other tetrapod lineages (birds and mammals), which could partly be driven by the expansion of genome size. Our study reveals novel evolutionary patterns of the MHC region in amphibians and provides a comprehensive framework for further studies on the MHC macroevolution across vertebrates.

Chromosome-level genome assembly of the endangered plant Tetraena mongolica.

Tetraena mongolica is an endangered xerophytic shrub with high ecological value for the restoration of desert vegetation because of its high tolerance to drought and heat stress. Here, we generated a high-quality chromosome-level reference genome of T. mongolica by combining PacBio HiFi data and Hi-C sequencing technologies, which was approximately 1.12 Gb (contig N50 of 25.5 Mb) in size and contained 61,888 protein-coding genes; repetitive sequences comprised 44.8% of the genome. This genome of T. mongolica is the first published genome sequence of a member of the order Zygophyllales. Genome analysis showed that T. mongolica has undergone a recent whole genome duplication event, and a recent burst of long terminal repeat insertions afterward, which may be responsible for its genome size expansion and drought adaptation. We also conducted searches for gene homologues and identified terpene synthase (TPS) gene families and candidate genes involved in triacylglycerol biosynthesis. The T. mongolica genome sequence could aid future studies aimed at functional gene identification, germplasm resource management, molecular breeding efforts, as well as evolutionary studies of Fabids and angiosperm taxa.

Whole-genome sequencing of a biocontrol Myxococcus xanthus R31 isolate and comparative genomic analysis

Tomato bacterial wilt (TBW) caused by Ralstonia solanacearum is one of the most destructive soil-borne diseases. Myxococcus xanthus R31, isolated from healthy tomato rhizosphere soil using the R. solanacearum baiting method, exhibiting good biocontrol efficacy against TBW. However, the genomic information and evolutionary features of R31 are largely unclear. Here, the high-quality genome assembly of R31 was presented. Using Nanopore sequencing technology, we assembled the 9.25 Mb complete genome of R31 and identified several extracellular enzyme proteins, including carbohydrate-active enzymes (CAZymes) and peptidases. We also performed a comparative genome analysis of R31 and 17 other strains of M. xanthus with genome sequences in the NCBI database to gain insights into myxobacteria predation and genome size expansion. Average nucleotide identity and digital DNA-DNA hybridization calculation and phylogenetic analysis indicated that R31 was closely related to the species M. xanthus. Further comparative genomics analysis suggested that, in addition to characteristics of predatory microorganisms, R31 contains many strain-specific genes, which may provide a genetic basis for its proficient predatory ability. This study provides new insights into R31 and other closely related species and facilitates studies using genetic approaches to further elucidate the predation mechanism of myxobacteria.

The genome of Lyophyllum shimeji provides insight into the initial evolution of ectomycorrhizal fungal genomes

Mycorrhizae are one of the most fundamental symbioses between plants and fungi, with ectomycorrhizae being the most widespread in boreal forest ecosystems. Ectomycorrhizal fungi are hypothesized to have evolved convergently from saprotrophic ancestors in several fungal clades, especially members of the subdivision Agaricomycotina. Studies on fungal genomes have identified several typical characteristics of mycorrhizal fungi, such as genome size expansion and decreases in plant cell-wall degrading enzymes (PCWDEs). However, genomic changes concerning the evolutionary transition to the ectomycorrhizal lifestyle are largely unknown. In this study, we sequenced the genome of Lyophyllum shimeji, an ectomycorrhizal fungus that is phylogenetically related to saprotrophic species and retains some saprotroph-like traits. We found that the genome of Ly. shimeji strain AT787 lacks both incremental increases in genome size and reduced numbers of PCWDEs. Our findings suggest that the previously reported common genomic traits of mycorrhizal fungi are not essential for the ectomycorrhizal lifestyle, but are a result of abolishing saprotrophic activity. Since Ly. shimeji is commercially consumed as an edible mushroom, the newly available genomic information may also impact research designed to enhance the cultivation of this mushroom.

The Evolution of tRNA Copy Number and Repertoire in Cellular Life.

tRNAs are universal decoders that bridge the gap between transcriptome and proteome. They can also be processed into small RNA fragments with regulatory functions. In this work, we show that tRNA copy number is largely controlled by genome size in all cellular organisms, in contrast to what is observed for protein-coding genes that stop expanding between ~20,000 and ~35,000 loci per haploid genome in eukaryotes, regardless of genome size. Our analyses indicate that after the bacteria/archaea split, the tRNA gene pool experienced the evolution of increased anticodon diversity in the archaeal lineage, along with a tRNA gene size increase and mature tRNA size decrease. The evolution and diversification of eukaryotes from archaeal ancestors involved further expansion of the tRNA anticodon repertoire, additional increase in tRNA gene size and decrease in mature tRNA length, along with an explosion of the tRNA gene copy number that emerged coupled with accelerated genome size expansion. Our findings support the notion that macroscopic eukaryotes with a high diversity of cell types, such as land plants and vertebrates, independently evolved a high diversity of tRNA anticodons along with high gene redundancy caused by the expansion of the tRNA copy number. The results presented here suggest that the evolution of tRNA genes played important roles in the early split between bacteria and archaea, and in eukaryogenesis and the later emergence of complex eukaryotes, with potential implications in protein translation and gene regulation through tRNA-derived RNA fragments.

Haplotype-phased and chromosome-level genome assembly of Puccinia polysora, a giga-scale fungal pathogen causing southern corn rust.

Rust fungi are characterized by large genomes with high repeat content and have two haploid nuclei in most life stages, which makes achieving high-quality genome assemblies challenging. Here, we described a pipeline using HiFi reads and Hi-C data to assemble a gigabase-sized fungal pathogen, Puccinia polysora f.sp. zeae, to haplotype-phased and chromosome-scale. The final assembled genome is 1.71 Gbp, with ~850 Mbp and 18 chromosomes in each haplotype, being currently one of the two giga-scale fungi assembled to chromosome level. Transcript-based annotation identified 47,512 genes for the dikaryotic genome with a similar number for each haplotype. A high level of interhaplotype variation was found with 10% haplotype-specific BUSCO genes, 5.8 SNPs/kbp, and structural variation accounting for 3% of the genome size. The P.polysora genome displayed over 85% repeat contents, with genome-size expansion and copy number increasing of species-specific orthogroups. Interestingly, these features did not affect overall synteny with other Puccinia species having smaller genomes. Fine-time-point transcriptomics revealed seven clusters of coexpressed secreted proteins that are conserved between two haplotypes. The fact that candidate effectors interspersed with all genes indicated the absence of a "two-speed genome" evolution in P.polysora. Genome resequencing of 79 additional isolates revealed a clonal population structure of P.polysora in China with low geographic differentiation. Nevertheless, a minor population differentiated from the major population by having mutations on secreted proteins including AvrRppC, indicating the ongoing virulence to evade recognition by RppC, a major resistance gene in Chinese corn cultivars. The high-quality assembly provides valuable genomic resources for future studies on disease management and the evolution of P.polysora.

A comparative analysis of mitochondrial ORFs provides new insights on expansion of mitochondrial genome size in Arcidae

BackgroundArcidae, comprising about 260 species of ark shells, is an ecologically and economically important lineage of bivalve mollusks. Interestingly, mitochondrial genomes of several Arcidae species are 2–3 times larger than those of most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. The large mitochondrial genome size is mainly due to expansion of unassigned regions (regions that are functionally unassigned). Previous work on unassigned regions of Arcidae mtDNA genomes has focused on nucleotide-level analyses to observe sequence characteristics, however the origin of expansion remains unclear.ResultsWe assembled six new mitogenomes and sequenced six transcriptomes of Scapharca broughtonii to identify conserved functional ORFs that are transcribed in unassigned regions. Sixteen lineage-specific ORFs with different copy numbers were identified from seven Arcidae species, and 11 of 16 ORFs were expressed and likely biologically active. Unassigned regions of 32 Arcidae mitogenomes were compared to verify the presence of these novel mitochondrial ORFs and their distribution. Strikingly, multiple structural analyses and functional prediction suggested that these additional mtDNA-encoded proteins have potential functional significance. In addition, our results also revealed that the ORFs have a strong connection to the expansion of Arcidae mitochondrial genomes and their large-scale duplication play an important role in multiple expansion events. We discussed the possible origin of ORFs and hypothesized that these ORFs may originate from duplication of mitochondrial genes.ConclusionsThe presence of lineage-specific mitochondrial ORFs with transcriptional activity and potential functional significance supports novel features for Arcidae mitochondrial genomes. Given our observation and analyses, these ORFs may be products of mitochondrial gene duplication. These findings shed light on the origin and function of novel mitochondrial genes in bivalves and provide new insights into evolution of mitochondrial genome size in metazoans.

Repeat sequences limit the effectiveness of lateral gene transfer and favored the evolution of meiotic sex in early eukaryotes

The transition from prokaryotic lateral gene transfer to eukaryotic meiotic sex is poorly understood. Phylogenetic evidence suggests that it was tightly linked to eukaryogenesis, which involved an unprecedented rise in both genome size and the density of genetic repeats. Expansion of genome size raised the severity of Muller's ratchet, while limiting the effectiveness of lateral gene transfer (LGT) at purging deleterious mutations. In principle, an increase in recombination length combined with higher rates of LGT could solve this problem. Here, we show using a computational model that this solution fails in the presence of genetic repeats prevalent in early eukaryotes. The model demonstrates that dispersed repeat sequences allow ectopic recombination, which leads to the loss of genetic information and curtails the capacity of LGT to prevent mutation accumulation. Increasing recombination length in the presence of repeat sequences exacerbates the problem. Mutational decay can only be resisted with homology along extended sequences of DNA. We conclude that the transition to homologous pairing along linear chromosomes was a key innovation in meiotic sex, which was instrumental in the expansion of eukaryotic genomes and morphological complexity.

The ecology of palm genomes: repeat-associated genome size expansion is constrained by aridity.

Genome size varies 2400-fold across plants, influencing their evolution through changes in cell size and cell division rates which impact plants' environmental stress tolerance. Repetitive element expansion explains much genome size diversity, and the processes structuring repeat 'communities' are analogous to those structuring ecological communities. However, which environmental stressors influence repeat community dynamics has not yet been examined from an ecological perspective. We measured genome size and leveraged climatic data for 91% of genera within the ecologically diverse palm family (Arecaceae). We then generated genomic repeat profiles for 141 palm species, and analysed repeats using phylogenetically informed linear models to explore relationships between repeat dynamics and environmental factors. We show that palm genome size and repeat 'community' composition are best explained by aridity. Specifically, Ty3-gypsy and TIR elements were more abundant in palm species from wetter environments, which generally had larger genomes, suggesting amplification. By contrast, Ty1-copia and LINE elements were more abundant in drier environments. Our results suggest that water stress inhibits repeat expansion through selection on upper genome size limits. However, elements that may associate with stress-response genes (e.g. Ty1-copia) have amplified in arid-adapted palm species. Overall, we provide novel evidence of climate influencing the assembly of repeat 'communities'.

Independent Size Expansions and Intron Proliferation in Red Algal Plastid and Mitochondrial Genomes.

Proliferation of selfish genetic elements has led to significant genome size expansion in plastid and mitochondrial genomes of various eukaryotic lineages. Within the red algae, such expansion events are only known in the plastid genomes of the Proteorhodophytina, a highly diverse group of mesophilic microalgae. By contrast, they have never been described in the much understudied red algal mitochondrial genomes. Therefore, it remains unclear how widespread such organellar genome expansion events are in this eukaryotic phylum. Here, we describe new mitochondrial and plastid genomes from 25 red algal species, thereby substantially expanding the amount of organellar sequence data available, especially for Proteorhodophytina, and show that genome expansions are common in this group. We confirm that large plastid genomes are limited to the classes Rhodellophyceae and Porphyridiophyceae, which, in part, are caused by lineage-specific expansion events. Independently expanded mitochondrial genomes—up to three times larger than typical red algal mitogenomes—occur across Proteorhodophytina classes and a large shift toward high GC content occurred in the Stylonematophyceae. Although intron proliferation is the main cause of plastid and mitochondrial genome expansion in red algae, we do not observe recent intron transfer between different organelles. Phylogenomic analyses of mitochondrial and plastid genes from our expanded taxon sampling yielded well-resolved phylogenies of red algae with strong support for the monophyly of Proteorhodophytina. Our work shows that organellar genomes followed different evolutionary dynamics across red algal lineages.

Chromosome-level genome assembly of a xerophytic plant, Haloxylon ammodendron.

Haloxylon ammodendron is a xerophytic perennial shrub or small tree that has a high ecological value in anti-desertification due to its high tolerance to drought and salt stress. Here, we report a high-quality, chromosome-level genome assembly of H. ammodendron by integrating PacBio’s high-fidelity sequencing and Hi-C technology. The assembled genome size was 685.4 Mb, of which 99.6% was assigned to nine pseudochromosomes with a contig N50 value of 23.6 Mb. Evolutionary analysis showed that both the recent substantial amplification of long terminal repeat retrotransposons and tandem gene duplication may have contributed to its genome size expansion and arid adaptation. An ample amount of low-GC genes was closely related to functions that may contribute to the desert adaptation of H. ammodendron. Gene family clustering together with gene expression analysis identified differentially expressed genes that may play important roles in the direct response of H. ammodendron to water-deficit stress. We also identified several genes possibly related to the degraded scaly leaves and well-developed root system of H. ammodendron. The reference-level genome assembly presented here will provide a valuable genomic resource for studying the genome evolution of xerophytic plants, as well as for further genetic breeding studies of H. ammodendron.

The comparative studies of complete chloroplast genomes in Actinidia (Actinidiaceae): novel insights into heterogenous variation, clpP gene annotation and phylogenetic relationships.

The genus Actinidia, also called kiwifruit, is characterized with abundant balanced nutritional metabolites, including exceptionally high vitamin C content. However, the traditional classification could not fully reflect the actual Actinidia species' relationships, which need further revision through more accurate approaches. Compared to the nuclear genome, the chloroplast genome has simple heredity characteristics, conserved genome structure and small size, suitable for deciphering complicated species' phylogenetic relationships. Here, the genome-wide comprehensive comparative analyses were performed over 29 independent chloroplast genomes' sequences derived from 25 Actinidia taxa. The average genome size is 156,673.38bp, with an average 37.20% GC content. The long repeat sequences rather than SSRs (simple sequence repeats) in Actinidia were revealed to be the causal agent leading to the chloroplast genome size expansion. The clpP gene sequences with exon merge and intron deletion were annotated in all the 29 chloroplast genomes tested, which has been previously reported to be lost in Actinidia species. Comprehensive sequence analyses indicated the distinct variation at the clpP gene locus was Actinidiaceae-specific, emerging after the Actinidiaceae-other Ericales species divergence. Four highly divergent sequences (i.e., rps16 ~ trnQ-UUG, rps4 ~ trnT-UGU, petA ~ psbJ, and rps12 ~ psbB) evolved in the LSC (large single-copy) and SSC (small single-copy) regions embodying rps12 ~ psbB (including clpP gene and its up/downstream noncoding sequence) were identified as variation hot spots in Actinidia species. Based on either LSC region alone, combined sequences of LSC and SSC or the whole chloroplast genome sequences, three identical phylogenetic trees of the 25 Actinidia taxa with relatively improved resolution were reconstructed, consistently supporting the reticulate evolutionary lineage in Actinidia. Our findings could help to better understand the evolution characteristics of chloroplast genomes and phylogenetic relationships among Actinidia species.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi.

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Rapid Genome Evolution and Adaptation of Thlaspi arvense Mediated by Recurrent RNA-Based and Tandem Gene Duplications.

Retrotransposons are the most abundant group of transposable elements (TEs) in plants, providing an extraordinarily versatile source of genetic variation. Thlaspi arvense, a close relative of the model plant Arabidopsis thaliana with worldwide distribution, thrives from sea level to above 4,000 m elevation in the Qinghai-Tibet Plateau (QTP), China. Its strong adaptability renders it an ideal model system for studying plant adaptation in extreme environments. However, how the retrotransposons affect the T. arvense genome evolution and adaptation is largely unknown. We report a high-quality chromosome-scale genome assembly of T. arvense with a scaffold N50 of 59.10 Mb. Long terminal repeat retrotransposons (LTR-RTs) account for 56.94% of the genome assembly, and the Gypsy superfamily is the most abundant TEs. The amplification of LTR-RTs in the last six million years primarily contributed to the genome size expansion in T. arvense. We identified 351 retrogenes and 303 genes flanked by LTRs, respectively. A comparative analysis showed that orthogroups containing those retrogenes and genes flanked by LTRs have a higher percentage of significantly expanded orthogroups (SEOs), and these SEOs possess more recent tandem duplicated genes. All present results indicate that RNA-based gene duplication (retroduplication) accelerated the subsequent tandem duplication of homologous genes resulting in family expansions, and these expanded gene families were implicated in plant growth, development, and stress responses, which were one of the pivotal factors for T. arvense’s adaptation to the harsh environment in the QTP regions. In conclusion, the high-quality assembly of the T. arvense genome provides insights into the retroduplication mediated mechanism of plant adaptation to extreme environments.