Background: Based on the conserved sequences of a known NBS resistance gene, a pair of degenerate primers was designed to amplify the NBS-LRR resistance gene from peanut using PCR and RACE methods.Results: Analyzing the amino acid sequence by BLAST on NCBI, which was deduced from the 1088bp-long gene named PnAG1-2, showed that it had a certain homology with some resistance proteins, among which Arachis cardenasii resistance protein gene had the highest homology (66%).Relative quantification PCR analysis indicated that PnAG1-2 gene expresses more in J11 (an A. flavus-resistant variety) than in JH1012 (an A. flavus-susceptible variety) when the harvest time was coming.Conclusions: In this study, the NBS-LRR resistance sequence was successfully cloned from peanut and prokaryotic expression was done on the gene, which provided a foundation for cultivating anti-A.flavus peanut varieties.
Abstract Background/purpose Cancer immunotherapy has revolutionized the clinical treatment of several tumors. Immune infiltration has been found to be closely related to clinical prognosis, but it shows limited activity in breast cancer (BC). Therefore, this study aimed to explore the infiltration pattern of immune cells in BC, and to find potential prognostic markers and new therapeutic targets.Patients and methods We downloaded the immune genome data of BC from the Cancer Genome Atlas (TCGA), and analyzed the tumor- infiltrating immune cells (TIICs) in BC for the first time using the CIBERSORT algorithm. The aim of this study was to assess the proportions of 22 immune cell subsets in BC and examine the correlation between each TIIC and overall survival (OS) as well as clinical characteristics.Results The results indicated that: (1) there was a significant difference between the immune infiltration spectrum of cancerous and adjacent tissues, with M2 macrophages, M0 macrophages, and CD4 + T cells being highly expressed in BC; (2) CD8 + T cells were positively correlated with activated CD4 + memory T cells and negatively correlated with M0 macrophages, and M2 macrophages was inversely correlated with M1 macrophages, T cells regulatory, T cells CD8; (3) T cells, macrophages and BC TNM stage, age, clinical stage were correlated (P < 0.05); and (4) high expression of M2 macrophage markers could be an independent biomarker of poor prognosis and a potential therapeutic target for BC.Conclusion This study provides a new research method for the systematic study of immune cells in the BC tumor microenvironment, and provides theoretical guidance for further experiments to verify M2 macrophages and T cell subsets as a potential target for immunotherapy and prognosis.
Evodiae Fructus (EF) is a common herbal medicine with thousands of years of medicinal history in China, which has been demonstrated with many promising pharmacological effects on cancer, cardiovascular diseases and Alzheimer's disease. However, there have been increasing reports of hepatotoxicity associated with EF consumption. Unfortunately, in a long term, many implicit constituents of EF as well as their toxic mechanisms remain poorly understood. Recently, metabolic activation of hepatotoxic compounds of EF to generate reactive metabolites (RMs) has been implicated. Herein, we capture metabolic reactions relevant to hepatotoxicity of these compounds. Initially, catalyzed by the hepatic cytochrome P450 enzymes (CYP450s), the hepatotoxic compounds of EF are oxidized to generate RMs. Subsequently, the highly electrophilic RMs could react with nucleophilic groups contained in biomolecules, such as hepatic proteins, enzymes, and nucleic acids to form conjugates and/or adducts, leading to a sequence of toxicological consequences. In addition, currently proposed biological pathogenesis, including oxidative stress, mitochondrial damage and dysfunction, endoplasmic reticulum (ER) stress, hepatic metabolism disorder, and cell apoptosis are represented. In short, this review updates the knowledge on the pathways of metabolic activation of seven hepatotoxic compounds of EF and provides considerable insights into the relevance of proposed molecular hepatotoxicity mechanisms from a biochemical standpoint, for the purpose of providing a theoretical guideline for the rational application of EF in clinics.
Objective of this paper is to excavate the siderophore synthesis gene from Brevibacillus brevis GZDF3 and verify its type and antibacterial effects. The method is using genome mining technology to analyze the siderophore synthesis genes and the phylogenetic tree of each synthesis gene was constructed separately. Iron free medium was utilized to induce the synthesis of siderophore and CAS liquid detection method was used for qualitative and quantitative analysis on siderophore. The type of siderophore was preliminaries identified by Arnow and its antibacterial effects were analyzed according to the agar punching method. The results show that a siderophore synthesis gene cluster with 83% similarity to Petrobactin was found in Brevibacillus brevis GZDF3 genome. Iron free medium could induce siderophore synthesis and the optimal incubation time cultured in iron free medium was 30 h and 48 h. Antagonistic strain GZDF3 had the capacity to synthesize catechol-type siderophore. Also, GZDF3 had a powerful antibacterial effect on pathogenic fungus Fusarium oxysporum of rotted root on Pinellia ternata. Therefore, Brevibacillus brevis GZDF3 can produce catechol-type siderophore in an iron-deficient culture medium, which was also a main antifungal active substance.
Peptides are a class of protein fragments with relatively high biological activity, strong specificity. Some of them have stimulated considerable interest because of their unique advantages in many diseases. Modern studies have shown that multiple anticoagulant protein play a crucial role in Shuxuetong injection (SXT). However, the extraordinary complexity of Chinese medicinal formulate and the lack of identification method are primary challenges for ingredients. In addition, infinitesimal peptides contents further hinder the identification and structural characterization of polypeptide by traditional means. In this paper, we described a strategy that LC-MS combined with molecular docking to illustrate the peptide components of SXT. The key to this strategy was use of gene sequencing to establish a SXT protein database to further achieve the separation and enrichment of chemical methods. Moreover, the ADRA2A, PAR4, DRD3 were precisely docked with the identified peptides. The result indicated that the 14 compounds had stable binding ability, which were speculated to be the latent bioactive monomers for the treatment of stroke. Additionally, 14 peptides were verified by cell-based experiment. The results showed that YLKTT could indeed protect astrocytes from OGD/R. The YLKTT showed values of higher activity than the others in vitro. It may be a completely new compound that has never been reported before and provides the basis for further research and a new paradigm for stroke.
Precision livestock management requires animal traceability and disease trajectory, for which discriminating between or re-identifying individual animals is of significant importance. Existing re-identification (re-ID) methods are mostly proposed for persons and vehicles, compared with which animals are extraordinarily more challenging to be re-identified because of subtle visual differences between individuals. In this paper, we focus on image-based re-ID of yaks (Bos grunniens), which are indispensable livestock in local animal husbandry economy in Qinghai-Tibet Plateau. We establish the first yak re-ID dataset (called YakReID-103) which contains 2, 247 images of 103 different yaks with bounding box, direction-based pose, and identity annotations. Moreover, according to the characteristics of yaks, we modifiy several person re-ID and animal re-ID methods as baselines for yak re-ID. Experimental results of the baselines on YakReID-103 demonstrate the challenges in yak re-ID. We expect that the proposed benchmark will promote the research of animal biometrics and extend the application scope of re-ID techniques.
Due to the rapidly increasing global interest in the use of herbs, phytomedicines and other natural products as medical or complementary remedies, concerns about the clinical medication safety have drawn much attention worldwide. Particularly, many natural ingredients exhibit inhibitory effects on cytochrome P450 (CYP) enzymes, which are the most important Phase I metabolism enzymes in liver. CYP2C9 is one of the most abundant CYP enzymes and responsible for the metabolism of over 15% clinical drugs, including oral sulfonylurea hypoglycemics, nonsteroidal anti-inflammatory agents, selective cyclooxygenase-2 inhibitors, antiepileptics, angiotensin II receptor inhibitors and anticoagulants. Diclofenac (4'-hydroxylase) and tolbutamide (methylhydroxylation) are widely used as probe substrates for CYP2C9. To date, numerous natural products have been reported to have the capabilities of inhibiting the catalytic activity of CYP2C9 and further influencing the pharmacokinetic and pharmacodynamic behaviors of drugs that are mainly metabolized by CYP2C9, leading to potential herb-drug interactions. Moreover, some fatal adverse interactions may occur for drugs with a narrow therapeutic window when they are coadministered with a CYP2C9 inhibitor, especially irreversible inactivators. For the purpose of better understanding the interactions of natural products with CYP2C9, we comprehensively reviewed the characteristics of CYP2C9, the natural ingredients that inhibit CYP2C9, the related research approaches and strategies, the types of inhibition and the underlying mechanisms.
TP53 mutations have been observed in diffuse large B-cell lymphoma (DLBCL), with a mean frequency of ~20%. Studies on TP53 mutations as prognostic markers have historically been controversial, and the results have not been consistent across different studies on DLBCL. Considering the complex pathophysiological mechanisms involved in DLBCL, we wondered whether the interaction of TP53 with other genetic variants could further promote the development of DLBCL, and thus be more prognostically predictive. Moreover, whether the genetic interactions between TP53 and other oncogenic mutations could shape the discrepant immune landscape in DLBCL remains unknown, as these genetic alterations usually drive the malignant phenotype and directly or indirectly affect the tumor microenvironment (TME) and support tumor survival. In this study, we performed a comprehensive analysis of the genomic characteristics of TP53 through high-throughput sequencing in patients with de novo DLBCL. Patients' characteristics are reported in Table S1. Detailed methods are provided in the Supplementary Material. A total of 227 significantly mutated genes were identified (Table S2), of which TP53 was the second most frequently mutated gene, with a rate of 30% (53 of 176) and 62 sequence variants detected. Among these variants, 74% (n = 46/62) were missense mutations, and the remaining were inactivating frameshift indels (n = 7), nonsense mutations (n = 3), coding sequencing indels (n = 4), and splicing mutations (n = 2). Mutation patterns and distributions are shown in Figure S1 and Table S3. Importantly, most mutations (56/64, 87.5%) occurred in exons 5–8, which encoded the DNA-binding domain (DBD) region of TP53 (Figure S1C,E). Codons 175, 273, and 248 of the p53 protein had the highest mutation frequency, which are also the hot spots of TP53 mutation found in most human cancers (Figure S1D). Given that the DBD of TP53 is the functional central core domain and mutations in this region potentially have a strong impact on TP53 function, we mainly focused on mutations in this region. Patients were divided into TP53-MUT and TP53-WT groups according to TP53 mutation status in the DBD region. There were no differences in the number of small deletions/insertions and single nucleotide variants (SNVs) between the TP53-MUT and TP53-WT groups (Figure S2). Moreover, the tumor mutation burden (TMB) was similar between the two groups (Figure S2E). Clinical relevance between TP53 mutation and clinicopathological characteristics, such as age, sex, B symptoms, stage, number of extranodal sites, performance status, LDH level, and International Prognostic Index was not observed (Table S4). Note that TP53 mutations significantly enriched in the GCB subtype (p = .033, Table S4). Among the 176 patients, 155 patients having the complete follow-up data were enrolled in the survival analysis. Overall, patients with TP53-MUT tended to have inferior overall survival compared with patients with TP53-WT (median: 92.3 versus 110.8 months, respectively, p = .17, Figure 1A), but it did not reach statistical significance. A subgroup analysis showed that the potential predictive value of TP53 mutations was mainly attributed to the GCB subtype (Figure S3). We next recognized the genomic variants that co-occur or are mutually exclusive with TP53. We observed that DDX3X, MYLK2, and FUT6 mutations co-occurred with TP53 mutations, and CD58 mutations were mutually exclusive with TP53 mutations (Figure 1B and Table S5). Specifically, patients were divided into four groups based on the mutation status of these genes. No significant difference was observed in survival among the three groups according to the combination of co-occurring mutation genes with TP53 (p = .37 for DDX3X; p = .11 for MYLK2; p = .54 for FUT6; Figure S4). However, we found that the combination of TP53 and CD58 mutations could significantly distinguish the prognosis of patients with DLBCL (p = .033, Figure 1C). Patients with both wild-type TP53 and CD58 had a better prognosis than patients with either of the two mutually exclusive modes of CD58 and TP53 mutations. Unexpectedly, patients with TP53 wild-type and CD58 mutations (TP53WT&CD58MUT) had worse survival than those with TP53 mutations and CD58 wild-type (TP53MUT&CD58WT). Because TP53 and CD58 mutations were mutually exclusive, only one patient harbored both TP53 and CD58 mutations, and the patient still alive at the last follow-up. The predictive value of TP53WT&CD58MUT group was also observed in patients with GCB-DLBCL (Figure S5). Moreover, the relationship of a mutually exclusive mutant between CD58 and TP53 and the prognostic significance of this interaction were validated using publicly available data from 1001 patients with DLBCL from the Duke University's cohort1 (Figure S6). We then explored whether the cooperation of the mutually exclusive mutations between TP53 and CD58 may profoundly influence the microenvironment in DLBCL. We found that the overall TMB was significantly higher in the TP53WT&CD58MUT group than in the TP53MUT&CD58WT group (p = .0177, Figure 1D), while there was no difference in TMB when dividing patients only according to TP53 mutation status (p = .5348, Figure S2E). In addition, the ESTIMATE immune scores in the TP53WT&CD58MUT group were significantly higher than that in the TP53MUT&CD58WT groups (p = .0047) (Figure 1D). Moreover, the exhausted T cell, macrophage cell, NK cell, and Th1 cell enriched in the TP53WT&CD58MUT group (Figure 1E). The difference between the two groups was mainly due to the combined influence of the mutation pattern "TP53WT&CD58MUT," but rather only affected by the CD58 mutations, given the immune cell infiltration was similar when dividing patients just according to CD58 mutation status (Figure S7). Furthermore, the co-inhibitory receptors such as PD-1, TIM3, and LAG3 were preferentially expressed in the TP53WT&CD58MUT group (Figure 1F and Table S6). However, there was no difference in the expression of co-stimulatory molecules (Table S6). Inhibitory immunomodulators were also significantly upregulated in the TP53WT&CD58MUT group when comparing with the TP53WT&CD58WT group (Figure S8), suggesting the unique immune phenotype in the TP53WT&CD58MUT group. The findings that high immune scores and abundant infiltrating exhausted T cells in the TP53WT&CD58MUT group were validated in an independent external cohort from the REMoDL-B trail (N = 400)2 (Figure S9 and Table S7). Finally, we investigated the differentially biological pathways between the TP53WT&CD58MUT and the TP53MUT&CD58WT groups. Five hundred differentially expressed genes were identified with a false discovery rate less than 0.05 and |log2foldchange| > 1. One hundred genes were significantly upregulated and 400 genes were significantly downregulated in the TP53WT&CD58MUT group (Figure 1G). Figure S10 presented the enriched gene ontology terms in the TP53WT&CD58MUT group, including cytokine and chemokine production, binding and activity, and interferon-γ pathways. GSEA showed significantly activated interferon-α and interferon-γ responses and IL-6/JAK/STAT3 signaling in the TP53WT&CD58MUT group (Figure 1H). The immune landscape of patients with TP53WT&CD58MUT harbored was summarized and conceptualized in Figure 1I. The value of TP53 gene alterations in predicting survival in DLBCL remains controversial, even in the era of sequencing. In the L.M. Staudt' study, four prominent genetic subtypes were identified based on the molecular classifications. Notably, TP53 was not significantly enriched in one of these subtypes, although TP53 was the much frequently mutated gene (25.2%). On the basis of this classification, L.M. Staudt and colleague further distinguished ST2, A53, and mixed subtypes, and the survival of A53, characterized by inactivation of TP53, was intermediate in this model.3 The C2 molecular subtype identified in Chapuy's study corresponding to A53 also showed a tendency for a poor prognosis.4 A subsequent study demonstrated that the impact of TP53 mutations on survival was relied on the genetic context of the lymphoma, conferring no effect in the SOCS1/SGK1 clusters and NOTCH2 subtype and inferior prognosis in the MYD88 subtype.5 These results demonstrated that TP53 alterations have limited ability to identify a subset of patients at high risk. DLBCL is highly molecular heterogeneity and there are still a proportion of patients who could not be accurately classified according to the existing molecular classifications. In this study, we found that patients with TP53WT&CD58MUT had the worst outcomes. Interestingly, this subtype of patients harbored enhanced immune escape capacity, giving the abundant infiltration of exhausted T cells and multiple upregulated inhibitory immunomodulatory molecules. Persistent interferon signaling could augment the expression of T-cell inhibitory immune checkpoints such as PD-1, TIM-3, and LAG-3 through JAK/STAT pathway.6 We found that interferon responses and JAK/STAT pathway enriched in the TP53WT&CD58MUT group. Moreover, inhibitory immunomodulatory molecules were preferentially expressed in this group. It suggests that interferon/JAK/STAT pathway-mediated up-regulated expression of inhibitory immunomodulatory molecules could be one potential mechanism through which the inactivation of CD58 facilitated immune evasion and accelerated tumor growth in DLBCL. Despite many of such immune dysregulations had an adverse prognosis; they provided new opportunities for anti-tumor immunotherapy of the subset of DLBCL patients with TP53WT&CD58MUT. Historically, the response rate to anti-PD-1/PD-L1 therapy in unselected DLBCL patients was generally low. Consequently, patients with TP53WT&CD58MUT may be optimal candidates for novel immunotherapy in clinical trials. In conclusion, our results suggest that TP53 mutation alone is insufficient to effectively differentiate the risk of DLBCL. The mutually exclusive patterns between TP53 and CD58 mutations accurately stratified patients with DLBCL to permit the optional immunotherapy. This study was supported by Natural Science Foundation of Tianjin grants (19JCYBJC26500 and 18JCZDJC45100), National Natural Science Foundation of China grants (81770213), Clinical Oncology Research Fund of CSCO grants (Y-XD2019-162 and Y-Roche20192-0097), The Science and Technology Research Program of Tianjin Education Commission (2019KJ191), National Key New Drug Creation Special Programs grants (2018ZX09201015), and National Human Genetic Resources Sharing Service Platform/Cancer Biobank of Tianjin Medical University Cancer Institute and Hospital grant (2005DKA21300). The authors thank the Marvel Medical Laboratory, Tianjin Marvelbio Technology Co., Ltd., for providing the assistance in next-generation sequencing and bioinformatics analysis. The authors declare no conflict of interest. Xianhuo Wang conceived and designed the study; Xianhuo Wang and Huilai Zhang supervised all aspects of research project and interpreted data; Tingting Zhang, Yaxiao Lu, and Xia Liu performed the research and statistical and bioinformatics analyses; Mengmeng Zhao performed the next-generation sequencing and bioinformatics analysis; Jin He, Xia Liu, Lanfang Li, Lihua Qiu, Zhengzi Qian, and Shiyong Zhou collected samples and clinical information; Bin Meng and Qiongli Zhai reviewed the diagnosis of DLBCL; Xianhuo Wang, Huilai Zhang, and Xiubao Ren provided the clinical samples and material support; Tingting Zhang wrote the manuscript and finalized the figures; Xianhuo Wang. and Huilai Zhang reviewed the manuscript. All authors read and approved the final version of the manuscript. DNA and RNA sequencing data have been submitted to the CNGB Sequence Archive of China National GeneBank DataBase (https://db.cngb.org/cnsa/) under the accession numbers CNP0001322 and CNP0001327, respectively. The mutation data of 1001 DLBCL patients from the Duke University's cohort were downloaded from cBioPortal (http://www.cbioportal.org/). The mutation and RNA expression data of the patients from the REMoDL-B trail were downloaded from the Supplementary Material from the publication. DNA and RNA sequencing data have been submitted to the CNGB Sequence Archive of China National GeneBank DataBase (https://db.cngb.org/cnsa/) under the accession numbers CNP0001322 and CNP0001327, respectively. The mutation data of 1001 DLBCL patients from the Duke University's cohort were downloaded from cBioPortal (http://www.cbioportal.org/). The mutation and RNA expression data of the patients from the REMoDL-B trail were downloaded from the Supplementary Material from the publication. Appendix S1 Supporting Information Figure S1. Mutation profile of the TP53 gene in diffuse large B-cell lymphoma. (A) Proportions of TP53 mutations according to the effect on the protein sequence. (B) Proportions of classified point mutations. (C) Distribution of mutation numbers according to TP53 exons. (D) Codon distribution of TP53 mutations. E. Mapping of the TP53 mutation sites Figure S2. TP53 mutation status and genomic instability at the individual nucleotide level. (A) The distribution of small deletion, small insertion, single nucleotide variant (SNV), and tumor mutation burden (TMB) in patients with TP53-MUT and TP53-WT. (B–E) Comparison of the number of small deletions, small insertions, SNVs, and TMB in patients with TP53-MUT and TP53-WT, respectively. Significance threshold was set at p < .05 Figure S3. Survival analysis by TP53 mutations in the molecular subtypes of diffuse large B-cell lymphoma (DLBCL). (A) Overall survival of patients with GCB-DLBCL with TP53 mutations. (B) Overall survival of patients with non-GCB-DLBCL with TP53 mutations Figure S4. Survival analysis by genomic alterations correlated with TP53 mutations in diffuse large B-cell lymphoma (DLBCL) patients. (A) Overall survival of DLBCL patients based on the mutation status of TP53 and DDX3X, MYLK2 (B) and FUT6 (C) Figure S5. Survival analysis by TP53 combined with CD58 in the molecular subtypes of diffuse large B-cell lymphoma (DLBCL). (A) Overall survival of DLBCL patients based on the mutation status of CD58 and TP53 in GCB subtype. (B) Overall survival of DLBCL patients based on the mutation status of CD58 and TP53 in non-GCB subtype Figure S6. Validation of mutually exclusive mutational patterns of TP53 and CD58 using a public data set from 1001 diffuse large B-cell lymphoma (DLBCL) patients. (A) CD58 mutations were mutually exclusive with TP53 mutations according to Fisher's exact test. (B) Overall survival in DLBCL patients according to the different mutational patterns of TP53 and CD58. (C) Overall survival in DLBCL patients, excluding the TP53MUT&CD58MUT mutational pattern due to the lower occurrence of co-mutations between TP53 and CD58 Figure S7. Comparison of immune infiltrating cells between the CD58 mutation group and CD58 wild-type group. *p < .05; NS, p > .05 Figure S8. Comparison of immunomodulatory molecule expression between the TP53WT&CD58MUT group and the TP53WT&CD58WT group Figure S9. Enhanced immune escape in the TP53WT&CD58MUT group in the external validation cohort from the REMoDL-B trail. (A) Comparison of ESTIMATE immune scores, tumor purity, and stromal scores between the TP53MUT&CD58WT and the TP53WT&CD58MUT groups. (B) Comparison of the exhausted CD8+ T cell abundance between the TP53MUT&CD58WT and the TP53WT&CD58MUT groups. (C) Comparison of the inhibitory immunomodulatory molecule expression between the TP53MUT&CD58WT and the TP53WT&CD58MUT groups. The adjusted p values were .076 for PDCD1, 0.009 for LAG3, 0.009 for HAVCR2, and 0.023 for IDO1, respectively Figure S10. Gene ontology analysis based on the significantly upregulated genes in the TP53WT&CD58MUT group. Selected and significantly enriched (FDR-adjusted p-value < .05) gene ontology annotations for biological processes and molecular functions are represented as dots Table S1. Clinicopathologic features of the 176 diffuse large B cell lymphoma patients enrolled in the study Table S2. List of the significantly mutated genes identified in our study Table S3. List of identified somatic non-silent mutations of TP53 in 53 DLBCL cases Table S4. Association of TP53 mutation statues and clinicopathologic parameters Table S5. Genes that are significantly co-mutated or mutually exclusive with TP53 mutations Table S6. Comparison of immunomodulatory molecule expression in the TP53WT&CD58MUT group versus the TP53MUT&CD58WT group Table S7. Comparison of immunomodulatory molecule expression in the TP53WT&CD58MUT group versus the TP53MUT&CD58WT group in the external cohort from the REMoDL-B trail Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.
Abstract Background Berberidis Radix (Sankezhen), a typical multi-origin Chinese medicinal material, originates from the dried roots of the Berberis genus plants. Previous studies have shown that there were significant differences in chemical compositions in different Berberis species. Therefore, it was exceedingly important to accurately identify different species of Berberis . Results In this study, for the first time, we systematically compared the complete chloroplast genome sequences of the six Berberis species ( B. julianae , B. tsienii , B. pruinose , B. thunbergii , B. poiretii and B. wilsoniae ), which commonly were used as medicinal herb Berberidis Radix. The ndhD-ccsA as highly divergent region was found and taken as a potential marker for species identification. Subsequently, the barcode was applied to the Chinese patent medicines containing Berberidis Radix (Sankezhen) combined with DNA metabarcoding technology. The results showed that the six complete chloroplast genomes exhibited a typical quadripartite structure which ranging from 165,934 to 168,828 bp in length. A total of 147 unique genes were identified in each chloroplast genome, comprising 101 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. Comparative genome analysis demonstrated that the six chloroplast genomes were highly conserved in genome size, gene organization and GC contents. The phylogenetic relationships of six Berberis plants were revealed and the results showed that Mahonia was supported as separate clade in the Berberis genera tree, which was coincident with previous studies. The nucleotide diversity analysis revealed seven variable loci in protein coding regions, and four variable loci in gene spacer regions, respectively. The primer pair 1508F-1864R on ndhD-ccsA region was proven to precisely discriminate the six studied Berberis species and recovered the biodiversity of Berberis species in Chinese patent medicines. Conclusions In general, this study provides meaningful genetic information for Berberis plants, and establishes a method to realize the identification of Berberidis Radix as multi-origin Chinese medicinal materials, which can be applied to Chinese patent medicines containing Berberidis Radix.
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)
