A chromosome-level genome assembly of the Knoxia roxburghii (Rubiaceae)

Knoxia roxburghii is a well-known medicinal plant that is widely distributed in southern China and Southeast Asia. Its dried roots, known as hongdaji in traditional Chinese medicine, are used to treat a range of diseases, including cancers, carbuncles, and ascites. In this study, we report a de novo chromosome-level genome sequence for this diploid plant, which has a length of approximately 446.30 Mb with a contig N50 size of 42.26 Mb and scaffold N50 size of 44.38 Mb. Approximately 99.78% of the assembled sequences were anchored to 10 pseudochromosomes and 3 gapless assembled chromosomes were included in this assembly. A total of 24,507 genes were annotated, along with 68.92% of repetitive elements. Overall, our results will facilitate further active component biosynthesis for K. roxburghii and provide insights for future functional genomic studies and DNA-informed breeding.

Allelic frequencies of 22 short tandem repeats loci and tri-allelic patterns of Penta D and TPOX identified in Gabonese population

Short tandem repeats (STRs) are repeating DNA sequences used in forensic human identity testing and the diagnosis of aneuploidies. Many STRs like Penta D and TPOX are used routinely for paternity tests, but these tests are not widely used in sub-Saharan Africa. In this study we recruited individuals from Gabonese families seeking a paternity test. After DNA extraction from buccal swabs, we genotyped samples using a panel of 22 STRs. A total of 115 unrelated subjects from 39 families were included. Allele frequencies of the 22 STR loci were determined in unrelated Gabonese subjects. The most polymorphic loci were D21S11 (16 alleles) and FGA (17 alleles), while D3S1358 and TH01 loci were less polymorphic, with five alleles each. Deviation from Hardy-Weinberg equilibrium was observed for TPOX, D3S1358, CSFPO and D7S820 loci. We reported tri-allelic patterns that indicate aneuploidies at a combined frequency of 4% (4/115) with 3% for Penta D (1/35) and 3% for TPOX (3/102). Furthermore, we identified a new tri-allelic genotype 5-8-16 for the Penta D locus located on chromosome 21 in a healthy subject. In addition, we observed three tri-allelic variants of TPOX, located on chromosome 2, in healthy subjects, namely 8-10-11, 8-9-10, and 8-8-10. Our study revealed unsuspected polymorphic variations in Penta D and TPOX for the first time in Gabon, raising several questions about chromosomal disorders. Further population genetics studies are needed in Gabon to better characterize these variations, both qualitatively and quantitative.

Metagenomics of African Empogona and Tricalysia (Rubiaceae) reveals the presence of leaf endophytes

BACKGROUND

Leaf symbiosis is a phenomenon in which host plants of Rubiaceae interact with bacterial endophytes within their leaves. To date, it has been found in around 650 species belonging to eight genera in four tribes; however, the true extent in Rubiaceae remains unknown. Our aim is to investigate the possible occurrence of leaf endophytes in the African plant genera Empogona and Tricalysia and, if present, to establish their identity.

METHODS

Total DNA was extracted from the leaves of four species of the Coffeeae tribe (Empogona congesta, Tricalysia hensii, T. lasiodelphys, and T. semidecidua) and sequenced. Bacterial reads were filtered out and assembled. Phylogenetic analysis of the endophytes was used to reveal their identity and their relationship with known symbionts.

RESULTS

All four species have non-nodulated leaf endophytes, which are identified as Caballeronia. The endophytes are distinct from each other but related to other nodulated and non-nodulated endophytes. An apparent phylogenetic or geographic pattern appears to be absent in endophytes or host plants. Caballeronia endophytes are present in the leaves of Empogona and Tricalysia, two genera not previously implicated in leaf symbiosis. This interaction is likely to be more widespread, and future discoveries are inevitable.

The Complete Chloroplast Genome of An Ophiorrhiza baviensis Drake Species Reveals Its Molecular Structure, Comparative, and Phylogenetic Relationships

Ophiorrhiza baviensis Drake, a flowering medical plant in the Rubiaceae, exists uncertainly within the Ophiorrhiza genus' evolutionary relationships. For the first time, the whole chloroplast (cp) genome of an O. baviensis Drake species was sequenced and annotated. Our findings demonstrate that the complete cp genome of O. baviensis is 154,770 bp in size, encoding a total of 128 genes, including 87 protein-coding genes, 8 rRNAs, and 33 tRNAs. A total of 59 SSRs were screened in the studied cp genome, along with six highly variable loci, which can be applied to generate significant molecular markers for the Ophiorrhiza genus. The comparative analysis of the O. baviensis cp genome with two published others of the Ophiorrhiza genus revealed a high similarity; however, there were some notable gene rearrangements in the O. densa plastome. The maximum likelihood phylogenetic trees were constructed based on the concatenation of the rps16 gene and the trnL-trnF intergenic spacer sequence, indicating a close relationship between the studied O. baviensis and other Ophiorrhiza. This study will provide a theoretical molecular basis for identifying O. baviensis Drake, as well as species of the Ophiorrhiza genus, and contribute to shedding light on the chloroplast genome evolution of Rubiaceae.

OpNAC1 transcription factor regulates the biosynthesis of the anticancer drug camptothecin by targeting loganic acid O-methyltransferase in Ophiorrhiza pumila

Camptothecin (CPT) is an anticancer pentacyclic quinoline alkaloid widely used to treat cancer patients worldwide. However, the biosynthetic pathway and transcriptional regulation of camptothecin are largely unknown. Ophiorrhiza pumila, the herbaceous plant from the Rubiaceae family, has emerged as a model plant for studying camptothecin biosynthesis and regulation. In this study, a high-quality reference genome of O. pumila with estimated size of ~456.90 Mb was reported, and the accumulation level of camptothecin in roots was higher than that in stems and leaves. Based on its spatial distribution in the plant, we examined gene functions and expression by combining genomics with transcriptomic analysis. Two loganic acid O-methyltransferase (OpLAMTs) were identified in strictosidine-producing plant O. pumila, and enzyme catalysis assays showed that OpLAMT1 and not OpLAMT2 could convert loganic acid into loganin. Further knock-out of OpLAMT1 expression led to the elimination of loganin and camptothecin accumulation in O. pumila hairy roots. Four key residues were identified in OpLAMT1 protein crucial for the catalytic activity of loganic acid to loganin. By co-expression network, we identified a NAC transcription factor, OpNAC1, as a candidate gene for regulating camptothecin biosynthesis. Transgenic hairy roots and biochemical assays demonstrated that OpNAC1 suppressed OpLAMT1 expression. Here, we reported on two camptothecin metabolic engineering strategies paving the road for industrial-scale production of camptothecin in CPT-producing plants.

Phylogeny of the Madagascar-centred tribe Danaideae (Rubiaceae) as a precursor to taxonomic revision: insights into its generic and species limits, affinities and distribution

BACKGROUND AND AIMS

The tribe Danaideae (Rubiaceae) is almost exclusively endemic to the Western Indian Ocean Region (WIOR), and encompasses the genera Danais, Payera and Schismatoclada that occur in humid, sub-humid and mountain and mountain bio climate zones. Much of the species diversity is endemic to restricted, remote and/or mountainous areas of Madagascar and recent field work on the island indicates substantial unknown diversity of the Danaideae. Furthermore, the monophyly of the Malagasy genera Payera and Schismatoclada has been questioned in previous work, species delimitations and phylogenetic relationships within the genera are poorly understood, and the distribution and evolution of gross morphological features have not been assessed.

METHODS

We conducted morphological investigations, and produced robust phylogenies of Danaideae based on nuclear and plastid sequence data from 193 terminals. Ample plant material has been newly collected in the WIOR for the purpose of the present study, including potentially new species unknown to science. We performed Bayesian non-clock and relaxed-clock analyses employing three alternative clock models of a dataset with a dense sample of taxa from the entire geographical ranges of Danaideae. Based on the results, we discuss species diversity and distribution, relationships, and morphology in Danaideae.

KEY RESULTS

Our results demonstrate the monophyly of Danaideae, its three genera and 42 species. Nine species are resolved as non-monophyletic. Many geographically distinct but morphologically heterogeneous lineages were identified, and morphological features traditionally considered diagnostic of subgroups of the genera, used for example in species identification keys, are not clade-specific.

CONCLUSIONS

Our results demonstrate that Madagascar contains ample previously undocumented morphological and species diversity of Danaideae. Our novel approach to molecular phylogenetic analyses as a precursor to taxonomic revisions provides numerous benefits for the latter. There are tentative indications of parallel northward diversification in Payera and Schismatoclada on Madagascar, and of geographical phylogenetic clustering despite the anemochorous condition of Danaideae.

