Luke R. Tembrock

Are differences in genomic data sets due to true biological variants or errors in genome assembly: an example from two chloroplast genomes.

DNA sequencing has been revolutionized by the development of high-throughput sequencing technologies. Plummeting costs and the massive throughput capacities of second and third generation sequencing platforms have transformed many fields of biological research. Concurrently, new data processing pipelines made rapid de novo genome assemblies possible. However, high quality data are critically important for all investigations in the genomic era. We used chloroplast genomes of one Oryza species (O. australiensis) to compare differences in sequence quality: one genome (GU592209) was obtained through Illumina sequencing and reference-guided assembly and the other genome (KJ830774) was obtained via target enrichment libraries and shotgun sequencing. Based on the whole genome alignment, GU592209 was more similar to the reference genome (O. sativa: AY522330) with 99.2% sequence identity (SI value) compared with the 98.8% SI values in the KJ830774 genome; whereas the opposite result was obtained when the SI values in coding and noncoding regions of GU592209 and KJ830774 were compared. Additionally, the junctions of two single copies and repeat copies in the chloroplast genome exhibited differences. Phylogenetic analyses were conducted using these sequences, and the different data sets yielded dissimilar topologies: phylogenetic replacements of the two individuals were remarkably different based on whole genome sequencing or SNP data and insertions and deletions (indels) data. Thus, we concluded that the genomic composition of GU592209 was heterogeneous in coding and non-coding regions. These findings should impel biologists to carefully consider the quality of sequencing and assembly when working with next-generation data.

The Complete Plastid Genome of Lagerstroemia fauriei and Loss of rpl2 Intron from Lagerstroemia (Lythraceae)

Lagerstroemia (crape myrtle) is an important plant genus used in ornamental horticulture in temperate regions worldwide. As such, numerous hybrids have been developed. However, DNA sequence resources and genome information for Lagerstroemia are limited, hindering evolutionary inferences regarding interspecific relationships. We report the complete plastid genome of Lagerstroemia fauriei. To our knowledge, this is the first reported whole plastid genome within Lythraceae. This genome is 152,440 bp in length with 38% GC content and consists of two single-copy regions separated by a pair of 25,793 bp inverted repeats. The large single copy and the small single copy regions span 83,921 bp and 16,933 bp, respectively. The genome contains 129 genes, including 17 located in each inverted repeat. Phylogenetic analysis of genera sampled from Geraniaceae, Myrtaceae, and Onagraceae corroborated the sister relationship between Lythraceae and Onagraceae. The plastid genomes of L. fauriei and several other Lythraceae species lack the rpl2 intron, which indicating an early loss of this intron within the Lythraceae lineage. The plastid genome of L. fauriei provides a much needed genetic resource for further phylogenetic research in Lagerstroemia and Lythraceae. Highly variable markers were identified for application in phylogenetic, barcoding and conservation genetic applications.

Multilocus analysis of nucleotide variation and speciation in three closely related Populus (Salicaceae) species

Historical tectonism and climate oscillations can isolate and contract the geographical distributions of many plant species, and they are even known to trigger species divergence and ultimately speciation. Here, we estimated the nucleotide variation and speciation in three closely related Populus species, Populus tremuloides, P. tremula and P. davidiana, distributed in North America and Eurasia. We analysed the sequence variation in six single‐copy nuclear loci and three chloroplast (cpDNA) fragments in 497 individuals sampled from 33 populations of these three species across their geographic distributions. These three Populus species harboured relatively high levels of nucleotide diversity and showed high levels of nucleotide differentiation. Phylogenetic analysis revealed that P. tremuloides diverged earlier than the other two species. The cpDNA haplotype network result clearly illustrated the dispersal route from North America to eastern Asia and then into Europe. Molecular dating results confirmed that the divergence of these three species coincided with the sundering of the Bering land bridge in the late Miocene and a rapid uplift of the Qinghai‐Tibetan Plateau around the Miocene/Pliocene boundary. Vicariance‐driven successful allopatric speciation resulting from historical tectonism and climate oscillations most likely played roles in the formation of the disjunct distributions and divergence of these three Populus species.

A Multiplex Real-Time PCR Assay to Diagnose and Separate Helicoverpa armigera and H. zea (Lepidoptera: Noctuidae) in the New World.

