Periodontal disease increases the severity of chronic obstructive pulmonary disease: a Mendelian randomization study

BACKGROUND

Recent research suggests that periodontitis can increase the risk of chronic obstructive pulmonary disease (COPD). In this study, we performed two-sample Mendelian randomization (MR) and investigated the causal effect of periodontitis (PD) on the genetic prediction of COPD. The study aimed to estimate how exposures affected outcomes.

METHODS

Published data from the Gene-Lifestyle Interaction in the Dental Endpoints (GLIDE) Consortium's genome-wide association studies (GWAS) for periodontitis (17,353 cases and 28,210 controls) and COPD (16,488 cases and 169,688 controls) from European ancestry were utilized. This study employed a two-sample MR analysis approach and applied several complementary methods, including weighted median, inverse variance weighted (IVW), and MR-Egger regression. Multivariable Mendelian randomization (MVMR) analysis was further conducted to mitigate the influence of smoking on COPD.

RESULTS

We chose five single-nucleotide polymorphisms (SNPs) as instrumental variables for periodontitis. A strong genetically predicted causal link between periodontitis and COPD, that is, periodontitis as an independent risk factor for COPD was detected. PD (OR = 1.102951, 95% CI: 1.005-1.211, p = 0.039) MR-Egger regression and weighted median analysis results were coincident with those of the IVW method. According to the sensitivity analysis, horizontal pleiotropy's effect on causal estimations seemed unlikely. However, reverse MR analysis revealed no significant genetic causal association between COPD and periodontitis. IVW (OR = 1.048 > 1, 95%CI: 0.973-1.128, p = 0.2082) MR Egger (OR = 0.826, 95%CI:0.658-1.037, p = 0.1104) and weighted median (OR = 1.043, 95%CI: 0.941-1.156, p = 0.4239). The results of multivariable Mendelian randomization (MVMR) analysis, after adjusting for the confounding effect of smoking, suggest a potential causal relationship between periodontitis and COPD (P = 0.035).

CONCLUSION

In this study, periodontitis was found to be independent of COPD and a significant risk factor, providing new insights into periodontitis-mediated mechanisms underlying COPD development.

The evolution of ephemeral flora in Xinjiang, China: insights from plastid phylogenomic analyses of Brassicaceae

BACKGROUND

The ephemeral flora of northern Xinjiang, China, plays an important role in the desert ecosystems. However, the evolutionary history of this flora remains unclear. To gain new insights into its origin and evolutionary dynamics, we comprehensively sampled ephemeral plants of Brassicaceae, one of the essential plant groups of the ephemeral flora.

RESULTS

We reconstructed a phylogenetic tree using plastid genomes and estimated their divergence times. Our results indicate that ephemeral species began to colonize the arid areas in north Xinjiang during the Early Miocene and there was a greater dispersal of ephemeral species from the surrounding areas into the ephemeral community of north Xinjiang during the Middle and Late Miocene, in contrast to the Early Miocene or Pliocene periods.

CONCLUSIONS

Our findings, together with previous studies, suggest that the ephemeral flora originated in the Early Miocene, and species assembly became rapid from the Middle Miocene onwards, possibly attributable to global climate changes and regional geological events.

A DNA barcode library for woody plants in tropical and subtropical China

The application of DNA barcoding has been significantly limited by the scarcity of reliable specimens and inadequate coverage and replication across all species. The deficiency of DNA barcode reference coverage is particularly striking for highly biodiverse subtropical and tropical regions. In this study, we present a comprehensive barcode library for woody plants in tropical and subtropical China. Our dataset includes a standard barcode library comprising the four most widely used barcodes (rbcL, matK, ITS, and ITS2) for 2,520 species from 4,654 samples across 49 orders, 144 families, and 693 genera, along with 79 samples identified at the genus level. This dataset also provides a super-barcode library consisting of 1,239 samples from 1,139 species, 411 genera, 113 families, and 40 orders. This newly developed library will serve as a valuable resource for DNA barcoding research in tropical and subtropical China and bordering countries, enable more accurate species identification, and contribute to the conservation and management of tropical and subtropical forests.

SHED-exos promote saliva secretion by suppressing p-ERK1/2-mediated apoptosis in glandular cells

OBJECTIVES

Confirm that stem cells from human exfoliated deciduous teeth-derived exosomes (SHED-exos) can limit inflammation-triggered epithelial cell apoptosis and explore the molecular mechanism.

METHODS

SHED-exos were injected into the submandibular glands (SMGs) of non-obese diabetic (NOD) mice, an animal model of Sjögren's syndrome (SS). Cell death was evaluated by western blotting and terminal deoxynucleotidyl transferase-mediated dUTP nick-end labelling staining.

RESULTS

SHED-exos treatment promoted the saliva flow rates of NOD mice, accompanied by decreased cleaved caspase-3 levels and apoptotic cell numbers in SMGs. SHED-exos inhibited autophagy, pyroptosis, NETosis, ferroptosis, necroptosis and oxeiptosis marker expression in SS-damaged glands. Mechanistically, Kyoto Encyclopedia of Genes and Genomes analysis of exosomal miRNAs suggested that the rat sarcoma virus (RAS)/mitogen-activated protein kinase kinase (MEK)/extracellular signal-regulated kinase (ERK) pathway might play an important role. In vivo, the expression of Kirsten RAS, Harvey RAS, MEK1/2 and p-ERK1/2 was upregulated in SMGs, and this change was blocked by SHED-exos treatment. In vitro, SHED-exos suppressed p-ERK1/2 activation and increased cleaved caspase-3 and apoptotic cell numbers, which were induced by IFN-γ.

CONCLUSION

SHED-exos suppress epithelial cell death, which is responsible for promoting salivary secretion. SHED-exos inhibited inflammation-triggered epithelial cell apoptosis by suppressing p-ERK1/2 activation, which is involved in these effects.

The phased telomere-to-telomere reference genome of Musa acuminata, a main contributor to banana cultivars

Musa acuminata is a main wild contributor to banana cultivars. Here, we reported a haplotype-resolved and telomere-to-telomere reference genome of M. acuminata by incorporating PacBio HiFi reads, Nanopore ultra-long reads, and Hi-C data. The genome size of the two haploid assemblies was estimated to be 469.83 Mb and 470.21 Mb, respectively. Multiple assessments confirmed the contiguity (contig N50: 16.53 Mb and 18.58 Mb; LAI: 20.18 and 19.48), completeness (BUSCOs: 98.57% and 98.57%), and correctness (QV: 45.97 and 46.12) of the genome. The repetitive sequences accounted for about half of the genome size. In total, 40,889 and 38,269 protein-coding genes were annotated in the two haploid assemblies, respectively, of which 9.56% and 3.37% were newly predicted. Genome comparison identified a large reciprocal translocation involving 3 Mb and 10 Mb from chromosomes 01 and 04 within M. acuminata. This reference genome of M. acuminata provides a valuable resource for further understanding of subgenome evolution of Musa species, and precise genetic improvement of banana.

Telomere-to-telomere haplotype-resolved reference genome reveals subgenome divergence and disease resistance in triploid Cavendish banana

Banana is one of the most important crops of the world. Cavendish-type bananas, which have a monospecific Musa acuminata origin (AAA), account for around half of the global banana production, thereby are of great significance for human societies. However, until now, the high-quality haplotype-resolved reference genome was still undecoded for banana cultivars. Here, we reported the telomere-to-telomere (T2T) and haplotype-resolved reference genome of 'Baxijiao' (Cavendish) consisting of three haploid assemblies. The sizes of the three haploid assemblies were estimated to be 477.16 Mb, 477.18 Mb, and 469.57 Mb, respectively. Although with monospecific origins, the three haploid assemblies showed great differences with low levels of sequence collinearity. Several large reciprocal translocations were identified among chromosomes 1, 4, and 7. An expansion of gene families that might affect fruit quality and aroma was detected, such as those belonging to sucrose/disaccharide/oligosaccharide catabolic processes, sucrose metabolic process, starch metabolic process, and aromatic compound biosynthetic process. Besides, an expansion of gene families related to anther and pollen development was observed, which could be associated with parthenocarpy and sterility of the Cavendish cultivar. Finally, much fewer resistance genes were identified in 'Baxijiao' than in M. acuminata, particularly in the gene clusters in chromosomes 3 and 10, providing potential targets to explore for molecular analysis of disease resistance in banana. This T2T haplotype-resolved reference genome will thus be a valuable genetic resource for biological studies, molecular breeding, and genetic improvement of banana.