Physiological and Transcriptomic Responses of Growth in Neolamarckia cadamba Stimulated by Exogenous Gibberellins

(1) The phytohormones gibberellins (GAs) play a crucial role in plant growth and development, such as seed germination, flowering, fruiting, and stem elongation. Although many biological roles of GAs have been studied intensively, the molecular mechanisms of GAs in woody plants are still unclear. (2) In this study, we investigated the effects of exogenous application of GAs on Neolamarckia cadamba. (3) The height and biomass of N. cadamba increased after 7 days of GA treatment, especially on the second internode. Transcriptome analysis showed that although the majority of genes involved in the GA signaling pathway were up-regulated, the expression of GA20 oxidase (GA20ox) and GA3 oxidase (GA3ox) was down-regulated in the 3 days GA-treated group compared to the CK group. The expression of the cell elongation-related basic helix-loop-helix genes bHLH74 and bHLH49 was up-regulated in the GA-treated group compared with the CK group. Transcriptional expression levels of transcription factors involved in hormone signaling were changed, mainly including bHLH, ethylene response factor (ERF), and WRKY families. In addition, the transcriptional expression level of the key enzymes engaged in the phenylalanine pathway was downregulated after GA treatment. (4) In brief, our findings reveal the physiological and molecular mechanisms of exogenous GA treatment stimulation in N. cadamba.

The Bacterial and Fungal Microbiota of the Mexican Rubiaceae Family Medicinal Plant Bouvardia ternifolia

Bouvardia ternifolia is a medicinal plant considered a source of therapeutic compounds, like the antitumoral cyclohexapeptide bouvardin. It is known that large number of secondary metabolites produced by plants results from the interaction of the host and adjacent or embedded microorganisms. Using high-throughput DNA sequencing of V3-16S and V5-18S ribosomal gene libraries, we characterized the endophytic, endophytic + epiphyte bacterial, and fungal communities associated to flowers, leaves, stems, and roots, as well as the rhizosphere. The Proteobacteria (average 80.7%) and Actinobacteria (average 14.7%) were the most abundant bacterial phyla, while Leotiomycetes (average 54.8%) and Dothideomycetes (average 27.4%) were the most abundant fungal classes. Differential abundance for the bacterial endophyte group showed a predominance of Erwinia, Propionibacterium, and Microbacterium genera, while Sclerotinia, Coccomyces, and Calycina genera predominated for fungi. The predictive metagenome analysis for bacteria showed significative abundance of pathways for secondary metabolite production, while a FUNguild analysis revealed the presence of pathotroph, symbiotroph, and saprotrophs in the fungal community. Intra and inter copresence and mutual exclusion interactions were identified for bacterial and fungal kingdoms in the endophyte communities. This work provides a description of the diversity and composition of bacterial and fungal microorganisms living in flowers, leaves, stems, roots, and the rhizosphere of this medicinal plant; thus, it paves the way towards an integral understanding in the production of therapeutic metabolites.

Phylogenomic and comparative analyses of Coffeeae alliance (Rubiaceae): deep insights into phylogenetic relationships and plastome evolution

BACKGROUND

The large and diverse Coffeeae alliance clade of subfamily Ixoroideae (Rubiaceae) consists of 10 tribes, > 90 genera, and > 2000 species. Previous molecular phylogenetics using limited numbers of markers were often unable to fully resolve the phylogenetic relationships at tribal and generic levels. Also, the structural variations of plastomes (PSVs) within the Coffeeae alliance tribes have been poorly investigated in previous studies. To fully understand the phylogenetic relationships and PSVs within the clade, highly reliable and sufficient sampling with superior next-generation analysis techniques is required. In this study, 71 plastomes (40 newly sequenced and assembled and the rest from the GenBank) were comparatively analyzed to decipher the PSVs and resolve the phylogenetic relationships of the Coffeeae alliance using four molecular data matrices.

RESULTS

All plastomes are typically quadripartite with the size ranging from 153,055 to 155,908 bp and contained 111 unique genes. The inverted repeat (IR) regions experienced multiple contraction and expansion; five repeat types were detected but the most abundant was SSR. The size of the Coffeeae alliance clade plastomes and its elements are affected by the IR boundary shifts and the repeat types. However, the emerging PSVs had no taxonomic and phylogenetic implications. Eight highly divergent regions were identified within the plastome regions ndhF, ccsA, ndhD, ndhA, ndhH, ycf1, rps16-trnQ-UUG, and psbM-trnD. These highly variable regions may be potential molecular markers for further species delimitation and population genetic analyses for the clade. Our plastome phylogenomic analyses yielded a well-resolved phylogeny tree with well-support at the tribal and generic levels within the Coffeeae alliance.

CONCLUSIONS

Plastome data could be indispensable in resolving the phylogenetic relationships of the Coffeeae alliance tribes. Therefore, this study provides deep insights into the PSVs and phylogenetic relationships of the Coffeeae alliance and the Rubiaceae family as a whole.

Chromosome-level assembly of the Neolamarckia cadamba genome provides insights into the evolution of cadambine biosynthesis

Neolamarckia cadamba (Roxb.), a close relative of Coffea canephora and Ophiorrhiza pumila, is an important traditional medicine in Southeast Asia. Three major glycosidic monoterpenoid indole alkaloids (MIAs), cadambine and its derivatives 3β-isodihydrocadambine and 3β-dihydrocadambine, accumulate in the bark and leaves, and exhibit antimalarial, antiproliferative, antioxidant, anticancer and anti-inflammatory activities. Here, we report a chromosome-scale N. cadamba genome, with 744.5 Mb assembled into 22 pseudochromosomes with contig N50 and scaffold N50 of 824.14 Kb and 29.20 Mb, respectively. Comparative genomic analysis of N. cadamba with Co. canephora revealed that N. cadamba underwent a relatively recent whole-genome duplication (WGD) event after diverging from Co. canephora, which contributed to the evolution of the MIA biosynthetic pathway. We determined the key intermediates of the cadambine biosynthetic pathway and further showed that NcSTR1 catalyzed the synthesis of strictosidine in N. cadamba. A new component, epoxystrictosidine (C27H34N2O10, m/z 547.2285), was identified in the cadambine biosynthetic pathway. Combining genome-wide association study (GWAS), population analysis, multi-omics analysis and metabolic gene cluster prediction, this study will shed light on the evolution of MIA biosynthetic pathway genes. This N. cadamba reference sequence will accelerate the understanding of the evolutionary history of specific metabolic pathways and facilitate the development of tools for enhancing bioactive productivity by metabolic engineering in microbes or by molecular breeding in plants.

Unlocking the Complete Chloroplast Genome of a Native Tree Species from the Amazon Basin, Capirona (Calycophyllum Spruceanum, Rubiaceae), and Its Comparative Analysis with Other Ixoroideae Species

Capirona (Calycophyllum spruceanum Benth.) belongs to subfamily Ixoroideae, one of the major lineages in the Rubiaceae family, and is an important timber tree. It originated in the Amazon Basin and has widespread distribution in Bolivia, Peru, Colombia, and Brazil. In this study, we obtained the first complete chloroplast (cp) genome of capirona from the department of Madre de Dios located in the Peruvian Amazon. High-quality genomic DNA was used to construct libraries. Pair-end clean reads were obtained by PE 150 library and the Illumina HiSeq 2500 platform. The complete cp genome of C. spruceanum has a 154,480 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (84,813 bp) and a small single-copy (SSC) region (18,101 bp), separated by two inverted repeat (IR) regions (25,783 bp). The annotation of C. spruceanum cp genome predicted 87 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and one pseudogene. A total of 41 simple sequence repeats (SSR) of this cp genome were divided into mononucleotides (29), dinucleotides (5), trinucleotides (3), and tetranucleotides (4). Most of these repeats were distributed in the noncoding regions. Whole chloroplast genome comparison with the other six Ixoroideae species revealed that the small single copy and large single copy regions showed more divergence than inverted regions. Finally, phylogenetic analyses resolved that C. spruceanum is a sister species to Emmenopterys henryi and confirms its position within the subfamily Ixoroideae. This study reports for the first time the genome organization, gene content, and structural features of the chloroplast genome of C. spruceanum, providing valuable information for genetic and evolutionary studies in the genus Calycophyllum and beyond.

Comparative anatomy and genetic bases of fruit development in selected Rubiaceae (Gentianales)

PREMISE

