Divergence times, historical biogeography, and shifts in speciation rates of Myrtales

Disconnecting trapdoors: Phylogenomic analyses reveal evolutionary contrasts in trapdoor spiders with intercontinental distribution (Idiopidae, Idiopinae)

Among Mygalomorphae spiders, the family Idiopidae is the second most diverse, consisting exclusively of trapdoor spiders and is divided into three subfamilies: Arbanitinae, Genysinae, and Idiopinae. The subfamily Idiopinae, distinguished mainly by anterior lateral eyes that project forward, includes 153 species across seven genera, distributed throughout South America, Africa, and parts of Asia. Within this subfamily, the genus Idiops includes the greatest diversity and is the only genus recorded in both the New and Old Worlds. Utilizing a taxon set from the Neotropical and Afrotropical regions, with specimens collected from 1947 to 2021, our study provides the first phylogenomic analysis of the family based on Ultraconserved Elements (UCEs). To assess the monophyly of Idiopinae and Idiops, as well as the relationships among genera within the subfamily, we conducted phylogenetic analyses employing maximum likelihood, Bayesian inference, and coalescent-based methods. The phylogenetic trees reveal that Idiopinae forms a monophyletic lineage, split into two geographically distinct groups: one with African species and the other with Neotropical species. We did not recover monophyly of the genus Idiops; Neotropical Idiops form a monophyletic lineage, while African Idiops species group with Titanidiops, forming a sister lineage to the remaining African Idiopinae. The relationship between the phylogenetic results obtained and the main morphological differences observed among the genera is discussed. Finally, our study challenges the monophyly of Idiopidae by including Neocteniza, which is found to be an independent lineage sister to Ctenizidae and the rest of Idiopidae.

New insights into rhizosphere bacterial community shaped by lettuce genotypes for divergent degradation efficiencies of phthalates

Rhizosphere dissipation of organic pollutants benefits safe utilization of the polluted agricultural soil. Nevertheless, dissipation variation of phthalates (PAEs) in rhizosphere among different vegetable genotypes and the related microbial mechanisms remain unknown. Here, twelve lettuce cultivars with different genetic relationships identified by 18S rRNA gene sequencing were cultivated in soil spiked with di-(2-ethylhexyl) phthalate (DEHP). Bacterial communities and function genes in rhizosphere of lettuce were analyzed by 16S rRNA gene and metagenomic sequencing. Results showed significant variations in DEHP concentrations of roots (2.8-15.3 mg/kg) and shoots (0.70-1.8 mg/kg) among 12 cultivars. Notably, cultivars L11 and L12 showed the lowest DEHP accumulation in roots and shoots, being lower by 82 % and 58 % than the highest accumulators (cultivars L5 and L6), respectively. This accumulation variation was closely connected with their genetic relationships and exhibited genotype-dependent trait. The significantly different bacterial community diversities and structures were recorded in rhizosphere among 12 cultivars. Especially, bacterial communities in rhizosphere of cultivars L11 and L12 (low-DEHP accumulators with high DEHP dissipation) strengthened their adaptation by enriching pollutant-resistant taxa, increasing extracellular polymeric substance contents and biofilm formation, as well as constructing complex ecological networks under DEHP pollution. Moreover, PAE-degrading bacteria and genes (e.g., hydrolase65, phtAb, and pcaI) in rhizosphere were enriched by low-DEHP accumulators, which benefited DEHP removal and subsequently safe agricultural products. This study provides new insights into microbial mechanisms on rhizosphere DEHP degradation and its correlation with accumulation variation among different crop genotypes.

Continental-scale interactions of Australian showy mistletoes and their hosts

PREMISE

Showy mistletoes are obligate hemiparasites of woody plants. Host specificity is therefore a fundamental determinant of mistletoe diversity, persistence, geographic distribution, and abundance. Investigations of host specificity in Australian Loranthaceae have focused mostly on host range (taxon counts), but additional insights into specificity are gained by quantifying mistletoe prevalence on taxa in their host range and by exploring specificity in a phylogenetic context.

METHODS

We estimated measures of host specificity to characterize mistletoe-host interactions at a continental scale by using occurrence records in the Atlas of Living Australia. We calculated host taxon richness, mistletoe prevalence, and phylogenetic diversity, and used rarefaction curves to evaluate sampling coverage.

RESULTS

Many mistletoe taxa were represented by few records that listed the host, which often was identified to genus only. Mistletoe genera were recorded on 29 orders and 80 families, and no association was observed between host richness and number of records per genus. Rarefaction curves suggested that additional host orders and families remain to be discovered for Amylotheca, Decaisnina, Dendrophthoe, and Muellerina. Four mistletoe genera were most prevalent on Myrtales, one on Fabales, and one on Laurales. Rosids were most often the recorded hosts (84.3% of all records). We found evidence of significant phylogenetic clustering in host use by Amyema, Amylotheca, and Decasinina.

CONCLUSIONS

Our results, particularly the high prevalence on rosids, suggest that relationships of mistletoes with rainforest lineages may have been established early in the history of Australian Loranthaceae and that some lineages co-diversified with their hosts in arid regions.

The complete sequence of chloroplast genome of Baeckea frutescens Linaeus 1753 (Myrtoideae), a traditional folk medicinal plant

Baeckea frutescens Linaeus 1753, as a traditional folk medicine in South East Asia, possesses sesquiterpenes, phloroglucinols, chromones, and essential oil, and is utilized for traditional Chinese medicinal purposes. The genetic diversity of the plant must be better understood, considering its significance. The complete chloroplast (cp) genome of B. frutescens was sequenced and assembled by using Illumina paired-end data, marking a significant advancement toward comprehending its genetic composition. The complete cp genome is 158,939 bp in length and contains 128 genes, consisting of 83 protein-coding genes, 8 ribosomal RNA genes, and 37 transfer RNA genes. Phylogenetic analyses indicated that B. frutescens and other the 13 were clustered to the family of Myrtaceae. These findings are crucial for the conservation and utilization of this important plant species. Additionally, they underscore the potential for future research on the evolution and preservation of B. frutescens, which could be advantageous in pharmaceutical applications.

YOLO-Granada: a lightweight attentioned Yolo for pomegranates fruit detection

Pomegranate is an important fruit crop that is usually managed manually through experience. Intelligent management systems for pomegranate orchards can improve yields and address labor shortages. Fast and accurate detection of pomegranates is one of the key technologies of this management system, crucial for yield and scientific management. Currently, most solutions use deep learning to achieve pomegranate detection, but deep learning is not effective in detecting small targets and large parameters, and the computation speed is slow; therefore, there is room for improving the pomegranate detection task. Based on the improved You Only Look Once version 5 (YOLOv5) algorithm, a lightweight pomegranate growth period detection algorithm YOLO-Granada is proposed. A lightweight ShuffleNetv2 network is used as the backbone to extract pomegranate features. Using grouped convolution reduces the computational effort of ordinary convolution, and using channel shuffle increases the interaction between different channels. In addition, the attention mechanism can help the neural network suppress less significant features in the channels or space, and the Convolutional Block Attention Module attention mechanism can improve the effect of attention and optimize the object detection accuracy by using the contribution factor of weights. The average accuracy of the improved network reaches 0.922. It is only less than 1% lower than the original YOLOv5s model (0.929) but brings a speed increase and a compression of the model size. and the detection speed is 17.3% faster than the original network. The parameters, floating-point operations, and model size of this network are compressed to 54.7%, 51.3%, and 56.3% of the original network, respectively. In addition, the algorithm detects 8.66 images per second, achieving real-time results. In this study, the Nihui convolutional neural network framework was further utilized to develop an Android-based application for real-time pomegranate detection. The method provides a more accurate and lightweight solution for intelligent management devices in pomegranate orchards, which can provide a reference for the design of neural networks in agricultural applications.

Comparative genomics revealed new insights into the plastome evolution of Ludwigia (Onagraceae, Myrtales)

The primrose-willow (Ludwigia L.), a well-defined genus of the Onagraceae family, comprises 87 species widely distributed worldwide. In this study, we sequenced and characterized the complete chloroplast (cp) genomes of three species in the genus, including Ludwigia adscendens, Ludwigia hyssopifolia, and Ludwigia prostrata. Three Ludwigia cp genomes ranged from 158,354 to 159,592 bp in size, and each contained 113 genes, including 79 unique protein-coding genes (PCGs), four rRNA genes, and 30 tRNA genes. A comparison of the Ludwigia cp genomes revealed that they were highly conserved in gene composition, gene orientation, and GC content. Moreover, we compared the structure of cp genomes and reconstructed phylogenetic relationships with related species in the Onagraceae family. Regarding contraction/expansion of inverted repeat (IR) region, two kinds of expansion IR region structures were found in Oenothera, Chamaenerion, and Epilobium genera, with primitive IR structures in Ludwigia and Circeae genera. The regions clpP, ycf2, and ycf1 genes possessed highly divergent nucleotides among all available cp genomes of the Onagraceae family. The phylogenetic reconstruction using 79 PCGs from 39 Onagraceae cp genomes inferred that Ludwigia (including L. adscendens, L. hyssopifolia, L. prostrata, and Ludwigia octovalvis) clade was monophyletic and well-supported by the bootstrap and posterior probability values. This study provides the reference cp genomes of three Ludwigia species, which can be used for species identification and phylogenetic reconstruction of Ludwigia and Onagraceae taxa.

