Data mining of transcriptional biomarkers at different cotton fiber developmental stages.

Advancement of the gene expression study provides comprehensive information on pivotal genes at different cotton fiber development stages. For the betterment of cotton fiber yield and their quality, genetic improvement is a major target point for the cotton community. Therefore, various studies were carried out to understand the transcriptional machinery of fiber leading to the detailed integrative as well as innovative study. Through data mining and statistical approaches, we identified and validated the transcriptional biomarkers for staged specific differentiation of fiber. With the unique mapping read matrix of ~ 200 cotton transcriptome data and sequential statistical analysis, we identified several important genes that have a deciding and specific role in fiber cell commitment, initiation and elongation, or secondary cell wall synthesis stage. Based on the importance score and validation analysis, IQ domain 26, Aquaporin, Gibberellin regulated protein, methionine gamma lyase, alpha/beta hydrolases, and HAD-like superfamily have shown the specific and determining role for fiber developmental stages. These genes are represented as transcriptional biomarkers that provide a base for molecular characterization for cotton fiber development which will ultimately determine the high yield.

Dynamicity of histone H3K27ac and H3K27me3 modifications regulate the cold-responsive gene expression in Oryza sativa L. ssp. indica.

Cultivated in tropical and subtropical regions, Oryza sativa L. ssp. indica is largely affected by cold-stress, especially at the seedling stage. The present model of the stress-responsive regulatory network in plants entails the role of genetic and epigenetic factors in stress-responsive gene expression. Despite extensive transcriptomic studies, the regulation of various epigenetic factors in plants cold-stress response is less explored. The present study addresses the effect of genome-wide changes of H3K27 modifications on gene expression in IR64 rice, during cold-stress. Our results suggest a positive correlation between the changes in H3K27 modifications and stress-responsive gene activation in indica rice. Cold-induced enrichment of H3K27 acetylation promotes nucleosomal rearrangement, thereby facilitating the accessibility of the transcriptional machinery at the stress-responsive loci for transcription activation. Although H3K27ac exhibits uniform distribution throughout the loci of enriched genes; occupancy of H3K27me3 is biased to intergenic regions. Integration of the ChIP-seq data with transcriptome indicated that upregulation of stress-responsive TFs, photosynthesis-TCA-related, water-deficit genes, redox and JA signalling components, was associated with differential changes of H3K27ac and H3K27me3 levels. Furthermore, cold-induced upregulation of histone acetyltransferases and downregulation of DNA methyltransferases was noted through the antagonistic switch of H3K27ac and H3K27me3. Moreover, motif analysis of H3K27ac and H3K27me3 enriched regions are associated with putative stress responsive transcription factors binding sites, GAGA element and histone H3K27demethylase. Collectively our analysis suggests that differential expression of various chromatin and DNA modifiers coupled with increased H3K27ac and depleted H3K27me3 increases DNA accessibility, thereby promoting transcription of the cold-responsive genes in indica rice.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits.

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

Identification of genome-wide targets and DNA recognition sequence of the Arabidopsis HMG-box protein AtHMGB15 during cold stress response.

AtHMGB15 belongs to a group of ARID-HMG proteins which are plant specific. The presence of two known DNA binding domains: AT rich interacting domain (ARID) and High Mobility Group (HMG)-box, in one polypeptide, makes this protein intriguing. Although proteins containing individual HMG and ARID domains have been characterized, not much is known about the role of ARID-HMG proteins. Promoter analysis of AtHMGB15 showed the presence of various stress responsive cis regulatory elements along with MADS-box containing transcription factors. Our result shows that the expression of AtHMGB15 increased significantly upon application of cold stress. Using ChIP-chip approach, we have identified 6128 and 4689 significantly enriched loci having AtHMGB15 occupancy under control and cold stressed condition respectively. GO analysis shows genes belonging to abiotic stress response, cold response and root development were AtHMGB15 targets during cold stress. DNA binding and footprinting assays further identified A(A/C)--ATA---(A/T)(A/T) as AtHMGB15 binding motif. The enriched probe distribution in both control and cold condition shows a bias of AtHMGB15 binding towards the transcribed (gene body) region. Further, the expression of cold stress responsive genes decreased in athmgb15 knockout plants compared to wild-type. Taken together, binding enrichment of AtHMGB15 to the promoter and upstream to stress loci suggest an unexplored role of the protein in stress induced transcription regulation.

PLANET-SNP pipeline: PLants based ANnotation and Establishment of True SNP pipeline.