The Rubiaceae are ideal for studying the diversity of fruits that develop from flowers with inferior ovary. We aimed to identify morpho-anatomical changes during fruit development that distinguish those derived from the carpel versus the extra-carpellary tissues. In addition, we present the fruit genetic core regulatory network in selected Rubiaceae species and compare it in terms of copy number and expression patterns to model core eudicots in the Brassicaceae and the Solanaceae.

METHODS

We used light microscopy to follow morphoanatomical changes in four selected species with different fruit types. We generated reference transcriptomes for seven selected Rubiaceae species and isolated homologs of major transcription factors involved in fruit development histogenesis, assessed their homology, identified conserved and new protein motifs, and evaluated their expression in three species with different fruit types.

RESULTS

Our studies revealed ovary-derived pericarp tissues versus floral-cup-derived epicarp tissues. Gene evolution analyses of FRUITFULL, SHATTERPROOF, ALCATRAZ, INDEHISCENT and REPLUMLESS homologs suggest that the gene complement in Rubiaceae is simpler compared to that in Brassicaceae or Solanaceae. Expression patterns of targeted genes vary in response to the fruit type and the developmental stage evaluated.

CONCLUSIONS

Morphologically similar fruits can have different anatomies as a result of convergent tissues developed from the epicarps covering the anatomical changes from the pericarps. Expression analyses suggest that the fruit patterning regulatory network established in model core eudicots cannot be extrapolated to asterids with inferior ovaries.

Conservation biology of threatened Mediterranean chasmophytes: The case of Asperula naufraga endemic to Zakynthos island (Ionian islands, Greece)

Asperula naufraga is a rare and threatened obligate chasmophyte, endemic to Zakynthos island (Ionian islands, Greece). In this study, we provide a combined approach (including monitoring of demographic and reproductive parameters and study of genetic diversity) to assess the current conservation status of the species and to estimate its future extinction risk. The five subpopulations of A. naufraga were monitored for five years (2014-2018). Population size markedly fluctuated between 68-130 mature individuals during the monitoring period. The extent of occurrence (EOO) was estimated at 28.7 km2 and the area of occupancy (AOO) was 8 km2. Stage-structure recordings were similar for all subpopulations, characterized by high proportions of adult and senescent individuals, following a common pattern, which has been observed in other cliff-dwelling plants. Preliminary genetic analysis with SSRs markers revealed low heterozygosity within subpopulations and significant departure from H-W equilibrium, which combined with small population size suggest increased threat of genetic diversity loss. Our results indicate that the species should be placed in the Critically Endangered (CR) IUCN threat category, while according to Population Viability Analysis results its extinction risk increases to 47.8% in the next 50 years. The small population size combined with large fluctuations in its size, low recruitment and low genetic diversity, indicate the need of undertaking effective in situ and ex situ conservation measures.

Physiological, Biochemical, and Transcriptomic Responses of Neolamarckia cadamba to Aluminum Stress

Aluminum is the most abundant metal of the Earth's crust accounting for 7% of its mass, and release of toxic Al3+ in acid soils restricts plant growth. Neolamarckia cadamba, a fast-growing tree, only grows in tropical regions with acidic soils. In this study, N. cadamba was treated with high concentrations of aluminum under acidic condition (pH 4.5) to study its physiological, biochemical, and molecular response mechanisms against high aluminum stress. High aluminum concentration resulted in significant inhibition of root growth with time in N. cadamba. The concentration of Al3+ ions in the root tip increased significantly and the distribution of absorbed Al3+ was observed in the root tip after Al stress. Meanwhile, the concentration of Ca, Mg, Mn, and Fe was significantly decreased, but P concentration increased. Aluminum stress increased activities of antioxidant enzymes such as superoxide dismutase (SOD), catalase from micrococcus lysodeiktic (CAT), and peroxidase (POD) in the root tip, while the content of MDA was decreased. Transcriptome analysis showed 37,478 differential expression genes (DEGs) and 4096 GOs terms significantly associated with treatments. The expression of genes regulating aluminum transport and abscisic acid synthesis was significantly upregulated; however, the genes involved in auxin synthesis were downregulated. Of note, the transcripts of several key enzymes affecting lignin monomer synthesis in phenylalanine pathway were upregulated. Our results shed light on the physiological and molecular mechanisms of aluminum stress tolerance in N. cadamba.

Genome assembly of Chiococca alba uncovers key enzymes involved in the biosynthesis of unusual terpenoids

Chiococca alba (L.) Hitchc. (snowberry), a member of the Rubiaceae, has been used as a folk remedy for a range of health issues including inflammation and rheumatism and produces a wealth of specialized metabolites including terpenes, alkaloids, and flavonoids. We generated a 558 Mb draft genome assembly for snowberry which encodes 28,707 high-confidence genes. Comparative analyses with other angiosperm genomes revealed enrichment in snowberry of lineage-specific genes involved in specialized metabolism. Synteny between snowberry and Coffea canephora Pierre ex A. Froehner (coffee) was evident, including the chromosomal region encoding caffeine biosynthesis in coffee, albeit syntelogs of N-methyltransferase were absent in snowberry. A total of 27 putative terpene synthase genes were identified, including 10 that encode diterpene synthases. Functional validation of a subset of putative terpene synthases revealed that combinations of diterpene synthases yielded access to products of both general and specialized metabolism. Specifically, we identified plausible intermediates in the biosynthesis of merilactone and ribenone, structurally unique antimicrobial diterpene natural products. Access to the C. alba genome will enable additional characterization of biosynthetic pathways responsible for health-promoting compounds in this medicinal species.

Chloroplast genomes of Rubiaceae: Comparative genomics and molecular phylogeny in subfamily Ixoroideae

In Rubiaceae phylogenetics, the number of markers often proved a limitation with authors failing to provide well-supported trees at tribal and generic levels. A robust phylogeny is a prerequisite to study the evolutionary patterns of traits at different taxonomic levels. Advances in next-generation sequencing technologies have revolutionized biology by providing, at reduced cost, huge amounts of data for an increased number of species. Due to their highly conserved structure, generally recombination-free, and mostly uniparental inheritance, chloroplast DNA sequences have long been used as choice markers for plant phylogeny reconstruction. The main objectives of this study are: 1) to gain insight in chloroplast genome evolution in the Rubiaceae (Ixoroideae) through efficient methodology for de novo assembly of plastid genomes; and, 2) to test the efficiency of mining SNPs in the nuclear genome of Ixoroideae based on the use of a coffee reference genome to produce well-supported nuclear trees. We assembled whole chloroplast genome sequences for 27 species of the Rubiaceae subfamily Ixoroideae using next-generation sequences. Analysis of the plastid genome structure reveals a relatively good conservation of gene content and order. Generally, low variation was observed between taxa in the boundary regions with the exception of the inverted repeat at both the large and short single copy junctions for some taxa. An average of 79% of the SNP determined in the Coffea genus are transferable to Ixoroideae, with variation ranging from 35% to 96%. In general, the plastid and the nuclear genome phylogenies are congruent with each other. They are well-resolved with well-supported branches. Generally, the tribes form well-identified clades but the tribe Sherbournieae is shown to be polyphyletic. The results are discussed relative to the methodology used and the chloroplast genome features in Rubiaceae and compared to previous Rubiaceae phylogenies.

Phylogeny, character evolution and spatiotemporal diversification of the species-rich and world-wide distributed tribe Rubieae (Rubiaceae)

The Rubiaceae tribe Rubieae has a world-wide distribution with up to 1,000 species. These collectively exhibit an enormous ecological and morphological diversity, making Rubieae an excellent group for macro- and microevolutionary studies. Previous molecular phylogenetic analyses used only a limited sampling within the tribe or missed lineages crucial for understanding character evolution in this group. Here, we analyze sequences from two plastid spacer regions as well as morphological and biogeographic data from an extensive and evenly distributed sampling to establish a sound phylogenetic framework. This framework serves as a basis for our investigation of the evolution of important morphological characters and the biogeographic history of the Rubieae. The tribe includes three major clades, the Kelloggiinae Clade (Kelloggia), the Rubiinae Clade (Didymaea, Rubia) and the most species-rich Galiinae Clade (Asperula, Callipeltis, Crucianella, Cruciata, Galium, Mericarpaea, Phuopsis, Sherardia, Valantia). Within the Galiinae Clade, the largest genera Galium and Asperula are para- and polyphyletic, respectively. Smaller clades, however, usually correspond to currently recognized taxa (small genera or sections within genera), which may be used as starting points for a refined classification in this clade. Life-form (perennial versus annual), flower shape (long versus short corolla tube) and fruit characters (dry versus fleshy, with or without uncinate hairs) are highly homoplasious and have changed multiple times independently. Inference on the evolution of leaf whorls, a characteristic feature of the tribe, is sensitive to model choice. Multi-parted leaf whorls appear to have originated from opposite leaves with two small interpetiolar stipules that are subsequently enlarged and increased in number. Early diversification of Rubieae probably started during the Miocene in western Eurasia. Disjunctions between the Old and the New World possibly are due to connections via a North Atlantic land bridge. Diversification of the Galiineae Clade started later in the Miocene, probably in the Mediterranean, from where lineages reached, often multiple times, Africa, eastern Asia and further on the Americas and Australia.