Estimating ancestral ranges and biogeographical processes in early hominins

Historical biogeography provides crucial insights into understanding the evolutionary history of hominins. We applied maximum-likelihood and biogeographical stochastic mapping to infer the ancestral ranges of hominins and estimate the frequency of biogeographical events. These events were inferred using two time-calibrated phylogenetic trees that differ in the position of Australopithecus sediba. Results suggest that regardless of which phylogeny was selected, Northcentral Africa was the preferred ancestral region for the ancestor of the Homo-Pan clade, as well as the ancestor of Sahelanthropus and later hominins. The northern and middle part of eastern Africa was the preferred ancestral region for several clades originating at subsequent deep nodes of the trees (∼5-4 Ma). The choice of tree topology had one important effect on results: whether hominin ancestors appearing after ∼4 Ma were widespread or endemic. These different patterns highlight the biogeographic significance of the phylogenetic relationships of A. sediba. Overall, the results showed that dispersal, local extinction, and sympatry played vital roles in creating the hominin distribution, whereas vicariance and jump dispersal were not as common. The results suggested symmetry in the directionality of dispersals. Distance probably influenced how rapidly taxa colonized a new region, and dispersals often followed the closest path. These findings are potentially impacted by the imperfection of the fossil record, suggesting that the results should be interpreted cautiously.

De novo genome assembly and population genomics of a shrub tree Barthea barthei (Hance) krass provide insights into the adaptive color variations

Flower color is a classic example of an ecologically important trait under selection in plants. Understanding the genetic mechanisms underlying shifts in flower color can provide key insights into ecological speciation. In this study, we investigated the genetic basis of flower color divergence in Barthea barthei, a shrub tree species exhibiting natural variation in flower color. We assembled a high-quality genome assembly for B. barthei with a contig N50 of 2.39 Mb and a scaffold N50 of 16.21 Mb. The assembly was annotated with 46,430 protein-coding genes and 1,560 non-coding RNAs. Genome synteny analysis revealed two recent tetraploidization events in B. barthei, estimated to have occurred at approximately 17 and 63 million years ago. These tetraploidization events resulted in massive duplicated gene content, with over 70% of genes retained in collinear blocks. Gene family members of the core regulators of the MBW complex were significantly expanded in B. barthei compared to Arabidopsis, suggesting that these duplications may have provided raw genetic material for the evolution of novel regulatory interactions and the diversification of anthocyanin pigmentation. Transcriptome profiling of B. barthei flowers revealed differential expression of 9 transcription factors related to anthocyanin biosynthesis between the two ecotypes. Six of these differentially expressed transcription factors were identified as high-confidence candidates for adaptive evolution based on positive selection signals. This study provides insights into the genetic basis of flower color divergence and the evolutionary mechanisms underlying ecological adaptation in plants.

Pangenome of water caltrop reveals structural variations and asymmetric subgenome divergence after allopolyploidization

Water caltrop (Trapa spp., Lythraceae) is a traditional but currently underutilized non-cereal crop. Here, we generated chromosome-level genome assemblies for the two diploid progenitors of allotetraploid Trapa. natans (4x, AABB), i.e., diploid T. natans (2x, AA) and Trapa incisa (2x, BB). In conjunction with four published (sub)genomes of Trapa, we used gene-based and graph-based pangenomic approaches and a pangenomic transposable element (TE) library to develop Trapa genomic resources. The pangenome displayed substantial gene-content variation with dispensable and private gene clusters occupying a large proportion (51.95%) of the total cluster sets in the six (sub)genomes. Genotyping of presence-absence variation (PAVs) identified 40 453 PAVs associated with 2570 genes specific to A- or B-lineages, of which 1428 were differentially expressed, and were enriched in organ development process, organic substance metabolic process and response to stimulus. Comparative genome analyses showed that the allotetraploid T. natans underwent asymmetric subgenome divergence, with the B-subgenome being more dominant than the A-subgenome. Multiple factors, including PAVs, asymmetrical amplification of TEs, homeologous exchanges (HEs), and homeolog expression divergence, together affected genome evolution after polyploidization. Overall, this study sheds lights on the genome architecture and evolution of Trapa, and facilitates its functional genomic studies and breeding program.

A comprehensive genus-level phylogeny and biogeographical history of the Lythraceae based on whole plastome sequences

BACKGROUND AND AIMS

The Lythraceae are a mainly subtropical to tropical family of the order Myrtales with 28 currently accepted genera and approximately 600 species. There is currently no well-supported phylogenetic and biogeographical hypothesis of the Lythraceae incorporating all currently accepted genera, which we sought to provide.

METHODS

Plastomes of representative species of 18 distinct Lythraceae genera were sequenced and annotated. Together with existing sequences, plastomes of all 28 currently accepted genera in the Lythraceae were brought together for the first time. The plastomes were aligned and a Bayesian phylogenetic hypothesis was produced. We then conducted a time-calibrated Bayesian analysis and a biogeographical analysis.

KEY RESULTS

Plastome-based Bayesian and maximum-likelihood phylogenetic trees are generally congruent with recent nuclear phylogenomic data and resolve two deeply branching major clades in the Lythraceae. One major clade concentrates shrubby and arboreal South American and African genera that inhabit seasonally dry environments, with larger, often winged seeds, adapted to dispersal by the wind. The second major clade concentrates North American, Asian, African and several near-cosmopolitan herbaceous, shrubby and arboreal genera, often inhabiting humid or aquatic environments, with smaller seeds possessing structures that facilitate dispersal by water.

CONCLUSIONS

We hypothesize that the Lythraceae dispersed early in the Late Cretaceous from South American to North American continents, with subsequent expansion in the Late Cretaceous of a North American lineage through Laurasia to Africa via a boreotropical route. Two later expansions of South American clades to Africa in the Palaeocene and Eocene, respectively, are also hypothesized. Transoceanic dispersal in the family is possibly facilitated by adaptations to aquatic environments that are common to many extant genera of the Lythraceae, where long-distance dispersal and vicariance may be invoked to explain several remarkable disjunct distributions in Lythraceae clades.

A dataset of pomegranate growth stages for machine learning-based monitoring and analysis

Machine learning and deep learning have grown very rapidly in recent years and are widely used in agriculture. Neat and clean datasets are a major requirement for building accurate and robust machine learning models and minimizing misclassification in real-time environments. To achieve this goal, we created a dataset of images of pomegranate growth stages. These images of pomegranate growth stages were taken from May to September from an orchard inside the Henan Institute of Science and Technology in China. The dataset contains 5857 images of pomegranates at different growth stages, which are labeled and classified into five periods: bud, flower, early-fruit, mid-growth and ripe. The dataset consists of four folders, which respectively store the images, two formats of annotation files, and the record files for the division of training, validation, and test sets. The authors have confirmed the usability of this dataset through previous research. The dataset may help researchers develop computer applications using machine learning and computer vision algorithms.

High-quality genome of a pioneer mangrove Laguncularia racemosa explains its advantages for intertidal zone reforestation

Ecological restoration of mangrove ecosystems that became susceptible to recent habitat perturbations is crucial for tropical coast conservation. The white mangrove Laguncularia racemosa, a pioneer species inhabiting intertidal environments of the Atlantic East Pacific (AEP) region, has been used for reforestation in China for decades. However, the molecular mechanisms underlying its fast growth and high adaptive potential remain unknown. Using PacBio single-molecule real-time sequencing, we completed a high-quality L. racemosa genome assembly covering 1105 Mb with scaffold N50 of 3.46 Mb. Genomic phylogeny shows that L. racemosa invaded intertidal zones during a period of global warming. Multi-level genomic convergence analyses between L. racemosa and three native dominant mangrove clades show that they experienced convergent changes in genes involved in nutrient absorption and high salinity tolerance. This may explain successful L. racemosa adaptation to stressful intertidal environments after introduction. Without recent whole-genome duplications or activated transposable elements, L. racemosa has retained many tandem gene duplications. Some of them are involved in auxin biosynthesis, intense light stress and cold stress response pathways, associated with L. racemosa's ability to grow fast under high light or cold conditions when used for reforestation. In summary, our study identifies shared mechanisms of intertidal environmental adaptation and unique genetic changes underlying fast growth in mangrove-unfavourable conditions and sheds light on the molecular mechanisms of the white mangrove utility in ecological restoration.

Assembly and analysis of the first complete mitochondrial genome of Punica granatum and the gene transfer from chloroplast genome

Pomegranate (Punica granatum L.) is one of the oldest fruits with edible, medicinal and ornamental values. However, there is no report on the mitochondrial genome of pomegranate. In this study, the mitochondrial genome of P. granatum was sequenced, assembled and analyzed in detail, while the chloroplast genome was assembled using the same set of data. The results showed that the P. granatum mitogenome had a multi branched structure, using BGI + Nanopore mixed assembly strategy. The total genome length was 404,807 bp, with the GC content of 46.09%, and there were 37 protein coding genes, 20 tRNA genes and three rRNA genes. In the whole genome, 146 SSRs were identified. Besides, 400 pairs of dispersed repeats were detected, including 179 palindromic, 220 forward and one reverse. In the P. granatum mitochondrial genome, 14 homologous fragments of chloroplast genome were found, accounting for 0.54% of the total length. Phylogenetic analysis showed that among the published mitochondrial genomes of related genera, P. granatum had the closest genetic relationship with Lagerstroemia indica of Lythraceae. The 580 and 432 RNA editing sites were predicted on 37 protein coding genes of mitochondrial genome using BEDTools software and online website PREPACT respectively, but all were from C to U, of which ccmB and nad4 gene were most frequently edited, with 47 sites. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of pomegranate germplasm resources.

