Origin of angiosperms and the puzzle of the Jurassic gap

Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

Complete Plastid Genomes of Nine Species of Ranunculeae (Ranunculaceae) and Their Phylogenetic Inferences

The tribe Ranunculeae, Ranunculaceae, comprising 19 genera widely distributed all over the world. Although a large number of Sanger sequencing-based molecular phylogenetic studies have been published, very few studies have been performed on using genomic data to infer phylogenetic relationships within Ranunculeae. In this study, the complete plastid genomes of nine species (eleven samples) from Ceratocephala, Halerpestes, and Ranunculus were de novo assembled using a next-generation sequencing method. Previously published plastomes of Oxygraphis and other related genera of the family were downloaded from GenBank for comparative analysis. The complete plastome of each Ranunculeae species has 112 genes in total, including 78 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The plastome structure of Ranunculeae samples is conserved in gene order and arrangement. There are no inverted repeat (IR) region expansions and only one IR contraction was found in the tested samples. This study also compared plastome sequences across all the samples in gene collinearity, codon usage, RNA editing sites, nucleotide variability, simple sequence repeats, and positive selection sites. Phylogeny of the available Ranunculeae species was inferred by the plastome data using maximum-likelihood and Bayesian inference methods, and data partitioning strategies were tested. The phylogenetic relationships were better resolved compared to previous studies based on Sanger sequencing methods, showing the potential value of the plastome data in inferring the phylogeny of the tribe.

Geographical and life-history traits associated with low and high species richness across angiosperm families

Introduction

The phenomenal expansion of angiosperms has prompted many investigations into the factors driving their diversification, but there remain significant gaps in our understanding of flowering plant species diversity.

Methods

Using the crown age of families from five studies, we used a maximum likelihood approach to classify families as having poor, predicted or high species richness (SR) using strict consensus criteria. Using these categories, we looked for associations between family SR and i) the presence of an inferred familial ancestral polyploidization event, ii) 23 life history and floral traits compiled from previously published datasets and papers, and iii) sexual system (dioecy) or genetically determined self-incompatibility (SI) mating system using an updated version of our own database and iv) geographic distribution using a new database describing the global distribution of plant species/families across realms and biomes and inferred range.

Results

We find that more than a third of angiosperm families (65%) had predicted SR, a large proportion (30.2%) were species poor, while few (4.8%) had high SR. Families with poor SR were less likely to have undergone an ancestral polyploidization event, exhibited deficits in diverse traits, and were more likely to have unknown breeding systems and to be found in only one or few biomes and realms, especially the Afrotropics or Australasia. On the other hand, families with high SR were more likely to have animal mediated pollination or dispersal, are enriched for epiphytes and taxa with an annual life history, and were more likely to harbour sporophytic SI systems. Mapping the global distribution of georeferenced taxa by their family DR, we find evidence of regions dominated by taxa from lineages with high vs low SR.

Discussion

These results are discussed within the context of the literature describing "depauperons" and the factors contributing to low and high biodiversity in angiosperm clades.

New Insights into Phylogenetic Relationship of Hydrocotyle (Araliaceae) Based on Plastid Genomes

Hydrocotyle, belonging to the Hydrocotyloideae of Araliaceae, consists of 95 perennial and 35 annual species. Due to the lack of stable diagnostic morphological characteristics and high-resolution molecular markers, the phylogenetic relationships of Hydrocotyle need to be further investigated. In this study, we newly sequenced and assembled 13 whole plastid genomes of Hydrocotyle and performed comparative plastid genomic analyses with four previously published Hydrocotyle plastomes and phylogenomic analyses within Araliaceae. The plastid genomes of Hydrocotyle exhibited typical quadripartite structures with lengths from 152,659 bp to 153,669 bp, comprising a large single-copy (LSC) region (83,958-84,792 bp), a small single-copy (SSC) region (18,585-18,768 bp), and a pair of inverted repeats (IRs) (25,058-25,145 bp). Each plastome encoded 113 unique genes, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Comparative analyses showed that the IR boundaries of Hydrocotyle plastomes were highly similar, and the coding and IR regions exhibited more conserved than non-coding and single-copy (SC) regions. A total of 2932 simple sequence repeats and 520 long sequence repeats were identified, with specificity in the number and distribution of repeat sequences. Six hypervariable regions were screened from the SC region, including four intergenic spacers (IGS) (ycf3-trnS, trnS-rps4, petA-psbJ, and ndhF-rpl32) and two coding genes (rpl16 and ycf1). Three protein-coding genes (atpE, rpl16, and ycf2) were subjected to positive selection only in a few species, implying that most protein-coding genes were relatively conserved during the plastid evolutionary process. Plastid phylogenomic analyses supported the treatment of Hydrocotyle from Apiaceae to Araliaceae, and topologies with a high resolution indicated that plastome data can be further used in the comprehensive phylogenetic research of Hydrocotyle. The diagnostic characteristics currently used in Hydrocotyle may not accurately reflect the phylogenetic relationships of this genus, and new taxonomic characteristics may need to be evaluated and selected in combination with more comprehensive molecular phylogenetic results.

First record of the genus Camptotheca (Nyssaceae) in Vietnam and the lectotypification of C.acuminata

As a primary source of anticancer camptothecin, Camptotheca (Nyssaceae) is an economically valuable genus and has long been recorded as endemic to China. Here, Camptotheca is reported as a new record to the flora of Vietnam with the discovery of a wild population of C.acuminata from Lai Chau Province of northern Vietnam. Based on the consultation of relevant literature and type specimens of C.acuminata, a lectotype of the species is designated. Photographic illustrations, morphological description and a distribution map of C.acuminata is provided, and a key to all known species of Nyssaceae in Vietnam is presented, too.

Different roles of host and habitat in determining the microbial communities of plant-feeding true bugs

BACKGROUND

The true bugs (Heteroptera) occupy nearly all of the known ecological niches of insects. Among them, as a group containing more than 30,000 species, the phytophagous true bugs are making increasing impacts on agricultural and forestry ecosystems. Previous studies proved that symbiotic bacteria play important roles in these insects in fitting various habitats. However, it is still obscure about the evolutionary and ecological patterns of the microorganisms of phytophagous true bugs as a whole with comprehensive taxon sampling.

RESULTS

Here, in order to explore the symbiotic patterns between plant-feeding true bugs and their symbiotic microorganisms, 209 species belonging to 32 families of 9 superfamilies had been sampled, which covered all the major phytophagous families of true bugs. The symbiotic microbial communities were surveyed by full-length 16S rRNA gene and ITS amplicons respectively for bacteria and fungi using the PacBio platform. We revealed that hosts mainly affect the dominant bacteria of symbiotic microbial communities, while habitats generally influence the subordinate ones. Thereafter, we carried out the ancestral state reconstruction of the dominant bacteria and found that dramatic replacements of dominant bacteria occurred in the early Cretaceous and formed newly stable symbiotic relationships accompanying the radiation of insect families. In contrast, the symbiotic fungi were revealed to be horizontally transmitted, which makes fungal communities distinctive in different habitats but not significantly related to hosts.

CONCLUSIONS

Host and habitat determine microbial communities of plant-feeding true bugs in different roles. The symbiotic bacterial communities are both shaped by host and habitat but in different ways. Nevertheless, the symbiotic fungal communities are mainly influenced by habitat but not host. These findings shed light on a general framework for future microbiome research of phytophagous insects. Video Abstract.

Reconciling fossils with phylogenies reveals the origin and macroevolutionary processes explaining the global cycad biodiversity

The determinants of biodiversity patterns can be understood using macroevolutionary analyses. The integration of fossils into phylogenies offers a deeper understanding of processes underlying biodiversity patterns in deep time. Cycadales are considered a relict of a once more diverse and globally distributed group but are restricted to low latitudes today. We still know little about their origin and geographic range evolution. Combining molecular data for extant species and leaf morphological data for extant and fossil species, we study the origin of cycad global biodiversity patterns through Bayesian total-evidence dating analyses. We assess the ancestral geographic origin and trace the historical biogeography of cycads with a time-stratified process-based model. Cycads originated in the Carboniferous on the Laurasian landmass and expanded in Gondwana in the Jurassic. Through now-vanished continental connections, Antarctica and Greenland were crucial biogeographic crossroads for cycad biogeography. Vicariance is an essential speciation mode in the deep and recent past. Their latitudinal span increased in the Jurassic and restrained toward subtropical latitudes in the Neogene in line with biogeographic inferences of high-latitude extirpations. We show the benefits of integrating fossils into phylogenies to estimate ancestral areas of origin and to study evolutionary processes explaining the global distribution of present-day relict groups.

A comprehensive evaluation of flowering plant diversity and conservation priority for national park planning in China