Comparative Chloroplast Genome Analyses of Species in Gentiana section Cruciata (Gentianaceae) and the Development of Authentication Markers

Gentiana section Cruciata is widely distributed across Eurasia at high altitudes, and some species in this section are used as traditional Chinese medicine. Accurate identification of these species is important for their utilization and conservation. Due to similar morphological and chemical characteristics, correct discrimination of these species still remains problematic. Here, we sequenced three complete chloroplast (cp) genomes (G. dahurica, G. siphonantha and G. officinalis). We further compared them with the previously published plastomes from sect. Cruciata and developed highly polymorphic molecular markers for species authentication. The eight cp genomes shared the highly conserved structure and contained 112 unique genes arranged in the same order, including 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. We analyzed the repeats and nucleotide substitutions in these plastomes and detected several highly variable regions. We found that four genes (accD, clpP, matK and ycf1) were subject to positive selection, and sixteen InDel-variable loci with high discriminatory powers were selected as candidate barcodes. Our phylogenetic analyses based on plastomes further confirmed the monophyly of sect. Cruciata and primarily elucidated the phylogeny of Gentianales. This study indicated that cp genomes can provide more integrated information for better elucidating the phylogenetic pattern and improving discriminatory power during species authentication.

Kupeantha (Coffeeae, Rubiaceae), a new genus from Cameroon and Equatorial Guinea

Two new coffee relatives (tribe Coffeeae, Rubiaceae), discovered during botanical expeditions to Cameroon, are examined for generic placement, and the placement of three previously known species (Argocoffeopsis fosimondi, A. spathulata and Calycosiphonia pentamera) is reinvestigated using plastid sequence (accD-psa1, rpl16, trnL-F) and morphological data. Seed biochemistry of the new species and pollen micromorphology (only one of the two species) are also studied. Based on the plastid sequence data, the new taxa are nested in a well-supported monophyletic group that includes Argocoffeopsis and Calycosiphonia. Within this clade, three well-supported subclades are recovered that are morphologically easy to diagnose: (1) Calycosiphonia (excluding C. pentamera), (2) Argocoffeopsis (excluding A. fosimondi and A. spathulata), and (3) a clade including the above excluded species, in addition to the new species. Based on the results, Kupeantha, a new genus of five species, is described, including two new Critically Endangered taxa from the Highlands of Cameroon: Kupeantha ebo and K. kupensis. Phytochemical analysis of Kupeantha seeds reveals compounds assigned as hydroxycinnamic acid derivatives, amino acids and ent-kaurane diterpenoids; caffeine was not detected. Kupeantha is the first new genus described in tribe Coffeeae in 40 years.

Selection and Validation of Reference Genes for mRNA Expression by Quantitative Real-Time PCR Analysis in Neolamarckia cadamba

Neolamarckia cadamba is an economically-important fast-growing tree species in South China and Southeast Asia. As a prerequisite first step for future gene expression studies, we have identified and characterized a series of stable reference genes that can be used as controls for quantitative real time PCR (qRT-PCR) expression analysis in this study. The expression stability of 15 candidate reference genes in various tissues and mature leaves under different conditions was evaluated using four different algorithms, i.e., geNorm, NormFinder, BestKeeper and RefFinder. Our results showed that SAMDC was the most stable of the selected reference genes across the set of all samples, mature leaves at different photosynthetic cycles and under drought stress, whereas RPL10A had the most stable expression in various tissues. PGK and RPS25 were considered the most suitable reference for mature leaves at different developmental stages and under cold treatment, respectively. Additionally, the gene expression profiles of sucrose transporter 4 (NcSUT4), and 9-cis-epoxycarotenoid dioxygenase 3 (NcNCED3) were used to confirm the validity of candidate reference genes. Collectively, our study is the first report to validate the optimal reference genes for normalization under various conditions in N. cadamba and will benefit the future discovery of gene function in this species.

Plant reference genes for development and stress response studies

Many reference genes are used by different laboratories for gene expression analyses to indicate the relative amount of input RNA/DNA in the experiment. These reference genes are supposed to show least variation among the treatments and with the control sets in a given experiment. However, expression of reference genes varies significantly from one set of experiment to the other. Thus, selection of reference genes depends on the experimental conditions. Sometimes the average expression of two or three reference genes is taken as standard. This review consolidated the details of about 120 genes attempted for normalization during comparative expression analysis in 16 different plants. Plant species included in this review are Arabidopsis thaliana, cotton (Gossypium hirsutum), tobacco (Nicotiana benthamiana and N. tabacum), soybean (Glycine max), rice (Oryza sativa), blueberry (Vaccinium corymbosum), tomato (Solanum lycopersicum), wheat (Triticum aestivum), potato (Solanum tuberosum), sugar cane (Saccharum sp.), carrot (Daucus carota), coffee (Coffea arabica), cucumber (Cucumis sativus), kiwi (Actinidia deliciosa) and grape (Vitis vinifera). The list includes model and cultivated crop plants from both monocot and dicot classes. We have categorized plant-wise the reference genes that have been used for expression analyses in any or all of the four different conditions such as biotic stress, abiotic stress, developmental stages and various organs and tissues, reported till date. This review serves as a guide during the reference gene hunt for gene expression analysis studies.

Role of the Qinghai-Tibetan Plateau uplift in the Northern Hemisphere disjunction: evidence from two herbaceous genera of Rubiaceae

To assess the role of the Qinghai-Tibetan Plateau uplift in shaping the intercontinental disjunction in Northern Hemisphere, we analyzed the origin and diversification within a geological timeframe for two relict herbaceous genera, Theligonum and Kelloggia (Rubiaceae). Phylogenetic relationships within and between Theligonum and Kelloggia as well as their relatives were inferred using five chloroplast markers with parsimony, Bayesian and maximum-likelihood approaches. Migration routes and evolution of these taxa were reconstructed using Bayesian relaxed molecular clock and ancestral area reconstruction. Our results suggest the monophyly of each Theligonum and Kelloggia. Eastern Asian and North American species of Kelloggia diverged at ca.18.52 Mya and the Mediterranean species of Theligonum diverged from eastern Asian taxa at ca.13.73 Mya. Both Kelloggia and Theligonum are Tethyan flora relicts, and their ancestors might have been occurred in warm tropical to subtropical environments along the Tethys coast. The Qinghai-Tibetan Plateau separated the eastern and western Tethyan area may contribute significantly to the disjunct distributions of Theligonum, and the North Atlantic migration appears to be the most likely pathway of expansion of Kelloggia to North America. Our results highlight the importance role of the QTP uplift together with corresponding geological and climatic events in shaping biodiversity and biogeographic distribution in the Northern Hemisphere.

Conflicting results from mitochondrial genomic data challenge current views of Rubiaceae phylogeny

PREMISE OF THE STUDY

Reconstruction of plant phylogeny has heavily relied on single-gene or multigene plastid data. New sequencing methods have led to an increasing number of studies based on data from the entire plastid, but the mitochondrion has rarely been used to infer plant phylogeny because of an assumed information poverty and demonstrated lateral transfer of mitochondrial gene regions between distantly related species.

METHODS

We explored phylogenetic information from the plant mitochondrion using 57 representatives of the species-rich coffee family as study system and assessed consistency with previous results based (mostly) on plastid data.

KEY RESULTS

We showed that the mitochondrial genome can provide structured and statistically significant information on plant phylogeny. While most of our results are consistent with those based on plastid data, some surprising and statistically significant conflicts emerge, and our study demonstrates with striking clarity that the phylogeny of Rubiaceae is far from resolved.

CONCLUSIONS

It appears unlikely that conflicts between results retrieved from the different genomic compartments would be restricted to Rubiaceae. Rather, they are probably a general phenomenon and an important factor behind longstanding "difficult" phylogenetic questions. The biological processes responsible for the conflicting results detected here are unclear, but some conflicts are likely caused by hybridization events that occurred tens of millions of years ago. Whether such ancient events can be reconstructed based on molecular data from extant plants remains to be seen, but future studies of the nuclear genome may provide a way forward.

AFLP diversity and spatial structure of Calycophyllum candidissimum (Rubiaceae), a dominant tree species of Nicaragua's critically endangered seasonally dry forest