Acute prediction of SNPs (Single Nucleotide Polymorphisms) from high throughput sequencing data is a challenging problem, having potential to explore possible variation within plants species. For the extraction of profitable information from bulk of data, machine learning (ML) could lead to development of accurate model based on the learning of prior information. We performed state of art, in-depth learning on six different plant species. Comparative evaluation of five different algorithms showed that Random Forest substantially outperformed in selection of potential SNPs, with markedly improved prediction accuracy via 10-fold cross validation technique and integrated in system known as PLANET-SNP. We present the accurate method to extract the potential SNPs with user specific customizable parameters. It will facilitate the identification of efficient and functional SNPs in most easy and intuitive way. PLANET-SNP pipeline is very flexible in terms of data input and output formats. PLANET-SNP Pipeline is available at http://www.ncgd.nbri.res.in/PLANET-SNP-Pipeline.aspx.

In Silico identification of SNP diversity in cultivated and wild tomato species: insight from molecular simulations.

Single Nucleotide Polymorphisms (SNPs), an important source of genetic variations, are often used in crop improvement programme. The present study represented comprehensive In silico analysis of nucleotide polymorphisms in wild (Solanum habrochaites) and cultivated (Solanum lycopersicum) species of tomato to explore the consequence of substitutions both at sequence and structure level. A total of 8978 SNPs having Ts/Tv (Transition/Transversion) ratio 1.75 were identified from the Expressed Sequence Tag (EST) and Next Generation Sequence (NGS) data of both the species available in public databases. Out of these, 1838 SNPs were non-synonymous and distributed in 988 protein coding genes. Among these, 23 genes containing 96 SNPs were involved in traits markedly different between the two species. Furthermore, there were 28 deleterious SNPs distributed in 27 genes and a few of these genes were involved in plant pathogen interaction and plant hormone pathways. Molecular docking and simulations of several selected proteins showed the effect of SNPs in terms of compactness, conformation and interaction ability. Observed SNPs exhibited various types of motif binding effects due to nucleotide changes. SNPs that provide the evidence of differential motif binding and interaction behaviour could be effectively used for the crop improvement program.

PanCoreGen - Profiling, detecting, annotating protein-coding genes in microbial genomes.

A large amount of genomic data, especially from multiple isolates of a single species, has opened new vistas for microbial genomics analysis. Analyzing the pan-genome (i.e. the sum of genetic repertoire) of microbial species is crucial in understanding the dynamics of molecular evolution, where virulence evolution is of major interest. Here we present PanCoreGen - a standalone application for pan- and core-genomic profiling of microbial protein-coding genes. PanCoreGen overcomes key limitations of the existing pan-genomic analysis tools, and develops an integrated annotation-structure for a species-specific pan-genomic profile. It provides important new features for annotating draft genomes/contigs and detecting unidentified genes in annotated genomes. It also generates user-defined group-specific datasets within the pan-genome. Interestingly, analyzing an example-set of Salmonella genomes, we detect potential footprints of adaptive convergence of horizontally transferred genes in two human-restricted pathogenic serovars - Typhi and Paratyphi A. Overall, PanCoreGen represents a state-of-the-art tool for microbial phylogenomics and pathogenomics study.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp.

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

Comparative interactions of withanolides and sterols with two members of sterol glycosyltransferases from Withania somnifera.

BACKGROUND: Sterol glycosyltransferases (SGTs) are ubiquitous but one of the most diverse group of enzymes of glycosyltransferases family. Members of this family modulate physical and chemical properties of secondary plant products important for various physiological processes. The role of SGTs has been demonstrated in the biosynthesis of pharmaceutically important molecules of medicinal plants like Withania somnifera. RESULTS: Analysis suggested conserved behaviour and high similarity in active sites of WsSGTs with other plant GTs. Substrate specificity of WsSGTs were analysed through docking performance of WsSGTs with different substrates (sterols and withanolides). Best docking results of WsSGTL1 in the form of stable enzyme-substrate complex having lowest binding energies were obtained with brassicasterol, transandrosteron and WsSGTL4 with solasodine, stigmasterol and 24-methylene cholesterol. CONCLUSION: This study reveals topological characters and conserved nature of two SGTs from W. somnifera (WsSGTs) i.e. WsSGTL1 and WsSGTL4. However, besides being ubiquitous in nature and with broad substrate specificity, difference between WsSGTL1 and WsSGTL4 is briefly described by difference in stability (binding energy) of enzyme-substrate complexes through comparative docking.

Comparative Transcriptional Profiling of Contrasting Rice Genotypes Shows Expression Differences during Arsenic Stress.