Establishment of a national park protection system in China, including the latest target proposed to protect at least 30% of the land area, calls for a comprehensive exploration of conservation priorities incorporating multiple diversity facets. We herein evaluate the spatial distribution of Chinese flowering plants from the perspectives of richness, uniqueness, vulnerability, and evolutionary history, by integrating three mega-phylogenies and comprehensive distribution data. We detect significantly high consistency among hotspots of different diversity measures for Chinese flowering plants, suggesting that multiple facets of evolutionary diversity are concentrically distributed in China. Affording legal protection to these areas is expected to maximize positive conservation outcomes. We propose two integrative diversity indices by incorporating three richness-based and three phylogeny-based measures, respectively. Both methods identify areas with high species richness, but the integrative phylogeny-based index also locates key areas with ancient and unique evolutionary histories (e.g., Ailao-Wuliang Mts, Dabie Mts, Hainan rainforest, Karst area of Yunnan-Guizhou-Guangxi, Nanling Mts, and southeast coastal regions). Of all the diversity indices explored, phylogenetic endemism maximizes the incidental protection of other indices in most cases, emphasizing its significance for conservation planning. Finally, 42 priority areas are identified by combining the 5%-criterion hotspots of two integrative indices and the minimum area to protect all threatened species analyzed. These priorities cover only 13.3% of China's land area, but host 97.1% of species richness (23,394/24,095), 96.5% of endemic species (11,841/12,274), 100% of threatened species (2,613/2,613), and 99.3% of phylogenetic diversity for flowering plants involved in this study. These frameworks provide a solid scientific basis for national park planning in China.

Pseudanthia in angiosperms: a review

BACKGROUND

Pseudanthia or 'false flowers' are multiflowered units that resemble solitary flowers in form and function. Over the last century the term 'pseudanthium' has been applied to a wide array of morphologically divergent blossoms, ranging from those with easily noticeable florets to derived, reduced units in which individual flowers become almost indistinguishable. Although initially admired mostly by botanists, the diversity and widespread distribution of pseudanthia across angiosperms has already made them a fascinating topic for evolutionary and developmental comparative studies.

SCOPE

This review synthesizes historical and current concepts on the biology of pseudanthia. Our first aim is to establish a clear, operational definition of pseudanthium and disentangle common terminological misconceptions surrounding that term. Our second aim is to summarize knowledge of the morphological and developmental diversity of pseudanthia and embed it within a modern phylogenetic framework. Lastly, we want to provide a comprehensive overview on the evolution and ecological importance of pseudanthia and outline perspectives for future studies.

CONCLUSIONS

The understanding of pseudanthia has changed multiple times and reflects three different interpretations of their 'flower-like' qualities: developmental (similarity in structure), figural (similarity in form and function) and phylogenetic (homology between angiosperm flowers and monoecious reproductive shoots in gymnosperms). Here, we propose to narrow the term pseudanthium to multiflowered blossoms resembling zoophilous flowers in form, i.e. in being structurally subdivided in a showy periphery and a reproductive centre. According to this definition, pseudanthia sensu stricto evolved independently in at least 41 angiosperm families. The recurrent acquisition of pseudanthia sensu stricto in all major lineages of flowering plants indicates repeated interactions between developmental constraints (smallness of flowers, meristematic conditions) and selective pressures, such as demands of pollinators and/or environmental conditions.

Multiple configurations of the plastid and mitochondrial genomes of Caragana spinosa

MAIN CONCLUSION

In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes. Caragana spinosa belongs to the Papilionoidea subfamily and has significant pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled and analyzed the plastome and mitogenome of C. spinosa using the Illumina and Nanopore DNA sequencing data. The plastome of C. spinosa was 129,995 bp belonging to the inverted repeat lacking clade (IRLC), which contained 77 protein-coding genes, 29 tRNA genes, and four rRNA genes. The mitogenome was 378,373 bp long and encoded 54 unique genes, including 33 protein-coding, three ribosomal RNA (rRNA), and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, alternative conformations mediated by one and four repetitive sequences in the plastome and mitogenome were identified and validated, respectively. The inverted repeat (PDR12, the 12th dispersed repeat sequence in C. spinosa plastome) of plastome mediating recombinant was conserved in the genus Caragana. Furthermore, we identified 14 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. A phylogenetic analysis of protein-coding genes extracted from the plastid and mitochondrial genomes revealed congruent topologies. Analyses of sequence divergence found one intergenic region, trnN-GUU-ycf1, exhibiting a high degree of variation, which can be used to develop novel molecular markers to distinguish the nine Caragana species accurately. This plastome and mitogenome of C. spinosa could provide critical information for the molecular breeding of C. spinosa and be used as a reference genome for other species of Caragana. In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes.

Assembly and comparative analysis of the complete mitochondrial genome of Viburnum chinshanense

BACKGROUND

Viburnum chinshanense is an endemic species found exclusively in the North-Central and South-Central regions of China. This species is a lush garden ornamental tree and is extensively utilized for vegetation restoration in rocky desertification areas.

RESULTS

In this study, we obtained 13.96 Gb of Oxford Nanopore data for the whole genome, and subsequently, by combining Illumina short-reads, we successfully assembled the complete mitochondrial genome (mitogenome) of the V. chinshanense using a hybrid assembly strategy. The assembled genome can be described as a circular genome. The total length of the V. chinshanense mitogenome measures 643,971 bp, with a GC content of 46.18%. Our annotation efforts have revealed a total of 39 protein-coding genes (PCGs), 28 tRNA genes, and 3 rRNA genes within the V. chinshanense mitogenome. The analysis of repeated elements has identified 212 SSRs, 19 long tandem repeat elements, and 325 pairs of dispersed repeats in the V. chinshanense mitogenome. Additionally, we have investigated mitochondrial plastid DNAs (MTPTs) and identified 21 MTPTs within the mitogenome and plastidial genome. These MTPTs collectively span a length of 9,902 bp, accounting for 1.54% of the mitogenome. Moreover, employing Deepred-mt, we have confidently predicted 623 C to U RNA editing sites across the 39 protein-coding genes. Furthermore, extensive genomic rearrangements have been observed between V. chinshanense and the mitogenomes of related species. Interestingly, we have also identified a bacterial-derived tRNA gene (trnC-GCA) in the V. chinshanense mitogenome. Lastly, we have inferred the phylogenetic relationships of V. chinshanense with other angiosperms based on mitochondrial PCGs.

CONCLUSIONS

This study marks the first report of a mitogenome from the Viburnum genus, offering a valuable genomic resource for exploring the evolution of mitogenomes within the Dipsacales order.

Extensive characterization of 28 complete chloroplast genomes of Hydrangea species: A perspective view of their organization and phylogenetic and evolutionary relationships

The tribe Hydrangeeae displays a unique, distinctive disjunct distribution encompassing East Asia, North America and Hawaii. Despite its complex trait variations and polyphyletic nature, comprehensive phylogenomic and biogeographical studies on this tribe have been lacking. To address this gap, we sequenced and characterized 28 plastomes of Hydrangeeae. Our study highlights the highly conserved nature of Hydrangeaceae chloroplast (cp) genomes in terms of gene content and arrangement. Notably, synapomorphic characteristics of tandem repeats in the conserved domain of accD were observed in the Macrophyllae, Chinenses, and Dichroa sections within the Hydrangeeae tribe. Additionally, we found lower expression of accD in these sections using structure prediction and quantitative real-time PCR analysis. Phylogenomic analyses revealed the subdivision of the Hydrangeeae tribe into two clades with robust support values. Consistent with polyphyletic relationships, sect. Broussaisia was identified as the basal group in the tribe Hydrangeeae. Our study also provides insights into the phylogenetic relationships of Hydrangea petiolaris in the Jeju and Ulleung Island populations, suggesting the need for further studies with more samples and molecular data. Divergence time estimation and biogeographical analyses suggested that the common ancestors of the tribe Hydrangeeae likely originated from North America and East Asia during the Paleocene period via the Bering Land Bridge, potentially facilitating migration within the tribe between these regions. In conclusion, this study enhances our understanding of the evolutionary history and biogeography of the tribe Hydrangeeae, shedding light on the dispersal patterns and origins of this intriguing plant group with its unique disjunct distribution.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

The complete chloroplast genome of Diplodiscus trichospermus and phylogenetic position of Brownlowioideae within Malvaceae

BACKGROUND

Malvaceae is an economically important plant family of 4,225 species in nine subfamilies. Phylogenetic relationships among the nine subfamilies have always been controversial, especially for Brownlowioideae, whose phylogenetic position remains largely unknown due to the lack of samples in previous analysis datasets. To greatly clarify the phylogenetic relationship of Malvaceae, we newly sequenced and assembled the plastome of Diplodiscus trichospermus taxonomically located in Brownlowioideae, and downloaded the allied genomes from public database to build a dataset covering all subfamily members of Malvaceae.

RESULTS