Historical biogeography and diversification of ringless Amanita (section Vaginatae) support an African origin and suggest niche conservatism in the Americas

Ectomycorrhizal fungi (ECM) sustain nutrient recycling in most terrestrial ecosystems, yet we know little about what major biogeographical events gave rise to present-day diversity and distribution patterns. Given the strict relationship between some ECM lineages and their hosts, geographically well-sampled phylogenies are central to understanding major evolutionary processes of fungal biodiversity patterns. Here, we focus on Amanita sect. Vaginatae to address global diversity and distribution patterns. Ancestral-state-reconstruction based on a 4-gene timetree with over 200 species supports an African origin between the late Paleocene and the early Eocene (ca. 56 Ma). Major biogeographic "out-of-Africa" events include multiple dispersal events to Southeast Asia (ca. 45-21 Ma), Madagascar (ca. 18 Ma), and the current Amazonian basin (ca. 45-36 Ma), the last two likely trans-oceanic. Later events originating in Southeast Asia involve Nearctic dispersal to North America (ca. 20-5 Ma), Oceania (Australia and New Zealand; ca. 15 Ma), and Europe (ca. 10-5 Ma). Subsequent dispersals were also inferred from Southeast Asia to East Asia (ca. 4 Ma); from North America to East Asia (ca. 11-8 Ma), Southeast Asia (ca. 19-2 Ma), Northern Andes (ca. 15 Ma), and Europe (ca. 15-2 Ma), respectively; and from the Amazon to the Caribbean region (ca. 25-20 Ma). Finally, we detected a significant increase in the net diversification rates in the branch leading to most northern temperate species in addition to higher state-dependent diversification rates in temperate lineages, consistent with previous findings. These results suggest that species of sect. Vaginatae likely have higher dispersal ability and higher adaptability to new environments, in particular compared to those of its sister clade, sect. Caesareae. Overall, the much wider distribution of A. sect. Vaginatae, from pan-tropical to pan-arctic, provides a unique window to understanding niche conservatism across a species-rich clade of ECM fungi.

Genomic survey of TCP transcription factors in plants: Phylogenomics, evolution and their biology

The TEOSINTE BRANCHED1 (TBI1), CYCLOIDEA (CYC), and PROLIFERATING CELL NUCLEAR ANTIGEN FACTORS (PCF1 and PCF2) proteins truncated as TCP transcription factors carry conserved basic-helix-loop-helix (bHLH) structure, related to DNA binding functions. Evolutionary history of the TCP genes has shown their presence in early land plants. In this paper, we performed a comparative discussion on the current knowledge of the TCP Transcription Factors in lower and higher plants: their evolutionary history based on the phylogenetics of 849 TCP proteins from 37 plant species, duplication events, and biochemical roles in some of the plants species. Phylogenetics investigations confirmed the classification of TCP TFs into Class I (the PCF1/2), and Class II (the C- clade) factors; the Class II factors were further divided into the CIN- and CYC/TB1- subclade. A trace in the evolution of the TCP Factors revealed an absence of the CYC/TB1subclade in lower plants, and an independent evolution of the CYC/TB1subclade in both eudicot and monocot species. 54% of the total duplication events analyzed were biased towards the dispersed duplication, and we concluded that dispersed duplication events contributed to the expansion of the TCP gene family. Analysis in the TCP factors functional roles confirmed their involvement in various biochemical processes which mainly included promoting cell proliferation in leaves in Class I TCPs, and cell division during plant development in Class II TCP Factors. Apart from growth and development, the TCP Factors were also shown to regulate hormonal and stress response pathways. Although this paper does not exhaust the present knowledge of the TCP Transcription Factors, it provides a base for further exploration of the gene family.

Reference quality genome sequence of Indian pomegranate cv. 'Bhagawa' (Punica granatum L.)

Pomegranate is an important fruit crop for ensuring livelihood and nutrition security in fragile semi-arid regions of the globe having limited irrigation resources. This is a high-value, nutritionally rich, and export-oriented agri-commodity that ensures high returns on investment to growers across the world. Although it is a valuable fruit crop, it has received only a limited genomics research outcome. To fast-track the pomegranate improvement program, de novo whole-genome sequencing of the main Indian cultivar 'Bhagawa' was initiated by the Indian Council of Agricultural Research-National Research Center on Pomegranate (ICAR-NRCP). We have demonstrated that a combination of commercially available technologies from Illumina, PacBio, 10X Genomics, and BioNano Genomics could be used efficiently for sequencing and reference-grade de novo assembly of the pomegranate genome. The research led to a final reference-quality genome assembly for 'Bhagawa' of 346.08 Mb in 342 scaffolds and an average N50 of 16.12 Mb and N90 of 1088.62 Kb. This assembly covered more than 98% of the estimated pomegranate genome size, 352.54 Mb. The LTR assembly index (LAI) value of 10 and 93.68% Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score over the 1,440 ortholog genes of the completed pomegranate genome indicates the quality of the assembled pomegranate genome. Furthermore, 29,435 gene models were discovered with a mean transcript length of 2,954 bp and a mean coding sequence length 1,090 bp. Four transcript data samples of pomegranate tissues were mapped over the assembled 'Bhagawa' genome up to 95% significant matches, indicating the high quality of the assembled genome. We have compared the 'Bhagawa' genome with the genomes of the pomegranate cultivars 'Dabenzi' and 'Taishanhong.' We have also performed whole-genome phylogenetic analysis using Computational Analysis of Gene Family Evolution (CAFE) and found that Eucalyptus grandis and pomegranate diverged 64 (60-70) million years ago. About 1,573 protein-coding resistance genes identified in the 'Bhagawa' genome were classified into 32 domains. In all, 314 copies of miRNA belonging to 26 different families were identified in the 'Bhagawa' genome. The reference-quality genome assembly of 'Bhagawa' is certainly a significant genomic resource for accelerated pomegranate improvement.

Chromosome-specific potential intron polymorphism markers for large-scale genotyping applications in pomegranate

Despite the availability of whole genome assemblies, the identification and utilization of gene-based marker systems has been limited in pomegranate. In the present study, we performed a genome-wide survey of intron length (IL) markers in the 36,524 annotated genes of the Tunisia genome. We identified and designed a total of 8,812 potential intron polymorphism (PIP) markers specific to 3,445 (13.40%) gene models that span 8 Tunisia chromosomes. The ePCR validation of all these PIP markers on the Tunisia genome revealed single-locus amplification for 1,233 (14%) markers corresponding to 958 (27.80%) genes. The markers yielding single amplicons were then mapped onto Tunisia chromosomes to develop a saturated linkage map. The functional categorization of 958 genes revealed them to be a part of the nucleus and the cytoplasm having protein binding and catalytic activity, and these genes are mainly involved in the metabolic process, including photosynthesis. Further, through ePCR, 1,233 PIP markers were assayed on multiple genomes, which resulted in the identification of 886 polymorphic markers with an average PIC value of 0.62. In silico comparative mapping based on physically mapped PIP markers indicates a higher synteny of Tunisia with the Dabenzi and Taishanhong genomes (>98%) in comparison with the AG2017 genome (95%). We then performed experimental validation of a subset of 100 PIP primers on eight pomegranate genotypes and identified 76 polymorphic markers, with 15 having PIC values ≥0.50. We demonstrated the potential utility of the developed markers by analyzing the genetic diversity of 31 pomegranate genotypes using 24 PIP markers. This study reports for the first time large-scale development of gene-based and chromosome-specific PIP markers, which would serve as a rich marker resource for genetic variation studies, functional gene discovery, and genomics-assisted breeding of pomegranate.

The clove (Syzygium aromaticum) genome provides insights into the eugenol biosynthesis pathway

The clove (Syzygium aromaticum) is an important tropical spice crop in global trade. Evolving environmental pressures necessitate modern characterization and selection techniques that are currently inaccessible to clove growers owing to the scarcity of genomic and genetic information. Here, we present a 370-Mb high-quality chromosome-scale genome assembly for clove. Comparative genomic analysis between S. aromaticum and Eucalyptus grandis-both species of the Myrtaceae family-reveals good genome structure conservation and intrachromosomal rearrangements on seven of the eleven chromosomes. We report genes that belong to families involved in the biosynthesis of eugenol, the major bioactive component of clove products. On the basis of our transcriptomic and metabolomic findings, we propose a hypothetical scenario in which eugenol acetate plays a key role in high eugenol accumulation in clove leaves and buds. The clove genome is a new contribution to omics resources for the Myrtaceae family and an important tool for clove research.

Genome sequencing and transcriptome analyses provide insights into the origin and domestication of water caltrop (Trapa spp., Lythraceae)