Genetic Predisposition to Periodontitis and Risk of Migraine: A Two-Sample Mendelian Randomization Study

INTRODUCTION

Previous observational studies have associated periodontitis (PD) with migraine; however, the results are inconclusive and the causality of the association between PD and migraine remains unclear. This two-sample Mendelian randomization (MR) study was performed to explore the bi-directional causal relationship between PD and migraine.

METHODS

To investigate the relationship between PD (17,353 cases; 28,210 controls) and migraine (1072 cases; 360,122 controls), we used genetic tools from the largest available genome-wide association study of European descent. Inverse variance-weighted (IVW) and a series of sensitivity analyses were used to explore the association between migraine and PD. We performed an MR study using seven SNPs (single nucleotide polymorphisms) as instrumental variables for PD to investigate the causal relationship between migraine and PD.

RESULTS

We found no significant causal relationship between PD and migraine (odds ratio, OR = 1.000; 95% confidence interval, CI = 0.99-1.00; p = 0.65). Similarly, no evidence supported a causal relationship between migraine and PD (OR = 0.07; CI = 2.04 × 10-9-2.65 × 106; p = 0.77). A sensitivity analysis revealed that no potential polymorphic effect (p = 0.356) and heterogeneity (p = 0.652) exists for the variants used in constructing the genetic instrument.

CONCLUSIONS

Based on the results of our MR study, there is no causal relationship between PD and migraines or migraines and PD.

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.

Flora diversity survey and establishment of a plant DNA barcode database of Lomas ecosystems in Peru

Lomas formations or "fog oases" are islands of vegetation in the desert belt of the west coast of South America, with a unique vegetation composition among the world's deserts. However, plant diversity and conservation studies have long been neglected, and there exists a severe gap in plant DNA sequence information. To address the lack of DNA information, we conducted field collections and laboratory DNA sequencing to establish a DNA barcode reference library of Lomas plants from Peru. This database provides 1,207 plant specimens and 3,129 DNA barcodes data corresponding with collections from 16 Lomas locations in Peru, during 2017 and 2018. This database will facilitate both rapid species identification and basic studies on plant diversity, thereby enhancing our understanding of Lomas flora's composition and temporal variation, and providing valuable resources for conserving plant diversity and maintaining the stability of the fragile Lomas ecosystems.

Plant DNA-barcode library for native flowering plants in the arid region of northwestern China

DNA barcoding is a well-established tool for rapid species identification and biodiversity monitoring. A reliable and traceable DNA barcode reference library with extensive coverage is necessary but unavailable for many geographic regions. The arid region in northwestern China, a vast area of about 2.5 million km² , is ecologically fragile and often overlooked in biodiversity studies. In particular, DNA barcode data from the arid region in China are lacking. We develop and evaluate the efficacy of an extensive DNA-barcode library for native flowering plants in the arid region of northwestern China. Plant specimens were collected, identified, and vouchered for this purpose. The database utilized four DNA barcode markers, namely rbcL, matK, ITS, and ITS2, for 1,816 accessions (representing 890 species from 386 genera and 72 families) and consisted of 5,196 barcode sequences. Individual barcodes varied in resolution rates; species- and genus-level rates for rbcL, matK, ITS, and ITS2 were 79.9-51.1%/76.1%, 79.9-67.2%/88.9%, 85.0-72.0%/88.2%, and 81.0-67.4%/84.9%, respectively. The three-barcode combination of rbcL + matK + ITS (RMI) revealed a higher species- and genus-level resolution (75.5%/92.1%, respectively). A total of 110 plastomes were newly generated as super-barcodes to increase species resolution for seven species-rich genera, namely Astragalus, Caragana, Lactuca, Lappula, Lepidium, Silene, and Zygophyllum. Plastomes revealed higher species resolution compared to standard DNA barcodes and their combination. We suggest future databases include super-barcodes, especially for species-rich and complex genera. The plant DNA barcode library in the current study provides a valuable resource for future biological investigations in the arid regions of China.

A plastid phylogenomic framework for the palm family (Arecaceae)

BACKGROUND

Over the past decade, phylogenomics has greatly advanced our knowledge of angiosperm evolution. However, phylogenomic studies of large angiosperm families with complete species or genus-level sampling are still lacking. The palms, Arecaceae, are a large family with ca. 181 genera and 2600 species and are important components of tropical rainforests bearing great cultural and economic significance. Taxonomy and phylogeny of the family have been extensively investigated by a series of molecular phylogenetic studies in the last two decades. Nevertheless, some phylogenetic relationships within the family are not yet well-resolved, especially at the tribal and generic levels, with consequent impacts for downstream research.

RESULTS

Plastomes of 182 palm species representing 111 genera were newly sequenced. Combining these with previously published plastid DNA data, we were able to sample 98% of palm genera and conduct a plastid phylogenomic investigation of the family. Maximum likelihood analyses yielded a robustly supported phylogenetic hypothesis. Phylogenetic relationships among all five palm subfamilies and 28 tribes were well-resolved, and most inter-generic phylogenetic relationships were also resolved with strong support.

CONCLUSIONS

The inclusion of nearly complete generic-level sampling coupled with nearly complete plastid genomes strengthened our understanding of plastid-based relationships of the palms. This comprehensive plastid genome dataset complements a growing body of nuclear genomic data. Together, these datasets form a novel phylogenomic baseline for the palms and an increasingly robust framework for future comparative biological studies of this exceptionally important plant family.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

[Comparison of dentinal microcracks caused by new generation of single reciprocating files after root canal preparation]

PURPOSE

To compare the incidence of dentinal microcracks after root canal preparation by new generation of nickel-titanium instrument WaveOne Gold, Reciproc Blue with previous WaveOne and Reciproc.

METHODS

Ninety extracted single-rooted mandibular premolars were randomly divided into 6 groups(n=15). The root canals were instrumented by using Hand K files, WaveOne, Reciproc, WaveOne Gold and Reciproc Blue. Fifteen teeth were left unprepared and served as negative controls. The root canals were all prepared to 25#. The roots were then sectioned at 3 mm, 6 mm and 9 mm from the apical orifice using a hard tissue slicer. The slices were observed under stereoscopic microscope at ×25 magnification. SPSS 17.0 software package was used for statistical analysis.

RESULTS

No dentinal microcrack was found in the hand K files group and negative control group. The reciprocating single files WaveOne, WaveOne Gold, Reciproc and Reciproc Blue all produced dentinal microcracks after root canal preparation. The WaveOne generated the most dentinal microcracks than the hand K files(P＜0.05), and the microcracks were mainly concentrated in the middle part of the root. The number of dentinal microcracks caused by Reciproc and Reciproc Blue was the same, with no significant difference(P＞0.05).

CONCLUSIONS

The new generation of reciprocating files of WaveOne Gold and Reciproc Blue may not increase the incidence of dentinal microcracks after root canal preparation.

Comparative analyses of three complete Primula mitogenomes with insights into mitogenome size variation in Ericales

BACKGROUND