The annotation results showed that the plastome of Diplodiscus trichospermus has a typical quadripartite structure, comprising 112 unique genes, namely 78 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The total length was 158,570 bp with 37.2% GC content. Based on the maximum likelihood method and Bayesian inference, a robust phylogenetic backbone of Malvaceae was reconstructed. The topology showed that Malvaceae was divided distinctly into two major branches which were previously recognized as Byttneriina and Malvadendrina. In the Malvadendrina clade, Malvoideae and Bombacoideae formed, as always, a close sister clade named as Malvatheca. Subfamily Helicteroideae occupied the most basal position and was followed by Sterculioideae which was sister to the alliance of Malvatheca, Brownlowioideae, Dombeyoideae, and Tilioideae. Brownlowioideae together with the clade comprising Dombeyoideae and Tilioideae formed a sister clade to Malvatheca. In addition, one specific conservation SSR and three specific palindrome sequences were observed in Brownlowioideae.

CONCLUSIONS

In this study, the phylogenetic framework of subfamilies in Malvaceae has been resolved clearly based on plastomes, which may contribute to a better understanding of the classification and plastome evolution for Malvaceae.

Convergent resistance to GABA receptor neurotoxins through plant-insect coevolution

The molecular mechanisms of coevolution between plants and insects remain elusive. GABA receptors are targets of many neurotoxic terpenoids, which represent the most diverse array of natural products known. Over deep evolutionary time, as plant terpene synthases diversified in plants, so did plant terpenoid defence repertoires. Here we show that herbivorous insects and their predators evolved convergent amino acid changing substitutions in duplicated copies of the Resistance to dieldrin (Rdl) gene that encodes the GABA receptor, and that the evolution of duplicated Rdl and terpenoid-resistant GABA receptors is associated with the diversification of moths and butterflies. These same substitutions also evolved in pests exposed to synthetic insecticides that target the GABA receptor. We used in vivo genome editing in Drosophila melanogaster to evaluate the fitness effects of each putative resistance mutation and found that pleiotropy both facilitates and constrains the evolution of GABA receptor resistance. The same genetic changes that confer resistance to terpenoids across 300 Myr of insect evolution have re-evolved in response to synthetic analogues over one human lifespan.

Comparative Plastomes of Curcuma alismatifolia (Zingiberaceae) Reveal Diversified Patterns among 56 Different Cut-Flower Cultivars

Curcuma alismatifolia (Zingiberaceae) is an ornamental species with high economic value due to its recent rise in popularity among floriculturists. Cultivars within this species have mixed genetic backgrounds from multiple hybridization events and can be difficult to distinguish via morphological and histological methods alone. Given the need to improve identification resources, we carried out the first systematic study using plastomic data wherein genomic evolution and phylogenetic relationships from 56 accessions of C. alismatifolia were analyzed. The newly assembled plastomes were highly conserved and ranged from 162,139 bp to 164,111 bp, including 79 genes that code for proteins, 30 tRNA genes, and 4 rRNA genes. The A/T motif was the most common of SSRs in the assembled genomes. The Ka/Ks values of most genes were less than 1, and only two genes had Ka/Ks values above 1, which were rps15 (1.15), and ndhl (1.13) with petA equal to 1. The sequence divergence between different varieties of C. alismatifolia was large, and the percentage of variation in coding regions was lower than that in the non-coding regions. Such data will improve cultivar identification, marker assisted breeding, and preservation of germplasm resources.

Phylogenetic Relationships and Next-Generation Barcodes in the Genus Torreya Reveal a High Proportion of Misidentified Cultivated Plants

Accurate species identification is key to conservation and phylogenetic inference. Living plant collections from botanical gardens/arboretum are important resources for the purpose of scientific research, but the proportion of cultivated plant misidentification are un-tested using DNA barcodes. Here, we assembled the next-generation barcode (complete plastid genome and complete nrDNA cistron) and mitochondrial genes from genome skimming data of Torreya species with multiple accessions for each species to test the species discrimination and the misidentification proportion of cultivated plants used in Torreya studies. A total of 38 accessions were included for analyses, representing all nine recognized species of genus Torreya. The plastid phylogeny showed that all 21 wild samples formed species-specific clades, except T. jiulongshanensis. Disregarding this putative hybrid, seven recognized species sampled here were successfully discriminated by the plastid genome. Only the T. nucifera accessions grouped into two grades. The species identification rate of the nrDNA cistron was 62.5%. The Skmer analysis based on nuclear reads from genome skims showed promise for species identification with seven species discriminated. The proportion of misidentified cultivated plants from arboreta/botanical gardens was relatively high with four accessions (23.5%) representing three species. Interspecific relationships within Torreya were fully resolved with maximum support by plastomes, where Torreya jackii was on the earliest diverging branch, though sister to T. grandis in the nrDNA cistron tree, suggesting that this is likely a hybrid species between T. grandis and an extinct Torreya ancestor lineage. The findings here provide quantitative insights into the usage of cultivated samples for phylogenetic study.

Profile of a flower: How rates of morphological evolution drive floral diversification in Ericales and angiosperms

PREMISE

Recent studies of floral disparity in the asterid order Ericales have shown that flowers vary strongly among families and that disparity is unequally distributed between the three flower modules (perianth, androecium, gynoecium). However, it remains unknown whether these patterns are driven by heterogeneous rates of morphological evolution or other factors.

METHODS

Here, we compiled a data set of 33 floral characters scored for 414 species of Ericales sampled from 346 genera and all 22 families. We conducted ancestral state reconstructions using an equal-rates Markov model for each character. We estimated rates of morphological evolution for Ericales and for a separate angiosperm-wide data set of 19 characters and 792 species, creating "rate profiles" for Ericales, angiosperms, and major angiosperm subclades. We compared morphological rates among flower modules within each data set separately and between data sets, and we compared rates among angiosperm subclades using the angiosperm data set.

RESULTS

The androecium exhibits the highest evolutionary rates across most characters, whereas most perianth and gynoecium characters evolve more slowly in both Ericales and angiosperms. Both high and low rates of morphological evolution can result in high floral disparity in Ericales. Analyses of an angiosperm-wide floral data set reveal that this pattern appears to be conserved across most major angiosperm clades.

CONCLUSIONS

Elevated rates of morphological evolution in the androecium of Ericales may explain the higher disparity reported for this floral module. Comparing rates of morphological evolution through rate profiles proves to be a powerful tool in understanding floral evolution.

Intra-gastric phytoliths provide evidence for folivory in basal avialans of the Early Cretaceous Jehol Biota

Angiosperms became the dominant plant group in early to middle Cretaceous terrestrial ecosystems, coincident with the timing of the earliest pulse of bird diversification. While living birds and angiosperms exhibit strong interactions across pollination/nectivory, seed dispersal/frugivory, and folivory, documentation of the evolutionary origins and construction of that ecological complexity remains scarce in the Mesozoic. Through the first study of preserved in situ dietary derived phytoliths in a nearly complete skeleton of the early diverging avialan clade Jeholornithidae, we provide direct dietary evidence that Jeholornis consumed leaves likely from the magnoliid angiosperm clade, and these results lend further support for early ecological connections among the earliest birds and angiosperms. The broad diet of the early diverging avialan Jeholornis including at least fruits and leaves marks a clear transition in the early evolution of birds in the establishment of an arboreal (angiosperm) herbivore niche in the Early Cretaceous occupied largely by birds today. Morphometric reanalysis of the lower jaw of Jeholornis further supports a generalized morphology shared with other herbivorous birds, including an extant avian folivore, the hoatzin.

Plastome evolution and phylogenomic insights into the evolution of Lysimachia (Primulaceae: Myrsinoideae)

BACKGROUND

Lysimachia L., the second largest genus within the subfamily Myrsinoideae of Primulaceae, comprises approximately 250 species worldwide. China is the species diversity center of Lysimachia, containing approximately 150 species. Despite advances in the backbone phylogeny of Lysimachia, species-level relationships remain poorly understood due to limited genomic information. This study analyzed 50 complete plastomes for 46 Lysimachia species. We aimed to identify the plastome structure features and hypervariable loci of Lysimachia. Additionally, the phylogenetic relationships and phylogenetic conflict signals in Lysimachia were examined.

RESULTS

These fifty plastomes within Lysimachia had the typical quadripartite structure, with lengths varying from 152,691 to 155,784 bp. Plastome size was positively correlated with IR and intron length. Thirteen highly variable regions in Lysimachia plastomes were identified. Additionally, ndhB, petB and ycf2 were found to be under positive selection. Plastid ML trees and species tree strongly supported that L. maritima as sister to subg. Palladia + subg. Lysimachia (Christinae clade), while the nrDNA ML tree clearly placed L. maritima and subg. Palladia as a sister group.

CONCLUSIONS

The structures of these plastomes of Lysimachia were generally conserved, but potential plastid markers and signatures of positive selection were detected. These genomic data provided new insights into the interspecific relationships of Lysimachia, including the cytonuclear discordance of the position of L. maritima, which may be the result of ghost introgression in the past. Our findings have established a basis for further exploration of the taxonomy, phylogeny and evolutionary history within Lysimachia.

