P. Moolhuijzen

Deep Sequencing of Plant and Animal DNA Contained within Traditional Chinese Medicines Reveals Legality Issues and Health Safety Concerns

Traditional Chinese medicine (TCM) has been practiced for thousands of years, but only within the last few decades has its use become more widespread outside of Asia. Concerns continue to be raised about the efficacy, legality, and safety of many popular complementary alternative medicines, including TCMs. Ingredients of some TCMs are known to include derivatives of endangered, trade-restricted species of plants and animals, and therefore contravene the Convention on International Trade in Endangered Species (CITES) legislation. Chromatographic studies have detected the presence of heavy metals and plant toxins within some TCMs, and there are numerous cases of adverse reactions. It is in the interests of both biodiversity conservation and public safety that techniques are developed to screen medicinals like TCMs. Targeting both the p-loop region of the plastid trnL gene and the mitochondrial 16S ribosomal RNA gene, over 49,000 amplicon sequence reads were generated from 15 TCM samples presented in the form of powders, tablets, capsules, bile flakes, and herbal teas. Here we show that second-generation, high-throughput sequencing (HTS) of DNA represents an effective means to genetically audit organic ingredients within complex TCMs. Comparison of DNA sequence data to reference databases revealed the presence of 68 different plant families and included genera, such as Ephedra and Asarum, that are potentially toxic. Similarly, animal families were identified that include genera that are classified as vulnerable, endangered, or critically endangered, including Asiatic black bear (Ursus thibetanus) and Saiga antelope (Saiga tatarica). Bovidae, Cervidae, and Bufonidae DNA were also detected in many of the TCM samples and were rarely declared on the product packaging. This study demonstrates that deep sequencing via HTS is an efficient and cost-effective way to audit highly processed TCM products and will assist in monitoring their legality and safety especially when plant reference databases become better established.

Genome sequence of the pathogenic intestinal spirochete brachyspira hyodysenteriae reveals adaptations to its lifestyle in the porcine large intestine.

Brachyspira hyodysenteriae is an anaerobic intestinal spirochete that colonizes the large intestine of pigs and causes swine dysentery, a disease of significant economic importance. The genome sequence of B. hyodysenteriae strain WA1 was determined, making it the first representative of the genus Brachyspira to be sequenced, and the seventeenth spirochete genome to be reported. The genome consisted of a circular 3,000,694 base pair (bp) chromosome, and a 35,940 bp circular plasmid that has not previously been described. The spirochete had 2,122 protein-coding sequences. Of the predicted proteins, more had similarities to proteins of the enteric Escherichia coli and Clostridium species than they did to proteins of other spirochetes. Many of these genes were associated with transport and metabolism, and they may have been gradually acquired through horizontal gene transfer in the environment of the large intestine. A reconstruction of central metabolic pathways identified a complete set of coding sequences for glycolysis, gluconeogenesis, a non-oxidative pentose phosphate pathway, nucleotide metabolism, lipooligosaccharide biosynthesis, and a respiratory electron transport chain. A notable finding was the presence on the plasmid of the genes involved in rhamnose biosynthesis. Potential virulence genes included those for 15 proteases and six hemolysins. Other adaptations to an enteric lifestyle included the presence of large numbers of genes associated with chemotaxis and motility. B. hyodysenteriae has diverged from other spirochetes in the process of accommodating to its habitat in the porcine large intestine.

Immunological Profiles of Bos taurus and Bos indicus Cattle Infested with the Cattle Tick, Rhipicephalus (Boophilus) microplus

ABSTRACT The cattle tick, Rhipicephalus (Boophilus) microplus, is a major threat to the improvement of cattle production in tropical and subtropical countries worldwide. Bos indicus cattle are naturally more resistant to infestation with the cattle tick than are Bos taurus breeds, although considerable variation in resistance occurs within and between breeds. It is not known which genes contribute to the resistant phenotype, nor have immune parameters involved in resistance to R. microplus been fully described for the bovine host. This study was undertaken to determine whether selected cellular and antibody parameters of the peripheral circulation differed between tick-resistant Bos indicus and tick-susceptible Bos taurus cattle following a period of tick infestations. This study demonstrated significant differences between the two breeds with respect to the percentage of cellular subsets comprising the peripheral blood mononuclear cell population, cytokine expression by peripheral blood leukocytes, and levels of tick-specific immunoglobulin G1 (IgG1) antibodies measured in the peripheral circulation. In addition to these parameters, the Affymetrix bovine genome microarray was used to analyze gene expression by peripheral blood leukocytes of these animals. The results demonstrate that the Bos indicus cattle developed a stabilized T-cell-mediated response to tick infestation evidenced by their cellular profile and leukocyte cytokine spectrum. The Bos taurus cattle demonstrated cellular and gene expression profiles consistent with a sustained innate, inflammatory response to infestation, although high tick-specific IgG1 titers suggest that these animals have also developed a T-cell response to infestation.