Humans have domesticated diverse species from across the plant kingdom; however, our current understanding of plant domestication is largely founded on major cereal crops. Here, we examine the evolutionary processes and genetic basis underlying the domestication of water caltrop (Trapa spp., Lythraceae), a traditional, yet presently underutilized non-cereal crop that sustained early Chinese agriculturalists. We generated a chromosome-level genome assembly of tetraploid T. natans, and then divided the allotetraploid genome into two subgenomes. Based on resequencing data from 57 accessions, representing cultivated diploid T. natans, wild T. natans (2x and 4x) and diploid T. incisa, we showed that water caltrop was likely first domesticated in the Yangtze River Valley as early as 6300 yr BP, and experienced a second improvement c. 800 years ago. We also provided strong support for an allotetraploid origin of T. natans within the past 230 000-310 000 years. By integrating selective sweep and transcriptome profiling analyses, we identified a number of genes potentially selected and/or differentially expressed during domestication, some of which likely contributed not only to larger fruit sizes but also to a more vigorous root system, facilitating nutrient uptake, environmental stress response and underwater photosynthesis. Our results shed light on the evolutionary and domestication history of water caltrop, one of the earliest domesticated crops in China. This study has implications for genomic-assisted breeding of this presently underutilized aquatic plant, and improves our general understanding of plant domestication.

Phylogenetic congruence between Neotropical primates and plants is driven by frugivory

Seed dispersal benefits plants and frugivores, and potentially drives co-evolution, with consequences to diversification evidenced for, e.g., primates. Evidence for macro-coevolutionary patterns in multi-specific, plant-animal mutualisms is scarce, and the mechanisms driving them remain unexplored. We tested for phylogenetic congruences in primate-plant interactions and showed strong co-phylogenetic signals across Neotropical forests, suggesting that both primates and plants share evolutionary history. Phylogenetic congruence between Platyrrhini and Angiosperms was driven by the most generalist primates, modulated by their functional traits, interacting with a wide-range of Angiosperms. Consistently similar eco-evolutionary dynamics seem to be operating irrespective of local assemblages, since co-phylogenetic signal emerged independently across three Neotropical regions. Our analysis supports the idea that macroevolutionary, coevolved patterns among interacting mutualistic partners are driven by super-generalist taxa. Trait convergence among multiple partners within multi-specific assemblages appears as a mechanism favouring these likely coevolved outcomes.

Food Plant Shifts Drive the Diversification of Sack-Bearer Moths

AbstractLepidoptera are a highly diverse group of herbivorous insects; however, some superfamilies have relatively few species. Two alternative hypotheses for drivers of Lepidoptera diversity are shifts in food plant use or shifts from concealed to external feeding as larvae. Many studies address the former hypothesis but with bias toward externally feeding taxa. One of the most striking examples of species disparity between sister lineages in Lepidoptera is between the concealed-feeding sack-bearer moths (Mimallonoidea), which contain about 300 species, and externally feeding Macroheterocera, which have over 74,000 species. We provide the first dated tree of Mimallonidae to understand the diversification dynamics of these moths in order to fill a knowledge gap pertaining to drivers of diversity within an important concealed-feeding clade. We find that Mimallonidae is an ancient Lepidoptera lineage that originated in the Cretaceous ∼105 million years ago and has had a close association with the plant order Myrtales for the past 40 million years. Diversification dynamics are tightly linked with food plant usage in this group. Reliance on Myrtales may have influenced diversification of Mimallonidae because clades that shifted away from the ancestral condition of feeding on Myrtales have the highest speciation rates in the family.

Ultraconserved elements-based systematics reveals evolutionary patterns of host-plant family shifts and phytophagy within the predominantly parasitoid braconid wasp subfamily Doryctinae

Phytophagy has promoted species diversification in many insect groups, including Hymenoptera, one of the most diverse animal orders on Earth. In the predominantly parasitoid family Braconidae, an association with insect-induced, plant galls in angiosperms have been reported in three subfamilies, but in particular in the Doryctinae, where it has been recorded to occur in species of ten genera. Allorhogas Gahan is the most species-rich of these genera, with its species having different phytophagous strategies. Here we conducted a comprehensive phylogenomic study for the doryctine gall-associated genera, with an emphasis on Allorhogas, using ultraconserved elements (UCEs). Based on this estimate of phylogeny we: (1) evaluated their taxonomic composition, (2) estimated the timing of origin of the gall-associated clade and divergence of its main subclades, and (3) performed ancestral state reconstruction analyses for life history traits related to their host-plant association. Our phylogenetic hypothesis confirmed Allorhogas as polyphyletic, with most of its members being nested in a main clade composed of various subclades, each comprising species with a particular host-plant family and herbivorous feeding habit. The origin of gall-association was estimated to have occurred during the late Oligocene to early Miocene, with a subsequent diversification of subclades during the middle to late Miocene and Pliocene. Overlap in divergence timing appears to occur between some taxa and their host-associated plant lineages. Evolution of the feeding strategies in the group shows "inquilinism-feeding" as the likely ancestral state, with gall-formation in different plant organs and seed predation having independently evolved on multiple occasions.

Two centuries of distribution data: detection of areas of endemism for the Brazilian angiosperms

Brazil has high levels of biodiversity and has received strong criticism for the increasing country-wide deforestation that threatens it. Although a significant percentage of land area in Brazil is protected, the areas are insufficient and unevenly distributed. Many studies have contributed to the biogeographical knowledge of Brazilian flora, but no endemicity analysis (EA) has been conducted including all endemic angiosperms. We investigated the spatial component, drawing on a huge and taxonomically diverse dataset based on 827 016 records collected over the last two centuries. We conducted an EA for 15 034 species from 173 families using an optimality criterion with 2° and 3° grid sizes, in order to search for distributional concordance, to identify the biogeographical units and discuss the implications for conservation. Six analyses were run for basal angiosperms, monocots and eudicots. The EA recovered 66 consensus areas (CAs). The concordance of CAs enabled the identification of five best-supported areas of endemism--three in the Atlantic and Parana Forest and two areas in the Cerrado province--supported by species of 120 families. The age of divergence for some genera that contributed significantly to the identification of areas recovered in the Cerrado coincides with the recent, <10 Ma, estimated age of that province. By contrast, the areas in the Atlantic and Parana Forest are supported by genera with earlier diversification >30 Ma, supporting an ancient origin. Most areas in the Atlantic Forest are partially superimposed with the limits of the protected areas, however, big gaps were identified in the Cerrado. Protecting Brazilian biomes was at the heart of Brazil's environmental policy. Regrettably, this scenario has radically changed based on misleading divergences in conservation policy. Areas of endemism are pivotal for biodiversity conservation due to the common evolutionary history shared by their endemic taxa. Thus, we hope that these congruent patterns of endemism support the establishment of biodiversity priorities.

Chloroplast Genome Sequencing, Comparative Analysis, and Discovery of Unique Cytoplasmic Variants in Pomegranate (Punica granatum L.)

Here we report on comprehensive chloroplast (cp) genome analysis of 16 pomegranate (Punica granatum L.) genotypes representing commercial cultivars, ornamental and wild types, through large-scale sequencing and assembling using next-generation sequencing (NGS) technology. Comparative genome analysis revealed that the size of cp genomes varied from 158,593 bp (in wild, “1201” and “1181”) to 158,662 bp (cultivar, “Gul-e-Shah Red”) among the genotypes, with characteristic quadripartite structures separated by a pair of inverted repeats (IRs). The higher conservation for the total number of coding and non-coding genes (rRNA and tRNA) and their sizes, and IRs (IR-A and IR-B) were observed across all the cp genomes. Interestingly, high variations were observed in sizes of large single copy (LSC, 88,976 to 89,044 bp) and small single copy (SSC, 18,682 to 18,684 bp) regions. Although, the structural organization of newly assembled cp genomes were comparable to that of previously reported cp genomes of pomegranate (“Helow,” “Tunisia,” and “Bhagawa”), the striking differences were observed with the Lagerstroemia lines, viz., Lagerstroemia intermedia (NC_0346620) and Lagerstroemia speciosa (NC_031414), which clearly confirmed previous findings. Furthermore, phylogenetic analysis also revealed that members outside the genus Punica were clubbed into a separate clade. The contraction and expansion analysis revealed that the structural variations in IRs, LSC, and SSC have significantly accounted for the evolution of cp genomes of Punica and L. intermedia over the periods. Microsatellite survey across cp genomes resulted in the identification of a total of 233 to 234 SSRs, with majority of them being mono- (A/T or C/G, 164–165 numbers), followed by di- (AT/AT or AG/CT, 54), tri- (6), tetra- (8), and pentanucleotides (1). Furthermore, the comparative structural variant analyses across cp genomes resulted in the identification of many varietal specific SNP/indel markers. In summary, our study has offered a successful development of large-scale cp genomics resources to leverage future genetic, taxonomical, and phylogenetic studies in pomegranate.

Low bee visitation rates explain pollinator shifts to vertebrates in tropical mountains

Evolutionary shifts from bee to vertebrate pollination are common in tropical mountains. Reduction in bee pollination efficiency under adverse montane weather conditions was proposed to drive these shifts. Although pollinator shifts are central to the evolution and diversification of angiosperms, we lack experimental evidence of the ecological processes underlying such shifts. Here, we combine phylogenetic and distributional data for 138 species of the Neotropical plant tribe Merianieae (Melastomataceae) with pollinator observations of 11 and field pollination experiments of six species to test whether the mountain environment may indeed drive such shifts. We demonstrate that shifts from bee to vertebrate pollination coincided with occurrence at high elevations. We show that vertebrates were highly efficient pollinators even under the harsh environmental conditions of tropical mountains, whereas bee pollination efficiency was lowered significantly through reductions in flower visitation rates. Furthermore, we show that pollinator shifts in Merianieae coincided with the final phases of the Andean uplift and were contingent on adaptive floral trait changes to alternative rewards and mechanisms facilitating pollen dispersal. Our results provide evidence that abiotic environmental conditions (i.e. mountain climate) may indeed reduce the efficiency of a plant clade's ancestral pollinator group and correlate with shifts to more efficient new pollinators.