Molecular phylogeny and inflorescence evolution of Prunus (Rosaceae) based on RAD-seq and genome skimming analyses

Prunus is an economically important genus widely distributed in the temperate Northern Hemisphere. Previous studies on the genus using a variety of loci yielded conflicting phylogenetic hypotheses. Here, we generated nuclear reduced representation sequencing data and plastid genomes for 36 Prunus individuals and two outgroups. Both nuclear and plastome data recovered a well-resolved phylogeny. The species were divided into three main clades corresponding to their inflorescence types, - the racemose group, the solitary-flower group and the corymbose group - with the latter two sister to one another. Prunus was inferred to have diversified initially in the Late Cretaceous around 67.32 million years ago. The diversification of the three major clades began between the Paleocene and Miocene, suggesting that paleoclimatic events were an important driving force for Prunus diversification. Ancestral state reconstructions revealed that the most recent common ancestor of Prunus had racemose inflorescences, and the solitary-flower and corymb inflorescence types were derived by reduction of flower number and suppression of the rachis, respectively. We also tested the hybrid origin hypothesis of the racemose group proposed in previous studies. Prunus has undergone extensive hybridization events, although it is difficult to identify conclusively specific instances of hybridization when using SNP data, especially deep in the phylogeny. Our study provides well-resolved nuclear and plastid phylogenies of Prunus, reveals substantial cytonuclear discord at shallow scales, and sheds new light on inflorescence evolution in this economically important lineage.

The complete chloroplast genome of Ulmus mianzhuensis with insights into structural variations, adaptive evolution, and phylogenetic relationships of Ulmus (Ulmaceae)

BACKGROUND

Ulmus mianzhuensis is an endemic tree species in China with high ornamental and economic value. Currently, little is known regarding its genomic architecture, phylogenetic position, or adaptive evolution. Here, we sequenced the complete chloroplast genome (cp genome) of U. mianzhuensis and further compared the variations in gene organization and structure within Ulmus species to define their genomic evolution, then reconstructed the phylogenomic relationship of 31 related Ulmus species to explore the systematic position of U. mianzhuensis and the utility of cp genome for resolving phylogenetics among Ulmus species.

RESULTS

Our results revealed that all the Ulmus species exhibited a typical quadripartite structure, with a large single copy (LSC) region of 87,170 - 88,408 bp, a small single copy (SSC) region of 18,650 - 19,038 bp and an inverted repeat (IR) region of 26,288 - 26,546 bp. Within Ulmus species, gene structure and content of cp genomes were highly conserved, although slight variations were found in the boundary of SC/IR regions. Moreover, genome-wide sliding window analysis uncovered the variability of ndhC-trnV-UAC, ndhF-rpl32, and psbI-trnS-GCU were higher among 31 Ulmus that may be useful for the population genetics and potential DNA barcodes. Two genes (rps15 and atpF) were further detected under a positive selection of Ulmus species. Comparative phylogenetic analysis based on the cp genome and protein-coding genes revealed consistent topology that U. mianzhuensis is a sister group to U. parvifolia (sect. Microptelea) with a relatively low-level nucleotide variation of the cp genome. Additionally, our analyses also found that the traditional taxonomic system of five sections in Ulmus is not supported by the current phylogenomic topology with a nested evolutionary relationship between sections.

CONCLUSIONS

Features of the cp genome length, GC content, organization, and gene order were highly conserved within Ulmus. Furthermore, molecular evidence from the low variation of the cp genome suggested that U. mianzhuensis should be merged into U. parvifolia and regarded as a subspecies of U. parvifolia. Overall, we demonstrated that the cp genome provides valuable information for understanding the genetic variation and phylogenetic relationship in Ulmus.

Interspecific variation and phylogenetic relationship between mangrove and non-mangrove species of a same family (Meliaceae)-insights from comparative analysis of complete chloroplast genome

The mahogany family, Meliaceae, contains 58 genera with only one mangrove genus: Xylocarpus. Two of the three species of the genus Xylocarpus are true mangroves (X. granatum and X. moluccensis), and one is a non-mangrove (X. rumphii). In order to resolve the phylogenetic relationship between the mangrove and non-mangrove species, we sequenced chloroplast genomes of these Xylocarpus species along with two non-mangrove species of the Meliaceae family (Carapa guianensis and Swietenia macrophylla) and compared the genome features and variations across the five species. The five Meliaceae species shared 130 genes (85 protein-coding genes, 37 tRNA, and eight rRNA) with identical direction and order, with a few variations in genes and intergenic spacers. The repetitive sequences identified in the rpl22 gene region only occurred in Xylocarpus, while the repetitive sequences in accD were found in X. moluccensis and X. rumphii. The TrnH-GUG and rpl32 gene regions and four non-coding gene regions showed high variabilities between X. granatum and the two non-mangrove species (S. macrophylla and C. guianensis). In addition, among the Xylocarpus species, only two genes (accD and clpP) showed positive selection. Carapa guianensis and S. macrophylla owned unique RNA editing sites. The above genes played an important role in acclimation to different stress factors like heat, low temperature, high UV light, and high salinity. Phylogenetic analysis with 22 species in the order Sapindales supported previous studies, which revealed that the non-mangrove species X. rumphii is closer to X. moluccensis than X. granatum. Overall, our results provided important insights into the variation of genetic structure and adaptation mechanism at interspecific (three Xylocarpus species) and intergeneric (mangrove and non-mangrove genera) levels.

Insights into phylogenetic relationships in Pinus inferred from a comparative analysis of complete chloroplast genomes

BACKGROUND

Pinus is the largest genus of Pinaceae and the most primitive group of modern genera. Pines have become the focus of many molecular evolution studies because of their wide use and ecological significance. However, due to the lack of complete chloroplast genome data, the evolutionary relationship and classification of pines are still controversial. With the development of new generation sequencing technology, sequence data of pines are becoming abundant. Here, we systematically analyzed and summarized the chloroplast genomes of 33 published pine species.

RESULTS

Generally, pines chloroplast genome structure showed strong conservation and high similarity. The chloroplast genome length ranged from 114,082 to 121,530 bp with similar positions and arrangements of all genes, while the GC content ranged from 38.45 to 39.00%. Reverse repeats showed a shrinking evolutionary trend, with IRa/IRb length ranging from 267 to 495 bp. A total of 3,205 microsatellite sequences and 5,436 repeats were detected in the studied species chloroplasts. Additionally, two hypervariable regions were assessed, providing potential molecular markers for future phylogenetic studies and population genetics. Through the phylogenetic analysis of complete chloroplast genomes, we offered novel opinions on the genus traditional evolutionary theory and classification.

CONCLUSION

We compared and analyzed the chloroplast genomes of 33 pine species, verified the traditional evolutionary theory and classification, and reclassified some controversial species classification. This study is helpful in analyzing the evolution, genetic structure, and the development of chloroplast DNA markers in Pinus.

Phylogenetic analysis and temporal diversification of the tribe Alsineae (Caryophyllaceae) with the description of three new genera, Hesperostellaria, Reniostellaria and Torreyostellaria

Alsineae are one of the most taxonomically difficult tribes in Caryophyllaceae and consist of over 500 species distributed in the northern temperate zone. Recent phylogenetic results have improved our understanding on the evolutionary relationships among Alsineae members. Nevertheless, there are still some unresolved taxonomic and phylogenetic problems at the generic level, and the evolutionary history of major clades within the tribe was unexplored to date. In this study, we carried out phylogenetic analyses and divergence time estimation of Alsineae using the nuclear ribosomal internal transcribed spacer (nrITS) and four plastid regions (matK, rbcL, rps16, trnL-F). The present analyses yielded a robustly supported phylogenetic hypothesis of the tribe. Our results showed that the monophyletic Alsineae are strongly supported to be the sister of Arenarieae, and the inter-generic relationships within Alsineae were mostly resolved with strong support. Both molecular phylogenetic and morphological evidence supported the Asian species Stellaria bistylata and the two North American species Pseudostellaria jamesiana and Stellaria americana all should be recognized as new monotypic genera respectively, and three new genera Reniostellaria, Torreyostellaria, and Hesperostellaria were thereby proposed here. Additionally, molecular and morphological evidence also supported the proposal of the new combination Schizotechium delavayi. Nineteen genera were accepted within Alsineae and a key to these genera was provided. Molecular dating analysis suggested that Alsineae splitted from its sister tribe at ca. 50.2 million-years ago (Ma) during the early Eocene and began to diverge at ca. 37.9 Ma during the late Eocene, and divergent events within Alsineae occurred mainly since the late Oligocene. Results from the present study provide insights into the historical assembly of herbaceous flora in northern temperate regions.

The genome of Acorus deciphers insights into early monocot evolution