Although knowledge of the sizes, contents, and forms of plant mitochondrial genomes (mitogenomes) is increasing, little is known about the mechanisms underlying their structural diversity. Evolutionary information on the mitogenomes of Primula, an important ornamental taxon, is more limited than the information on their nuclear and plastid counterparts, which has hindered the comprehensive understanding of Primula mitogenomic diversity and evolution. The present study reported and compared three Primula mitogenomes and discussed the size expansion of mitogenomes in Ericales.

RESULTS

Mitogenome master circles were sequenced and successfully assembled for three Primula taxa and were compared with publicly available Ericales mitogenomes. The three mitogenomes contained similar gene contents and varied primarily in their structures. The Primula mitogenomes possessed relatively high nucleotide diversity among all examined plant lineages. In addition, high nucleotide diversity was found among Primula species between the Mediterranean and Himalaya-Hengduan Mountains. Most predicted RNA editing sites appeared in the second amino acid codon, increasing the hydrophobic character of the protein. An early stop in atp6 caused by RNA editing was conserved across all examined Ericales species. The interfamilial relationships within Ericales and interspecific relationships within Primula could be well resolved based on mitochondrial data. Transfer of the two longest mitochondrial plastid sequences (MTPTs) occurred before the divergence of Primula and its close relatives, and multiple independent transfers could also occur in a single MTPT sequence. Foreign sequence [MTPTs and mitochondrial nuclear DNA sequences (NUMTs)] uptake and repeats were to some extent associated with changes in Ericales mitogenome size, although none of these relationships were significant overall.

CONCLUSIONS

The present study revealed relatively conserved gene contents, gene clusters, RNA editing, and MTPTs but considerable structural variation in Primula mitogenomes. Relatively high nucleotide diversity was found in the Primula mitogenomes. In addition, mitogenomic genes, collinear gene clusters, and locally collinear blocks (LCBs) all showed phylogenetic signals. The evolutionary history of MTPTs in Primula was complicated, even in a single MTPT sequence. Various reasons for the size variation observed in Ericales mitogenomes were found.

The banana genome hub: a community database for genomics in the Musaceae

The Banana Genome Hub provides centralized access for genome assemblies, annotations, and the extensive related omics resources available for bananas and banana relatives. A series of tools and unique interfaces are implemented to harness the potential of genomics in bananas, leveraging the power of comparative analysis, while recognizing the differences between datasets. Besides effective genomic tools like BLAST and the JBrowse genome browser, additional interfaces enable advanced gene search and gene family analyses including multiple alignments and phylogenies. A synteny viewer enables the comparison of genome structures between chromosome-scale assemblies. Interfaces for differential expression analyses, metabolic pathways and GO enrichment were also added. A catalogue of variants spanning the banana diversity is made available for exploration, filtering, and export to a wide variety of software. Furthermore, we implemented new ways to graphically explore gene presence-absence in pangenomes as well as genome ancestry mosaics for cultivated bananas. Besides, to guide the community in future sequencing efforts, we provide recommendations for nomenclature of locus tags and a curated list of public genomic resources (assemblies, resequencing, high density genotyping) and upcoming resources-planned, ongoing or not yet public. The Banana Genome Hub aims at supporting the banana scientific community for basic, translational, and applied research and can be accessed at https://banana-genome-hub.southgreen.fr.

Plastome structure, phylogenomics, and divergence times of tribe Cinnamomeae (Lauraceae)

BACKGROUND

Tribe Cinnamomeae is a species-rich and ecologically important group in tropical and subtropical forests. Previous studies explored its phylogenetic relationships and historical biogeography using limited loci, which might result in biased molecular dating due to insufficient parsimony-informative sites. Thus, 15 plastomes were newly sequenced and combined with published plastomes to study plastome structural variations, gene evolution, phylogenetic relationships, and divergence times of this tribe.

RESULTS

Among the 15 newly generated plastomes, 14 ranged from 152,551 bp to 152,847 bp, and the remaining one (Cinnamomum chartophyllum XTBGLQM0164) was 158,657 bp. The inverted repeat (IR) regions of XTBGLQM0164 contained complete ycf2, trnI^CAU, rpl32, and rpl2. Four hypervariable plastid loci (ycf1, ycf2, ndhF-rpl32-trnL^UAG, and petA-psbJ) were identified as candidate DNA barcodes. Divergence times based on a few loci were primarily determined by prior age constraints rather than by DNA data. In contrast, molecular dating using complete plastid protein-coding genes (PCGs) was determined by DNA data rather than by prior age constraints. Dating analyses using PCGs showed that Cinnamomum sect. Camphora diverged from C. sect. Cinnamomum in the late Oligocene (27.47 Ma).

CONCLUSIONS

This study reports the first case of drastic IR expansion in tribe Cinnamomeae, and indicates that plastomes have sufficient parsimony-informative sites for molecular dating. Besides, the dating analyses provide preliminary insights into the divergence time within tribe Cinnamomeae and can facilitate future studies on its historical biogeography.

Fingerprints of climatic changes through the late Cenozoic in southern Asian flora: Magnolia section Michelia (Magnoliaceae)

BACKGROUND AND AIMS

Ongoing global warming is a challenge for humankind. A series of drastic climatic changes have been proven to have occurred throughout the Cenozoic based on a variety of geological evidence, which helps to better understand our planet's future climate. Notably, extant biomes have recorded drastic environmental shifts. The climate in southern Asia, which hosts high biodiversity, is deeply impacted by the Asian monsoon. The origins and evolutionary dynamics of biomes occurring between the tropics and sub-tropics in southern Asia have probably been deeply impacted by climatic changes; however, these aspects remain poorly studied. We tested whether the evolutionary dynamics of the above biomes have recorded the drastic, late Cenozoic environmental shifts, by focusing on Magnolia section Michelia of the family Magnoliaceae.

METHODS

We established a fine time-calibrated phylogeny of M. section Michelia based on complete plastid genomes and inferred its ancestral ranges. Finally, we estimated the evolutionary dynamics of this section through time, determining its diversification rate and the dispersal events that occurred between tropical and sub-tropical areas.

KEY RESULTS

The tropical origin of M. section Michelia was dated to the late Oligocene; however, the diversification of its core group (i.e. M. section Michelia subsection Michelia) has occurred mainly from the late Miocene onward. Two key evolutionary shifts (dated approx. 8 and approx. 3 million years ago, respectively) were identified, each of them probably in response to drastic climatic changes.

CONCLUSION

Here, we inferred the underlying evolutionary dynamics of biomes in southern Asia, which probably reflect late Cenozoic climatic changes. The occurrence of modern Asian monsoons was probably fundamental for the origin of M. section Michelia; moreover, the occurrence of asymmetric dispersal events between the tropics and sub-tropics hint at an adaptation strategy of M. section Michelia to global cooling, in agreement with the tropical conservatism hypothesis.

The relative role of plasticity and demographic history in Capsella bursa-pastoris: a common garden experiment in Asia and Europe

The colonization success of a species depends on the interplay between its phenotypic plasticity, adaptive potential and demographic history. Assessing their relative contributions during the different phases of a species range expansion is challenging, and requires large-scale experiments. Here, we investigated the relative contributions of plasticity, performance and demographic history to the worldwide expansion of the shepherd's purse, Capsella bursa-pastoris. We installed two large common gardens of the shepherd's purse, a young, self-fertilizing, allopolyploid weed with a worldwide distribution. One common garden was located in Europe, the other in Asia. We used accessions from three distinct genetic clusters (Middle East, Europe and Asia) that reflect the demographic history of the species. Several life-history traits were measured. To explain the phenotypic variation between and within genetic clusters, we analysed the effects of (i) the genetic clusters, (ii) the phenotypic plasticity and its association to fitness and (iii) the distance in terms of bioclimatic variables between the sampling site of an accession and the common garden, i.e. the environmental distance. Our experiment showed that (i) the performance of C. bursa-pastoris is closely related to its high phenotypic plasticity; (ii) within a common garden, genetic cluster was a main determinant of phenotypic differences; and (iii) at the scale of the experiment, the effect of environmental distance to the common garden could not be distinguished from that of genetic clusters. Phenotypic plasticity and demographic history both play important role at different stages of range expansion. The success of the worldwide expansion of C. bursa-pastoris was undoubtedly influenced by its strong phenotypic plasticity.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae

BACKGROUND

Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales.

FINDINGS

Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11.

CONCLUSIONS

The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.

Comparative plastome analysis of Musaceae and new insights into phylogenetic relationships

Background

Musaceae is an economically important family consisting of 70-80 species. Elucidation of the interspecific relationships of this family is essential for a more efficient conservation and utilization of genetic resources for banana improvement. However, the scarcity of herbarium specimens and quality molecular markers have limited our understanding of the phylogenetic relationships in wild species of Musaceae. Aiming at improving the phylogenetic resolution of Musaceae, we analyzed a comprehensive set of 49 plastomes for 48 species/subspecies representing all three genera of this family.

Results

Musaceae plastomes have a relatively well-conserved genomic size and gene content, with a full length ranging from 166,782 bp to 172,514 bp. Variations in the IR borders were found to show phylogenetic signals to a certain extent in Musa. Codon usage bias analysis showed different preferences for the same codon between species and three genera and a common preference for A/T-ending codons. Among the two genes detected under positive selection (dN/dS > 1), ycf2 was indicated under an intensive positive selection. The divergent hotspot analysis allowed the identification of four regions (ndhF-trnL, ndhF, matK-rps16, and accD) as specific DNA barcodes for Musaceae species.

Bayesian and maximum likelihood phylogenetic analyses using full plastome resulted in nearly identical tree topologies with highly supported relationships between species. The monospecies genus Musella is sister to Ensete, and the genus Musa was divided into two large clades, which corresponded well to the basic number of n = x = 11 and n = x =10/9/7, respectively. Four subclades were divided within the genus Musa. A dating analysis covering the whole Zingiberales indicated that the divergence of Musaceae family originated in the Palaeocene (59.19 Ma), and the genus Musa diverged into two clades in the Eocene (50.70 Ma) and then started to diversify from the late Oligocene (29.92 Ma) to the late Miocene. Two lineages (Rhodochlamys and Australimusa) radiated recently in the Pliocene /Pleistocene periods.

Conclusions

The plastome sequences performed well in resolving the phylogenetic relationships of Musaceae and generated new insights into its evolution. Plastome sequences provided valuable resources for population genetics and phylogenetics at lower taxon.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08454-3.

Plastid phylogenomics of tribe Perseeae (Lauraceae) yields insights into the evolution of East Asian subtropical evergreen broad-leaved forests

BACKGROUND

The East Asian subtropical evergreen broad-leaved forests (EBLFs) harbor remarkable biodiversity. However, their historical assembly remains unclear. To gain new insights into the assembly of this biome, we generated a molecular phylogeny of one of its essential plant groups, the tribe Perseeae (Lauraceae).

RESULTS

Our plastid tree topologies were robust to analyses based on different plastid regions and different strategies for data partitioning, nucleotide substitution saturation, and gap handling. We found that tribe Perseeae comprised six major clades and began to colonize the subtropical EBLFs of East Asia in the early Miocene. The diversification rates of tribe Perseeae accelerated twice in the late Miocene.

CONCLUSIONS

Our findings suggest that the intensified precipitation in East Asia in the early Miocene may have facilitated range expansions of the subtropical EBLFs and establishment of tribe Perseeae within this biome. By the late Miocene, species assembly and diversification within the EBLFs had become rapid.

The complete chloroplast genome of the ephemeral plant Isatis minima (brassicaceae) of northwest China

The complete chloroplast genome of Isatis minima, a typical ephemeral plant of Brassicaceae in the Central Asia desert, was sequenced and characterized in this study. The genome 153,642 bp in size, contains a typical quadripartite genome organization including LSC and SSC regions of 83,423 bp and 17,709 bp, and two copies of the IR regions of 26,255 bp. It has 113 unique genes, including 79 protein-coding, 30 tRNA, and four rRNA genes. Phylogenetic analysis fully resolved I. minima in a monophyletic clade with I. tinctoria. This bioinformatic data contributes to the phylogenetics systematics and evolutionary history of Brassicaceae.

Plastome-based phylogeny improves community phylogenetics of subtropical forests in China

Phylogenetic trees have been extensively used in community ecology. However, how the phylogeny construction affects ecological inferences is poorly understood. In this study, we constructed three different types of phylogenetic trees (a synthetic-tree generated using V.PhyloMaker, a barcode-tree generated using rbcL+matK+trnH-psbA, and a plastome-tree generated from plastid genomes) that represented an increasing level of phylogenetic resolution among 580 woody plant species from six forest dynamic plots in subtropical evergreen broadleaved forests of China. We then evaluated the performance of each phylogeny in estimations of community phylogenetic structure, turnover and phylogenetic signal in functional traits. As expected, the plastome-tree was most resolved and most supported for relationships among species. For local phylogenetic structure, the three trees showed consistent results with Faith's PD and MPD; however, only the synthetic-tree produced significant clustering patterns using MNTD for some plots. For phylogenetic turnover, contrasting results between the molecular trees and the synthetic-tree occurred only with nearest neighbor distance. The barcode-tree agreed more with the plastome-tree than the synthetic-tree for both phylogenetic structure and turnover. For functional traits, both the barcode-tree and plastome-tree detected phylogenetic signal in maximum height, but only the plastome-tree detected signal in leaf width. This is the first study that uses plastid genomes in large-scale community phylogenetics. Our results highlight the improvement of plastome-trees over barcode-trees and synthetic-trees for the analyses studied here. Our results also point to the possibility of type I and II errors in estimation of phylogenetic structure and turnover and detection of phylogenetic signal when using synthetic-trees.

The complete chloroplast genome of Lithocarpus hancei (Benth.) Rehd (Fagaceae) from Zhejiang, China

Lithocarpus hancei (Benth.) Rehd is a widely distributed evergreen tree with broad-leaves that dominates the lower stories of the forest in China. Here, we sequenced and assembled the complete chloroplast genome of L. hancei. The genome is 161,304 bp with one large single copy (LSC: 90,585 bp), one small single copy (SSC: 18,959 bp), and two inverted repeat (IR) regions (IRa and IRb, each 25,880 bp). It contains 117 genes, including 80 protein-coding genes, 33 tRNA genes, and four rRNA genes. Phylogenetic analysis of 21 representative cp genomes of the Fagaceae suggests Lithocarpus is monophyletic with strong bootstrap support and also that L. hancei is closely related to L. polystachyus. The cp genome is important for constructing a robust phylogeny of Lithocarpus and Fagaceae for future study.

Draft genome of Puya raimondii (Bromeliaceae), the Queen of the Andes

Puya raimondii, the Queen of the Andes, is an endangered high Andean species in the Bromeliaceae family. Here, we report its first genome to promote its conservation and evolutionary study. Comparative genomics showed P. raimondii diverged from Ananas comosus about 14.8 million years ago, and the long terminal repeats were likely to contribute to the genus diversification in last 3.5 million years. The gene families related to plant reproductive development and stress responses significantly expanded in the genome. At the same time, gene families involved in disease defense, photosynthesis and carbohydrate metabolism significantly contracted, which may be an evolutionary strategy to adapt to the harsh conditions in high Andes. The demographic history analysis revealed the P. raimondii population size sharply declined in the Pleistocene and then increased in the Holocene. We also designed and tested 46 pairs of universal primers for amplifying orthologous single-copy nuclear genes in Puya species.

Comparative Analyses of the Chloroplast Genomes of Patchouli Plants and Their Relatives in Pogostemon (Lamiaceae)