A nuclear phylogenomic study of the angiosperm order Myrtales, exploring the potential and limitations of the universal Angiosperms353 probe set

PREMISE

To further advance the understanding of the species-rich, economically and ecologically important angiosperm order Myrtales in the rosid clade, comprising nine families, approximately 400 genera and almost 14,000 species occurring on all continents (except Antarctica), we tested the Angiosperms353 probe kit.

METHODS

We combined high-throughput sequencing and target enrichment with the Angiosperms353 probe kit to evaluate a sample of 485 species across 305 genera (76% of all genera in the order).

RESULTS

Results provide the most comprehensive phylogenetic hypothesis for the order to date. Relationships at all ranks, such as the relationship of the early-diverging families, often reflect previous studies, but gene conflict is evident, and relationships previously found to be uncertain often remain so. Technical considerations for processing HTS data are also discussed.

CONCLUSIONS

High-throughput sequencing and the Angiosperms353 probe kit are powerful tools for phylogenomic analysis, but better understanding of the genetic data available is required to identify genes and gene trees that account for likely incomplete lineage sorting and/or hybridization events.

Contrasting Phylogeographic Patterns in Lumnitzera Mangroves Across the Indo-West Pacific

Mangroves are ecologically important forest communities in tropical and subtropical coasts, the effective management of which requires understanding of their phylogeographic patterns. However, these patterns often vary among different species, even among ecologically similar taxa or congeneric species. Here, we investigated the levels and patterns of genetic variation within Lumnitzera consisting of two species (L. racemosa and L. littorea) with nearly sympatric ranges across the Indo-West Pacific (IWP) region by sequencing three chloroplast DNA regions (for both species) and genotyping 11 nuclear microsatellite loci (for L. littorea). Consistent with findings in studies on other mangrove species, we found that both L. racemosa and L. littorea showed relatively high genetic variation among populations but low genetic variation within populations. Haplotype network and genetic clustering analyses indicated two well-differentiated clades in both L. racemosa and L. littorea. The relationship between geographic and genetic distances and divergence time estimates of the haplotypes indicated that limited dispersal ability of the propagules, emergence of land barriers during ancient sea-level changes, and contemporary oceanic circulation pattern in the IWP influenced the current population structure of the two species. However, the position of genetic break was found to vary between the two species: in L. racemosa, strong divergence was observed between populations from the Indian Ocean and the Pacific Ocean possibly due to land barrier effect of the Malay Peninsula; in L. littorea, the phylogeographic pattern was created by a more eastward genetic break along the biogeographic barrier identified as the Huxley's line. Overall, our findings strongly supported previous hypothesis of mangrove species divergence and revealed that the two Lumnitzera species have different phylogeographic patterns despite their close genetic relationship and similar current geographic distribution. The findings also provided references for the management of Lumnitzera mangroves, especially for the threatened L. littorea.

Chloroplast phylogenomics and divergence times of Lagerstroemia (Lythraceae)

Background

Crape myrtles, belonging to the genus Lagerstroemia L., have beautiful paniculate inflorescences and are cultivated as important ornamental tree species for landscaping and gardening. However, the phylogenetic relationships within Lagerstroemia have remained unresolved likely caused by limited sampling and the insufficient number of informative sites used in previous studies.

Results

In this study, we sequenced 20 Lagerstroemia chloroplast genomes and combined with 15 existing chloroplast genomes from the genus to investigate the phylogenetic relationships and divergence times within Lagerstroemia. The phylogenetic results indicated that this genus is a monophyletic group containing four clades. Our dating analysis suggested that Lagerstroemia originated in the late Paleocene (~ 60 Ma) and started to diversify in the middle Miocene. The diversification of most species occurred during the Pleistocene. Four variable loci, trnD-trnY-trnE, rrn16-trnI, ndhF-rpl32-trnL and ycf1, were discovered in the Lagerstroemia chloroplast genomes.

Conclusions

The chloroplast genome information was successfully utilized for molecular characterization of diverse crape myrtle samples. Our results are valuable for the global genetic diversity assessment, conservation and utilization of Lagerstroemia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07769-x.

Comparative analysis of chloroplast genome structure and molecular dating in Myrtales

BACKGROUND

Myrtales is a species rich branch of Rosidae, with many species having important economic, medicinal, and ornamental value. At present, although there are reports on the chloroplast structure of Myrtales, a comprehensive analysis of the chloroplast structure of Myrtales is lacking. Phylogenetic and divergence time estimates of Myrtales are mostly constructed by using chloroplast gene fragments, and the support for relationships is low. A more reliable method to reconstruct the species divergence time and phylogenetic relationships is by using whole chloroplast genomes. In this study, we comprehensively analyzed the structural characteristics of Myrtales chloroplasts, compared variation hotspots, and reconstructed the species differentiation time of Myrtales with four fossils and one secondary calibration point.

RESULTS

A total of 92 chloroplast sequences of Myrtales, representing six families, 16 subfamilies and 78 genera, were obtained including nine newly sequenced chloroplasts by whole genome sequencing. Structural analyses showed that the chloroplasts range in size between 152,214-171,315 bp and exhibit a typical four part structure. The IR region is between 23,901-36,747 bp, with the large single copy region spanning 83,691-91,249 bp and the small single copy region spanning 11,150-19,703 bp. In total, 123-133 genes are present in the chloroplasts including 77-81 protein coding genes, four rRNA genes and 30-31 tRNA genes. The GC content was 36.9-38.9%, with the average GC content being 37%. The GC content in the LSC, SSC and IR regions was 34.7-37.3%, 30.6-36.8% and 39.7-43.5%, respectively. By analyzing nucleotide polymorphism of the chloroplast, we propose 21 hypervariable regions as potential DNA barcode regions for Myrtales. Phylogenetic analyses showed that Myrtales and its corresponding families are monophyletic, with Combretaceae and the clade of Onagraceae + Lythraceae (BS = 100%, PP = 1) being sister groups. The results of molecular dating showed that the crown of Myrtales was most likely to be 104.90 Ma (95% HPD = 87.88-114.18 Ma), and differentiated from the Geraniales around 111.59 Ma (95% HPD = 95.50-118.62 Ma).

CONCLUSIONS

The chloroplast genome structure of Myrtales is similar to other angiosperms and has a typical four part structure. Due to the expansion and contraction of the IR region, the chloroplast genome sizes in this group are slightly different. The variation of noncoding regions of the chloroplast genome is larger than those of coding regions. Phylogenetic analysis showed that Combretaceae and Onagraceae + Lythraceae were well supported as sister groups. Molecular dating indicates that the Myrtales crown most likely originated during the Albian age of the Lower Cretaceous. These chloroplast genomes contribute to the study of genetic diversity and species evolution of Myrtales, while providing useful information for taxonomic and phylogenetic studies of Myrtales.

Phylogenomics and Historical Biogeography of Seahorses, Dragonets, Goatfishes, and Allies (Teleostei: Syngnatharia): Assessing Factors Driving Uncertainty in Biogeographic Inferences

The charismatic trumpetfishes, goatfishes, dragonets, flying gurnards, seahorses, and pipefishes encompass a recently defined yet extraordinarily diverse clade of percomorph fishes-the series Syngnatharia. This group is widely distributed in tropical and warm-temperate regions, with a great proportion of its extant diversity occurring in the Indo-Pacific. Because most syngnatharians feature long-range dispersal capabilities, tracing their biogeographic origins is challenging. Here, we applied an integrative phylogenomic approach to elucidate the evolutionary biogeography of syngnatharians. We built upon a recently published phylogenomic study that examined ultraconserved elements by adding 62 species (total 169 species) and one family (Draconettidae), to cover ca. 25% of the species diversity and all 10 families in the group. We inferred a set of time-calibrated trees and conducted ancestral range estimations. We also examined the sensitivity of these analyses to phylogenetic uncertainty (estimated from multiple genomic subsets), area delimitation, and biogeographic models that include or exclude the jump-dispersal parameter (j). Of the three factors examined, we found that the j parameter has the strongest effect in ancestral range estimates, followed by number of areas defined, and tree topology and divergence times. After accounting for these uncertainties, our results reveal that syngnatharians originated in the ancient Tethys Sea ca. 87 Ma (84-94 Ma; Late Cretaceous) and subsequently occupied the Indo-Pacific. Throughout syngnatharian history, multiple independent lineages colonized the eastern Pacific (6-8 times) and the Atlantic (6-14 times) from their center of origin, with most events taking place following an east-to-west route prior to the closure of the Tethys Seaway ca. 12-18 Ma. Ultimately, our study highlights the importance of accounting for different factors generating uncertainty in macroevolutionary and biogeographic inferences.

Austropuccinia psidii, causing myrtle rust, has a gigabase-sized genome shaped by transposable elements