Acorales is the sister lineage to all the other extant monocot plants. Genomic resource enhancement of this genus can help to reveal early monocot genomic architecture and evolution. Here, we assemble the genome of Acorus gramineus and reveal that it has ~45% fewer genes than the majority of monocots, although they have similar genome size. Phylogenetic analyses based on both chloroplast and nuclear genes consistently support that A. gramineus is the sister to the remaining monocots. In addition, we assemble a 2.2 Mb mitochondrial genome and observe many genes exhibit higher mutation rates than that of most angiosperms, which could be the reason leading to the controversies of nuclear genes- and mitochondrial genes-based phylogenetic trees existing in the literature. Further, Acorales did not experience tau (τ) whole-genome duplication, unlike majority of monocot clades, and no large-scale gene expansion is observed. Moreover, we identify gene contractions and expansions likely linking to plant architecture, stress resistance, light harvesting, and essential oil metabolism. These findings shed light on the evolution of early monocots and genomic footprints of wetland plant adaptations.

Deconstructing age estimates for angiosperms

Estimates of the age of angiosperms from molecular phylogenies vary considerably. As in all estimates of evolutionary timescales from phylogenies, generating these estimates requires assumptions about the rate that molecular sequences are evolving (using clock models) and the time duration of the branches in a phylogeny (using fossil calibrations and branching processes). Often, it is difficult to demonstrate that these assumptions reflect current knowledge of molecular evolution or the fossil record. In this study we re-estimate the age of angiosperms using a minimal set of assumptions, therefore avoiding many of the assumptions inherent to other methods. The age estimates we generate are similar for each of the four datasets analysed, ranging from 130 to 400 Ma, but are far less precise than in previous studies. We demonstrate that this reduction in precision results from making less stringent assumptions about both rate and time, and that the analysed molecular dataset has very little effect on age estimates.

Comparitive Analysis of the Chloroplast Genomes of Three Houpoea Plants

The genus Houpoea belongs to the family Magnoliaceae, and the species in this genus have important medicinal values. However, the investigation of the correlation between the evolution of the genus and its phylogeny has been severely hampered by the unknown range of species within the genus and the paucity of research on its chloroplast genome. Thus, we selected three species of Houpoea: Houpoea officinalis var officinalis (OO), Houpoea officinalis var. biloba (OB), and Houpoea rostrata (R). With lengths of 160,153 bp (OO), 160,011 bp (OB), and 160,070 bp (R), respectively, the whole chloroplast genomes (CPGs) of these three Houpoea plants were acquired via Illumina sequencing technology, and the findings were annotated and evaluated. These three chloroplast genomes were revealed by the annotation findings to be typical tetrads. A total of 131, 132, and 120 different genes were annotated. The CPGs of the three species had 52, 47, and 56 repeat sequences, which were primarily found in the ycf2 gene. A useful tool for identifying species is the approximately 170 simple sequence repeats (SSRs) that have been found. The border area of the reverse repetition region (IR) was studied, and it was shown that across the three Houpoea plants, it is highly conservative, with only changes between H. rostrata and the other two plants observed. Numerous highly variable areas (rps3-rps19, rpl32-trnL, ycf1, ccsA, etc.) have the potential to serve as the barcode label for Houpoea, according to an examination of mVISTA and nucleotide diversity (Pi). Phylogenetic relation indicates that Houpoea is a monophyletic taxon, and its genus range and systematic position are consistent with the Magnoliaceae system of Sima Yongkang-Lu Shugang, including five species and varieties of H. officinalis var. officinalis, H. rostrata, H. officinalis var. biloba, Houpoea obovate, and Houpoea tripetala, which evolved and differentiated from the ancestors of Houpoea to the present Houpoea in the above order. This study provides valuable information on the genus Houpoea, enriches the CPG information on Houpoea genus, and provides genetic resources for the further classification of and phylogenetic research on Houpoea.

A new remarkable Early Cretaceous nelumbonaceous fossil bridges the gap between herbaceous aquatic and woody protealeans

Dating back to the late Early Cretaceous, the macrofossil record of the iconic lotus family (Nelumbonaceae) is one of the oldest of flowering plants and suggests that their unmistakable leaves and nutlets embedded in large pitted receptacular fruits evolved relatively little in the 100 million years since their first known appearance. Here we describe a new fossil from the late Barremian/Aptian Crato Formation flora (NE Brazil) with both vegetative and reproductive structures, Notocyamus hydrophobus gen. nov. et sp. nov., which is now the oldest and most complete fossil record of Nelumbonaceae. In addition, it displays a unique mosaic of ancestral and derived macro- and micromorphological traits that has never been documented before in this family. This new Brazilian fossil-species also provides a rare illustration of the potential morphological and anatomical transitions experienced by Nelumbonaceae prior to a long period of relative stasis. Its potential plesiomorphic and apomorphic features shared with Proteaceae and Platanaceae not only fill a major morphological gap within Proteales but also provide new support for their unexpected relationships first suggested by molecular phylogenies.

A complete chloroplast genome of a traditional Chinese medicine herb, Rubia podantha, and phylogenomics of Rubiaceae

Rubia podantha Diels is endemic to southwestern China and belongs to the family Rubiaceae. It is used in traditional Chinese medicines. To enrich the genetic data and resolve Rubiaceae's phylogeny, we assembled a complete chloroplast (cp) genome of R. podantha using Illumina HiSeq reads. The whole length of the cp genome was 154,866 bp. Annotation using PGA software found 113 genes, including 79 protein coding genes, 30 tRNA genes, and four rRNA genes. The large single-copy region was 84,717 bp, the inverted repeat B (IRa) region was 26,516 bp, the small single copy was 17,117 bp, and the inverted repeats B (IRb) region was 26,516 bp. Moreover, 64 SSRs were identified. Phylogenomic analysis using cp genomes of 109 Rubiaceae species found that R. podantha is closely related to R. cordifola. Rubiaceae was separated into three subfamilies: Ixoroideae, Cinchonoideae, and Rubiodeae. The genus Saprosma was not imbedded within the Spermacoceae alliance as reported in previous studies. Instead, it was imbedded within the Psychotrieae alliance. Divergence time estimation indicated that R. podantha split from its relative R. cordifolia around 1.25 million years ago. The assembled chloroplast genome in this study provided useful molecular information about the evolution of R. podantha and was a basis for phylogenetic analyses of Rubiaceae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01302-y.

A global phylogeny of butterflies reveals their evolutionary history, ancestral hosts and biogeographic origins

Butterflies are a diverse and charismatic insect group that are thought to have evolved with plants and dispersed throughout the world in response to key geological events. However, these hypotheses have not been extensively tested because a comprehensive phylogenetic framework and datasets for butterfly larval hosts and global distributions are lacking. We sequenced 391 genes from nearly 2,300 butterfly species, sampled from 90 countries and 28 specimen collections, to reconstruct a new phylogenomic tree of butterflies representing 92% of all genera. Our phylogeny has strong support for nearly all nodes and demonstrates that at least 36 butterfly tribes require reclassification. Divergence time analyses imply an origin ~100 million years ago for butterflies and indicate that all but one family were present before the K/Pg extinction event. We aggregated larval host datasets and global distribution records and found that butterflies are likely to have first fed on Fabaceae and originated in what is now the Americas. Soon after the Cretaceous Thermal Maximum, butterflies crossed Beringia and diversified in the Palaeotropics. Our results also reveal that most butterfly species are specialists that feed on only one larval host plant family. However, generalist butterflies that consume two or more plant families usually feed on closely related plants.

Chloroplast genome characteristics and phylogeny of the sinodielsia clade (apiaceae: apioideae)

BACKGROUND

The Sinodielsia clade of the subfamily Apioideae (Apiacieae) was established in 2008, and it is composed of 37 species from 17 genera. Its circumscription is still poorly delimited and unstable, and interspecific relationships in the clade lack comprehensive analysis. Chloroplast (cp.) genomes provide valuable and informative data sources for evolutionary biology and have been widely used in studies on plant phylogeny. To infer the phylogenetic history of the Sinodielsia clade, we assembled complete cp. genomes of 39 species and then performed phylogenetic analysis based on these cp. genome sequence data combined with 66 published cp. genomes from 16 genera relative to the Sinodielsia clade.

RESULTS

These 39 newly assembled genomes had a typical quadripartite structure with two inverted repeat regions (IRs: 17,599-31,486 bp) separated by a large single-copy region (LSC: 82,048-94,046 bp) and a small single-copy region (SSC: 16,343-17,917 bp). The phylogenetic analysis showed that 19 species were clustered into the Sinodielsia clade, and they were divided into two subclades. Six mutation hotspot regions were detected from the whole cp. genomes among the Sinodielsia clade, namely, rbcL-accD, ycf4-cemA, petA-psbJ, ycf1-ndhF, ndhF-rpl32 and ycf1, and it was found that ndhF-rpl32 and ycf1 were highly variable in the 105 sampled cp. genomes.