De novo assembly of Euphorbia fischeriana root transcriptome identifies prostratin pathway related genes.

BackgroundEuphorbia fischeriana is an important medicinal plant found in Northeast China. The plant roots contain many medicinal compounds including 12-deoxyphorbol-13-acetate, commonly known as prostratin that is a phorbol ester from the tigliane diterpene series. Prostratin is a protein kinase C activator and is effective in the treatment of Human Immunodeficiency Virus (HIV) by acting as a latent HIV activator. Latent HIV is currently the biggest limitation for viral eradication. The aim of this study was to sequence, assemble and annotate the E. fischeriana transcriptome to better understand the potential biochemical pathways leading to the synthesis of prostratin and other related diterpene compounds.ResultsIn this study we conducted a high throughput RNA-seq approach to sequence the root transcriptome of E. fischeriana. We assembled 18,180 transcripts, of these the majority encoded protein-coding genes and only 17 transcripts corresponded to known RNA genes. Interestingly, we identified 5,956 protein-coding transcripts with high similarity (> = 75%) to Ricinus communis, a close relative to E. fischeriana. We also evaluated the conservation of E. fischeriana genes against EST datasets from the Euphorbeacea family, which included R. communis, Hevea brasiliensis and Euphorbia esula. We identified a core set of 1,145 gene clusters conserved in all four species and 1,487 E. fischeriana paralogous genes. Furthermore, we screened E. fischeriana transcripts against an in-house reference database for genes implicated in the biosynthesis of upstream precursors to prostratin. This identified 24 and 9 candidate transcripts involved in the terpenoid and diterpenoid biosyntehsis pathways, respectively. The majority of the candidate genes in these pathways presented relatively low expression levels except for 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase (HDS) and isopentenyl diphosphate/dimethylallyl diphosphate synthase (IDS), which are required for multiple downstream pathways including synthesis of casbene, a proposed precursor to prostratin.ConclusionThe resources generated in this study provide new insights into the upstream pathways to the synthesis of prostratin and will likely facilitate functional studies aiming to produce larger quantities of this compound for HIV research and/or treatment of patients.

Aligning a New Reference Genetic Map of Lupinus angustifolius with the Genome Sequence of the Model Legume, Lotus japonicus

We have developed a dense reference genetic map of Lupinus angustifolius (2n = 40) based on a set of 106 publicly available recombinant inbred lines derived from a cross between domesticated and wild parental lines. The map comprised 1090 loci in 20 linkage groups and three small clusters, drawing together data from several previous mapping publications plus almost 200 new markers, of which 63 were gene-based markers. A total of 171 mainly gene-based, sequence-tagged site loci served as bridging points for comparing the Lu. angustifolius genome with the genome sequence of the model legume, Lotus japonicus via BLASTn homology searching. Comparative analysis indicated that the genomes of Lu. angustifolius and Lo. japonicus are highly diverged structurally but with significant regions of conserved synteny including the region of the Lu. angustifolius genome containing the pod-shatter resistance gene, lentus. We discuss the potential of synteny analysis for identifying candidate genes for domestication traits in Lu. angustifolius and in improving our understanding of Fabaceae genome evolution.

The complete genome sequence of the pathogenic intestinal spirochete Brachyspira pilosicoli and comparison with other Brachyspira genomes