Austropuccinia psidii, originating in South America, is a globally invasive fungal plant pathogen that causes rust disease on Myrtaceae. Several biotypes are recognized, with the most widely distributed pandemic biotype spreading throughout the Asia-Pacific and Oceania regions over the last decade. Austropuccinia psidii has a broad host range with more than 480 myrtaceous species. Since first detected in Australia in 2010, the pathogen has caused the near extinction of at least three species and negatively affected commercial production of several Myrtaceae. To enable molecular and evolutionary studies into A. psidii pathogenicity, we assembled a highly contiguous genome for the pandemic biotype. With an estimated haploid genome size of just over 1 Gb (gigabases), it is the largest assembled fungal genome to date. The genome has undergone massive expansion via distinct transposable element (TE) bursts. Over 90% of the genome is covered by TEs predominantly belonging to the Gypsy superfamily. These TE bursts have likely been followed by deamination events of methylated cytosines to silence the repetitive elements. This in turn led to the depletion of CpG sites in TEs and a very low overall GC content of 33.8%. Compared to other Pucciniales, the intergenic distances are increased by an order of magnitude indicating a general insertion of TEs between genes. Overall, we show how TEs shaped the genome evolution of A. psidii and provide a greatly needed resource for strategic approaches to combat disease spread.

Evolution and biogeography of Memecylon

PREMISE

The woody plant group Memecylon (Melastomataceae) is a large clade occupying diverse forest habitats in the Old World tropics and exhibiting high regional endemism. Its phylogenetic relationships have been previously studied using ribosomal DNA with extensive sampling from Africa and Madagascar. However, divergence times, biogeography, and character evolution of Memecylon remain uninvestigated. We present a phylogenomic analysis of Memecylon to provide a broad evolutionary perspective of this clade.

METHODS

One hundred supercontigs of 67 Memecylon taxa were harvested from target enrichment. The data were subjected to coalescent and concatenated phylogenetic analyses. A timeline was provided for Memecylon evolution using fossils and secondary calibration. The calibrated Memecylon phylogeny was used to elucidate its biogeography and ancestral character states.

RESULTS

Relationships recovered by the phylogenomic analyses are strongly supported in both maximum likelihood and coalescent-based species trees. Memecylon is inferred to have originated in Africa in the Eocene and subsequently dispersed predominantly eastward via long-distance dispersal (LDD), although a reverse dispersal from South Asia westward to the Seychelles was postulated. Morphological data exhibited high levels of homoplasy, but also showed that several vegetative and reproductive characters were phylogenetically informative.

CONCLUSIONS

The current distribution of Memecylon appears to be the result of multiple ancestral LDD events. Our results demonstrate the importance of the combined effect of geographic and paleoclimatic factors in shaping the distribution of this group in the Old World tropics. Memecylon includes a number of evolutionarily derived morphological features that contribute to diversity within the clade.

Uncertainty in Divergence Time Estimation

Understanding and representing uncertainty is crucial in academic research because it enables studies to build on the conclusions of previous studies, leading to robust advances in a particular field. Here, we evaluate the nature of uncertainty and the manner by which it is represented in divergence time estimation, a field that is fundamental to many aspects of macroevolutionary research, and where there is evidence that uncertainty has been seriously underestimated. We address this issue in the context of methods used in divergence time estimation, and with respect to the manner by which time-calibrated phylogenies are interpreted. With respect to methods, we discuss how the assumptions underlying different methods may not adequately reflect uncertainty about molecular evolution, the fossil record, or diversification rates. Therefore, divergence time estimates may not adequately reflect uncertainty and may be directly contradicted by subsequent findings. For the interpretation of time-calibrated phylogenies, we discuss how the use of time-calibrated phylogenies for reconstructing general evolutionary timescales leads to inferences about macroevolution that are highly sensitive to methodological limitations in how uncertainty is accounted for. By contrast, we discuss how the use of time-calibrated phylogenies to test specific hypotheses leads to inferences about macroevolution that are less sensitive to methodological limitations. Given that many biologists wish to use time-calibrated phylogenies to reconstruct general evolutionary timescales, we conclude that the development of methods of divergence time estimation that adequately account for uncertainty is necessary. [Divergence time estimation; macroevolution; uncertainty.].

Chloroplast Genome Sequences and Comparative Analyses of Combretaceae Mangroves with Related Species

In the Combretaceae family, only two species of Lumnitzera and one species of Laguncularia belong to mangroves. Among them, Lumnitzera littorea (Jack) Voigt. is an endangered mangrove plant in China for the limited occurrence and seed abortion. In contrast, Lumnitzera racemosa Willd. is known as the most widespread mangrove plant in China. Laguncularia racemosa C. F. Gaertn., an exotic mangrove in China, has the fast growth and high adaptation ability. To better understand the phylogenetic positions of these mangroves in Combretaceae and in Myrtales and to provide information for studies on evolutionary adaptation for intertidal habitat, the complete chloroplast (cp) genomes of Lu. racemosa and La. racemosa were sequenced. Furthermore, we present here the results from the assembly and annotation of the two cp genomes, which were further subjected to the comparative analysis with Lu. littorea cp genomes we published before and other eleven closely related species within Myrtales. The chloroplast genomes of the three Combretaceae mangrove species: Lu. littorea, Lu. racemosa, and La. racemosa are 159,687 bp, 159,473 bp, and 158,311 bp in size. All three cp genomes host 130 genes including 85 protein-coding genes, 37 tRNAs, and 4 rRNAs. A comparative analysis of those three genomes revealed the high similarity of genes in coding-regions and conserved gene order in the IR and LSC/SSC regions. The differences between Lumnitzera and Laguncularia cp genomes are the locations of rps19 and rpl2 genes in the IR/SC boundary regions. Investigating the effects of selection events on shared protein-coding genes showed a relaxed selection had acted on the ycf2, ycf1, and matK genes of Combretaceae mangroves compared to the nonmangrove species Eucalyptus aromaphloia. The phylogenetic analysis based on the whole chloroplast genome sequence with one outgroup species strongly supported three Combretaceae mangroves together with other two Combretaceae species formed a cluster in Combretaceae. This study is the first report on the comparative analysis of three Combretaceae mangrove chloroplast genomes, which will provide the significant information for understanding photosynthesis and evolution in Combretaceae mangrove plants.

The complete chloroplast genome sequence of Trapa bicornis Osbeck (Lythraceae)

Trapa (Lythraceae) is an aquatic plant genus widely distributed in the old world. Although Trapa species have great edible and medical value, studies related to species identification and utilization are still lacking. Here, we reported the complete chloroplast genome sequence of a cultivated species, T. bicornis. The chloroplast genome size of T. bicornis was 155,539 bp, consisting of a pair of inverted repeat (IR) regions (24,386 bp), separated by a large single copy (LSC) region (88,493 bp) and a small single copy (SSC) region (18,274 bp). A total of 130 genes were annotated, including 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The phylogenomic analysis supported the monophyly of Trapa, and a sister relationship between T. bicornis and T. natans.

Historical biogeography of Vochysiaceae reveals an unexpected perspective of plant evolution in the Neotropics

PREMISE

Despite the fast pace of exploration of the patterns and processes influencing Neotropical plant hyperdiversity, the taxa explored are mostly from large groups that are widely distributed, morphologically diverse, or economically important. Vochysiaceae is an example of an undersampled taxon, providing an excellent system for investigating Neotropical biogeography. We present a phylogenomics-based hypothesis of species relationships in Vochysiaceae to investigate its evolutionary history through space and time.

METHODS

We inferred a phylogeny for 122 species from Vochysiaceae and seven other families of Myrtales. Fossils from four myrtalean families were used to estimate the divergence times within Vochysiaceae. Historical biogeography was estimated using ancestral range probabilities and stochastic mapping.

RESULTS

Monophyly of all genera was supported except for Qualea, which was split by Ruizterania into two clades. Vochysiaceae originated ~100 mya, splitting into an Afrotropical and a Neotropical lineage ~50 mya, and its ancestral range is in the area currently occupied by the Cerrado.

CONCLUSIONS

The most recent common ancestor of Vochysiaceae + Myrtaceae had a West Gondwanan distribution, supporting a South American + African ancestral range of Vochysiaceae. On a global scale, geographic range reduction was the principal biogeographic event. At a finer scale, initial range reduction was also important and the Cerrado region was the most ancestral area with multiple colonization events to the Amazon, Central America, and the Atlantic Forest. Colonization events occurred from open areas to forest vegetation, an unusual finding regarding the evolution of plants in the Neotropics.

Complete Chloroplast Genomes of 14 Mangroves: Phylogenetic and Comparative Genomic Analyses

Mangroves are a group of plant species that occupy the coastal intertidal zone and are major components of this ecologically important ecosystem. Mangroves belong to about twenty diverse families. Here, we sequenced and assembled chloroplast genomes of 14 mangrove species from eight families spanning five rosid orders and one asterid order: Fabales (Pongamia pinnata), Lamiales (Avicennia marina), Malpighiales (Excoecaria agallocha, Bruguiera sexangula, Kandelia obovata, Rhizophora stylosa, and Ceriops tagal), Malvales (Hibiscus tiliaceus, Heritiera littoralis, and Thespesia populnea), Myrtales (Laguncularia racemosa, Sonneratia ovata, and Pemphis acidula), and Sapindales (Xylocarpus moluccensis). These chloroplast genomes range from 149 kb to 168 kb in length. A conserved structure of two inverted repeats (IRa and IRb, ~25.8 kb), one large single-copy region (LSC, ~89.0 kb), and one short single-copy region (SSC, ~18.9 kb) as well as ~130 genes (85 protein-coding, 37 tRNAs, and 8 rRNAs) was observed. We found the lowest divergence in the IR regions among the four regions. We also identified simple sequence repeats (SSRs), which were found to be variable in numbers. Most chloroplast genes are highly conserved, with only four genes under positive selection or relaxed pressure. Combined with publicly available chloroplast genomes, we carried out phylogenetic analysis and confirmed the previously reported phylogeny within rosids, including the positioning of obscure families in Malpighiales. Our study reports 14 mangrove chloroplast genomes and illustrates their genome features and evolution.