Pogostemon Desf., the largest genus of the tribe Pogostemoneae (Lamiaceae), consists of ca. 80 species distributed mainly from South and Southeast Asia to China. The genus contains many patchouli plants, which are of great economic importance but taxonomically difficult. Therefore, it is necessary to characterize more chloroplast (cp) genomes for infrageneric phylogeny analyses and species identification of Pogostemon, especially for patchouli plants. In this study, we newly generated four cp genomes for three patchouli plants (i.e., Pogostemon plectranthoides Desf., P. septentrionalis C. Y. Wu et Y. C. Huang, and two cultivars of P. cablin (Blanoco) Benth.). Comparison of all samples (including online available cp genomes of P. yatabeanus (Makino) Press and P. stellatus (Lour.) Kuntze) suggested that Pogostemon cp genomes are highly conserved in terms of genome size and gene content, with a typical quadripartite circle structure. Interspecific divergence of cp genomes has been maintained at a relatively low level, though seven divergence hotspot regions were identified by stepwise window analysis. The nucleotide diversity (Pi) value was correlated significantly with gap proportion (indels), but significantly negative with GC content. Our phylogenetic analyses based on 80 protein-coding genes yielded high-resolution backbone topologies for the Lamiaceae and Pogostemon. For the overall mean substitution rates, the synonymous (d_S) and nonsynonymous (d_N) substitution rate values of protein-coding genes varied approximately threefold, while the d_N values among different functional gene groups showed a wider variation range. Overall, the cp genomes of Pogostemon will be useful for phylogenetic reconstruction, species delimitation and identification in the future.

Conflicting phylogenetic signals in plastomes of the tribe Laureae (Lauraceae)

Background

Gene tree discordance is common in phylogenetic analyses. Many phylogenetic studies have excluded non-coding regions of the plastome without evaluating their impact on tree topology. In general, plastid loci have often been treated as a single unit, and tree discordance among these loci has seldom been examined. Using samples of Laureae (Lauraceae) plastomes, we explored plastome variation among the tribe, examined the influence of non-coding regions on tree topology, and quantified intra-plastome conflict.

Results

We found that the plastomes of Laureae have low inter-specific variation and are highly similar in structure, size, and gene content. Laureae was divided into three groups, subclades I, II and III. The inclusion of non-coding regions changed the phylogenetic relationship among the three subclades. Topologies based on coding and non-coding regions were largely congruent except for the relationship among subclades I, II and III. By measuring the distribution of phylogenetic signal across loci that supported different topologies, we found that nine loci (two coding regions, two introns and five intergenic spacers) played a critical role at the contentious node.

Conclusions

Our results suggest that subclade III and subclade II are successively sister to subclade I. Conflicting phylogenetic signals exist between coding and non-coding regions of Laureae plastomes. Our study highlights the importance of evaluating the influence of non-coding regions on tree topology and emphasizes the necessity of examining discordance among different plastid loci in phylogenetic studies.

The utility of DNA barcodes to confirm the identification of palm collections in botanical gardens

The palm family (Arecaceae) is of high ecological and economic value, yet identification in the family remains a challenge for both taxonomists and horticulturalists. The family consists of approximately 2600 species across 181 genera and DNA barcoding may be a useful tool for species identification within the group. However, there have been few systematic evaluations of DNA barcodes for the palm family. In the present study, five DNA barcodes (rbcL, matK, trnH-psbA, ITS, ITS2) were evaluated for species identification ability across 669 samples representing 314 species and 100 genera in the Arecaceae, employing four analytical methods. The ITS gene region was found to not be a suitable barcode for the palm family, due in part, to low recovery rates and paralogous gene copies. Among the four analyses used, species resolution for ITS2 was much higher than that achieved with the plastid barcodes alone (rbcL, matK, trnH-psbA), and the barcode combination ITS2 + matK + rbcL gave the highest resolution among all single barcodes and their combinations, followed by ITS2 + matK. Among 669 palm samples analyzed, 110 samples (16.3%) were found to be misidentified. The 2992 DNA barcode sequences generated in this study greatly enriches the existing identification toolbox available to plant taxonomists that are interested in researching genetic relationships among palm taxa as well as for horticulturalists that need to confirm palm collections for botanical garden curation and horticultural applications. Our results indicate that the use of the ITS2 DNA barcode gene region provides a useful and cost-effective tool to confirm the identity of taxa in the Palm family.

Complete plastome sequencing resolves taxonomic relationships among species of Calligonum L. (Polygonaceae) in China

BACKGROUND

Calligonum (Polygonaceae) is distributed from southern Europe through northern Africa to central Asia, and is typically found in arid, desert regions. Previous studies have revealed that standard DNA barcodes fail to discriminate Calligonum species. In this study, the complete plastid genomes (plastome) for 32 accessions of 21 Calligonum species is sequenced to not only generate the first complete plastome sequence for the genus Calligonum but to also 1) Assess the ability of the complete plastome sequence to discern species within the group, and 2) screen the plastome sequence for a cost-effective DNA barcode that can be used in future studies to resolve taxonomic relationships within the group.

RESULTS

The whole plastomes of Calligonum species possess a typical quadripartite structure. The size of the Calligonum plastome is approximately 161 kilobase pairs (kbp), and encodes 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Based on ML phylogenetic tree analyses, the complete plastome has higher species identification (78%) than combinations of standard DNA barcodes (rbcL + matK + nrITS, 56%). Five newly screened gene regions (ndhF, trnS-G, trnC-petN, ndhF-rpl32, rpl32-trnL) had high species resolution, where ndhF and trnS-G were able to distinguish the highest proportion of Calligonum species (56%).

CONCLUSIONS

The entire plastid genome was the most effective barcode for the genus Calligonum, although other gene regions showed great potential as taxon-specific barcodes for species identification in Calligonum.

Ancient Polyploidy and Genome Evolution in Palms

Mechanisms of genome evolution are fundamental to our understanding of adaptation and the generation and maintenance of biodiversity, yet genome dynamics are still poorly characterized in many clades. Strong correlations between variation in genomic attributes and species diversity across the plant tree of life suggest that polyploidy or other mechanisms of genome size change confer selective advantages due to the introduction of genomic novelty. Palms (order Arecales, family Arecaceae) are diverse, widespread, and dominant in tropical ecosystems, yet little is known about genome evolution in this ecologically and economically important clade. Here, we take a phylogenetic comparative approach to investigate palm genome dynamics using genomic and transcriptomic data in combination with a recent, densely sampled, phylogenetic tree. We find conclusive evidence of a paleopolyploid event shared by the ancestor of palms but not with the sister clade, Dasypogonales. We find evidence of incremental chromosome number change in the palms as opposed to one of recurrent polyploidy. We find strong phylogenetic signal in chromosome number, but no signal in genome size, and further no correlation between the two when correcting for phylogenetic relationships. Palms thus add to a growing number of diverse, ecologically successful clades with evidence of whole-genome duplication, sister to a species-poor clade with no evidence of such an event. Disentangling the causes of genome size variation in palms moves us closer to understanding the genomic conditions facilitating adaptive radiation and ecological dominance in an evolutionarily successful, emblematic tropical clade.

Homeologue-specific expression divergence in the recently formed tetraploid Capsella bursa-pastoris (Brassicaceae)

Following allopolyploid formation, extensive genome evolution occurs, with the eventual loss of many homeologous gene copies. Although this process of diploidization has occurred many times independently, the evolutionary forces determining the probability and rate of gene loss remain poorly understood. Here, we conduct genome and transcriptome sequencing in a broad sample of Chinese accessions of Capsella bursa-pastoris, a recently formed allotetraploid. Our whole genome data reveal three groups of these accessions: an Eastern group from low-altitude regions, a Western group from high-altitude regions, and a much more differentiated Northwestern group. Population differentiation in total expression was limited among closely related populations; by contrast, the relative expression of the two homeologous copies closely mirrors the genome-wide SNP divergence. Consistent with this, we observe a negative correlation between expression changes in the two homeologues. However, genes showing population genomic evidence for adaptive evolution do not show an enrichment for expression divergence between homeologues, providing no clear evidence for adaptive shifts in relative gene expression. Overall, these patterns suggest that neutral drift may contribute to the population differentiation in the expression of the homeologues, and drive eventual gene loss over longer periods of time.