CONCLUSION

The Sinodielsia clade was subdivided into two subclades relevant to geographical distributions, except for cultivated and introduced species. Six mutation hotspot regions, especially ndhF-rpl32 and ycf1, could be used as potential DNA markers in the identification and phylogenetic analyses of the Sinodielsia clade and Apioideae. Our study provided new insights into the phylogeny of the Sinodielsia clade and valuable information on cp. genome evolution in Apioideae.

Phylogeny, Taxonomy and Evolutionary Trade-Offs in Reproductive Traits of Gomphoid Fungi (Gomphaceae, Gomphales)

Although functional ecology is a well-established field, our understanding of the evolutionary and ecological significance of the reproductive traits in macrofungi is still limited. Here, we reconstructed a phylogeny tree of gomphoid fungi in the narrower sense, including the species of the genera Gomphus and Turbinellus and used it to uncover the evolution of reproductive traits. Our analyses indicated that fungal fruit bodies and spores did not enlarge at a steady rate over time. Early gomphoid fungi essentially maintained their fruit body size, spore size and spore shape through the Mesozoic. In the Cenozoic, gomphoid fungi acquired significantly larger and more spherical spores by simultaneously expanding in length and width, with the fruit body size first decreasing and then enlarging. We argue that these trade-offs were driven by the effect of biological extinction and the dramatic climate changes of the Cenozoic. Gomphoid fungi initially increased in spore size and fruit body number as extinction survivors filled vacant niches. Both fruit bodies and spores eventually became larger as ecosystems saturated and competition intensified. One new species of Gomphus and nine new species of Turbinellus are described.

Genome-Wide Identification and Comparative Analysis of RALF Gene Family in Legume and Non-Legume Species

Rapid alkalinization factor (RALF) are small secreted peptide hormones that can induce rapid alkalinization in a medium. They act as signaling molecules in plants, playing a critical role in plant development and growth, especially in plant immunity. Although the function of RALF peptides has been comprehensively analyzed, the evolutionary mechanism of RALFs in symbiosis has not been studied. In this study, 41, 24, 17 and 12 RALFs were identified in Arabidopsis, soybean, Lotus and Medicago, respectively. A comparative analysis including the molecular characteristics and conserved motifs suggested that the RALF pre-peptides in soybean represented a higher value of isoelectric point and more conservative motifs/residues composition than other species. All 94 RALFs were divided into two clades according to the phylogenetic analysis. Chromosome distribution and synteny analysis suggested that the expansion of the RALF gene family in Arabidopsis mainly depended on tandem duplication, while segment duplication played a dominant role in legume species. The expression levels of most RALFs in soybean were significantly affected by the treatment of rhizobia. Seven GmRALFs are potentially involved in the release of rhizobia in the cortex cells. Overall, our research provides novel insights into the understanding of the role of the RALF gene family in nodule symbiosis.

Comparative plastome genomics, taxonomic delimitation and evolutionary divergences of Tetraena hamiensis var. qatarensis and Tetraena simplex (Zygophyllaceae)

The Zygophyllum and Tetraena genera are intriguingly important ecologically and medicinally. Based on morphological characteristics, T. hamiensis var. qatarensis, and T. simplex were transferred from Zygophyllum to Tetraena with the least genomic datasets available. Hence, we sequenced the T. hamiensis and T. simplex and performed in-depth comparative genomics, phylogenetic analysis, and estimated time divergences. The complete plastomes ranged between 106,720 and 106,446 bp-typically smaller than angiosperms plastomes. The plastome circular genomes are divided into large single-copy regions (~ 80,964 bp), small single-copy regions (~ 17,416 bp), and two inverted repeats regions (~ 4170 bp) in both Tetraena species. An unusual shrinkage of IR regions 16-24 kb was identified. This resulted in the loss of 16 genes, including 11 ndh genes which encode the NADH dehydrogenase subunits, and a significant size reduction of Tetraena plastomes compared to other angiosperms. The inter-species variations and similarities were identified using genome-wide comparisons. Phylogenetic trees generated by analyzing the whole plastomes, protein-coding genes, matK, rbcL, and cssA genes exhibited identical topologies, indicating that both species are sisters to the genus Tetraena and may not belong to Zygophyllum. Similarly, based on the entire plastome and proteins coding genes datasets, the time divergence of Zygophyllum and Tetraena was 36.6 Ma and 34.4 Ma, respectively. Tetraena stem ages were 31.7 and 18.2 Ma based on full plastome and protein-coding genes. The current study presents the plastome as a distinguishing and identification feature among the closely related Tetraena and Zygophyllum species. It can be potentially used as a universal super-barcode for identifying plants.

Cretaceous pollen cone with three-dimensional preservation sheds light on the morphological evolution of cycads in deep time

The Cycadales are an ancient and charismatic group of seed plants. However, their morphological evolution in deep time is poorly understood. While molecular divergence time analyses estimate a Cretaceous origin for most major living cycad clades, much of the extant diversity is inferred to be a result of Neogene diversifications. This leads to long branches throughout the cycadalean phylogeny that, with few exceptions, have yet to be rectified by unequivocal fossil cycads. We report a permineralized pollen cone from the Campanian Holz Shale located in Silverado Canyon, CA, USA (c. 80 million yr ago). This fossil was studied via serial sectioning, SEM, 3D reconstruction and phylogenetic analyses. Microsporophyll and pollen morphology indicate this cone is assignable to Skyttegaardia, a recently described genus based on disarticulated lignitized microsporophylls from the Early Cretaceous of Denmark. Data from this new species, including a simple cone architecture, anatomical details and vasculature organization, indicate cycadalean affinities for Skyttegaardia. Phylogenetic analyses support this assignment and recover Skyttegaardia as crown-group Cycadales, nested within Zamiaceae. Our findings support a Cretaceous diversification for crown-group Zamiaceae, which included the evolution of morphological divergent extinct taxa with unique traits that have yet to be widely identified in the fossil record.

Advances in the understanding of Blattodea evolution: Insights from phylotranscriptomics and spermathecae

Cockroaches, an ancient and diverse group of insects on earth that originated in the Carboniferous, displays a wide array of morphology or biology diversity. The spermatheca is an organ of the insect reproductive system; the diversity of spermathecae might be the adaption to different mating and sperm storage strategies. Yet a consensus about the phylogenetic relationships among the main lineages of Blattodea and the evolution of spermatheca has not been reached until now. Here we added the transcriptome data of Anaplectidae for the first time and supplemented other family level groups (such as Blaberidae, Corydiidae) to address the pending issues. Our results showed that Blattoidea was recovered as sister to Corydioidea, which was strongly supported by molecular evidence. In Blattoidea, (Lamproblattidae + Anaplectidae) + (Cryptocercidae + Termitoidae) was strongly supported by our molecular data. In Blaberoidea, Pseudophyllodromiidae and Blaberidae were recovered to be monophyletic, while Blattellidae was found to be paraphyletic with respect to Malaccina. Ectobius sylvestris + Malaccina discoidalis formed the sister group to other Blaberoidea; Blattellidae (except Malaccina discoidalis) + Nyctiboridae was found as the sister of Blaberidae. Corydiidae was recovered to be non-monophyletic due to the embedding of Nocticola sp. Our ASR analysis of spermatheca suggested that primary spermathecae were present in the common ancestor, and it transformed at least six times during the evolutionary history of Blattodea. The evolution of spermatheca could be described as a unidirectional trend: the increased size to accommodate more sperm. Furthermore, major splits within the existing genera of cockroaches occurred in the Upper Paleogene to Neogene. Our study provides strong support for the relationship among three superfamilies and offers some new insights into the phylogeny of cockroaches. Meanwhile, this study also provides basic knowledge on the evolution of spermathecae and reproductive patterns.

Out of Madagascar, repeatedly: The phylogenetics and biogeography of Dombeyoideae (Malvaceae s.l.)

