A low-cost platform suitable for sequencing-based recovery of natural variation in understudied plants.

Genetic characterization of wild and cultivated plants provides valuable knowledge for conservation and agriculture. DNA sequencing technologies are improving, and costs are dropping. Yet analysis of many species is hindered because they grow in regions that lack infrastructure for advanced molecular biology. The authors developed and adapted low-cost methods that address these issues. Tissue was collected and stored in silica gel, avoiding the need for liquid nitrogen and freezers. The authors optimized low-cost, homemade DNA extraction to increase yields, reduce costs and produce DNA suitable for next-generation sequencing. The authors describe how to build a gel documentation system for DNA quantification. As a proof of principle, the authors used these methods to evaluate wild Berberis darwinii, native to Southern Chile.

Inheritance of Virulence and Linkages of Virulence Genes in an Ethiopian Isolate of the Wheat Stripe Rust Pathogen (Puccinia striiformis f. sp. tritici) Determined Through Sexual Recombination on Berberis holstii.

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is one of the most devastating wheat diseases in Ethiopia. To study virulence genetics of the pathogen, 117 progeny isolates were produced through sexual reproduction of an Ethiopian isolate of the stripe rust pathogen on Berberis holstii plants under controlled conditions. The parental and progeny isolates were characterized by phenotyping on wheat lines carrying single Yr genes for resistance and genotyped using 10 polymorphic simple sequence repeated (SSR) markers. The progeny isolates were classified into 37 virulence phenotypes and 75 multilocus genotypes. The parental isolate and progeny isolates were all avirulent to resistance genes Yr5, Yr10, Yr15, Yr24, Yr32, YrTr1, YrSP, and Yr76 but virulent to Yr1 and Yr2, indicating that the parental isolate was homozygous avirulent or homozygous virulent at these loci. The progeny isolates segregated for virulence to 12 Yr genes. Virulence phenotypes to Yr6, Yr28, Yr43, and Yr44 were controlled by a single dominant gene; those to Yr7, Yr9, Yr17, Yr27, Yr25, Yr31, and YrExp2 were each controlled by two dominant genes; and the virulence phenotype to Yr8 was controlled by two complementary dominant genes. A linkage map was constructed with seven SSR markers, and 16 virulence loci corresponding to 11 Yr resistance genes were mapped with some loci linked to each other. These results are useful in understanding host-pathogen interactions and selecting resistance genes to develop wheat cultivars with highly effective resistance to stripe rust.

Mapping non-host resistance to the stem rust pathogen in an interspecific barberry hybrid.

BACKGROUND: Non-host resistance (NHR) presents a compelling long-term plant protection strategy for global food security, yet the genetic basis of NHR remains poorly understood. For many diseases, including stem rust of wheat [causal organism Puccinia graminis (Pg)], NHR is largely unexplored due to the inherent challenge of developing a genetically tractable system within which the resistance segregates. The present study turns to the pathogen's alternate host, barberry (Berberis spp.), to overcome this challenge. RESULTS: In this study, an interspecific mapping population derived from a cross between Pg-resistant Berberis thunbergii (Bt) and Pg-susceptible B. vulgaris was developed to investigate the Pg-NHR exhibited by Bt. To facilitate QTL analysis and subsequent trait dissection, the first genetic linkage maps for the two parental species were constructed and a chromosome-scale reference genome for Bt was assembled (PacBio + Hi-C). QTL analysis resulted in the identification of a single 13 cM region (~ 5.1 Mbp spanning 13 physical contigs) on the short arm of Bt chromosome 3. Differential gene expression analysis, combined with sequence variation analysis between the two parental species, led to the prioritization of several candidate genes within the QTL region, some of which belong to gene families previously implicated in disease resistance. CONCLUSIONS: Foundational genetic and genomic resources developed for Berberis spp. enabled the identification and annotation of a QTL associated with Pg-NHR. Although subsequent validation and fine mapping studies are needed, this study demonstrates the feasibility of and lays the groundwork for dissecting Pg-NHR in the alternate host of one of agriculture's most devastating pathogens.