High cryptic species diversity is revealed by genome-wide polymorphisms in a wild relative of banana, Musa itinerans, and implications for its conservation in subtropical China

BACKGROUND

Species delimitation is a challenging but essential task in conservation biology. Morphologically similar species are sometimes difficult to recognize even after examination by experienced taxonomists. With the advent of molecular approaches in species delimitation, this hidden diversity has received much recent attention. In addition to DNA barcoding approaches, analytical tools based on the multi-species coalescence model (MSC) have been developed for species delimitation. Musa itinerans is widely distributed in subtropical Asia, and at least six varieties have been documented. However, the number of evolutionarily distinct lineages remains unknown.

RESULTS

Using genome resequencing data of five populations (making up four varieties), we examined genome-wide variation and found four varieties that were evolutionary significant units. A Bayesian Phylogenetics and Phylogeography (BP&P) analysis using 123 single copy nuclear genes support three speciation events of M. itinerans varieties with robust posterior speciation probabilities; However, a Bayes factor delimitation of species with genomic data (BFD*) analysis using 1201 unlinked single nucleotide polymorphisms gave decisive support for a five-lineage model. When reconciling divergence time estimates with a speciation time scale, a modified three-lineage model was consistent with that of BP&P, in which the speciation time of two varieties (M. itinerans var. itinerans and M. itinerans var. lechangensis) were dated to 26.2 kya and 10.7 kya, respectively. In contrast, other two varieties (M. itinerans var. chinensis and M. itinerans var. guangdongensis) diverged only 3.8 kya in the Anthropocene; this may be a consequence of genetic drift rather than a speciation event.

CONCLUSION

Our results showed that the M. itinerans species complex harbours high cryptic species diversity. We recommend that M. itinerans var. itinerans and M. itinerans var. lechangensis be elevated to subspecies status, and the extremely rare latter subspecies be given priority for conservation. We also recommend that the very recently diverged M. itinerans var. chinensis and M. itinerans var. guangdongensis should be merged under the subspecies M. itinerans var. chinensis. Finally, we speculate that species delimitation of recently diverged lineages may be more effective using genome-wide bi-allelic SNP markers with BFD* than by using unlinked loci and BP&P.

What explains high plant richness in East Asia? Time and diversification in the tribe Lysimachieae (Primulaceae)

What causes the disparity in biodiversity among regions is a fundamental question in biogeography, ecology, and evolutionary biology. Evolutionary and biogeographic processes (speciation, extinction, dispersal) directly determine species richness patterns, and can be studied using integrative phylogenetic approaches. However, the strikingly high richness of East Asia relative to other Northern Hemisphere regions remains poorly understood from this perspective. Here, for the first time, we test two general hypotheses (older colonization time, faster diversification rate) to explain this pattern, using the plant tribe Lysimachieae (Primulaceae) as a model system. We generated a new time-calibrated phylogeny for Lysimachieae (13 genes, 126 species), to estimate colonization times and diversification rates for each region and to test the relative importance of these two factors for explaining regional richness patterns. We find that neither time nor diversification rates alone explain richness patterns among regions in Lysimachieae. Instead, a new index that combines both factors explains global richness patterns in the group and their high East Asian biodiversity. Based on our results from Lysimachieae, we suggest that the high richness of plants in East Asia may be explained by a combination of older colonization times and faster diversification rates in this region.

[Establishment of rat model with diabetes mellitus and concomitant periodontitis and the carotid artery lesions in the model rats]

Objectives: To establish SD rat model with type 2 diabetes mellitus (DM) and concomitant chronic periodontitis (CP) and to evaluate the influence of periodontitis on the vascular lesions of type 2 diabetes rats. Methods: Totally 241 clean level SD rats were randomly divided into four groups, group A (normal control, NC, n=27), group B (DM, n=34), group C (CP, n=90) and group D (DM+CP, n=90). The rats of DM group were fed with high-fat and high-sugar diet for 8 to 10 weeks, and then were multiply injected with small dose streptozotocin under the condition of ice bath. Blood sugar levels after the injection were dynamically monitored at 72 h, 1 week, 2 weeks and 4 weeks, respectively. The CP model was established by means of ligation. Bilateral maxillary first and second molars were selected and ligated using 0.2 mm orthodontic wires binding with 4-0 surgical suture soaked with Porphyromonas gingivalis (Pg) suspension. After a period of 14 weeks, all the rats were put to death. Maxillary samples were subjected to methylene blue staining to observe alveolar bone loss. Bilateral carotid artery specimens were collected. The left carotid artery specimens were used to detect the prevalence of Pg using quantitative real-time PCR. The right carotid artery specimens were used to observe pathological changes. Results: Blood sugar levels of rats in group B and D increased and changed sharply after Streptozotocin injection with in 1 week. Symptoms of 'more drink, more food and body weight loss' appeared. The fasting blood glucose (FBG) was more than 7.8 mmol/L and (or) the random blood glucose (RBG) was more than 17.8 mmol/L. Both FBG and RBG became stable after 2 to 3 weeks. Levels of HbA1C in group B and D ([7.32±0.45]%, [9.41±0.45]%) were significantly higher than that of group A ([4.02±0.45]%) (P<0.01). Rats of group D were observed the most severe bone loss showing wider interdental space and furcation involvement. Pathological results of carotid artery tissues of group D showed the worst lesions including thinning and calcification of vessel walls, and breaking down or disappearance of elastic fibers. The prevalences of DNA of Pg in groups of A, B, C and D were 3/7, 3/7, 6/7 and 7/7, respectively. The bacteria numbers detected by quantitative real-time PCR in groups C and D were significantly higher than that of groups A and B (P<0.01). Conclusions: Rat model of type 2 DM with periodontitis was successfully established in the present study. Carotid artery specimens from DM+CP model rats showed typical vascular lesions such as calcification and fiber disorders. Pg was found in all carotid specimens and the highest bacteria numbers were detected in the composite model rats. The Pg might play a role in the progress of diabetes vascular lesions.

Whole genome sequencing of a banana wild relative Musa itinerans provides insights into lineage-specific diversification of the Musa genus

Crop wild relatives are valuable resources for future genetic improvement. Here, we report the de novo genome assembly of Musa itinerans, a disease-resistant wild banana relative in subtropical China. The assembled genome size was 462.1 Mb, covering 75.2% of the genome (615.2Mb) and containing 32, 456 predicted protein-coding genes. Since the approximate divergence around 5.8 million years ago, the genomes of Musa itinerans and Musa acuminata have shown conserved collinearity. Gene family expansions and contractions enrichment analysis revealed that some pathways were associated with phenotypic or physiological innovations. These include a transition from wood to herbaceous in the ancestral Musaceae, intensification of cold and drought tolerances, and reduced diseases resistance genes for subtropical marginally distributed Musa species. Prevalent purifying selection and transposed duplications were found to facilitate the diversification of NBS-encoding gene families for two Musa species. The population genome history analysis of M. itinerans revealed that the fluctuated population sizes were caused by the Pleistocene climate oscillations, and that the formation of Qiongzhou Strait might facilitate the population downsizing on the isolated Hainan Island about 10.3 Kya. The qualified assembly of the M. itinerans genome provides deep insights into the lineage-specific diversification and also valuable resources for future banana breeding.

Characterization of polymorphic microsatellite markers for Primula sikkimensis (Primulaceae) using a 454 sequencing approach

PREMISE OF THE STUDY