Dispersals have been shown to be critical to the evolution of the long isolated but megadiverse flora of Madagascar and the surrounding islands of the western Indian Ocean, but we are just beginning to understand the directionality of these dispersals. With more than half of its species occurring in the western Indian Ocean region (WIOR), the paleotropical subfamily Dombeyoideae provides a particularly useful case study through which to better understand the biogeography of the WIOR, and yet its biogeography is poorly understood. Here we sampled six molecular markers from all 20 genera in the Dombeyoideae to reconstruct the most complete phylogeny to date for the subfamily. From this, divergence times, calibrated with three fossils (two dombeyoid, one malvoid), and ancestral range estimations were hypothesized. Biogeographic stochastic mapping (BSM) analyses on the maximum clade credibility tree were completed and compared to BSM analyses on 1,000 trees randomly sampled from the posterior distribution of trees resulting from the dating analysis. We found the Dombeyoideae crown node diverged ca. 53 million years ago out of a broad ancestral range involving all three major areas of its distribution: Madagascar, Africa, and Asia. The majority of diversification and dispersals in the subfamily occurred within the last ca. 10 million years, mostly from the Pliocene onwards. There were roughly five dispersals from Madagascar to Africa (and only one in reverse), at least six from Madagascar to surrounding islands of the WIOR (Mascarenes and Comoros), and one dispersal from Madagascar to Asia (and ca. 1 in reverse). Other long-distance dispersals included one from Africa to St. Helena and one from Africa to Australasia, both from within the most widespread clade, the Cheirolaena & allies clade, and one dispersal from Asia to Africa. Critically, the Dombeyoideae provide strong evidence for considering the island of Madagascar as a source for the colonization of continents, as well as the surrounding islands of the WIOR. Furthermore, narrow sympatry was a key process in the evolution of the subfamily, particularly in Madagascar and the Mascarenes.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Comprehensive phylogenetic analyses of Orchidaceae using nuclear genes and evolutionary insights into epiphytism

Orchidaceae (with >28,000 orchid species) are one of the two largest plant families, with economically and ecologically important species, and occupy global and diverse niches with primary distribution in rainforests. Among orchids, 70% grow on other plants as epiphytes; epiphytes contribute up to ~50% of the plant diversity in rainforests and provide food and shelter for diverse animals and microbes, thereby contributing to the health of these ecosystems. Orchids account for over two-thirds of vascular epiphytes and provide an excellent model for studying evolution of epiphytism. Extensive phylogenetic studies of Orchidaceae and subgroups have ;been crucial for understanding relationships among many orchid lineages, although some uncertainties remain. For example, in the largest subfamily Epidendroideae with nearly all epiphytic orchids, relationships among some tribes and many subtribes are still controversial, hampering evolutionary analyses of epiphytism. Here we obtained 1,450 low-copy nuclear genes from 610 orchid species, including 431 with newly generated transcriptomes, and used them for the reconstruction of robust Orchidaceae phylogenetic trees with highly supported placements of tribes and subtribes. We also provide generally well-supported phylogenetic placements of 131 genera and 437 species that were not sampled by previous plastid and nuclear phylogenomic studies. Molecular clock analyses estimated the Orchidaceae origin at ~132 million years ago (Ma) and divergences of most subtribes from 52 to 29 Ma. Character reconstruction supports at least 14 parallel origins of epiphytism; one such origin was placed at the most recent common ancestor of ~95% of epiphytic orchids and linked to modern rainforests. Ten occurrences of rapid increase in the diversification rate were detected within Epidendroideae near and after the K-Pg boundary, contributing to ~80% of the Orchidaceae diversity. This study provides a robust and the largest family-wide Orchidaceae nuclear phylogenetic tree thus far and new insights into the evolution of epiphytism in vascular plants.

Revisiting the phylogeny of the family Miridae (Heteroptera: Cimicomorpha), with updated insights into its origin and life history evolution

Heteroptera is one of the most successfully adapted groups on Earth and can be observed in almost every environment. Within the evolution of heteropteran insects, Miridae show remarkable diversity (>11,700 spp.), accounting for a quarter of all Heteroptera. However, their phylogeny is still unclear, and no plausible theory for the driving force of their diversification has been established. In this work, we provide new suggestions for the phylogeny of Miridae using a larger dataset than previous studies. In addition, we suggest an alternative evolutionary history based on newly calibrated divergence dates for Miridae and its subordinate groups, and present probable factors of the family's success in terms of species diversity. The entire dataset comprises 16 outgroups and 188 ingroup taxa including all seven known subfamilies and 37 out of 45 known tribes. Each species is aligned as 3,577bp with six molecular loci (COI, 16S rRNA, 18S rRNA, 28S rRNA D3 region, H2A, and H3A). Among the molecular markers, we are the first to test histone genes (H2A, H3A) in Miridae. Our results raise the following points about phylogenetic relationships: i) The earliest group to diverge from Miridae was Monaloniini (Bryocorinae). ii) Bryocorinae and Cylapinae are polyphyletic, Deraeocorinae and Orthotylinae also rendered as non-monophyletic group. iii) Termatophylini and Coridromiini separated from Deraeocorinae and Orthotylinae respectively. iv) Four large tribes, Orthotylini, Phylini, Deraeocorini and Mirini are non-monophyletic. The results from our ancestral state reconstruction and divergence date estimation suggest the following: i) Miridae first diverged during the Late Jurassic (approx. 163.4 Mya), and the divergence dates of most subfamilies and tribes overlap with angiosperm radiation, which perhaps synergized their diversification. ii) Ancestral reconstruction results for Miridae reveal it to be predominantly phytophagous and diverge to oligophagy mainly in plant-tissue habitats, which could have allowed the mirids to select optimal tactics as plant-dwellers. iii) The common ancestor of Miridae originated among plant-dwellers mainly on Eudicots, and that tendency was largely maintained, but sporadic host shifts also occurred.

Insect pollination in deep time

Inferring insect pollination from compression fossils and amber inclusions is difficult because of a lack of consensus on defining an insect pollinator and the challenge of recognizing this ecological relationship in deep time. We propose a conceptual definition for such insects and an operational classification into pollinator or presumed pollinator. Using this approach, we identified 15 insect families that include fossil pollinators and show that pollination relationships have existed since at least the Upper Jurassic (~163 Ma). Insects prior to this can only be classified as presumed pollinators. This gives a more nuanced insight into the origin and evolution of an ecological relationship that is vital to the establishment, composition and conservation of modern terrestrial ecosystems.

Is Santaniella a ranunculid? Reassessment of this enigmatic fossil angiosperm from the Lower Cretaceous (Aptian, Crato Konservat-Lagerstätte, Brazil) provides a new interpretation

PREMISE

The Lower Cretaceous Crato Konservat-Lagerstätte (CKL) preserves a rich flora that includes early angiosperms from northern Gondwana. From this area, the recently described fossil genus Santaniella was interpreted as a ranunculid (presumably Ranunculaceae). However, based on our examination of an additional specimen and a new phylogenetic analysis, we offer an alternative interpretation.

METHODS

The new fossil was collected from an active quarry for paving stones in the state of Ceará, northeastern Brazil. We assessed support for alternative phylogenetic hypotheses using a combined analysis of morphological data and DNA sequence data using Bayesian inference. We used a consensus network to visualize the posterior distribution of trees, and we used RoguePlot to illustrate the support for alternative positions on a scaffold tree.

RESULTS

The new material includes a flower-like structure not present in the original material and also includes follicles preserved at early stages of development. The flower-like structure is a compact terminal cluster of elliptical sterile laminar organs surrounding internal filamentous structures that occur on flexuous axes. Phylogenetic analyses did not support the fossil placement among eudicots. Instead, Santaniella appears to belong in the magnoliid clade.

CONCLUSIONS

The presence of seeds in a marginal-linear placentation and enclosed in a follicle supports the fossil as an angiosperm. However, even though most characters are clearly recognizable, its combination of characters does not provide strong support for a close relationship to any extant order of flowering plants. Its position in the magnoliid clade is intriguing and, based on plicate carpels, it is definitely a mesangiosperm.

The use of Anchored Hybrid Enrichment data to resolve higher-level phylogenetic relationships: A proof-of-concept applied to Asterales (Eudicotyledoneae; Angiosperms)

Anchored Hybrid Enrichment (AHE) is a tool for capturing orthologous regions of the nuclear genome shared in low or single copy across lineages. Despite the increasing number of studies using this method, its usefulness to estimate relationships at deeper taxonomic levels in plants has not been fully explored. Here we present a proof of concept about the performance of nuclear loci obtained with AHE to infer phylogenetic relationships and explore the use of gene sampling schemes to estimate divergence times in Asterales. We recovered low-copy nuclear loci using the AHE method from herbarium material and silica-preserved samples. Maximum likelihood, Bayesian inference, and coalescence approaches were used to reconstruct phylogenomic relationships. Dating analyses were conducted under a multispecies coalescent approach by jointly inferring species tree and divergence times with random gene sampling schemes and multiple calibrations. We recovered 403 low-copy nuclear loci for 63 species representing nine out of eleven families of Asterales. Phylogenetic hypotheses were congruent among the applied methods and previously published results. Analyses with concatenated datasets were strongly supported, but coalescence-based analyses showed low support for the phylogenetic position of families Argophyllaceae and Alseuosmiaceae. Estimated family ages were congruent among gene sampling schemes, with the mean age for Asterales around 130 Myr. Our study documents the usefulness of AHE for resolving phylogenetic relationships at deep phylogenetic levels in Asterales. Observed phylogenetic inconsistencies were possibly due to the non-inclusion of families Phellinceae and Pentaphragmataceae. Random gene sampling schemes produced consistent age estimates with coalescence and species tree relaxed clock approaches.

A complete genus-level phylogeny reveals the Cretaceous biogeographic diversification of the poppy family