Accumulation of arsenic (As) in rice (Oryza sativa L.) grain is a serious concern worldwide. Long-term exposure to As affects nutritional status in rice grain and is associated with higher rates of skin, bladder, and lung cancers, and heart disease. Genotypic variations in rice for As accumulation or tolerance are prevalent and are regulated by genetic and environmental factors. To understand molecular networks involved in As accumulation, genome-wide expression analysis was performed in roots of low- and high-As accumulating rice genotypes (LARGs and HARGs). Six rice genotypes with contrasting As accumulation potential and tolerance were used in this study. Genome-wide expression analysis suggested their differential response against As stress. This study suggests up- and downregulation of a number of unique genes involved in various pathways and biological processes in response to As stress in rice genotypes. A comparison of gene expression profiles, principal component analysis, and K-means clustering suggests that an independent pathway is operating during As stress tolerance or accumulation in contrasting genotypes. It was also observed that the differential behavior of aus genotype, Nayanmoni, from other LARGs might be due to its different genetic background. Cis-motif profiling of As-induced coexpressed genes in diverse rice genotypes led to the identification of unique cis-motifs present in differentially expressed genes. This study suggests that the genetic mechanism regulating the differential As accumulation in different genotypes may not be dependent on gene expression at the transcriptional level. However, many genes identified in this study can be analyzed and used for marker-trait associations related to As accumulation in diverse genotypes around the world.

Transcriptome analysis of ripe and unripe fruit tissue of banana identifies major metabolic networks involved in fruit ripening process.

BACKGROUND: Banana is one of the most important crop plants grown in the tropics and sub-tropics. It is a climacteric fruit and undergoes ethylene dependent ripening. Once ripening is initiated, it proceeds at a fast rate making postharvest life short, which can result in heavy economic losses. During the fruit ripening process a number of physiological and biochemical changes take place and thousands of genes from various metabolic pathways are recruited to produce a ripe and edible fruit. To better understand the underlying mechanism of ripening, we undertook a study to evaluate global changes in the transcriptome of the fruit during the ripening process. RESULTS: We sequenced the transcriptomes of the unripe and ripe stages of banana (Musa accuminata; Dwarf Cavendish) fruit. The transcriptomes were sequenced using a 454 GSFLX-Titanium platform that resulted in more than 7,00,000 high quality (HQ) reads. The assembly of the reads resulted in 19,410 contigs and 92,823 singletons. A large number of the differentially expressed genes identified were linked to ripening dependent processes including ethylene biosynthesis, perception and signalling, cell wall degradation and production of aromatic volatiles. In the banana fruit transcriptomes, we found transcripts included in 120 pathways described in the KEGG database for rice. The members of the expansin and xyloglucan transglycosylase/hydrolase (XTH) gene families were highly up-regulated during ripening, which suggests that they might play important roles in the softening of the fruit. Several genes involved in the synthesis of aromatic volatiles and members of transcription factor families previously reported to be involved in ripening were also identified. CONCLUSIONS: A large number of differentially regulated genes were identified during banana fruit ripening. Many of these are associated with cell wall degradation and synthesis of aromatic volatiles. A large number of differentially expressed genes did not align with any of the databases and might be novel genes in banana. These genes can be good candidates for future studies to establish their role in banana fruit ripening. The datasets developed in this study will help in developing strategies to manipulate banana fruit ripening and reduce post harvest losses.

Genome-wide identification and expression analysis of the mitogen-activated protein kinase gene family from banana suggest involvement of specific members in different stages of fruit ripening.

Mitogen-activated protein kinases (MAPKs) are important components of the tripartite mitogen-activated protein kinase signaling cascade and play an important role in plant growth and development. Although members of the MAPK gene family have been identified in model plants, little information is available regarding this gene family in fruit crops. In this study, we carried out a computational analysis using the Musa Genome database to identify members of the MAPK gene family in banana, an economically important crop and the most popular fruit worldwide. Our analysis identified 25 members of the MAP kinase (MAPK or MPK) gene family. Phylogenetic analyses of MPKs in Arabidopsis, Oryza, and Populus have classified these MPKs into four subgroups. The presence of conserved domains in the deduced amino acid sequences, phylogeny, and genomic organization strongly support their identity as members of the MPK gene family. Expression analysis during ethylene-induced banana fruit ripening suggests the involvement of several MPKs in the ethylene signal transduction pathway that are necessary for banana fruit ripening. Analysis of the cis-regulatory elements in the promoter regions and the involvement of the identified MPKs in various cellular processes, as analyzed using Pathway Studio, suggest a role for the banana MPK gene family in diverse functions related to growth, development, and the stress response. This report is the first concerning the identification of members of a gene family and the elucidation of their role in various processes using the Musa Genome database.