Microsatellite markers from Primula sikkimensis (Primulaceae) were developed for testing deep lineage divergence and speciation events.

METHODS AND RESULTS

A total of 3112 microsatellites were identified from 61,755 unique reads though 454 pyrosequencing technology. Twenty-nine microsatellite loci were selected for PCR amplification and polymorphic analyses. Among the 29 tested markers, 17 microsatellite loci were further used for genotyping in three wild P. sikkimensis populations. The number of alleles varied from one to eight, and the observed heterozygosity ranged from 0.111 to 1.000. Ten simple sequence repeat loci could be successfully cross-amplified in two Primula species. The transferability values were 76.5% in P. florindae and 58.8% in P. alpicola, respectively.

CONCLUSIONS

These microsatellite markers will be valuable for testing the hypothesis of lineage divergence, genetic introgression, and cryptic speciation events between P. sikkimensis and its closely related taxa.

Complete plastid genome sequence of Primula sinensis (Primulaceae): structure comparison, sequence variation and evidence for accD transfer to nucleus

Species-rich genus Primula L. is a typical plant group with which to understand genetic variance between species in different levels of relationships. Chloroplast genome sequences are used to be the information resource for quantifying this difference and reconstructing evolutionary history. In this study, we reported the complete chloroplast genome sequence of Primula sinensis and compared it with other related species. This genome of chloroplast showed a typical circular quadripartite structure with 150,859 bp in sequence length consisting of 37.2% GC base. Two inverted repeated regions (25,535 bp) were separated by a large single-copy region (82,064 bp) and a small single-copy region (17,725 bp). The genome consists of 112 genes, including 78 protein-coding genes, 30 tRNA genes and four rRNA genes. Among them, seven coding genes, seven tRNA genes and four rRNA genes have two copies due to their locations in the IR regions. The accD and infA genes lacking intact open reading frames (ORF) were identified as pseudogenes. SSR and sequence variation analyses were also performed on the plastome of Primula sinensis, comparing with another available plastome of P. poissonii. The four most variable regions, rpl36–rps8, rps16–trnQ, trnH–psbA and ndhC–trnV, were identified. Phylogenetic relationship estimates using three sub-datasets extracted from a matrix of 57 protein-coding gene sequences showed the identical result that was consistent with previous studies. A transcript found from P. sinensis transcriptome showed a high similarity to plastid accD functional region and was identified as a putative plastid transit peptide at the N-terminal region. The result strongly suggested that plastid accD has been functionally transferred to the nucleus in P. sinensis.

Characterization of Novel Microsatellite Loci for Primula poissonii (Primulaceae) Using High-Throughput Sequencing Technology

Primula poissonii (Primulaceae) is a perennial herb, widely distributed in the Hengduan Mountain region of Southwest China. In this study, Roche 454 pyrosequencing was used to isolate microsatellite markers. A total of 4528 unique sequences were identified from 68,070 unique reads. Of these, eighty-seven microsatellite loci were screened for utility using two criteria: successful PCR amplification and variation of these loci within three wild P. poissonii populations. Twenty loci were successfully amplified and exhibited polymorphic alleles. The number of observed alleles ranged from 1 to 9 with an average of 3.5. The observed and expected heterozygosities ranged from 0.087 to 1.000 and from 0.124 to 0.828, respectively. Among these SSR loci, only the P69 locus could not be cross-amplified successfully in two closely related species P. wilsonii and P. anisodora. The microsatellite loci developed in this study will be useful for studying genetic diversity and speciation events between P. poissonii and closely related Primula species.

The Use of DNA Barcoding on Recently Diverged Species in the Genus Gentiana (Gentianaceae) in China

DNA barcoding of plants poses particular challenges, especially in differentiating, recently diverged taxa. The genus Gentiana (Gentianaceae) is a species-rich plant group which rapidly radiated in the Himalaya-Hengduan Mountains in China. In this study, we tested the core plant barcode (rbcL + matK) and three promising complementary barcodes (trnH-psbA, ITS and ITS2) in 30 Gentiana species across 6 sections using three methods (the genetic distance-based method, Best Close Match and tree-based method). rbcL had the highest PCR efficiency and sequencing success (100%), while the lowest sequence recoverability was from ITS (68.35%). The presence of indels and inversions in trnH-psbA in Gentiana led to difficulties in sequence alignment. When using a single region for analysis, ITS exhibited the highest discriminatory power (60%-74.42%). Of the combinations, matK + ITS provided the highest discrimination success (71.43%-88.24%) and is recommended as the DNA barcode for the genus Gentiana. DNA barcoding proved effective in assigning most species to sections, though it performed poorly in some closely related species in sect. Cruciata because of hybridization events. Our analysis suggests that the status of G. pseudosquarrosa needs to be studied further. The utility of DNA barcoding was also verified in authenticating 'Qin-Jiao' Gentiana medicinal plants (G. macrophylla, G. crassicaulis, G. straminea, and G. dahurica), which can help ensure safe and correct usage of these well-known Chinese traditional medicinal herbs.

Phylotranscriptomic Analysis Based on Coalescence was Less Influenced by the Evolving Rates and the Number of Genes: A Case Study in Ericales

Advances in high-throughput sequencing have generated a vast amount of transcriptomic data that are being increasingly used in phylogenetic reconstruction. However, processing the vast datasets for a huge number of genes and even identifying optimal analytical methodology are challenging. Through de novo sequenced and retrieved data from public databases, we identified 221 orthologous protein-coding genes to reconstruct the phylogeny of Ericales, an order characterized by rapid ancient radiation. Seven species representing different families in Ericales were used as in-groups. Both concatenation and coalescence methods yielded the same well-supported topology as previous studies, with only two nodes conflicting with previously reported relationships. The results revealed that a partitioning strategy could improve the traditional concatenation methodology. Rapidly evolving genes negatively affected the concatenation analysis, while slowly evolving genes slightly affected the coalescence analysis. The coalescence methods usually accommodated rate heterogeneity better and required fewer genes to yield well-supported topologies than the concatenation methods with both real and simulated data.

Cross-Species, Amplifiable EST-SSR Markers for Amentotaxus Species Obtained by Next-Generation Sequencing

Amentotaxus, a genus of Taxaceae, is an ancient lineage with six relic and endangered species. Four Amentotaxus species, namely A. argotaenia, A. formosana, A. yunnanensis, and A. poilanei, are considered a species complex because of their morphological similarities. Small populations of these species are allopatrically distributed in Asian forests. However, only a few codominant markers have been developed and applied to study population genetic structure of these endangered species. In this study, we developed and characterized polymorphic expressed sequence tag-simple sequence repeats (EST-SSRs) from the transcriptome of A. formosana. We identified 4955 putative EST-SSRs from 68,281 unigenes as potential molecular markers. Twenty-six EST-SSRs were selected for estimating polymorphism and transferability among Amentotaxus species, of which 23 EST-SSRs were polymorphic within Amentotaxus species. Among these, the number of alleles ranged from 1-4, the polymorphism information content ranged from 0.000-0.692, and the observed and expected heterozygosity were 0.000-1.000 and 0.080-0.740, respectively. Population genetic structure analyses confirmed that A. argotaenia and A. formosana were separate species and A. yunnanensis and A. poilanei were the same species. These novel EST-SSRs can facilitate further population genetic structure research of Amentotaxus species.

DNA barcoding evaluation and its taxonomic implications in the species-rich genus Primula L. in China