A two-tier bioinformatic pipeline to develop probes for target capture of nuclear loci with applications in Melastomataceae

PREMISE

Putatively single-copy nuclear (SCN) loci, which are identified using genomic resources of closely related species, are ideal for phylogenomic inference. However, suitable genomic resources are not available for many clades, including Melastomataceae. We introduce a versatile approach to identify SCN loci for clades with few genomic resources and use it to develop probes for target enrichment in the distantly related Memecylon and Tibouchina (Melastomataceae).

METHODS

We present a two-tiered pipeline. First, we identified putatively SCN loci using MarkerMiner and transcriptomes from distantly related species in Melastomataceae. Published loci and genes of functional significance were then added (384 total loci). Second, using HybPiper, we retrieved 689 homologous template sequences for these loci using genome-skimming data from within the focal clades.

RESULTS

We sequenced 193 loci common to Memecylon and Tibouchina. Probes designed from 56 template sequences successfully targeted sequences in both clades. Probes designed from genome-skimming data within a focal clade were more successful than probes designed from other sources.

DISCUSSION

Our pipeline successfully identified and targeted SCN loci in Memecylon and Tibouchina, enabling phylogenomic studies in both clades and potentially across Melastomataceae. This pipeline could be easily applied to other clades with few genomic resources.

Flower morphology is correlated with distribution and phylogeny in Bertolonia (Melastomataceae), an herbaceous genus endemic to the Atlantic Forest

Several centers of endemism have been proposed for Melastomataceae, particularly in Amazonia and Atlantic Forest. Despite the high degree of human-caused degradation in the last 500 years, the Atlantic Forest still presents some of the largest levels of diversity and endemism across all angiosperms. With several recently described species in the last decade, the knowledge on Bertolonia's distribution and morphological characterization has changed, with most new species found in northern Atlantic Forest and with different flower color patterns than the species from southern Atlantic Forest. We first tested the monophyly of the genus sampling over 85% of its species to generate a reliable phylogenetic hypothesis. Afterwards, we used Bertolonia as a model group to study distribution patterns and morphological evolution of lineages in the Atlantic Forest. Bertolonia is particularly interesting to address such questions because it is endemic to this domain, with species distributed either in the southern, central or northern portions of the Atlantic Forest. The second step of our work aimed to respond (1) Do closely related species endemic to Atlantic Forest tend to have similar distributions and share similar morphological traits? and (2) Are floral traits more conserved within clades than vegetative characters? We hypothesize that both (1) and (2) are true due to the low dispersal ability and consequent microendemic distribution of most species in the genus. Our results confirm the veracity of the proposed hypotheses. Three major groups were recovered in our analysis: marmorata, formosa and nymphaeifolia clades. Most species that occur in northern Atlantic Forest were recovered in the marmorata clade. Most species distributed in central Atlantic Forest were recovered in the formosa clade, and the ones from the southern Atlantic Forest in the nymphaeifolia clade. A similar pattern was recovered with reproductive characters reconstructed across the phylogeny. Generally, species with pink petals and cream-colored anthers with an extrorse pore occur in northern Atlantic Forest, and species with white petals and yellow anthers with an introrse or apical pore tend to occur in southern Atlantic Forest. Some vegetative characters also have evolutionary congruence and are restricted to one or few lineages, while other characters, such as the texture of the leaf blade surface, have a strong taxonomic value and are useful to identify species, but are not homologous. Our analysis also indicates that the division between southern and northern Atlantic Forest could also be related to ancient events, not only linked with recent phylogeographic patterns. Moreover, we suggest that the orientation of the anther pore in Bertolonia could be related with diversity of species and stability of its populations. In summary, we corroborate, based on the evolutionary history of Bertolonia, that closely related species endemic to Atlantic Forest tend to have similar distributions and share similar morphological floral traits.

Is dispersal mode a driver of diversification and geographical distribution in the tropical plant family Melastomataceae?

Species of plants with different life history strategies may differ in their seed dispersal mechanisms, impacting their distribution and diversification patterns. Shorter or longer distance dispersal is favored by different dispersal modes, facilitating (or constraining) population isolation, which can, in turn, impact speciation and species range sizes. While these associations are intuitive, few studies have explicitly tested these hypotheses for large clades of angiosperms. The plant family Melastomataceae is found on disparate habitats with different dispersal modes, representing a good model to address these questions. In this study, we reconstruct the phylogeny of Melastomataceae and gather data on their dispersal mode and range size to test the impact of dispersal mode on diversification and range size evolution. We found that abiotic dispersal is ancestral in the family, while biotic dispersal evolved multiple times. Species richness distribution is very similar across dispersal modes, although abiotically dispersed species tend to be relatively more diverse in seasonal environments. Range sizes across dispersal modes are not significantly different, although biotically dispersed species have slightly wider distributions. Model comparisons indicate that factors other than dispersal mode might have driven diversification heterogeneity. We did not find evidence for the role of dispersal mode driving diversification rates or range size in the Melastomataceae, suggesting a complex macroevolutionary scenario for this diverse angiosperm family. The bulk of changes to biotic dispersal coinciding with an increase in passerine diversification suggests a possible "past" key innovation in Melastomataceae. Future studies should investigate the role of other diversification drivers in the family and the relatively higher diversity of abiotically dispersed species in open habitats.

Modularity increases rate of floral evolution and adaptive success for functionally specialized pollination systems

Angiosperm flowers have diversified in adaptation to pollinators, but are also shaped by developmental and genetic histories. The relative importance of these factors in structuring floral diversity remains unknown. We assess the effects of development, function and evolutionary history by testing competing hypotheses on floral modularity and shape evolution in Merianieae (Melastomataceae). Merianieae are characterized by different pollinator selection regimes and a developmental constraint: tubular anthers adapted to specialized buzz-pollination. Our analyses of tomography-based 3-dimensional flower models show that pollinators selected for functional modules across developmental units and that patterns of floral modularity changed during pollinator shifts. Further, we show that modularity was crucial for Merianieae to overcome the constraint of their tubular anthers through increased rates of evolution in other flower parts. We conclude that modularity may be key to the adaptive success of functionally specialized pollination systems by making flowers flexible (evolvable) for adaptation to changing selection regimes.

Dellinger et al. report the characterization of different pollinator selection regimes on 3D flower models. They show that pollinators selected for functional floral modules, and that this modularity allowed certain species to overcome morphological constraints by increased evolutionary rates in other flower parts.

Analyses of Plastome Sequences Improve Phylogenetic Resolution and Provide New Insight Into the Evolutionary History of Asian Sonerileae/Dissochaeteae

Sonerileae/Dissochaeteae (Melastomataceae) comprises ca. 50 genera, two thirds of which occur in Southeast Asia. Phylogenetic relationships within this clade remain largely unclear, which hampers our understanding of its origin, evolution, and biogeography. Here, we explored the use of chloroplast genomes in phylogenetic reconstruction of Sonerileae/Dissochaeteae, by sampling 138 species and 23 genera in this clade. A total of 151 complete plastid genomes were assembled for this study. Plastid genomic data provided better support for the backbone of the Sonerileae/Dissochaeteae phylogeny, and also for relationships among most closely related species, but failed to resolve the short internodes likely resulted from rapid radiation. Trees inferred from plastid genome and nrITS sequences were largely congruent regarding the major lineages of Sonerileae/Dissochaeteae. The present analyses recovered 15 major lineages well recognized in both nrITS and plastid phylogeny. Molecular dating and biogeographical analyses indicated a South American origin for Sonerileae/Dissochaeteae during late Eocene (stem age: 34.78 Mya). Two dispersal events from South America to the Old World were detected in late Eocene (33.96 Mya) and Mid Oligocene (28.33 Mya) respectively. The core Asian clade began to diversify around early Miocene in Indo-Burma and dispersed subsequently to Malesia and Sino-Japanese regions, possibly promoted by global temperature changes and East Asian monsoon activity. Our analyses supported previous hypothesis that Medinilla reached Madagascar by transoceanic dispersal in Miocene. In addition, generic limits of some genera concerned were discussed.

Incongruence between gene trees and species trees and phylogenetic signal variation in plastid genes