Distinct role of core promoter architecture in regulation of light-mediated responses in plant genes.

In the present study, we selected four distinct classes of light-regulated promoters. The light-regulated promoters can be distinctly grouped into either TATA-box-containing or TATA-less (initiator-containing) promoters. Further, using either native promoters or their swapped versions of core promoter elements, we established that TATA-box and Inr (Initiator) elements have distinct mechanisms which are involved in light-mediated regulation, and these elements are not swappable. We identified that mutations in either functional TATA-box or Inr elements lead to the formation of nucleosomal structure. The nucleotide diversity in either the TATA-box or Inr element in Arabidopsis ecotypes proposes that the nucleotide variation in core promoters can alter the gene expression. We show that motif overrepresentation in light-activated promoters encompasses different specific regulatory motifs present downstream of TSS (transcription start site), and this might serve as a key factor in regulating light promoters which are parallel with these elements. Finally, we conclude that the TATA-box or Inr element does not act in isolation, but our results clearly suggests the probable involvement of other distinct core promoter elements in concurrence with the TATA-box or Inr element to impart selectivity to light-mediated transcription.

Insights from the draft genome of Paenibacillus lentimorbus NRRL B-30488, a promising plant growth promoting bacterium.

Paenibacillus lentimorbus NRRL B-30488, a plant growth-promoting bacterium was isolated from Sahiwal cow's milk. The strain shows antagonism against phytopathogens, Fusarium oxysporum f. sp. ciceri and Alternaria solani. Its genome contains gene clusters involved in nonribosomal synthesis of secondary metabolites involved in antimicrobial activities. The genome sequence of P. lentimorbus NRRL B-30488 provides the genetic basis for application of this bacterial strain in plant growth promotion, plant protection and degradation of organic pollutants.

Large-scale resource development in Gossypium hirsutum L. by 454 sequencing of genic-enriched libraries from six diverse genotypes.

The sequence information has been proved to be an essential genomic resource in case of crop plants for their genetic improvement and better utilization by humans. To dissect the Gossypium hirsutum genome for large-scale development of genomic resources, we adopted hypomethylated restriction-based genomic enrichment strategy to sequence six diverse genotypes. Approximately 5.2-Gb data (more than 18.36 million reads) was generated which, after assembly, represents nearly 1.27-Gb genomic sequences. We predicted a total of 93,363 gene models (21,399 full length) and identified 35,923 gene models which were validated against already sequenced plant genomes. A total of 1,093 transcription factor-encoding genes, 3,135 promoter sequences and 78 miRNA (including 17 newly identified in Gossypium) were predicted. We identified significant no. of molecular markers including 47,093 novel simple sequence repeats and 66,364 novel single nucleotide polymorphisms. In addition, we developed NBRI-Comprehensive Cotton Genomics database, a web resource to provide access of cotton-related genomic resources developed at NBRI. This study contributes considerable amount of genomic resources and suggests a potential role of genic-enriched sequencing in genomic resource development for orphan crop plants.

Comparative transcriptome analysis of Gossypium hirsutum L. in response to sap sucking insects: aphid and whitefly.

BACKGROUND: Cotton (Gossypium hirsutum L.) is a major fiber crop that is grown worldwide; it faces extensive damage from sap-sucking insects, including aphids and whiteflies. Genome-wide transcriptome analysis was performed to understand the molecular details of interaction between Gossypium hirsutum L. and sap-sucking pests, namely Aphis gossypii (Aphid) and Bemisia tabacci (Whiteflies). Roche's GS-Titanium was used to sequence transcriptomes of cotton infested with aphids and whiteflies for 2 h and 24 h. RESULTS: A total of 100935 contigs were produced with an average length of 529 bp after an assembly in all five selected conditions. The Blastn of the non-redundant (nr) cotton EST database resulted in the identification of 580 novel contigs in the cotton plant. It should be noted that in spite of minimal physical damage caused by the sap-sucking insects, they can change the gene expression of plants in 2 h of infestation; further change in gene expression due to whiteflies is quicker than due to aphids. The impact of the whitefly 24 h after infestation was more or less similar to that of the aphid 2 h after infestation. Aphids and whiteflies affect many genes that are regulated by various phytohormones and in response to microbial infection, indicating the involvement of complex crosstalk between these pathways. The KOBAS analysis of differentially regulated transcripts in response to aphids and whiteflies indicated that both the insects induce the metabolism of amino acids biosynthesis specially in case of whiteflies infestation at later phase. Further we also observed that expression of transcript related to photosynthesis specially carbon fixation were significantly influenced by infestation of Aphids and Whiteflies. CONCLUSIONS: A comparison of different transcriptomes leads to the identification of differentially and temporally regulated transcripts in response to infestation by aphids and whiteflies. Most of these differentially expressed contigs were related to genes involved in biotic, abiotic stresses and enzymatic activities related to hydrolases, transferases, and kinases. The expression of some marker genes such as the overexpressors of cationic peroxidase 3, lipoxygenase I, TGA2, and non-specific lipase, which are involved in phytohormonal-mediated plant resistance development, was suppressed after infestation by aphids and whiteflies, indicating that insects suppressed plant resistance in order to facilitate their infestation. We also concluded that cotton shares several pathways such as phagosomes, RNA transport, and amino acid metabolism with Arabidopsis in response to the infestation by aphids and whiteflies.