The Central American seasonally dry tropical (SDT) forest biome is one of the worlds' most endangered ecosystems, yet little is known about the genetic consequences of its recent fragmentation. A prominent constituent of this biome is Calycophyllum candidissimum, an insect-pollinated and wind-dispersed canopy tree of high socio-economic importance, particularly in Nicaragua. Here, we surveyed amplified fragment length polymorphisms across 13 populations of this species in Nicaragua to elucidate the relative roles of contemporary vs historical factors in shaping its genetic variation. Genetic diversity was low in all investigated populations (mean HE=0.125), and negatively correlated with latitude. Overall population differentiation was moderate (ΦST=0.109, P<0.001), and Bayesian analysis of population structure revealed two major latitudinal clusters (I: 'Pacific North'+'Central Highland'; II: 'Pacific South'), along with a genetic cline between I and II. Population-based cluster analyses indicated a strong pattern of 'isolation by distance' as confirmed by Mantel's test. Our results suggest that (1) the low genetic diversity of these populations reflects biogeographic/population history (colonisation from South America, Pleistocene range contractions) rather than recent human impact; whereas (2) the underlying process of their isolation by distance pattern, which is best explained by 'isolation by dispersal limitation', implies contemporary gene flow between neighbouring populations as likely facilitated by the species' efficient seed dispersal capacity. Overall, these results underscore that even tree species from highly decimated forest regions may be genetically resilient to habitat fragmentation due to species-typical dispersal characteristics, the necessity of broad-scale measures for their conservation notwithstanding.

Historical biogeography and phylogeny of the pantropical Psychotrieae alliance (Rubiaceae), with particular emphasis on the Western Indian Ocean Region

PREMISE OF THE STUDY

The Western Indian Ocean Region (WIOR) is a biodiversity hotspot providing an ideal setting for exploring the origins of insular biodiversity and dynamics of island colonization. We aimed to investigate the origins of the WIOR Psychotrieae alliance (Rubiaceae) with typically small, probably mainly bird-dispersed drupes, and the timing and direction or sequence of its colonization events in the region.

METHODS

We used the program BEAST to estimate divergence times and Lagrange for biogeographic reconstruction.

KEY RESULTS

The alliance has reached the WIOR at least 14 times via dispersals from Africa along with Asia and the Pacific mostly during the last 10 My, with at least one back-colonization to Africa. We inferred the earliest dispersal to Madagascar from the Pacific or Asia in the Miocene and numerous out-of-Madagascar dispersals to the nearby archipelagos but no dispersal out of those archipelagos. Gynochthodes with multiple fruits reached Madagascar twice from the Pacific possibly via ocean drifting. Psychotria with dry fruits (schizocarps) colonized Madagascar from the Pacific or Asia before reaching the Comoros from Madagascar possibly via wind dispersal.

CONCLUSIONS

This study reinforces the pivotal role of dispersal in shaping the WIOR biodiversity and as the critical initiating step in the generation of endemic biodiversity on its islands. The WIOR alliance shows strong Asian and Pacific affinities despite the proximity of the region to Africa. Madagascar has served as a stepping-stone for subsequent dispersal to the rest of the region. The Afro-Malagasy-Seychelles genus Craterispermum and the Malagasy Puffia may represent relictual lineages.

Ixora (Rubiaceae) on the Philippines - crossroad or cradle?

BACKGROUND

The Philippine archipelago is globally one of the most important model island systems for studying evolutionary processes. However, most plant species on this archipelago have not yet been studied in sufficient detail. The main aim of this study is to unravel the evolutionary history and biogeographic relationships of the Philippine members of the pantropical genus Ixora.

RESULTS

The complex plastid and nuclear divergence patterns in Philippine Ixora, documented using tree and network approaches, reveal a highly dynamic evolution in Ixora, involving several phases of radiation and recolonization. Philippine Ixora comprises at least five lineages, of which one is most closely related to species from Wallacea, and the remaining four to species from Asia.

CONCLUSIONS

Our study highlights the importance of Philippine species for understanding phytogeographic patterns in the Indomalayan-Australasian eco-region. The overall genetic differentiation, as well as the incongruence between genealogies based on the biparentally inherited nucleome and the maternally inherited plastome in Ixora, reflect the complex tectonic history of the Philippine archipelago. The Ixora lineage related to Wallacean species supports the delimitation of different ecozones along Huxley's line, because it is absent from Palawan. The remaining four lineages are all allied with Asian taxa, reflecting several waves of colonization. Close relationships between some widespread Philippine species and locally adapted narrow endemics suggest that the widespread, genetically diverse species act as pools for the formation of new species in a process of ongoing speciation. Our results suggest that the species concepts of some of the more widespread taxa need to be revised.

Evolution of woody life form on tropical mountains in the tribe Spermacoceae (Rubiaceae)

PREMISE OF THE STUDY

Spermacoceae are mainly an herbaceous group in the Rubiaceae. However, a few lineages are woody and are found in a diverse range of habitat types. Three of the largest woody lineages (Arcytophyllum, Hedyotis, and Kadua) are characterized by their distribution in the moist tropical mountains and have disjunct distribution patterns with respect to their closest relatives. In this study, we explore the cases of derived woodiness in these three lineages and their diversification dynamics in the tropical mountains of Asia, the Pacific, and the Americas.

METHODS

By combining phylogenetic results with wood anatomical studies, we estimated timing of origin of the three woody groups, inferred their ancestral traits and ancestral distribution ranges, analyzed their associations with the tropical upland habitat, and elucidated their diversification across tropical mountains.

KEY RESULTS

The three woody clades originated and diversified from herbaceous ancestors in close association with the tropical upland habitat during the Miocene. The ancestral range for Asian-Pacific Hedyotis is Africa/Madagascar and continental Asia for Pacific Kadua. The complex geological history of tropical Asia allowed Hedyotis to diversify faster and create narrow endemics near oceans in the highlands of the Western Ghats in India, Sri Lanka, Southeast Asia including southeastern China, and New Guinea.

CONCLUSIONS

The three major woody clades in Spermacoceae have gained their woodiness independently from one another, subsequent to colonization by their ancestors from a different geographic environment. The evolution and diversification along the tropical mountain orogeny is strongly linked with the formation of woody habit and many narrow endemic species.

Discovery, structure, function, and applications of cyclotides: circular proteins from plants

Cyclotides are plant-derived cyclic peptides that have a head-to-tail cyclic backbone and three conserved disulphide bonds that form a cyclic cystine knot motif. They occur in plants from the Violaceae, Rubiaceae, Cucurbitaceae, Fabaceae, and Solanaceae families, typically with 10-100 cyclotides in a given plant species, in a wide range of tissues, including flowers, leaves, stems, and roots. Some cyclotides are expressed in large amounts (up to 1g kg(-1) wet plant weight) and their natural function appears to be to protect plants from pests or pathogens. This article provides a brief overview of their discovery, distribution in plants, and applications. In particular, their exceptional stability has led to their use as peptide-based scaffolds in drug design applications. They also have potential as natural 'ecofriendly' insecticides, and as protein engineering frameworks.

Transcriptomic Analysis of Multipurpose Timber Yielding Tree Neolamarckia cadamba during Xylogenesis Using RNA-Seq

Neolamarckia cadamba is a fast-growing tropical hardwood tree that is used extensively for plywood and pulp production, light furniture fabrication, building materials, and as a raw material for the preparation of certain indigenous medicines. Lack of genomic resources hampers progress in the molecular breeding and genetic improvement of this multipurpose tree species. In this study, transcriptome profiling of differentiating stems was performed to understand N. cadamba xylogenesis. The N. cadamba transcriptome was sequenced using Illumina paired-end sequencing technology. This generated 42.49 G of raw data that was then de novo assembled into 55,432 UniGenes with a mean length of 803.2bp. Approximately 47.8% of the UniGenes (26,487) were annotated against publically available protein databases, among which 21,699 and 7,754 UniGenes were assigned to Gene Ontology categories (GO) and Clusters of Orthologous Groups (COG), respectively. 5,589 UniGenes could be mapped onto 116 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Among 6,202 UniGenes exhibiting differential expression during xylogenesis, 1,634 showed significantly higher levels of expression in the basal and middle stem segments compared to the apical stem segment. These genes included NAC and MYB transcription factors related to secondary cell wall biosynthesis, genes related to most metabolic steps of lignin biosynthesis, and CesA genes involved in cellulose biosynthesis. This study lays the foundation for further screening of key genes associated with xylogenesis in N. cadamba as well as enhancing our understanding of the mechanism of xylogenesis in fast-growing trees.