Angiosperms, a trigger for the Cretaceous Terrestrial Revolution (KTR), underwent a rapid expansion and occupied all the environments during the Mid-Upper Cretaceous. Yet, Cretaceous biogeographic patterns and processes underlying the distribution of angiosperm diversity in the Northern Hemisphere are still poorly known. Here, we elucidated the biogeographic diversification of the angiosperm family Papaveraceae, an ancient Northern Hemisphere clade characterized by poor dispersal ability and high level of regional endemism. Based on both plastome and multi-locus datasets, we reconstructed a robust time-calibrated phylogeny that includes all currently recognized 45 genera of this family. Within the time-calibrated phylogenetic framework, we conducted 72 biogeographic analyses by testing the sensitivity of uncertainties of area delimitation, maxarea constraints, and the parameters of the model, i.e., j (describing jump-dispersal events) and w (modifying dispersal multiplier matrices), to ancestral range estimations. We also inferred ancestral habitat and ecological niches. Phylogenetic analyses strongly support Papaveraceae as monophyletic. Pteridophylloideae is strongly supported as sister to Hypecoideae-Fumarioideae. Our results indicate that the j parameter and number of predefined areas strongly affect ancestral range estimates, generating questionable ancestral ranges, whereas maxarea constraint and w parameter have no effect and improve model fit. After accounting for these uncertainties, our results indicate that Papaveraceae differentiated in Asian wet forests during the Lower Cretaceous and subsequently occupied the Asian and western North American arid and open areas. Three dispersals from Asia to western North America via the Bering land bridge occurred in the Mid-Upper Cretaceous, largely in agreement with the KTR. Habitat shift and ecological niche divergence resulted in the subsequent disjunctions between Asia and western North America. These findings suggest that the interplay of range expansion and niche divergence-driven vicariance might have shaped Cretaceous biogeographic patterns of angiosperms with Papaveraceae-like ecological requirements and dispersal abilities in the Northern Hemisphere, hence contributing to the knowledge on the geographic expansion of angiosperms during the KTR.

Atypical ‘long-tailed’ cockroaches arose during Cretaceous in response to angiosperm terrestrial revolution

Typical cockroaches are flat, broad, with large pronotum and wings covering the body. This conserved morphotype dates back to the Carboniferous, during which the ancestral cockroaches, or roachoids, originated. On the other hand, the ovipositor of cockroaches gradually reduced during the Mesozoic, coupled with a major shift of reproductive strategy. By the Cretaceous, long external ovipositors became rare, most cockroaches used very short or even hidden internal ovipositors to fabricate egg cases (oothecae), which is an innovation for egg protection. Here, we describe two cockroaches from mid-Cretaceous Myanmar amber: Ensiferoblatta oecanthoides gen. et sp. nov. (Ensiferoblattidae fam. nov.) and Proceroblatta colossea gen. et sp. nov. They are slim, elongate, fusiform, with longitudinal pronotum, and have long external ovipositors. The combination of these traits represents a unique morphotype, which resembles crickets and katydids (Ensifera) more than general cockroaches. Ensiferoblatta and Proceroblatta may be arboreal, feeding on and/or laying eggs into certain angiosperms that newly emerged. Their open habit causes latent impairment to viability, and may contribute to their extinction. These new taxa are the youngest members of the ancient, extinct group of cockroaches, namely Eoblattodea, which are characterized by long ovipositors. We speculate that the extinction of certain gymnosperm hosts almost ended the 200-My triumph of Eoblattodea. Despite an attempt to adapt to angiosperm hosts, Ensiferoblatta, Proceroblatta and suchlike cockroaches as an evolutionary dead end failed to save Eoblattodea from extinction. The lack of protection for eggs (maternal care in particular) might accelerate the extinction of Eoblattodea as a whole.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences

BACKGROUND

Plant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates.

RESULTS

The results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP.

CONCLUSION

The study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species.

Evolutionary history of two evergreen Rhododendron species as revealed by chromosome-level genome assembly

BACKGROUND

The genus Rhododendron (Ericaceae), a species-rich and widely distributed genus of woody plants, is distinguished for the beautiful and diverse flowers. Rhododendron delavayi Franch. and Rhododendron irroratum Franch., are highly attractive species widely distributed in south-west China and abundant new varieties have been selected from their genetic resources.

METHODS

We constructed chromosome-scale genome assemblies for Rhododendron delavayi and Rhododendron irroratum. Phylogenetic and whole-genome duplication analyses were performed to elucidate the evolutionary history of Rhododendron. Further, different types of gene duplications were identified and their contributions to gene family expansion were investigated. Finally, comprehensive characterization and evolutionary analysis of R2R3-MYB and NBS-encoding genes were conducted to explore their evolutionary patterns.

RESULTS

The phylogenetic analysis classified Rhododendron species into two sister clades, 'rhododendrons' and 'azaleas'. Whole-genome duplication (WGD) analysis unveiled only one WGD event that occurred in Rhododendron after the ancestral γ triplication. Gene duplication and gene family expansion analyses suggested that the younger tandem and proximal duplications contributed greatly to the expansion of gene families involved in secondary metabolite biosynthesis and stress response. The candidate R2R3-MYB genes likely regulating anthocyanin biosynthesis and stress tolerance in Rhododendron will facilitate the breeding for ornamental use. NBS-encoding genes had undergone significant expansion and experienced species-specific gain and loss events in Rhododendron plants.

CONCLUSIONS

The reference genomes presented here will provide important genetic resources for molecular breeding and genetic improvement of plants in this economically important Rhododendron genus.

New Fossil Evidence Suggests That Angiosperms Flourished in the Middle Jurassic

Angiosperms are a group of plants with the highest rate of evolution, the largest number of species, the widest distribution and the strongest adaptability. Needless to say, angiosperms are the most important group for the humans. The studies on the origin, evolution and systematics of angiosperms have been the major challenges in plant sciences. However, the origin and early history of angiosperms remains poorly understood and controversial among paleobotanists. Some paleobotanists insist that there were no angiosperms in the pre-Cretaceous age. However, this conclusion is facing increasing challenges from fossil evidence, especially Early Jurassic Nanjinganthus, which is based on over two hundred specimens of fossil flowers. Studying more fossil plants is the only reliable way to elucidate the origin and early evolution of angiosperms. Here, we document a new species of angiosperms, Qingganninginfructus formosa gen. et sp. nov, and provide the first detailed three-dimensional morphology of Qingganninginfructus gen. nov from the Middle Jurassic of Northwest China. A Micro-CT examination shows that the best-preserved fossil infructescence has eleven samaroid fruits, each with a single basal ovule. Since these fossils are distinct in morphology and organization from all organs of known gymnosperms and angiosperms (the latter are defined by their enclosed ovules), we interpret Qingganninginfructus as a new genus of angiosperms including a new species, Q. formosa gen. et sp. nov., and an unspecified species from the Middle Jurassic of Northwest China. The discovery of this new genus of angiosperms from the Middle Jurassic, in addition to the existing records, undermines the "no angiosperms until the Cretaceous" stereotype and updates the perspective on the origin and early history of angiosperms.

The phylogeny and evolutionary ecology of hoverflies (Diptera: Syrphidae) inferred from mitochondrial genomes

Hoverflies (Diptera: Syrphidae) are a diverse group of pollinators and a major research focus in ecology, but their phylogenetic relationships remain incompletely known. Using a genome skimming approach we generated mitochondrial genomes for 91 species, capturing a wide taxonomic diversity of the family. To reduce the required amount of input DNA and overall cost of the library construction, sequencing and assembly was conducted on mixtures of specimens, which raises the problem of chimera formation of mitogenomes. We present a novel chimera detection test based on gene tree incongruence, but identified only a single mitogenome of chimeric origin. Together with existing data for a final set of 127 taxa, phylogenetic analysis on nucleotide and amino acid sequences using Maximum Likelihood and Bayesian Inference revealed a basal split of Microdontinae from all other syrphids. The remainder consists of several deep clades assigned to the subfamily Eristalinae in the current classification, including a clade comprising the subfamily Syrphinae (plus Pipizinae). These findings call for a re-definition of subfamilies, but basal nodes had insufficient support to allow such action. Molecular-clock dating placed the origin of the Syrphidae crown group in the mid-Cretaceous while the Eristalinae-Syrphinae clade likely originated near the K/Pg boundary. Transformation of larval life history characters on the tree suggests that Syrphidae initially had sap feeding larvae, which diversified greatly in diet and habitat association during the Eocene and Oligocene, coinciding with the diversification of angiosperms and the evolution of various insect groups used as larval host, prey, or mimicry models. Mitogenomes proved to be a powerful phylogenetic marker for studies of Syrphidae at subfamily and tribe levels, allowing dense taxon sampling that provided insight into the great ecological diversity and rapid evolution of larval life history traits of the hoverflies.