The Old World bollworm, Helicoverpa armigera (Hübner), and the corn earworm, H. zea (Boddie), are two of the most important agricultural pests in the world. Diagnosing these two species is difficult—adults can only be separated with a complex dissection, and larvae cannot be identified to species using morphology, necessitating the use of geographic origin for identification in most instances. With the discovery of H. armigera in the New World, identification of immature Helicoverpa based on origin is no longer possible because H. zea also occurs in all of the geographic regions where H. armigera has been discovered. DNA barcoding and restriction fragment length polymorphism (RFLP) analyses have been reported in publications to distinguish these species, but these methods both require post-PCR processing (i.e., DNA sequencing or restriction digestion) to complete. We report the first real-time PCR assay to distinguish these pests based on two hydrolysis probes that bind to a segment of the internal transcribed spacer region 2 (ITS2) amplified using a single primer pair. One probe targets H. armigera, the second probe targets H. zea, and a third probe that targets a conserved segment of 18S rDNA is used as a control of DNA quality. The assay can be completed in 50 minutes when using isolated DNA and is successfully tested on larvae intercepted at ports of entry and adults captured during domestic surveys. We demonstrate that the assay can be run in triplex with no negative effects on sensitivity, can be run using alternative real-time PCR reagents and instruments, and does not cross react with other New World Heliothinae.

The complete plastid genome of the wild rice species Oryza brachyantha (Poaceae)

Abstract The whole plastid genome of wild rice (Oryza brachyantha) is characterized in this study. The genome is 134 604 bp in length and is arranged in a typical circular structure, including a pair of inverted repeats (IRs) of 20 832 bp in size separated by a large single-copy region (LSC) of 80 411 bp in length and a small single-copy region (SSC) of 12 529 bp in length. The overall GC content of the genome is 38.98%. One hundred and ten unique genes were annotated from the chloroplast genome, including 76 protein-coding genes, 4 ribosomal RNA genes and 30 tRNA genes. A total of 20 of these genes are duplicated in the IR regions, 13 genes contain 1 intron and 2 genes (rps12 and ycf3) have 2 introns.

The Complete Chloroplast Genome of Catha edulis: A Comparative Analysis of Genome Features with Related Species

Qat (Catha edulis, Celastraceae) is a woody evergreen species with great economic and cultural importance. It is cultivated for its stimulant alkaloids cathine and cathinone in East Africa and southwest Arabia. However, genome information, especially DNA sequence resources, for C. edulis are limited, hindering studies regarding interspecific and intraspecific relationships. Herein, the complete chloroplast (cp) genome of Catha edulis is reported. This genome is 157,960 bp in length with 37% GC content and is structurally arranged into two 26,577 bp inverted repeats and two single-copy areas. The size of the small single-copy and the large single-copy regions were 18,491 bp and 86,315 bp, respectively. The C. edulis cp genome consists of 129 coding genes including 37 transfer RNA (tRNA) genes, 8 ribosomal RNA (rRNA) genes, and 84 protein coding genes. For those genes, 112 are single copy genes and 17 genes are duplicated in two inverted regions with seven tRNAs, four rRNAs, and six protein coding genes. The phylogenetic relationships resolved from the cp genome of qat and 32 other species confirms the monophyly of Celastraceae. The cp genomes of C. edulis, Euonymus japonicus and seven Celastraceae species lack the rps16 intron, which indicates an intron loss took place among an ancestor of this family. The cp genome of C. edulis provides a highly valuable genetic resource for further phylogenomic research, barcoding and cp transformation in Celastraceae.

Characterization of the whole chloroplast genome of Chikusichloa mutica and its comparison with other rice tribe (Oryzeae) species

Chloroplast genomes are a significant genomic resource in plant species and have been used in many research areas. The complete genomic information from wild crop species could supply a valuable genetic reservoir for breeding. Chikusichloa mutica is one of the most important wild distant relatives of cultivated rice. In this study, we sequenced and characterized its complete chloroplast (cp) genome and compared it with other species in the same tribe. The whole cp genome sequence is 136,603 bp in size and exhibits a typical quadripartite structure with large and small single-copy regions (LSC, 82,327 bp; SSC, 12,598 bp) separated by a pair of 20,839-bp inverted repeats (IRA, B). A total of 110 unique genes are annotated, including 76 protein-coding genes, 4 ribosomal RNA genes and 30 tRNA genes. The genome structure, gene order, GC content, and other features are similar to those of other angiosperm cp genomes. When comparing the cp genomes between Oryzinae and Zizaniinae subtribes, the main differences were found between the junction regions and distribution of simple sequence repeats (SSRs). In comparing the two Chikusichloa species, the genomes were only 40 bp different in length and 108 polymorphic sites, including 83 single nucleotide substitutions (SNPs) and 25 insertion-deletions (Indels), were found between the whole cp genomes. The complete cp genome of C. mutica will be an important genetic tool for future breeding programs and understanding the evolution of wild rice relatives.