DNA Barcodes Confirm the Taxonomic and Conservation Status of a Species of Tree on the Brink of Extinction in the Pacific

The taxonomic status of a single island, narrow range endemic plant species from Palau, Micronesia (Timonius salsedoi) was assessed using DNA barcode markers, additional plastid loci, and morphology in order to verify its conservation status. DNA barcode loci distinguished T. salsedoi from all other Timonius species sampled from Palau, and were supported by sequence data from the atpB-rbcL intergenic spacer region. Timonius salsedoi was only known from two mature individual trees in 2012. Due to its extremely narrow range and population size, it had previously been recommended to be listed as Critically Endangered Status under three separate IUCN Criteria. In 2014 a second survey of the population following a typhoon revealed that the only two known trees had died suggesting that this species may now be extinct. Comprehensive follow up surveys of suitable habitat for this species are urgently required.

Evolutionary Relationships and Biogeography of the Ant-Epiphytic Genus Squamellaria (Rubiaceae: Psychotrieae) and Their Taxonomic Implications

Ecological research on ant/plant symbioses in Fiji, combined with molecular phylogenetics, has brought to light four new species of Squamellaria in the subtribe Hydnophytinae of the Rubiaceae tribe Psychotrieae and revealed that four other species, previously in Hydnophytum, need to be transferred to Squamellaria. The diagnoses of the new species are based on morphological and DNA traits, with further insights from microCT scanning of flowers and leaf δ¹³C ratios (associated with Crassulacean acid metabolism). Our field and phylogenetic work results in a new circumscription of the genus Squamellaria, which now contains 12 species (to which we also provide a taxonomic key), not 3 as in the last revision. A clock-dated phylogeny and a model-testing biogeographic framework were used to infer the broader geographic history of rubiaceous ant plants in the Pacific, specifically the successive expansion of Squamellaria to Vanuatu, the Solomon Islands, and Fiji. The colonization of Vanuatu may have occurred from Fiji, when these islands were still in the same insular arc, while the colonization of the Solomon islands may have occurred after the separation of this island from the Fiji/Vanuatu arc. Some of these ant-housing epiphytes must have dispersed with their specialized ants, for instance attached to floating timber. Others acquired new ant symbionts on different islands.

Reproductive isolation between sympatric sister species, Mussaenda kwangtungensis and M. pubescens var. alba

Reproductive isolation defines the biological species concept and plays a key role in the formation and maintenance of species. The relative contributions of different isolating stages has been suggested to be closely associated with phylogenetic relatedness. Few studies have focused on the relative contributions of pre- versus post-zygotic mechanisms, and even fewer have been conducted under strict phylogenetic frameworks. Pre- and post-zygotic reproductive isolation stages have been investigated in the sister species Mussaenda kwangtungensis and M. pubescens var. alba. The two species have partly overlapping distribution ranges and flowering times, while the principal pollinators differed strikingly for them, demonstrating strong pre-zygotic isolations. Natural hybrids were detected by simple sequence repeat markers and their maternal parents were identified based on chloroplast gene sequences. Five out of 81 individuals were suggested to be hybrids that fall into the categories F2, BC1, and BC2 by the NewHybrids analysis. Interspecific crossings resulted in significantly reduced fruit set and seed germination rates. Phylogenetic analysis revealed short Kimura-2-parameter distance between M. kwangtungensis and M. pubescens var. alba. These findings strongly supported the hypothesis that for species with a closer phylogenetic relationship, pre-zygotic isolation plays an important part in limiting gene exchange in sympatric areas.

Large distribution and high sequence identity of a Copia-type retrotransposon in angiosperm families

Retrotransposons are the main component of plant genomes. Recent studies have revealed the complexity of their evolutionary dynamics. Here, we have identified Copia25 in Coffea canephora, a new plant retrotransposon belonging to the Ty1-Copia superfamily. In the Coffea genomes analyzed, Copia25 is present in relatively low copy numbers and transcribed. Similarity sequence searches and PCR analyses show that this retrotransposon with LTRs (Long Terminal Repeats) is widely distributed among the Rubiaceae family and that it is also present in other distantly related species belonging to Asterids, Rosids and monocots. A particular situation is the high sequence identity found between the Copia25 sequences of Musa, a monocot, and Ixora, a dicot species (Rubiaceae). Our results reveal the complexity of the evolutionary dynamics of the ancient element Copia25 in angiosperm, involving several processes including sequence conservation, rapid turnover, stochastic losses and horizontal transfer.

Loganin and secologanin derived tryptamine-iridoid alkaloids from Palicourea crocea and Palicourea padifolia (Rubiaceae)