The genus Primula is extremely diverse in the east Himalaya-Hengduan Mountains (HHM) in China as a result of rapid radiation. In order to overcome the difficulty of morphological classification of this genus, we surveyed three plastid regions (rbcL, matK, and trnH-psbA) and two nuclear markers (ITS and ITS2) from 227 accessions representing 66 Primula species across 18 sections, to assess their discriminatory power as barcodes. We found that ITS alone or combined with plastid regions showed the best discrimination across different infrageneric ranks and at species level. We suggest rbcL + matK + ITS as the first choice at present to barcode Primula plants. Although the present barcoding combination performed poorly in many closely related species of Primula, it still provided many new insights into current Primula taxonomy, such as the underlying presence of cryptic species, and several potential improper taxonomic treatments. DNA barcoding is one useful technique in the integrative taxonomy of the genus Primula, but it still requires further efforts to improve its effectiveness in some taxonomically challenging groups.

Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae)

Young incipient species provide ideal materials for untangling the process of ecological speciation in the presence of gene flow. The Miscanthus floridulus/sinensis complex exhibits diverse phenotypic and ecological differences despite recent divergence (approximately 1.59 million years ago). To elucidate the process of genetic differentiation during early stages of ecological speciation, we analyzed genomic divergence in the Miscanthus complex using 72 randomly selected genes from a newly assembled transcriptome. In this study, rampant gene flow was detected between species, estimated as M = 3.36 × 10(-9) to 1.20 × 10(-6) , resulting in contradicting phylogenies across loci. Nevertheless, beast analyses revealed the species identity and the effects of extrinsic cohesive forces that counteracted the non-stop introgression. As expected, early in speciation with gene flow, only 3-13 loci were highly diverged; two to five outliers (approximately 2.78-6.94% of the genome) were characterized by strong linkage disequilibrium, and asymmetrically distributed among ecotypes, indicating footprints of diversifying selection. In conclusion, ecological speciation of incipient species of Miscanthus probably followed the parapatric model, whereas allopatric speciation cannot be completely ruled out, especially between the geographically isolated northern and southern M. sinensis, for which no significant gene flow across oceanic barriers was detected. Divergence between local ecotypes in early-stage speciation began at a few genomic regions under the influence of natural selection and divergence hitchhiking that overcame gene flow.

Origins and domestication of cultivated banana inferred from chloroplast and nuclear genes

BACKGROUND

Cultivated bananas are large, vegetatively-propagated members of the genus Musa. More than 1,000 cultivars are grown worldwide and they are major economic and food resources in numerous developing countries. It has been suggested that cultivated bananas originated from the islands of Southeast Asia (ISEA) and have been developed through complex geodomestication pathways. However, the maternal and parental donors of most cultivars are unknown, and the pattern of nucleotide diversity in domesticated banana has not been fully resolved.

METHODOLOGY/PRINCIPAL FINDINGS

We studied the genetics of 16 cultivated and 18 wild Musa accessions using two single-copy nuclear (granule-bound starch synthase I, GBSS I, also known as Waxy, and alcohol dehydrogenase 1, Adh1) and two chloroplast (maturase K, matK, and the trnL-F gene cluster) genes. The results of phylogenetic analyses showed that all A-genome haplotypes of cultivated bananas were grouped together with those of ISEA subspecies of M. acuminata (A-genome). Similarly, the B- and S-genome haplotypes of cultivated bananas clustered with the wild species M. balbisiana (B-genome) and M. schizocarpa (S-genome), respectively. Notably, it has been shown that distinct haplotypes of each cultivar (A-genome group) were nested together to different ISEA subspecies M. acuminata. Analyses of nucleotide polymorphism in the Waxy and Adh1 genes revealed that, in comparison to the wild relatives, cultivated banana exhibited slightly lower nucleotide diversity both across all sites and specifically at silent sites. However, dramatically reduced nucleotide diversity was found at nonsynonymous sites for cultivated bananas.

CONCLUSIONS/SIGNIFICANCE

Our study not only confirmed the origin of cultivated banana as arising from multiple intra- and inter-specific hybridization events, but also showed that cultivated banana may have not suffered a severe genetic bottleneck during the domestication process. Importantly, our findings suggested that multiple maternal origins and a reduction in nucleotide diversity at nonsynonymous sites are general attributes of cultivated bananas.

Isolation and characterization of microsatellite markers for Cycas debaoensis Y. C. Zhong et C. J. Chen (Cycadaceae)

Cycas debaoensis is one of the critically endangered cycad species endemic to China. In this study, we described the development of six microsatellite markers from the genome of C. debaoensis using the protocol of fast isolation by AFLP of sequences containing repeats (FIASCO) and two microsatellite markers derived from the database of expressed sequence tags (dbEST). Polymorphism of each locus was assessed in 60 adult individuals of the cycad. The average allele number of the microsatellites was 2.6 per locus, ranging from two to five. The observed and expected heterozygosities varied from 0.0833 to 0.7333 and from 0.0805 to 0.7188, respectively. Despite its rarity, only one locus (Y177) deviated from Hardy-Weinberg equilibrium due to the excessive homozygosity. The marker transferability of the eight primer pairs was tested on other four congeneric species that also occur in China.

Comparative transcriptome analysis and marker development of two closely related Primrose species (Primula poissonii and Primula wilsonii)

Background

Primula species are important early spring garden plants with a centre of diversity and speciation in the East Himalaya-Hengduan Mountains in Western China. Studies on population genetics, speciation and phylogeny of Primula have been impeded by a lack of genomic resources. In the present study, we sequenced the transcriptomes of two closely related primrose species, Primula poissonii and Primula wilsonii, using short reads on the Illumina Genome Analyzer platform.

Results

We obtained 55,284 and 55,011 contigs with N50 values of 938 and 1,085 for P. poissonii and P. wilsonii, respectively, and 6,654 pairs of putative orthologs were identified between the two species. Estimations of non-synonymous/synonymous substitution rate ratios for these orthologs indicated that 877 of the pairs may be under positive selection (Ka/Ks > 0.5), and functional enrichment analysis revealed that significant proportions of the orthologs were in the categories DNA repair, stress resistance, which may provide some hints as to how the two closely related Primula species adapted differentially to extreme environments, such as habitats characterized by aridity, high altitude and high levels of ionizing radiation. It was possible for the first time to estimate the divergence time between the radiated species pair, P. poissonii and P. wilsonii; this was found to be approximately 0.90 ± 0.57 Mya, which falls between the Donau and Gunz glaciation in the Middle Pleistocene. Primers based on 54 pairs of orthologous SSR-containing sequences between the two Primula species were designed and verified. About half of these pairs successfully amplified for both species. Of the 959 single copy nuclear genes shared by four model plants (known as APVO genes), 111 single copy nuclear genes were verified as being present in both Primula species and exon-anchored and intron-spanned primers were designed for use.

Conclusion

We characterized the transcriptomes for the two Primula species, and produced an unprecedented amount of genomic resources for these important garden plants. Evolutionary analysis of these two Primula species not only revealed a more precise divergence time, but also provided some novel insights into how differential adaptations occurred in extreme habitats. Furthermore, we developed two sets of genetic markers, single copy nuclear genes and nuclear microsatellites (EST-SSR). Both these sets of markers will facilitate studies on the genetic improvement, population genetics and phylogenetics of this rapidly adapting taxon.

Development and characterization of nine microsatellites for an endangered tree, Pinus wangii (Pinaceae)

PREMISE OF THE STUDY

Pinus wangii is an endemic and endangered species in southwestern China, and microsatellite primers were developed to characterize its genetic diversity and population structure. •

METHODS AND RESULTS

Using the Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO) protocol, nine sets of microsatellite primers were developed in P. wangii. One population with 26 individuals of P. wangii, as well as 11 individuals each for two congeneric species, P. taiwanensis and P. squamata, were used to test their polymorphism and transferability. The number of alleles per locus ranged from one to seven with an average of 3.7, and the observed heterozygosity and expected heterozygosity ranged from 0 to 0.91 and 0 to 0.75, respectively. •

CONCLUSIONS

We developed nine sets of polymorphic microsatellite loci that are suitable for investigating genetic diversity and population structure of P. wangii, and these markers may be useful for other Pinus species.