Two complete chloroplast genomes of white campion (Silene latifolia) from male and female individuals

Abstract In this study, we assembled two individuals’ complete chloroplast genomes with one male and one female from the dioecious plant species white campion (Silene latifolia). The two chloroplast genomes have an identical composition with each other as a circular molecule of 150 931 bp in length, with an overall GC content of 36.4%. The genomes consist of a pair of inverted repeats (IRs) of 25 503 bp, separated by a large single-copy (LSC) region and a small single-copy (SSC) region of 82 708 and 17 217 bp, respectively. The genomes contain 111 single copy genes, including 77 protein-coding genes, 4 ribosomal RNA genes and 30 transfer RNA genes. In these protein-coding genes, eight genes (rpl16, rpoC1, rps16, petD, petB, ndhB, ndhA and atpF) contain a single intron and three genes (rps12, clpP and ycf3) contain two introns. The two newly sequenced chloroplast genomes provide valuable information for detecting polymorphisms within species and between sexes.

Employing Two‐stage Derivatisation and GC–MS to Assay for Cathine and Related Stimulant Alkaloids across the Celastraceae

INTRODUCTION Catha edulis (qat, khat, mirra) is a woody plant species that is grown and consumed in East Africa and Yemen for its stimulant alkaloids cathinone, cathine and norephedrine. Two Celastraceae species, in addition to qat, have been noted for their stimulant properties in ethnobotanical literature. Recent phylogenetic reconstructions place four genera in a clade sister to Catha edulis, and these genera are primary candidates to search for cathine and related alkaloids. OBJECTIVE Determine if cathine or related alkaloids are present in species of Celastraceae other than Catha edulis. METHODS Leaf samples from 43 Celastraceae species were extracted in water followed by basification of the aqueous extract and partitioning with methyl-t-butyl ether to provide an alkaloid-enriched fraction. The extract was derivatised in a two-stage process and analysed using GC-MS for the presence of cathine. Related alkaloids and other metabolites in this alkaloid-enriched fraction were tentatively identified. RESULTS Cathinone, cathine and norephedrine were not detected in any of the 43 Celastraceae species assayed other than Catha edulis. However, the phenylalanine- or tyrosine-derived alkaloid phenylethylamine was identified in five species. Nine species were found to be enriched for numerous sterol- and terpene-like compounds. CONCLUSION These results indicate that cathine is unique to Catha edulis, and not the compound responsible for the stimulant properties reported in related Celastraceae species. Copyright

The Complete Chloroplast Genome of Heimia myrtifolia and Comparative Analysis within Myrtales

Heimia myrtifolia is an important medicinal plant with several pharmacologically active alkaloids and is also used as an ornamental landscape plant. The purpose of this study is to complete and characterize the chloroplast (cp) genome of H. myrtifolia and compare genomic features to other Myrtales species’ cp genomes. The analysis showed that H. myrtifolia has a total length of 159,219 bp with a typical quadripartite structure containing two identical inverted repeats (IRs) of 25,643 bp isolated by one large single copy (LSC) of 88,571 bp and one small single copy (SSC) of 18,822 bp. The H. myrtifolia cp genome contains 129 genes with eight ribosomal RNAs, 30 transfer RNAs, and 78 protein coding genes, in which 17 genes are duplicated in two IR regions. The genome organization including gene type and number and guanine-cytosine (GC) content is analyzed among the 12 cp genomes in this study. Approximately 255 simple sequence repeats (SSRs) and 16 forward, two reverses, and two palindromic repeats were identified in the H. myrtifolia cp genome. By comparing the whole H. myrtifolia cp genome with 11 other Myrtales species, the results showed that the sequence similarity was high between coding regions while sequence divergence was high between intergenic regions. By employing the full cp genomes for phylogenetic analysis, structural and sequence differences were characterized between H. myrtifolia and 11 Myrtales species illustrating what patterns are common in the evolution of cp genomes within the Myrtales. The first entire cp genome in the genus Heimia provides a valuable resource for further studies in these medicinally and ornamentally important taxa.