During comparative analysis on Palicourea species from Costa Rica, two unusual loganin derived tryptamine-iridoid alkaloids were isolated from an accession of Palicourea crocea. Besides the already known brachycerine (2), palicroceaine (1) features a novel hexacyclic backbone. A second provenance, however, yielded strictosidinic acid (3), belonging to the more common secologanin derived tryptamine-iridoid alkaloids, such as those found in Palicourea padifolia. From this species, strictosidine (4), lyaloside (5) and its derivative (E)-O-(6')-(4″-hydroxy-3″,5″-dimethoxy)-cinnamoyl lyaloside (6) could be isolated. A herbarium specimen-based screening was performed, indicating some degree of regional differentiation in alkaloid content and biosynthetic pathways within the widespread and variable Pal. crocea. It further shows its differentiation from the related strictosidine containing Palicourea croceoides. The occurrence of loganin derived tryptamine-iridoid alkaloids in Pal. crocea, Psychotria brachyceras and Psychotria brachypoda, all putatively unrelated members of the Palicourea s.l. clade, is a noteworthy exception within the genus, otherwise largely characterized by secologanin-derived tryptamine-iridoid alkaloids.

Reciprocal herkogamy promotes disassortative mating in a distylous species with intramorph compatibility

Mating patterns in heterostylous species with intramorph compatibility have the potential to deviate from symmetrical disassortative mating owing to ecological and reproductive factors influencing pollen dispersal. Here, we investigate potential and realized patterns of mating in distylous Luculia pinceana (Rubiaceae), a species with intramorph compatibility. Our analysis provides an opportunity to test Darwin's hypothesis that reciprocal herkogamy promotes disassortative pollen transfer. We combined measurements of sex-organ reciprocity and pollen production to predict potential pollen transfer and mating patterns in a population from SW China. Marker-based paternity analysis was then used to estimate realized patterns of disassortative and assortative mating at the individual and floral morph levels. Both potential and realized mating patterns indicated a significant component of disassortative mating, satisfying theoretical conditions for the maintenance of floral dimorphism. Levels of assortative mating (37.7%) were significantly lower than disassortative mating (62.3%), but numerous offspring resulting from intramorph mating were detected in the majority of maternal seed families in both floral morphs. Our results provide empirical support for Darwin's cross-promotion hypothesis on the function of reciprocal herkogamy, but indicate that in most heterostylous species strong diallelic incompatibility may be a general requirement for complete disassortative mating.

A revised time tree of the asterids: establishing a temporal framework for evolutionary studies of the coffee family (rubiaceae)

Divergence time analyses in the coffee family (Rubiaceae) have all relied on the same Gentianales crown group age estimate, reported by an earlier analysis of the asterids, for defining the upper age bound of the root node in their analyses. However, not only did the asterid analysis suffer from several analytical shortcomings, but the estimate itself has been used in highly inconsistent ways in these Rubiaceae analyses. Based on the original data, we here reanalyze the divergence times of the asterids using relaxed-clock models and 14 fossil-based minimum age constraints. We also expand the data set to include an additional 67 taxa from Rubiaceae sampled across all three subfamilies recognized in the family. Three analyses are conducted: a separate analysis of the asterids, which completely mirrors the original asterid analysis in terms of taxon sample and data; a separate analysis of the Gentianales, where the result from the first analysis is used for defining a secondary root calibration point; and a combined analysis where all taxa are analyzed simultaneously. Results are presented in the form of a time-calibrated phylogeny, and age estimates for asterid groups, Gentianales, and major groups of Rubiaceae are compared and discussed in relation to previously published estimates. Our updated age estimates for major groups of Rubiaceae provide a significant step forward towards the long term goal of establishing a robust temporal framework for the divergence of this biologically diverse and fascinating group of plants.

Genetic connectivity of the moth pollinated tree Glionnetia sericea in a highly fragmented habitat

Long-distance gene flow is thought to be one prerequisite for the persistence of plant species in fragmented environments. Human influences have led to severe fragmentation of native habitats in the Seychelles islands, with many species surviving only in small and isolated populations. The endangered Seychelles endemic tree Glionnetia sericea is restricted to altitudes between 450 m and 900 m where the native forest vegetation has been largely lost and replaced with exotic invasives over the last 200 years. This study explores the genetic and ecological consequences of population fragmentation in this species by analysing patterns of genetic diversity in a sample of adults, juveniles and seeds, and by using controlled pollination experiments. Our results show no decrease in genetic diversity and no increase in genetic structuring from adult to juvenile cohorts. Despite significant inbreeding in some populations, there is no evidence of higher inbreeding in juvenile cohorts relative to adults. A Bayesian structure analysis and a tentative paternity analysis indicate extensive historical and contemporary gene flow among remnant populations. Pollination experiments and a paternity analysis show that Glionnetia sericea is self-compatible. Nevertheless, outcrossing is present with 7% of mating events resulting from pollen transfer between populations. Artificial pollination provided no evidence for pollen limitation in isolated populations. The highly mobile and specialized hawkmoth pollinators (Agrius convolvuli and Cenophodes tamsi; Sphingidae) appear to promote extensive gene flow, thus mitigating the potential negative ecological and genetic effects of habitat fragmentation in this species. We conclude that contemporary gene flow is sufficient to maintain genetic connectivity in this rare and restricted Seychelles endemic, in contrast to other island endemic tree species with limited contemporary gene flow.

Pre- and post-zygotic reproductive isolation between co-occurring Mussaenda pubescens var. alba and M. shikokiana (Rubiaceae)

Reproductive isolation is a fundamental requirement for speciation and includes several sequential stages. Few studies have determined the relative contributions of pre- and post-zygotic reproductive isolation in plants, especially between relative species with clear differentiation in flower form. To investigate the mechanisms responsible for reproductive isolation in sympatric Mussaenda pubescens var. alba and Mussaenda shikokiana (Rubiaceae) in Guangxi Province, China, we made observations of flowering phenology, patterns of insect visitation, and conducted pollination experiments, including artificial hybridization. The two species had overlapping flowering times and were pollinated by overlapping pollinators; however, their relative importance differed significantly with M. pubescens visited more commonly by bees and M. shikokiana more frequently by butterflies. Using vegetative and floral characters and molecular evidence based on nuclear ribosomal internal and external transcribed spacer regions we detected seven naturally occurring hybrids among a sample of approximately 125 individuals. Hybrids were characterized by morphologies that most closely resembled their maternal parents based on chloroplast evidence. Studies of artificially synthesized and natural hybrids demonstrated that hybrid seed had very low germination rates and naturally occurring hybrids exhibited pollen sterility. Post-zygotic reproductive isolating mechanisms play a primary role in limiting gene exchange between co-occurring species and maintaining species integrity in areas of sympatry.

Phylogenetic lineages in Vanguerieae (Rubiaceae) associated with Burkholderia bacteria in sub-Saharan Africa

PREMISE OF THE STUDY

It is well known that mutualistic bacteria can provide substantial benefits to their host plants. However, 'how,' 'why,' and the possible applications of such an interaction are only second to the questions 'who is involved?', and 'where does it occur?'. In the coffee family (Rubiaceae), certain species closely interact with endophytic leaf bacteria that are freely distributed among the mesophyll cells. This type of interaction was recently discovered in South Africa. Our aim is to document the bacterial diversity associated with Rubiaceae ('who') and to establish the geographic range of the interaction ('where').

METHODS

Representatives of the Vanguerieae tribe in Rubiaceae were investigated for the presence of endophytes with special emphasis on the distributional range of the plant-bacteria association by collecting specimens from different African regions.

KEY RESULTS

The interaction is found in five genera and is restricted to three major host lineages. The endophytic bacteria belong to the genus Burkholderia and are part of the plant-associated beneficial and environmental group. Some endophytes are similar to B. caledonica, B. graminis, B. phenoliruptrix or B. phytofirmans, while others are classified in OTUs that show no similarity with any previously described Burkholderia species of bacteria.

CONCLUSIONS

The association is not obligate from the bacterial point of view and is considered a loose and recent interaction, which is demonstrated by the fact that there is no evidence for coevolution. The geographical distribution of the association is restricted by the distributional range of the host plants covering the whole of sub-Saharan Africa.

Evolutionary history of the Afro-Madagascan Ixora species (Rubiaceae): species diversification and distribution of key morphological traits inferred from dated molecular phylogenetic trees

BACKGROUND AND AIMS

Previous work on the pantropical genus Ixora has revealed an Afro-Madagascan clade, but as yet no study has focused in detail on the evolutionary history and morphological trends in this group. Here the evolutionary history of Afro-Madagascan Ixora spp. (a clade of approx. 80 taxa) is investigated and the phylogenetic trees compared with several key morphological traits in taxa occurring in Madagascar.

METHODS

Phylogenetic relationships of Afro-Madagascan Ixora are assessed using sequence data from four plastid regions (petD, rps16, rpoB-trnC and trnL-trnF) and nuclear ribosomal external transcribed spacer (ETS) and internal transcribed spacer (ITS) regions. The phylogenetic distribution of key morphological characters is assessed. Bayesian inference (implemented in BEAST) is used to estimate the temporal origin of Ixora based on fossil evidence.

KEY RESULTS

Two separate lineages of Madagascan taxa are recovered, one of which is nested in a group of East African taxa. Divergence in Ixora is estimated to have commenced during the mid Miocene, with extensive cladogenesis occurring in the Afro-Madagascan clade during the Pliocene onwards.

CONCLUSIONS

Both lineages of Madagascan Ixora exhibit morphological innovations that are rare throughout the rest of the genus, including a trend towards pauciflorous inflorescences and a trend towards extreme corolla tube length, suggesting that the same ecological and selective pressures are acting upon taxa from both Madagascan lineages. Novel ecological opportunities resulting from climate-induced habitat fragmentation and corolla tube length diversification are likely to have facilitated species radiation on Madagascar.

Major clades of Australasian Rutoideae (Rutaceae) based on rbcL and atpB sequences

Background

Rutaceae subfamily Rutoideae (46 genera, c. 660 species) is diverse in both rainforests and sclerophyll vegetation of Australasia. Australia and New Caledonia are centres of endemism with a number of genera and species distributed disjunctly between the two regions. Our aim was to generate a high-level molecular phylogeny for the Australasian Rutoideae and identify major clades as a framework for assessing morphological and biogeographic patterns and taxonomy.

Methodology/Principal Findings

Phylogenetic analyses were based on chloroplast genes, rbcL and atpB, for 108 samples (78 new here), including 38 of 46 Australasian genera. Results were integrated with those from other molecular studies to produce a supertree for Rutaceae worldwide, including 115 of 154 genera. Australasian clades are poorly matched with existing tribal classifications, and genera Philotheca and Boronia are not monophyletic. Major sclerophyll lineages in Australia belong to two separate clades, each with an early divergence between rainforest and sclerophyll taxa. Dehiscent fruits with seeds ejected at maturity (often associated with myrmecochory) are inferred as ancestral; derived states include woody capsules with winged seeds, samaras, fleshy drupes, and retention and display of seeds in dehisced fruits (the last two states adaptations to bird dispersal, with multiple origins among rainforest genera). Patterns of relationship and levels of sequence divergence in some taxa, mostly species, with bird-dispersed (Acronychia, Sarcomelicope, Halfordia and Melicope) or winged (Flindersia) seeds are consistent with recent long-distance dispersal between Australia and New Caledonia. Other deeper Australian/New Caledonian divergences, some involving ant-dispersed taxa (e.g., Neoschmidia), suggest older vicariance.

Conclusions/Significance

This comprehensive molecular phylogeny of the Australasian Rutoideae gives a broad overview of the group’s evolutionary and biogeographic history. Deficiencies of infrafamilial classifications of Rutoideae have long been recognised, and our results provide a basis for taxonomic revision and a necessary framework for more focused studies of genera and species.

Diversity and endemism in deglaciated areas: ploidy, relative genome size and niche differentiation in the Galium pusillum complex (Rubiaceae) in Northern and Central Europe

BACKGROUND AND AIMS

Plants endemic to areas covered by ice sheets during the last glaciation represent paradigmatic examples of rapid speciation in changing environments, yet very few systems outside the harsh arctic zone have been comprehensively investigated so far. The Galium pusillum aggregate (Rubiaceae) is a challenging species complex that exhibits a marked differentiation in boreal parts of Northern Europe. As a first step towards understanding its evolutionary history in deglaciated regions, this study assesses cytological variation and ecological preferences of the northern endemics and compares the results with corresponding data for species occurring in neighbouring unglaciated parts of Central and Western Europe.

METHODS

DNA flow cytometry was used together with confirmatory chromosome counts to determine ploidy levels and relative genome sizes in 1158 individuals from 181 populations. A formalized analysis of habitat preferences was applied to explore niche differentiation among species and ploidy levels.

KEY RESULTS

The G. pusillum complex evolved at diploid and tetraploid levels in Northern Europe, in contrast to the high-polyploid evolution of most other northern endemics. A high level of eco-geographic segregation was observed between different species (particularly along gradients of soil pH and competition) which is unusual for plants in deglaciated areas and most probably contributes to maintaining species integrity. Relative monoploid DNA contents of the species from previously glaciated regions were significantly lower than those of their counterparts from mostly unglaciated Central Europe, suggesting independent evolutionary histories.

CONCLUSIONS

The aggregate of G. pusillum in Northern Europe represents an exceptional case with a geographically vicariant and ecologically distinct diploid/tetraploid species endemic to formerly glaciated areas. The high level of interspecific differentiation substantially widens our perception of the evolutionary dynamics and speciation rates in the dramatically changing environments of Northern Europe.

Coupling deep transcriptome analysis with untargeted metabolic profiling in Ophiorrhiza pumila to further the understanding of the biosynthesis of the anti-cancer alkaloid camptothecin and anthraquinones

The Rubiaceae species, Ophiorrhiza pumila, accumulates camptothecin, an anti-cancer alkaloid with a potent DNA topoisomerase I inhibitory activity, as well as anthraquinones that are derived from the combination of the isochorismate and hemiterpenoid pathways. The biosynthesis of these secondary products is active in O. pumila hairy roots yet very low in cell suspension culture. Deep transcriptome analysis was conducted in O. pumila hairy roots and cell suspension cultures using the Illumina platform, yielding a total of 2 Gb of sequence for each sample. We generated a hybrid transcriptome assembly of O. pumila using the Illumina-derived short read sequences and conventional Sanger-derived expressed sequence tag clones derived from a full-length cDNA library constructed using RNA from hairy roots. Among 35,608 non-redundant unigenes, 3,649 were preferentially expressed in hairy roots compared with cell suspension culture. Candidate genes involved in the biosynthetic pathway for the monoterpenoid indole alkaloid camptothecin were identified; specifically, genes involved in post-strictosamide biosynthetic events and genes involved in the biosynthesis of anthraquinones and chlorogenic acid. Untargeted metabolomic analysis by Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) indicated that most of the proposed intermediates in the camptothecin biosynthetic pathway accumulated in hairy roots in a preferential manner compared with cell suspension culture. In addition, a number of anthraquinones and chlorogenic acid preferentially accumulated in hairy roots compared with cell suspension culture. These results suggest that deep transcriptome and metabolome data sets can facilitate the identification of genes and intermediates involved in the biosynthesis of secondary products including camptothecin in O. pumila.

