Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae.

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

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.

Forensic Implications of the Discrepancies Caused between NGS and CE Results by New Microvariant Allele at Penta E Microsatellite.

Examination of STR markers using the MPS technology is becoming more common in forensic genetics, but scientists still have insufficient experience in dealing with ambiguous results. However, it is always essential to resolve discordant data if we want to use the technology as an accredited method in routine forensic casework. During the internal laboratory validation of the Precision ID GlobalFiler NGS STR Panel v2 kit, we observed two discrepant genotypes at Penta E locus compared to the previous capillary electrophoresis results. Each NGS software that we applied (i.e., Converge, STRaitRazor and IGV) returned the same 12,14 and 12,16 genotypes in the two samples, respectively, instead of the 11.3,14 and 11.3,16 genotypes previously observed with CE (Capillary electrophoresis) typing. In the case of the length variant 11.3 alleles, traditional Sanger sequencing confirmed a complete twelve repeat unit structure in both samples. However, after sequencing was extended to the flanking regions of the variant alleles, sequence data revealed a two-bases GG deletion downstream of the last TCTTT repeat motif in the forward strand. The determined allele variant has not been previously reported in the scientific literature and highlights the need for a careful evaluation and thorough concordance studies before using NGS STR data in forensic cases.

Comprehensive analysis of OpHD-ZIP transcription factors related to the regulation of camptothecin biosynthesis in Ophiorrhiza pumila.

Ophiorrhiza pumila, as a folk herb belonging to the Rubiaceae family, has become a potential source of camptothecin (CPT), which is a monoterpenoid indole alkaloid with good antitumor property. However, the camptothecin content in this herb is low, and is far from meeting the increasing clinical demand. Understanding the transcriptional regulation of camptothecin biosynthesis provides an effective strategy for improvement of camptothecin yield. Previous studies have demonstrated several transcription factors that are related to camptothecin biosynthesis, while the functions of HD-ZIP members in O. pumila have not been investigated yet. In this study, 32 OpHD-ZIP transcription factor members were genome-wide identified. Phylogenetic tree showed that these OpHD-ZIP proteins are divided into four subfamilies. Based on the transcriptome data, nine OpHD-ZIP genes were shown to be predominantly expressed in O. pumila roots, which were in line with the camptothecin biosynthetic genes. Co-expression analysis showed that OpHD-ZIP7 and OpHD-ZIP20 were potentially related to the modulation of camptothecin biosynthesis. Dual-luciferase reporter assays (Dual-LUC) showed that both OpHD-ZIP7 and OpHD-ZIP20 could activate the expression of camptothecin biosynthetic genes OpIO and OpTDC. In conclusion, this study offered the promising data for exploring the roles of OpHD-ZIP transcription factors in regulating camptothecin biosynthesis.

Population Structure and Genetic Diversity in the Natural Distribution of Neolamarckia cadamba in China.