First report of bicolour FISH of Berberis diaphana and B. soulieana reveals interspecific differences and co-localization of (AGGGTTT)3 and rDNA 5S in B. diaphana.

BACKGROUND: Berberis consists of approximately 500 species and is the largest genus in Berberidaceae. Most Berberis species lack cytological data, and bicolour fluorescence in situ hybridization (FISH) has never been performed on Berberis. In this work, a karyotype of Berberis diaphana, an alpine Berberis species obtained from an altitude of 3600 m in Wolong National Nature Reserve, China, was analysed and compared with Berberis soulieana Schneid. via FISH using oligonucleotide telomere probes for (AGGGTTT)3 and 5S rDNA (41 bp) for the first time. RESULTS: Berberis diaphana belonged to cytotype 2A and had the karyotype formula 2n = 2x = 28 = 26 m + 2 sm (2SAT). The mitotic metaphase chromosome lengths ranged from 1.82 ± 0.04 µm to 2.75 ± 0.00 µm. Clear (AGGGTTT)3 signals were detected at two telomeres in every chromosome and were co-localized with 5S rDNA at the terminal regions of the long arms in the 6th pair of chromosomes. One pair of (AGGGTTT)3 sites was localized in the satellites of the 7th pair of chromosomes, which are the only submetacentric chromosomes in this species. Totally 28 chromosomes with one pair of satellited chromosomes were observed in B. soulieana. This species had four 5S rDNA signals with two weak signals at the end of long arms in the 5th pair of chromosomes and another two strong signals detected in the interstitial region close to the end of short arms in the 6th pair of chromosomes. Each large signal consisted of two smaller signals with secondary constrictions around them. CONCLUSIONS: FISH physical mapping of B. diaphana suggested that (AGGGTTT)3 and rDNA 5S co-localize at the 6th pair of chromosomes. The density, location and number difference of 5S rDNA loci indicated structural differences among the chromosomes between B. diaphana and B. soulieana. Our results provide information that may contribute to future studies on the physical assembly of the Berberis genome and the evolution of rDNA and telomere FISH patterns in Berberis.

Classification of barberry genotypes by multivariate analysis of biochemical constituents and HPLC profiles.