Post-Boreotropical dispersals explain the pantropical disjunction in Paederia (Rubiaceae)

BACKGROUND AND AIMS

Pantropical intercontinental disjunction is a common biogeographical pattern in flowering plants exhibiting a discontinuous distribution primarily in tropical Asia, Africa and the Americas. Only a few plant groups with this pattern have been investigated at the generic level with molecular phylogenetic and biogeographical methods. Paederia (Rubiaceae) is a pantropical genus of 31 species of woody lianas, with the greatest species diversity in continental Asia and Madagascar and only two species from tropical America. The aim of this study was to reconstruct the biogeographical history of Paederia based on phylogenetic analyses to explore how the genus attained its pantropical distribution.

METHODS

Maximum parsimony and Bayesian inference were used for phylogenetic analyses using sequences of five plastid markers (the rbcL gene, rps16 intron, trnT-F region, atpB-rbcL spacer and psbA-trnH spacer). Biogeographical inferences were based on a Bayesian uncorrelated lognormal relaxed molecular clock together with both Bayesian and likelihood ancestral area reconstructions.

KEY RESULTS

The data suggest an early diverged Asian lineage sister to the clade of the remaining species consisting of a predominantly Asian sub-clade and a primarily Malagasy sub-clade. Paederia is inferred to have originated in the Oligocene in tropical continental Asia. It then reached Africa in the early to middle Miocene, most probably via long-distance dispersal across the Indian Ocean. The two Neotropical species are inferred to have derived independently in the late Miocene from ancestors of Asia and East Africa, respectively.

CONCLUSIONS

The results demonstrate the importance of post-Boreotropical long-distance dispersals (across three major oceans) in shaping the global pantropical disjunction in some plants, such as Paederia, with small, winged diaspores adapted to long-distance dispersal by various agents including wind, ocean currents or birds. Overland migration is less likely to explain its palaeotropical disjunction between Asia and Africa.

Development of microsatellite markers from Mussaenda pubescens (Rubiaceae)

PREMISE OF THE STUDY

Microsatellite markers were developed and characterized in Mussaenda pubescens for further study of its levels of genetic diversity and changes in population genetic structure in reproductive character displacement and in shifts of sexual systems.

METHODS AND RESULTS

Nineteen microsatellite loci were amplified successfully in M. pubescens, 17 of which were polymorphic. A maximum of eight alleles were detected per locus in 68 individuals at population level. The observed and expected heterozygosities varied from 0 to 1.000 and 0 to 0.882, respectively.

CONCLUSIONS

These newly developed microsatellite markers will be useful in further investigations of genetic diversity and gene flow among populations of M. pubescens and its congeneric species.

Development and characterization of microsatellite markers for Emmenopterys henryi (Rubiaceae), a rare tree from China

PREMISE OF THE STUDY

Compound microsatellite primers were developed for Emmenopterys henryi, an endangered deciduous tree endemic to China, to assess its genetic diversity and population structure as well as its evolutionary history.

METHODS AND RESULTS

Using the compound microsatellite marker technique, 10 pairs of polymorphic microsatellite primers were isolated and characterized in E. henryi. Levels of polymorphism were tested across a total of 63 individuals from three natural populations. Allele numbers varied from 10 to 20 per locus, with an average of 14.50 alleles per locus. The observed heterozygosity per locus ranged from 0.125 to 0.962, and the expected heterozygosity ranged from 0.377 to 0.903.

CONCLUSIONS

The highly polymorphic markers developed and characterized in this study will facilitate evolutionary and population genetic studies in E. henryi.

Mechanism and control of Genipa americana seed germination

Genipa americana (Rubiaceae) is important for restoration of riparian forest in the Brazilian Cerrado. The objective was to characterize the mechanism and control of germination of G. americana to support uniform seedling production. Morphology and morphometrics of seeds, embryo and endosperm were assessed by light and scanning electron microscopy during germination. Imbibition and germination curves were generated and over the same time interval endosperm digestion and resistance were measured by puncture force analysis and activity assay of endo-β-mannanase (EBM) in water and in abscisic acid (ABA). The gene encoding for EBM was partially cloned and its expression monitored by quantitative real-time-polymerase chain reaction. Embryos displayed growth prior to radicle protrusion. A two-phase increase in EBM activity coincided with the two stages of weakening of the micropylar endosperm. The second stage also coincided with growth of the embryo prior to radicle protrusion. Enzyme activity was initiated in the micropylar endosperm but spread to the lateral endosperm. ABA completely inhibited germination by inhibiting embryo growth, the second stage of weakening and expression of the EBM gene, but EBM activity was not significantly inhibited. This suggests that a specific isoform of the enzyme is involved in endosperm weakening. EBM may cause a general 'softening' of micropylar endosperm cell walls, allowing the embryo to puncture the endosperm as the driving force of the decrease in puncture force.

A set of novel microsatellite markers developed for Luculia yunnanensis (Rubiaceae), an endangered plant endemic to Yunnan, China

The genus Luculia Sweet contains about five species of small trees or shrubs and is a member of the family Rubiaceae (tribe Cinchoneae). Luculia yunnanensis is an endangered ornamental shrub endemic to southwest China. Only two natural populations of L. yunnanensis exist in the wild according to our field investigation. It can be inferred that L. yunnanensis is facing a very high risk of extinction in the wild and an urgent conservation strategy is required. By using a modified biotin-sterptavidin capture method, 24 primer sets were identified in two wild populations. Of these primers, 11 displayed polymorphisms and 13 were monomorphic. The number of alleles per locus ranged from two to four, values for observed and expected heterozygosities ranged from 0.000 to 0.833 and from 0.431 to 0.771, with averages of 0.389 and 0.614, respectively. These markers will be useful for further investigation of conservation of resources, selecting parental types in cross-breeding, evolution of this species at the molecular level and related research in Luculia species.