Development and characterization of genomic and expressed SSRs for levant cotton (Gossypium herbaceum L.).

Four microsatellite-enriched genomic libraries for CA(15), GA(15), AAG(8) and ATG(8) repeats and transcriptome sequences of five cDNA libraries of Gossypium herbaceum were explored to develop simple sequence repeat (SSR) markers. A total of 428 unique clones from repeat enriched genomic libraries were mined for 584 genomic SSRs (gSSRs). In addition, 99,780 unigenes from transcriptome sequencing were explored for 8,900 SSR containing sequences with 12,471 expressed SSRs. The present study adds 1,970 expressed SSRs and 263 gSSRs to the public domain for the use of genetic studies of cotton. When 150 gSSRs and 50 expressed SSRs were tested on a panel of four species of cotton, 68 gSSRs and 12 expressed SSRs revealed polymorphism. These 200 SSRs were further deployed on 15 genotypes of levant cotton for the genetic diversity assessment. This is the first report on the successful use of repeat enriched genomic library and expressed sequence database for microsatellite markers development in G. herbaceum.

Indian genetic disease database.

Indians, representing about one-sixth of the world population, consist of several thousands of endogamous groups with strong potential for excess of recessive diseases. However, no database is available on Indian population with comprehensive information on the diseases common in the country. To address this issue, we present Indian Genetic Disease Database (IGDD) release 1.0 (http://www.igdd.iicb.res.in)--an integrated and curated repository of growing number of mutation data on common genetic diseases afflicting the Indian populations. Currently the database covers 52 diseases with information on 5760 individuals carrying the mutant alleles of causal genes. Information on locus heterogeneity, type of mutation, clinical and biochemical data, geographical location and common mutations are furnished based on published literature. The database is currently designed to work best with Internet Explorer 8 (optimal resolution 1440 × 900) and it can be searched based on disease of interest, causal gene, type of mutation and geographical location of the patients or carriers. Provisions have been made for deposition of new data and logistics for regular updation of the database. The IGDD web portal, planned to be made freely available, contains user-friendly interfaces and is expected to be highly useful to the geneticists, clinicians, biologists and patient support groups of various genetic diseases.

Transcriptomic and metabolomic shifts in rice roots in response to Cr (VI) stress.

BACKGROUND: Widespread use of chromium (Cr) contaminated fields due to careless and inappropriate management practices of effluent discharge, mostly from industries related to metallurgy, electroplating, production of paints and pigments, tanning, and wood preservation elevates its concentration in surface soil and eventually into rice plants and grains. In spite of many previous studies having been conducted on the effects of chromium stress, the precise molecular mechanisms related to both the effects of chromium phytotoxicity, the defense reactions of plants against chromium exposure as well as translocation and accumulation in rice remain poorly understood. RESULTS: Detailed analysis of genome-wide transcriptome profiling in rice root is reported here, following Cr-plant interaction. Such studies are important for the identification of genes responsible for tolerance, accumulation and defense response in plants with respect to Cr stress. Rice root metabolome analysis was also carried out to relate differential transcriptome data to biological processes affected by Cr (VI) stress in rice. To check whether the Cr-specific motifs were indeed significantly over represented in the promoter regions of Cr-responsive genes, occurrence of these motifs in whole genome sequence was carried out. In the background of whole genome, the lift value for these 14 and 13 motifs was significantly high in the test dataset. Though no functional role has been assigned to any of the motifs, but all of these are present as promoter motifs in the Database of orthologus promoters. CONCLUSION: These findings clearly suggest that a complex network of regulatory pathways modulates Cr-response of rice. The integrated matrix of both transcriptome and metabolome data after suitable normalization and initial calculations provided us a visual picture of the correlations between components. Predominance of different motifs in the subsets of genes suggests the involvement of motif-specific transcription modulating proteins in Cr stress response of rice.

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.