INTRODUCTION: Recently, there has been a growing interest in the use of edible barberry and their extracts as a source of natural antioxidants in food and pharmaceutical industries. OBJECTIVE: The aim of this study was to evaluate the biochemical constituents of 18 samples of barberry fruits and classification of barberry genotypes by multivariate analysis. METHODS: Total phenolic, total flavonoid, total anthocyanin, total tannin, total carbohydrate contents and antioxidant activity were determined using Folin-Ciocalteu, aluminum chloride, colorimetric, vanillin, anthron and DPPH (2,2'-diphenyl-1-picrylhydrazyl) assays, respectively. High-performance liquid chromatography (HPLC) system is used for quantitative determination of phytochemical constituents. The multivariate data analysis (principal component analysis and hierarchical cluster analysis) and heat map data visualisation techniques were performed to classify barberry genotypes using Minitab and GraphPad Prism software, respectively. RESULTS: The highest amounts of total phenolics and flavonoids were obtained in fruit extracts of G3 (Berberis vulgaris). The highest total anthocyanin content and antioxidant activity were observed in G8 (B. vulgaris) and G16 (B. vulgaris), respectively. HPLC analysis of phytochemicals (gallic acid, caffeic acid, chlorogenic acid, p-coumaric acid, cinnamic acid, rutin, apigenin, and quercetin) revealed that gallic acid and p-coumaric acid were found as the most abundant phytochemical compounds. Based on multivariate analysis and heat map visualisation techniques, Berberis genotypes were classified into three main clusters. CONCLUSIONS: These results showed that barberry species (especially B. vulgaris and B. carataegina) are promising sources of natural antioxidants and biochemical compounds beneficial to be used in the food industry and that the multivariate analysis was a suitable approach to classify the barberry samples.

Comparative Analysis of Genetic and Chemical Differences between Four Berberis Herbs Based on Molecular Phylogenetic and HPLC Methods.

In traditional Tibetan medicinal system, Berberis herbs mainly originate from the dried barks of Berberis kansuensis, Berberis dictyophylla, Berberis diaphana, and Berberis vernae. In this study, molecular phylogenetic method based on four markers (i.e., rbcL, internal transcribed spacer (ITS), ITS2, and psbA-trnH) and HPLC chemical analysis were used to evaluate the chemical and genetic differences between the four Berberis species. The results showed that the discriminatory power of ITS, ITS2 and psbA-trnH was low, but the rbcL marker was highly effective and reliable for the species differentiation. The four Berberis species can be successfully classified based on phylogenetic analysis of the rbcL sequences. Moreover, the results of chemical analysis showed that four main alkaloids (i.e., berberine, palmatine, magnoflorine, and jatrorrhizine) cannot be used as chemical markers for discrimination of the four Berberis species. These findings provide valuable information for distinguishing the four Berberis Tibetan herbs.

An interspecific barberry hybrid enables genetic dissection of non-host resistance to the stem rust pathogen Puccinia graminis.

Stem rust, caused by Puccinia graminis (Pg), remains a devastating disease of wheat, and the emergence of new Pg races virulent on deployed resistance genes fuels the ongoing search for sources of durable resistance. Despite its intrinsic durability, non-host resistance (NHR) is largely unexplored as a protection strategy against Pg, partly due to the inherent challenge of developing a genetically tractable system within which NHR segregates. Here, we demonstrate that Pg's far less studied ancestral host, barberry (Berberis spp.), provides such a unique pathosystem. Characterization of a natural population of B. ×ottawensis, an interspecific hybrid of Pg-susceptible B. vulgaris and Pg-resistant B. thunbergii (Bt), reveals that this uncommon nothospecies can be used to dissect the genetic mechanism(s) of Pg-NHR exhibited by Bt. Artificial inoculation of a natural population of B. ×ottawensis accessions, verified via genotyping by sequencing to be first-generation hybrids, revealed 51% susceptible, 33% resistant, and 16% intermediate phenotypes. Characterization of a B. ×ottawensis full sib family excluded the possibility of maternal inheritance of the resistance. By demonstrating segregation of Pg-NHR in a hybrid population, this study challenges the assumed irrelevance of Bt to Pg epidemiology and lays a novel foundation for the genetic dissection of NHR to one of agriculture's most studied pathogens.

Colonization in North American Arid Lands: The Journey of Agarito (Berberis trifoliolata) Revealed by Multilocus Molecular Data and Packrat Midden Fossil Remains.

Here we conduct research to understand the evolutionary history of a shrubby species known as Agarito (Berberis trifoliolata), an endemic species to the Chihuahuan Desert. We identify genetic signatures based on plastid DNA and AFLP markers and perform niche modelling and spatial connectivity analyses as well as niche modelling based on records in packrats to elucidate whether orogenic events such as mountain range uplift in the Miocene or the contraction/expansion dynamics of vegetation in response to climate oscillations in the Pliocene/Pleistocene had an effect on evolutionary processes in Agarito. Our results of current niche modelling and palaeomodelling showed that the area currently occupied by Berberis trifoliolata is substantially larger than it was during the Last Interglacial period and the Last Glacial Maximum. Agarito was probably confined to small areas in the Northeastern and gradually expanded its distribution just after the Last Glacial Maximum when the weather in the Chihuahuan Desert and adjacent regions became progressively warmer and drier. The most contracted range was predicted for the Interglacial period. Populations remained in stable areas during the Last Glacial Maximum and expanded at the beginning of the Holocene. Most genetic variation occured in populations from the Sierra Madre Oriental. Two groups of haplotypes were identified: the Mexican Plateau populations and certain Northeastern populations. Haplogroups were spatially connected during the Last Glacial Maximum and separated during interglacial periods. The most important prediction of packrat middens palaeomodelling lies in the Mexican Plateau, a finding congruent with current and past niche modelling predictions for agarito and genetic results. Our results corroborate that these climate changes in the Pliocene/Pleistocene affected the evolutionary history of agarito. The journey of agarito in the Chihuahuan Desert has been dynamic, expanding and contracting its distribution range and currently occupying the largest area in its history.

Paralogs are revealed by proportion of heterozygotes and deviations in read ratios in genotyping-by-sequencing data from natural populations.

Whole-genome duplications have occurred in the recent ancestors of many plants, fish, and amphibians, resulting in a pervasiveness of paralogous loci and the potential for both disomic and tetrasomic inheritance in the same genome. Paralogs can be difficult to reliably genotype and are often excluded from genotyping-by-sequencing (GBS) analyses; however, removal requires paralogs to be identified which is difficult without a reference genome. We present a method for identifying paralogs in natural populations by combining two properties of duplicated loci: (i) the expected frequency of heterozygotes exceeds that for singleton loci, and (ii) within heterozygotes, observed read ratios for each allele in GBS data will deviate from the 1:1 expected for singleton (diploid) loci. These deviations are often not apparent within individuals, particularly when sequence coverage is low; but, we postulated that summing allele reads for each locus over all heterozygous individuals in a population would provide sufficient power to detect deviations at those loci. We identified paralogous loci in three species: Chinook salmon (Oncorhynchus tshawytscha) which retains regions with ongoing residual tetrasomy on eight chromosome arms following a recent whole-genome duplication, mountain barberry (Berberis alpina) which has a large proportion of paralogs that arose through an unknown mechanism, and dusky parrotfish (Scarus niger) which has largely rediploidized following an ancient whole-genome duplication. Importantly, this approach only requires the genotype and allele-specific read counts for each individual, information which is readily obtained from most GBS analysis pipelines.

Male function for ensuring pollination and reproductive success in Berberis lycium Royle: A novel mechanism.

In Berberis lycium anthers on alternate stamens dehisce, thus prolonging the male function so that pollination is affected and reproduction is ensured. The large pollen sac of each bithecous anther after the appearance of longitudinal dehiscence slit moves away from the filament while remaining attached at the tip of the connective and then orients in such a way that pollen-laden surface faces the stigma. No pollen is available to receptive stigma as pollen grains remain stuck to the anther sac. They do not get dispersed even by wind. Pollination and consequently reproduction is ensured through the intervention of insect, which does not affect pollen transfer to the stigma directly but by touching the base of the staminal filament while foraging nectar secreted by nectaries at the base of corolla, thus leading to staminal movement. This makes the dehisced anthers stick to the stigma and deposit pollen there.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

[A new herbs traceability method based on DNA barcoding-origin-morphology analysis--an example from an adulterant of 'Heiguogouqi'].

The fruit of Lycium ruthenicum is a common folk medicine in China. Now it is popular for its antioxidative effect and other medical functions. The adulterants of the herb confuse consumers. In order to identify a new adulterant of L. ruthenicum, a research was performed based on NCBI Nucleotide Database ITS Sequence, combined analysis of the origin and morphology of the adulterant to traceable varieties. Total genomic DNA was isolated from the materials, and nuclear DNA ITS sequences were amplified and sequenced; DNA fragments were collated and matched by using ContingExpress. Similarity identification of BLAST analysis was performed. Besides, the distribution of plant origin and morphology were considered to further identification and verification. Families and genera were identified by molecular identification method. The adulterant was identified as plant belonging to Berberis. Origin analysis narrowed the range of sample identification. Seven different kinds of plants in Berberis were potential sources of the sample. Adulterants variety was traced by morphological analysis. The united molecular identification-origin-morphology research proves to be a preceding way to medical herbs traceability with time-saving and economic advantages and the results showed the new adulterant of L. ruthenicum was B. kaschgarica. The main differences between B. kaschgarica and L. ruthenicum are as follows: in terms of the traits, the surface of B. kaschgarica is smooth and crispy, and that of L. ruthenicum is shrinkage, solid and hard. In microscopic characteristics, epicarp cells of B. aschgarica thickening like a string of beads, stone cells as the rectangle, and the stone cell walls of L. ruthenicum is wavy, obvious grain layer. In molecular sequences, the length of ITS sequence of B. kaschgarica is 606 bp, L. ruthenicum is 654 bp, the similarity of the two sequences is 53.32%.

Development and characterization of microsatellite markers for Berberis thunbergii (Berberidaceae).

PREMISE OF THE STUDY: Microsatellite markers were isolated and characterized in Berberis thunbergii, an invasive and ornamental shrub in the eastern United States, to assess genetic diversity among populations and potentially identify horticultural cultivars. METHODS AND RESULTS: A total of 12 loci were identified for the species. Eight of the loci were polymorphic and were screened in 24 individuals from two native (Tochigi and Ibaraki prefectures, Japan) and one invasive (Connecticut, USA) population and 21 horticultural cultivars. The number of alleles per locus ranged from three to seven, and observed heterozygosity ranged from 0.048 to 0.636. CONCLUSIONS: These new markers will provide tools for examining genetic relatedness of B. thunbergii plants in the native and invasive range, including phylogeographic studies and assessment of rapid evolution in the invasive range. These markers may also provide tools for examining hybridization with other related species in the invasive range.

Heterologous expression of two FAD-dependent oxidases with (S)-tetrahydroprotoberberine oxidase activity from Arge mone mexicana and Berberis wilsoniae in insect cells.

Berberine, palmatine and dehydrocoreximine are end products of protoberberine biosynthesis. These quaternary protoberberines are elicitor inducible and, like other phytoalexins, are highly oxidized. The oxidative potential of these compounds is derived from a diverse array of biosynthetic steps involving hydroxylation, intra-molecular C-C coupling, methylenedioxy bridge formation and a dehydrogenation reaction as the final step in the biosynthesis. For the berberine biosynthetic pathway, the identification of the dehydrogenase gene is the last remaining uncharacterized step in the elucidation of the biosynthesis at the gene level. An enzyme able to catalyze these reactions, (S)-tetrahydroprotoberberine oxidase (STOX, EC 1.3.3.8), was originally purified in the 1980s from suspension cells of Berberis wilsoniae and identified as a flavoprotein (Amann et al. 1984). We report enzymatic activity from recombinant STOX expressed in Spodoptera frugiperda Sf9 insect cells. The coding sequence was derived successively from peptide sequences of purified STOX protein. Furthermore, a recombinant oxidase with protoberberine dehydrogenase activity was obtained from a cDNA library of Argemone mexicana, a traditional medicinal plant that contains protoberberine alkaloids. The relationship of the two enzymes is discussed regarding their enzymatic activity, phylogeny and the alkaloid occurrence in the plants. Potential substrate binding and STOX-specific amino acid residues were identified based on sequence analysis and homology modeling.

Universal plant DNA barcode loci may not work in complex groups: a case study with Indian berberis species.

BACKGROUND: The concept of DNA barcoding for species identification has gained considerable momentum in animals because of fairly successful species identification using cytochrome oxidase I (COI). In plants, matK and rbcL have been proposed as standard barcodes. However, barcoding in complex genera is a challenging task. METHODOLOGY AND PRINCIPAL FINDINGS: We investigated the species discriminatory power of four reportedly most promising plant DNA barcoding loci (one from nuclear genome--ITS, and three from plastid genome--trnH-psbA, rbcL and matK) in species of Indian Berberis L. (Berberidaceae) and two other genera, Ficus L. (Moraceae) and Gossypium L. (Malvaceae). Berberis species were delineated using morphological characters. These characters resulted in a well resolved species tree. Applying both nucleotide distance and nucleotide character-based approaches, we found that none of the loci, either singly or in combinations, could discriminate the species of Berberis. ITS resolved all the tested species of Ficus and Gossypium and trnH-psbA resolved 82% of the tested species in Ficus. The highly regarded matK and rbcL could not resolve all the species. Finally, we employed amplified fragment length polymorphism test in species of Berberis to determine their relationships. Using ten primer pair combinations in AFLP, the data demonstrated incomplete species resolution. Further, AFLP analysis showed that there was a tendency of the Berberis accessions to cluster according to their geographic origin rather than species affiliation. CONCLUSIONS/SIGNIFICANCE: We reconfirm the earlier reports that the concept of universal barcode in plants may not work in a number of genera. Our results also suggest that the matK and rbcL, recommended as universal barcode loci for plants, may not work in all the genera of land plants. Morphological, geographical and molecular data analyses of Indian species of Berberis suggest probable reticulate evolution and thus barcode markers may not work in this case.

Molecular cloning and characterization of CYP80G2, a cytochrome P450 that catalyzes an intramolecular C-C phenol coupling of (S)-reticuline in magnoflorine biosynthesis, from cultured Coptis japonica cells.

Cytochrome P450s (P450) play a key role in oxidative reactions in plant secondary metabolism. Some of them, which catalyze unique reactions other than the standard hydroxylation, increase the structural diversity of plant secondary metabolites. In isoquinoline alkaloid biosyntheses, several unique P450 reactions have been reported, such as methylenedioxy bridge formation, intramolecular C-C phenol-coupling and intermolecular C-O phenol-coupling reactions. We report here the isolation and characterization of a C-C phenol-coupling P450 cDNA (CYP80G2) from an expressed sequence tag library of cultured Coptis japonica cells. Structural analysis showed that CYP80G2 had high amino acid sequence similarity to Berberis stolonifera CYP80A1, an intermolecular C-O phenol-coupling P450 involved in berbamunine biosynthesis. Heterologous expression in yeast indicated that CYP80G2 had intramolecular C-C phenol-coupling activity to produce (S)-corytuberine (aporphine-type) from (S)-reticuline (benzylisoquinoline type). Despite this intriguing reaction, recombinant CYP80G2 showed typical P450 properties: its C-C phenol-coupling reaction required NADPH and oxygen and was inhibited by a typical P450 inhibitor. Based on a detailed substrate-specificity analysis, this unique reaction mechanism and substrate recognition were discussed. CYP80G2 may be involved in magnoflorine biosynthesis in C. japonica, based on the fact that recombinant C. japonica S-adenosyl-L-methionine:coclaurine N-methyltransferase could convert (S)-corytuberine to magnoflorine.

Taxonomic and phytogeographic implications from ITS phylogeny in Berberis (Berberidaceae).

A phylogeny based on the internal transcribed spacer (ITS) sequences from 79 taxa representing much of the diversity of Berberis L. (four major groups and 22 sections) was constructed for the first time. The phylogeny was basically congruent with the previous classification schemes at higher taxonomic levels, such as groups and subgroups. A notable exception is the non-monophyly of the group Occidentales of compound-leaved Berberis (previously separated as Mahonia). At lower levels, however, most of previous sections and subsections were not evident especially in simple-leaved Berberis. Possible relationship between section Horridae (group Occidentales) and the simple-leaved Berberis clade implies paraphyly of the compound-leaved Berberis. A well-known South America-Old World (mainly Asia) disjunctive distribution pattern of the simple-leaved Berberis is explained by a vicariance event occurring in the Cretaceous period. The ITS phylogeny also suggests that a possible connection between the Asian and South American groups through the North American species ( Berberis canadensis or B. fendleri) is highly unlikely.