Background The anaerobic spirochete Brachyspira pilosicoli colonizes the large intestine of various species of birds and mammals, including humans. It causes “intestinal spirochetosis”, a condition characterized by mild colitis, diarrhea and reduced growth. This study aimed to sequence and analyse the bacterial genome to investigate the genetic basis of its specialized ecology and virulence. Methodology/Principal Findings The genome of B. pilosicoli 95/1000 was sequenced, assembled and compared with that of the pathogenic Brachyspira hyodysenteriae and a near-complete sequence of Brachyspira murdochii. The B. pilosicoli genome was circular, composed of 2,586,443 bp with a 27.9 mol% G+C content, and encoded 2,338 genes. The three Brachyspira species shared 1,087 genes and showed evidence of extensive genome rearrangements. Despite minor differences in predicted protein functional groups, the species had many similar features including core metabolic pathways. Genes distinguishing B. pilosicoli from B. hyodysenteriae included those for a previously undescribed bacteriophage that may be useful for genetic manipulation, for a glycine reductase complex allowing use of glycine whilst protecting from oxidative stress, and for aconitase and related enzymes in the incomplete TCA cycle, allowing glutamate synthesis and function of the cycle during oxidative stress. B. pilosicoli had substantially fewer methyl-accepting chemotaxis genes than B. hyodysenteriae and hence these species are likely to have different chemotactic responses that may help to explain their different host range and colonization sites. B. pilosicoli lacked the gene for a new putative hemolysin identified in B. hyodysenteriae WA1. Both B. pilosicoli and B. murdochii lacked the rfbBADC gene cluster found on the B. hyodysenteriae plasmid, and hence were predicted to have different lipooligosaccharide structures. Overall, B. pilosicoli 95/1000 had a variety of genes potentially contributing to virulence. Conclusions/Significance The availability of the complete genome sequence of B. pilosicoli 95/1000 will facilitate functional genomics studies aimed at elucidating host-pathogen interactions and virulence.

Comparative microarray analysis of Rhipicephalus (Boophilus) microplus expression profiles of larvae pre-attachment and feeding adult female stages on Bos indicus and Bos taurus cattle

BackgroundRhipicephalus (Boophilus) microplus is an obligate blood feeder which is host specific to cattle. Existing knowledge pertaining to the host or host breed effects on tick transcript expression profiles during the tick - host interaction is poor.ResultsGlobal analysis of gene expression changes in whole R. microplus ticks during larval, pre-attachment and early adult stages feeding on Bos indicus and Bos taurus cattle were compared using gene expression microarray analysis. Among the 13,601 R. microplus transcripts from BmiGI Version 2 we identified 297 high and 17 low expressed transcripts that were significantly differentially expressed between R. microplus feeding on tick resistant cattle [Bos indicus (Brahman)] compared to R. microplus feeding on tick susceptible cattle [Bos taurus (Holstein-Friesian)] (p ≤ 0.001). These include genes encoding enzymes involved in primary metabolism, and genes related to stress, defence, cell wall modification, cellular signaling, receptor, and cuticle formation. Microarrays were validated by qRT-PCR analysis of selected transcripts using three housekeeping genes as normalization controls.ConclusionThe analysis of all tick stages under survey suggested a coordinated regulation of defence proteins, proteases and protease inhibitors to achieve successful attachment and survival of R. microplus on different host breeds, particularly Bos indicus cattle. R. microplus ticks demonstrate different transcript expression patterns when they encounter tick resistant and susceptible breeds of cattle. In this study we provide the first transcriptome evidence demonstrating the influence of tick resistant and susceptible cattle breeds on transcript expression patterns and the molecular physiology of ticks during host attachment and feeding.The microarray data used in this analysis have been submitted to NCBI GEO database under accession number GSE20605 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE20605.

Genomic analysis of Campylobacter fetus subspecies: identification of candidate virulence determinants and diagnostic assay targets