Neolamarckia cadamba (Roxb.) Bosser is a fast-growing deciduous tree species and belongs to the Neolamarckia genus of the Rubiaceae family. This species has great economic and medical values in addition to being an important timber species for multiple industrial purposes. However, few studies have examined the genetic diversity and population structure in the natural distribution of this species in China. Here, we applied both the haploid nrDNA ITS (619 bp for aligned sequences) and mtDNA (2 polymorphic loci) markers to investigate 10 natural populations (239 individuals in total) that covered most of the distribution of the species in China. The results showed that the nucleotide diversity was π = 0.1185 ± 0.0242 for the nrDNA ITS markers and π = 0.00038 ± 0.00052 for the mtDNA markers. The haplotype diversity for the mtDNA markers was h = 0.1952 ± 0.2532. The population genetic differentiation was small (Fstn = 0.0294) for the nrDNA ITS markers but large (Fstm = 0.6765) for the mtDNA markers. There were no significant effects of isolation by distance (IBD), by elevation, and by two climatic factors (annual average precipitation and tem perature). A geographic structure among populations (Nst

Frankia colletiae sp. nov., a nitrogen-fixing actinobacterium isolated from Colletia cruciata.

A nitrogen-fixing actinobacterium strain (Cc1.17T) isolated from a root nodule of Colletia cruciata was subjected to polyphasic taxonomic studies. The strain was characterized by the presence of meso-diaminopimelic acid in its peptidoglycan, galactose, glucose, mannose, rhamnose, ribose and xylose as cell-wall sugars, phosphatidylinositol, diphosphatidylglycerol, glycophospholipids, phosphatidylglycerol, glycophospholipid and uncharacterized lipids as its polar lipids, and C16 : 0, iso-C16 : 0, C17 : 1 ω9 and C18 : 1 ω9 as major fatty acids (>10 %). Strain Cc1.17T showed 16S rRNA gene sequence similarities of 97.4-99.8 % to validly named Frankia species. Phylogenetic trees based on 16S rRNA gene and genome sequences placed strain Cc1.17T in a new lineage within the genus Frankia. Digital DNA-DNA hybridization and average nucleotide identity values between strain Cc1.17T and its closest phylogenomic neighbours were well below the thresholds recommended for prokaryotic species delineation. Therefore, strain Cc1.17T (=DSM 43829T=CECT 9313T) merits recognition as the type strain of a new species for which the name Frankia colletiae sp. nov. is proposed.

Penta-Primer Amplification Refractory Mutation System (PARMS) with Direct PCR-Based SNP Marker-Assisted Selection (D-MAS).

The penta-primer amplification refractory mutation system (PARMS) is a high-throughput, low-cost, and automated genotyping assay system that utilizes competitive allele-specific polymerase chain reaction (AS-PCR) combined with a homogeneous fluorescence-based reporting system to detect genetic variation occurring at single-nucleotide polymorphism (SNP). It is flexible in terms of the number of SNPs and samples to be analyzed, and the whole process only needs standard liquid handling, thermal cycling instruments, and plate reading instruments, which are present in many labs. Its compatibility with DNA samples prepared from a variety of sources and extraction technologies, such as alkaline lysis, makes it suitable for a direct PCR-based SNP marker-assisted selection system (D-MAS), a simple, cost- and labor-saving, and robust SNP genotyping system. It combines rapid and high-throughput DNA extraction through modified alkaline lysis with PARMS to dramatically reduce the time of manual operation and result analysis in the molecular breeding of major crops.

Analysis of five rare alleles at the STR loci D1S1656, D12S391, D13S317, Penta D, and D2S441.

Short tandem repeat (STR) automatic typing technology is extensively used in forensic laboratories with commercial kits, in rare cases genotyping misinterpretations or mislabeling may occur due to unexpected rare alleles. This study refers to the investigation of several rare alleles observed from routine cases. Besides cross-kit verification with Goldeneye 25A (Beijing PeopleSpot Inc, China) and Huaxia platinum (Thermo Fisher Scientific, USA) kits, the next-generation sequencing technology by MiSeq FGx System (Illumina, USA) was applied to further validation. To solve the inconsistent outcomes reached by the above mentioned approaches at D2S441 locus, single gene amplification, gene cloning, and genetic sequencing was also performed. As a result, five rare alleles were detected. Two novel alleles of allele 3 at the D13S317 locus and allele 5 at the D2S441 locus were found; three previously reported alleles of allele 9 at D1S1656 locus, allele 19 at Penta D locus, and allele 28 at D12S391 locus in STRBase were initially supplemented with sequence information. We, therefore, propose that such uncommon observations with rare events should be carefully investigated and interpreted.

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.

Thymidine Kinase 2 and Mitochondrial Protein COX I in the Cerebellum of Patients with Spinocerebellar Ataxia Type 31 Caused by Penta-nucleotide Repeats (TTCCA)n.

Spinocerebellar ataxia type 31 (SCA31), an autosomal-dominant neurodegenerative disorder characterized by progressive cerebellar ataxia with Purkinje cell degeneration, is caused by a heterozygous 2.5-3.8 kilobase penta-nucleotide repeat of (TTCCA)n in intron 11 of the thymidine kinase 2 (TK2) gene. TK2 is an essential mitochondrial pyrimidine-deoxyribonucleoside kinase. Bi-allelic loss-of-function mutations of TK2 lead to mitochondrial DNA depletion syndrome (MDS) in humans through severe (~ 70%) reduction of mitochondrial electron-transport-chain activity, and tk2 knockout mice show Purkinje cell degeneration and ataxia through severe mitochondrial cytochrome-c oxidase subunit I (COX I) protein reduction. To clarify whether TK2 function is altered in SCA31, we investigated TK2 and COX I expression in human postmortem SCA31 cerebellum. We confirmed that canonical TK2 mRNA is transcribed from exons far upstream of the repeat site, and demonstrated that an extended version of TK2 mRNA ("TK2-EXT"), transcribed from exons spanning the repeat site, is expressed in human cerebellum. While canonical TK2 was conserved among vertebrates, TK2-EXT was specific to primates. Reverse transcription-PCR demonstrated that both TK2 mRNAs were preserved in SCA31 cerebella compared with control cerebella. The TK2 proteins, assessed with three different antibodies including our original polyclonal antibody against TK2-EXT, were detected as ~ 26 kilodalton proteins on western blot; their levels were similar in SCA31 and control cerebella. COX I protein level was preserved in SCA31 compared to nuclear DNA-encoded protein. We conclude that the expression and function of TK2 are preserved in SCA31, suggesting a mechanism distinct from that of MDS.

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.

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 complete organellar genomes of the entheogenic plant Psychotria viridis (Rubiaceae), a main component of the ayahuasca brew.

Psychotria viridis (Rubioideae: Rubiaceae), popularly known as chacrona, is commonly found as a shrub in the Amazon region and is well-known to produce psychoactive compounds, such as the N,N-dimethyltryptamine (DMT). Together with the liana Banisteropsis caapi, P. viridis is one of the main components of the Amerindian traditional, entheogenic beverage known as ayahuasca. In this work, we assembled and annotated the organellar genomes (ptDNA and mtDNA), presenting the first genomics resources for this species. The P. viridis ptDNA exhibits 154,106 bp, encoding all known ptDNA gene repertoire found in angiosperms. The Psychotria genus is a complex paraphyletic group, and according to phylogenomic analyses, P. viridis is nested in the Psychotrieae clade. Comparative ptDNA analyses indicate that most Rubiaceae plastomes present conserved ptDNA structures, often showing slight differences at the junction sites of the major four regions (LSC-IR-SSC). For the mitochondrion, assembly graph-based analysis supports a complex mtDNA organization, presenting at least two alternative and circular mitogenomes structures exhibiting two main repeats spanning 24 kb and 749 bp that may symmetrically isomerize the mitogenome into variable arrangements and isoforms. The circular mtDNA sequences (615,370 and 570,344 bp) encode almost all plant mitochondrial genes (except for the ccmC, rps7, rps10, rps14, rps19, rpl2 and rpl16 that appears as pseudogenes, and the absent genes sdh3, rps2, rsp4, rsp8, rps11, rpl6, and rpl10), showing slight variations related to exclusive regions, ptDNA integration, and relics of previous events of LTR-RT integration. The detection of two mitogenomes haplotypes is evidence of heteroplasmy as observed by the complex organization of the mitochondrial genome using graph-based analysis. Taken together, these results elicit the primary insights into the genome biology and evolutionary history of Psychotria viridis and may be used to aid strategies for conservation of this sacred, entheogenic species.

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.

The significance of recent diversification in the Northern Hemisphere in shaping the modern global flora revealed from the herbaceous tribe of Rubieae (Rubiaceae).

The global herbaceous flora is probably shaped by both ancient and/or recent diversification, companied with the impacts from geographic differences between the Northern and Southern Hemispheres. Therefore, its biogeographic pattern with respect to temporal and spatial divergence is far from full understanding. Tribe Rubieae, the largest herbaceous tribe in the woody-dominant Rubiaceae, provides an excellent opportunity for studying the macroevolution of worldwide colonization. Here, we aim to reconstruct the evolutionary history of Rubieae with regard to climate fluctuation and geological history in the Cenozoic. A total of 204 samples of Rubieae representing all the distribution areas of the tribe were used to infer its phylogenetic and biogeographic histories based on two nrDNA and six cpDNA regions. The ancestral area of Rubieae was reconstructed using a time-calibrated phylogeny in RASP and diversification rates were inferred using Bayesian analysis of macroevolutionary mixtures (BAMM). Our results show Rubieae probably originated in European region during the middle Oligocene, with the two subtribes separating at 26.8 million years ago (Ma). All the genera in Rubieae formed separate clades between 24.79 and 6.23 Ma. The ancestral area of the subtribe Rubiinae was the Madrean-Tethyan plant belt and the North Atlantic land bridge (NALB) provided passage between North America and Europe for Rubiinae. The subtribe Galiinae clade originated in Europe/central Asia during the late Oligocene. Two diversification shifts were detected within Rubieae in the late Neogene. Most extant Rubieae species diverged recently during the Neogene within clades that generally were established during the late Paleogene. The tribe shows complex migration/dispersal patterns within the North Hemisphere combined with multiple recent dispersals into Southern Hemisphere. Our results highlighted the important role of recent biogeographic diversification in the Northern Hemisphere in shaping the modern global herbaceous flora during the latest and rapid worldwide expansion in the Neogene.

Genome-wide survey of the GATA gene family in camptothecin-producing plant Ophiorrhiza pumila.

BACKGROUND: Ophiorrhiza pumila (Rubiaceae) is capable of producing camptothecin (CPT), one monoterpene indole alkaloid extensively employed in the treatment of multiple cancers. Transcription factors (TFs) GATA are a group of transcription regulators involved in plant development and metabolism, and show the feature of binding to the GATA motif within the promoters of target genes. However, GATA TFs have not been characterized in O. pumila. RESULT: In this study, a total of 18 GATA genes classified into four subfamilies were identified, which randomly distributed on 11 chromosomes of O. pumila. Synteny analysis of GATA genes between O. pumila and other plant species such as Arabidopsis thaliana, Oryza sativa, Glycine max, Solanum lycopersicum, Vitis vinifera, and Catharanthus roseus genomes were analyzed. Tissue expression pattern revealed that OpGATA1 and OpGATA18 were found to be correlated with ASA, MK, CPR and GPPS, which were highly expressed in leaves. OpGATA7, showed high expression in roots as most of the CPT biosynthetic pathway genes did, suggesting that these OpGATAs may be potential candidates regulating CPT biosynthesis in O. pumila. CONCLUSIONS: In this study, we systematically analyzed the OpGATA TFs, and provided insights into the involvement of OpGATA TFs from O. pumila in CPT biosynthesis.

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.

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.