The current classification of angiosperms is based primarily on concatenated plastid markers and maximum likelihood (ML) inference. This approach has been justified by the assumption that plastid DNA (ptDNA) is inherited as a single locus and that its individual genes produce congruent trees. However, structural and functional characteristics of ptDNA suggest that plastid genes may not evolve as a single locus and are experiencing different evolutionary forces. To examine this idea, we produced new complete plastid genome (plastome) sequences of 27 species and combined these data with publicly available sequences to produce a final dataset that includes 78 plastid genes for 89 species of rosids and five outgroups. We used four data matrices (i.e., gene, exon, codon-aligned, and amino acid) to infer species and gene trees using ML and multispecies coalescent (MSC) methods. Rosids include about one third of all angiosperms and their two major clades, fabids and malvids, were recovered in almost all analyses. However, we detected incongruence between species trees inferred with different matrices and methods and previously published plastid and nuclear phylogenies. We visualized and tested the significance of incongruence between gene trees and species trees. We then measured the distribution of phylogenetic signal across sites and genes supporting alternative placements of five controversial nodes at different taxonomic levels. Gene trees inferred with plastid data often disagree with species trees inferred using both ML (with unpartitioned or partitioned data) and MSC. Species trees inferred with both methods produced alternative topologies for a few taxa. Our results show that, in a phylogenetic context, plastid protein-coding genes may not be fully linked and behaving as a single locus. Furthermore, concatenated matrices may produce highly supported phylogenies that are discordant with individual gene trees. We also show that phylogenies inferred with MSC are accurate. We therefore emphasize the importance of considering variation in phylogenetic signal across plastid genes and the exploration of plastome data to increase accuracy of estimating relationships. We also support the use of MSC with plastome matrices in future phylogenomic investigations.

Tracking temporal shifts in area, biomes, and pollinators in the radiation of Salvia (sages) across continents: leveraging anchored hybrid enrichment and targeted sequence data

PREMISE OF THE STUDY

A key question in evolutionary biology is why some clades are more successful by being widespread geographically, biome diverse, or species-rich. To extend understanding of how shifts in area, biomes, and pollinators impact diversification in plants, we examined the relationships of these shifts to diversification across the mega-genus Salvia.

METHODS

A chronogram was developed from a supermatrix of anchored hybrid enrichment genomic data and targeted sequence data for over 500 of the nearly 1000 Salvia species. Ancestral areas and biomes were reconstructed using BioGeoBEARS. Pollinator guilds were scored, ancestral pollinators determined, shifts in pollinator guilds identified, and rates of pollinator switches compared.

KEY RESULTS

A well-resolved phylogenetic backbone of Salvia and updated subgeneric designations are presented. Salvia originated in Southwest Asia in the Oligocene and subsequently dispersed worldwide. Biome shifts are frequent from a likely ancestral lineage utilizing broadleaf and/or coniferous forests and/or arid shrublands. None of the four species diversification shifts are correlated to shifts in biomes. Shifts in pollination system are not correlated to species diversification shifts, except for one hummingbird shift that precedes a major shift in diversification near the crown of New World subgen. Calosphace. Multiple reversals back to bee pollination occurred within this hummingbird clade.

CONCLUSIONS

Salvia diversified extensively in different continents, biomes, and with both bee and bird pollinators. The lack of tight correlation of area, biome, and most pollinator shifts to the four documented species diversification shifts points to other important drivers of speciation in Salvia.

The plastid and mitochondrial genomes of Eucalyptus grandis

BACKGROUND

Land plant organellar genomes have significant impact on metabolism and adaptation, and as such, accurate assembly and annotation of plant organellar genomes is an important tool in understanding the evolutionary history and interactions between these genomes. Intracellular DNA transfer is ongoing between the nuclear and organellar genomes, and can lead to significant genomic variation between, and within, species that impacts downstream analysis of genomes and transcriptomes.

RESULTS

In order to facilitate further studies of cytonuclear interactions in Eucalyptus, we report an updated annotation of the E. grandis plastid genome, and the second sequenced and annotated mitochondrial genome of the Myrtales, that of E. grandis. The 478,813 bp mitochondrial genome shows the conserved protein coding regions and gene order rearrangements typical of land plants. There have been widespread insertions of organellar DNA into the E. grandis nuclear genome, which span 141 annotated nuclear genes. Further, we identify predicted editing sites to allow for the discrimination of RNA-sequencing reads between nuclear and organellar gene copies, finding that nuclear copies of organellar genes are not expressed in E. grandis.

CONCLUSIONS

The implications of organellar DNA transfer to the nucleus are often ignored, despite the insight they can give into the ongoing evolution of plant genomes, and the problems they can cause in many applications of genomics. Future comparisons of the transcription and regulation of organellar genes between Eucalyptus genotypes may provide insight to the cytonuclear interactions that impact economically important traits in this widely grown lignocellulosic crop species.

Seed size evolution and biogeography of Plukenetia (Euphorbiaceae), a pantropical genus with traditionally cultivated oilseed species

BACKGROUND

Plukenetia is a small pantropical genus of lianas and vines with variably sized edible oil-rich seeds that presents an ideal system to investigate neotropical and pantropical diversification patterns and seed size evolution. We assessed the biogeography and seed evolution of Plukenetia through phylogenetic analyses of a 5069 character molecular dataset comprising five nuclear and two plastid markers for 86 terminals in subtribe Plukenetiinae (representing 20 of ~ 23 Plukenetia species). Two nuclear genes, KEA1 and TEB, were used for phylogenetic reconstruction for the first time. Our goals were: (1) produce a robust, time-dependent evolutionary framework for Plukenetia using BEAST; (2) reconstruct its biogeographical history with ancestral range estimation in BIOGEOBEARS; (3) define seed size categories; (4) identify patterns of seed size evolution using ancestral state estimation; and (5) conduct regression analyses with putative drivers of seed size using the threshold model.

RESULTS

Plukenetia was resolved into two major groups, which we refer to as the pinnately- and palmately-veined clades. Our analyses suggest Plukenetia originated in the Amazon or Atlantic Forest of Brazil during the Oligocene (28.7 Mya) and migrated/dispersed between those regions and Central America/Mexico throughout the Miocene. Trans-oceanic dispersals explain the pantropical distribution of Plukenetia, including from the Amazon to Africa in the Early Miocene (17.4 Mya), followed by Africa to Madagascar and Africa to Southeast Asia in the Late Miocene (9.4 Mya) and Pliocene (4.5 Mya), respectively. We infer a single origin of large seeds in the ancestor of Plukenetia. Seed size fits a Brownian motion model of trait evolution and is moderately to strongly associated with plant size, fruit type/dispersal syndrome, and seedling ecology. Biome shifts were not drivers of seed size, although there was a weak association with a transition to fire prone semi-arid savannas.

CONCLUSIONS

The major relationships among the species of Plukenetia are now well-resolved. Our biogeographical analyses support growing evidence that many pantropical distributions developed by periodic trans-oceanic dispersals throughout the Miocene and Pliocene. Selection on a combination of traits contributed to seed size variation, while movement between forest edge/light gap and canopy niches likely contributed to the seed size extremes in Plukenetia.

The herophile species group of Calisto (Lepidoptera : Nymphalidae : Satyrinae), new taxa and historical biogeography

The genus Calisto is endemic tothe West Indiesand the only representative there of the Satyrinae. Here wereconstruct the evolutionary relationshipsof the herophile group and describe five new species from Cuba: Calisto gundlachi sp. nov., Calisto siguanensis sp. nov., Calisto disjunctus sp. nov., Calisto sharkeyae sp. nov. and Calisto lastrai sp. nov.We employ one mitochondrial and four nuclear markers to assess the phylogenetic position, Maximum Likelihood and Bayesian Inference approaches, of the new taxa. Our phylogenetic trees yielded two strongly supported main clades with four of the new species included within them and C. sharkeyae as sister group to the rest of the major main clade. We conduct time-divergence estimations and ancestral area reconstructions using BEAST and BioGeoBEARS. The group originated 12.15 million years ago during the middle Miocene in north-eastern Cuba, Nipe-Sagua-Baracoa Massif. After 6 million years of in situ evolution most lineages started to colonise other Cuban territories and the Bahamas. This scenario is consistent with key geological events, including the closure of the western Havana–Matanzas channel 8–6 million years ago, the uplift of the Sierra Maestra 6–5 million years ago, and the land connections among Cuban regions during the Miocene–Pleistocene sea level drops. Dispersal and vicariance processes may have occurred, with populations surviving floodings on the major and minor mountain ranges, which remained as ‘islands’. http://zoobank.org/urn:lsid:zoobank.org:act:03690F79-F938-42A0-B234-4A228D5C1913

Beyond buzz-pollination - departures from an adaptive plateau lead to new pollination syndromes

Pollination syndromes describe recurring adaptation to selection imposed by distinct pollinators. We tested for pollination syndromes in Merianieae (Melastomataceae), which contain bee- (buzz-), hummingbird-, flowerpiercer-, passerine-, bat- and rodent-pollinated species. Further, we explored trait changes correlated with the repeated shifts away from buzz-pollination, which represents an 'adaptive plateau' in Melastomataceae. We used random forest analyses to identify key traits associated with the different pollinators of 19 Merianieae species and estimated the pollination syndromes of 42 more species. We employed morphospace analyses to compare the morphological diversity (disparity) among syndromes. We identified three pollination syndromes ('buzz-bee', 'mixed-vertebrate' and 'passerine'), characterized by different pollen expulsion mechanisms and reward types, but not by traditional syndrome characters. Further, we found that 'efficiency' rather than 'attraction' traits were important for syndrome circumscription. Contrary to syndrome theory, our study supports the pooling of different pollinators (hummingbirds, bats, rodents and flowerpiercers) into the 'mixed-vertebrate' syndrome, and we found that disparity was highest in the 'buzz-bee' syndrome. We conclude that the highly adaptive buzz-pollination system may have prevented shifts towards classical pollination syndromes, but provided the starting point for the evolution of a novel set of distinct syndromes, all having retained multifunctional stamens that provide pollen expulsion, reward and attraction.