BackgroundCampylobacter fetus subspecies venerealis is the causative agent of bovine genital campylobacteriosis, asymptomatic in bulls the disease is spread to female cattle causing extensive reproductive loss. The microbiological and molecular differentiation of C. fetus subsp. venerealis from C. fetus subsp. fetus is extremely difficult. This study describes the analysis of the available C. fetus subsp. venerealis AZUL-94 strain genome (~75–80%) to identify elements exclusively found in C. fetus subsp. venerealis strains as potential diagnostic targets and the characterisation of subspecies virulence genes.ResultsEighty Kb of genomic sequence (22 contigs) was identified as unique to C. fetus subsp. venerealis AZUL-94 and consisted of type IV secretory pathway components, putative plasmid genes and hypothetical proteins. Of the 9 PCR assays developed to target C. fetus subsp. venerealis type IV secretion system genes, 4 of these were specific for C. fetus subsp. venerealis biovar venerealis and did not detect C. fetus subsp. venerealis biovar intermedius. Two assays were specific for C. fetus subsp. venerealis AZUL-94 strain, with a further single assay specific for the AZUL-94 strain and C. fetus subsp. venerealis biovar intermedius (and not the remaining C. fetus subsp. venerealis biovar venerealis strains tested). C. fetus subsp. fetus and C. fetus subsp. venerealis were found to share most common Campylobacter virulence factors such as SAP, chemotaxis, flagellar biosynthesis, 2-component systems and cytolethal distending toxin subunits (A, B, C). We did not however, identify in C. fetus the full complement of bacterial adherence candidates commonly found in other Campylobacter spp.ConclusionThe comparison of the available C. fetus subsp. venerealis genome sequence with the C. fetus subsp. fetus genome identified 80 kb of unique C. fetus subsp. venerealis AZUL94 sequence, with subsequent PCR confirmation demonstrating inconsistent amplification of these targets in all other C. fetus subsp. venerealis strains and biovars tested. The assays developed here highlight the complexity of targeting strain specific virulence genes for field studies for the molecular identification and epidemiology of C. fetus.

Complete genome sequence of sporisorium scitamineum and biotrophic interaction transcriptome with sugarcane

Sporisorium scitamineum is a biotrophic fungus responsible for the sugarcane smut, a worldwide spread disease. This study provides the complete sequence of individual chromosomes of S. scitamineum from telomere to telomere achieved by a combination of PacBio long reads and Illumina short reads sequence data, as well as a draft sequence of a second fungal strain. Comparative analysis to previous available sequences of another strain detected few polymorphisms among the three genomes. The novel complete sequence described herein allowed us to identify and annotate extended subtelomeric regions, repetitive elements and the mitochondrial DNA sequence. The genome comprises 19,979,571 bases, 6,677 genes encoding proteins, 111 tRNAs and 3 assembled copies of rDNA, out of our estimated number of copies as 130. Chromosomal reorganizations were detected when comparing to sequences of S. reilianum, the closest smut relative, potentially influenced by repeats of transposable elements. Repetitive elements may have also directed the linkage of the two mating-type loci. The fungal transcriptome profiling from in vitro and from interaction with sugarcane at two time points (early infection and whip emergence) revealed that 13.5% of the genes were differentially expressed in planta and particular to each developmental stage. Among them are plant cell wall degrading enzymes, proteases, lipases, chitin modification and lignin degradation enzymes, sugar transporters and transcriptional factors. The fungus also modulates transcription of genes related to surviving against reactive oxygen species and other toxic metabolites produced by the plant. Previously described effectors in smut/plant interactions were detected but some new candidates are proposed. Ten genomic islands harboring some of the candidate genes unique to S. scitamineum were expressed only in planta. RNAseq data was also used to reassure gene predictions.

The genetics of symbiotic nitrogen fixation: Comparative genomics of 14 rhizobia strains by resolution of protein clusters

The symbiotic relationship between legumes and nitrogen fixing bacteria is critical for agriculture, as it may have profound impacts on lowering costs for farmers, on land sustainability, on soil quality, and on mitigation of greenhouse gas emissions. However, despite the importance of the symbioses to the global nitrogen cycling balance, very few rhizobial genomes have been sequenced so far, although there are some ongoing efforts in sequencing elite strains. In this study, the genomes of fourteen selected strains of the order Rhizobiales, all previously fully sequenced and annotated, were compared to assess differences between the strains and to investigate the feasibility of defining a core ‘symbiome’—the essential genes required by all rhizobia for nodulation and nitrogen fixation. Comparison of these whole genomes has revealed valuable information, such as several events of lateral gene transfer, particularly in the symbiotic plasmids and genomic islands that have contributed to a better understanding of the evolution of contrasting symbioses. Unique genes were also identified, as well as omissions of symbiotic genes that were expected to be found. Protein comparisons have also allowed the identification of a variety of similarities and differences in several groups of genes, including those involved in nodulation, nitrogen fixation, production of exopolysaccharides, Type I to Type VI secretion systems, among others, and identifying some key genes that could be related to host specificity and/or a better saprophytic ability. However, while several significant differences in the type and number of proteins were observed, the evidence presented suggests no simple core symbiome exists. A more abstract systems biology concept of nitrogen fixing symbiosis may be required. The results have also highlighted that comparative genomics represents a valuable tool for capturing specificities and generalities of each genome.