Arun Jagannath

Comparative mapping of Brassica juncea and Arabidopsis thaliana using Intron Polymorphism (IP) markers: homoeologous relationships, diversification and evolution of the A, B and C Brassica genomes

BackgroundExtensive mapping efforts are currently underway for the establishment of comparative genomics between the model plant, Arabidopsis thaliana and various Brassica species. Most of these studies have deployed RFLP markers, the use of which is a laborious and time-consuming process. We therefore tested the efficacy of PCR-based Intron Polymorphism (IP) markers to analyze genome-wide synteny between the oilseed crop, Brassica juncea (AABB genome) and A. thaliana and analyzed the arrangement of 24 (previously described) genomic block segments in the A, B and C Brassica genomes to study the evolutionary events contributing to karyotype variations in the three diploid Brassica genomes.ResultsIP markers were highly efficient and generated easily discernable polymorphisms on agarose gels. Comparative analysis of the segmental organization of the A and B genomes of B. juncea (present study) with the A and B genomes of B. napus and B. nigra respectively (described earlier), revealed a high degree of colinearity suggesting minimal macro-level changes after polyploidization. The ancestral block arrangements that remained unaltered during evolution and the karyotype rearrangements that originated in the Oleracea lineage after its divergence from Rapa lineage were identified. Genomic rearrangements leading to the gain or loss of one chromosome each between the A-B and A-C lineages were deciphered. Complete homoeology in terms of block organization was found between three linkage groups (LG) each for the A-B and A-C genomes. Based on the homoeology shared between the A, B and C genomes, a new nomenclature for the B genome LGs was assigned to establish uniformity in the international Brassica LG nomenclature code.ConclusionIP markers were highly effective in generating comparative relationships between Arabidopsis and various Brassica species. Comparative genomics between the three Brassica lineages established the major rearrangements, translocations and fusions pivotal to karyotype diversification between the A, B and C genomes of Brassica species. The inter-relationships established between the Brassica lineages vis-à-vis Arabidopsis would facilitate the identification and isolation of candidate genes contributing to traits of agronomic value in crop Brassicas and the development of unified tools for Brassica genomics.

The use of a Spacer DNA fragment insulates the tissue-specific expression of a cytotoxic gene (barnase) and allows high-frequency generation of transgenic male sterile lines in Brassica juncea L.

Male-sterile lines were generated in oilseed mustard (Brassica juncea) with a cytotoxic gene (barnase) in conjunction with either of two tapetum-specific promoters, TA29 and A9. Several transformation vectors based on different promoter and marker gene combinations were developed and tested for their efficacy in generating agronomically viable male-sterile lines. Use of strong constitutive promoters (e.g. CaMV 35S or its double-enhancer variant) to express the marker gene (bar) in barnase constructs generated male-sterile plants at an extremely low frequency with most plants showing abnormalities in vegetative morphology, poor female fertility, low seed germination frequencies and/or distortion in segregation ratios of transgenes. Such abnormalities were considerably reduced on using weaker promoters (e.g. nos) to drive the marker gene (nptII) in barnase constructs and could therefore be attributed to leaky expression of the barnase gene under enhancing effects of strong constitutive promoters. We show that the use of a Spacer DNA fragment between the barnase gene (driven by a tapetum-specific promoter) and the CaMV 35S promoter-driven bar gene insulates tissue-specific expression of the barnase gene over all developmental stages of transgenic plants and significantly enhances recovery of agronomically viable male-sterile lines. All TA29-barnase male-sterile lines containing the Spacer DNA fragment exhibited normal morphology, growth and seed set on backcrossing as observed for wild-type plants. Around 75% of single-copy events tested further also showed proper segregation of the marker gene/male-sterile phenotype among backcross progeny. Constructs based on the use of Spacer DNA fragments as insulators could be successfully used to alleviate limitations associated with transformation of plant systems using cytotoxic genes for development of agronomically viable male-sterile lines in crop plants and for cell/tissue ablation studies in general.

Identification and characterization of miRNAome in root, stem, leaf and tuber developmental stages of potato (Solanum tuberosum L.) by high-throughput sequencing

BackgroundMicroRNAs (miRNAs) are ubiquitous components of endogenous plant transcriptome. miRNAs are small, single-stranded and ~21 nt long RNAs which regulate gene expression at the post-transcriptional level and are known to play essential roles in various aspects of plant development and growth. Previously, a number of miRNAs have been identified in potato through in silico analysis and deep sequencing approach. However, identification of miRNAs through deep sequencing approach was limited to a few tissue types and developmental stages. This study reports the identification and characterization of potato miRNAs in three different vegetative tissues and four stages of tuber development by high throughput sequencing.ResultsSmall RNA libraries were constructed from leaf, stem, root and four early developmental stages of tuberization and subjected to deep sequencing, followed by bioinformatics analysis. A total of 89 conserved miRNAs (belonging to 33 families), 147 potato-specific miRNAs (with star sequence) and 112 candidate potato-specific miRNAs (without star sequence) were identified. The digital expression profiling based on TPM (Transcripts Per Million) and qRT-PCR analysis of conserved and potato-specific miRNAs revealed that some of the miRNAs showed tissue specific expression (leaf, stem and root) while a few demonstrated tuberization stage-specific expressions. Targets were predicted for identified conserved and potato-specific miRNAs, and predicted targets of four conserved miRNAs, miR160, miR164, miR172 and miR171, which are ARF16 (Auxin Response Factor 16), NAM (NO APICAL MERISTEM), RAP1 (Relative to APETALA2 1) and HAM (HAIRY MERISTEM) respectively, were experimentally validated using 5′ RLM-RACE (RNA ligase mediated rapid amplification of cDNA ends). Gene ontology (GO) analysis for potato-specific miRNAs was also performed to predict their potential biological functions.ConclusionsWe report a comprehensive study of potato miRNAs at genome-wide level by high-throughput sequencing and demonstrate that these miRNAs have tissue and/or developmental stage-specific expression profile. Also, predicted targets of conserved miRNAs were experimentally confirmed for the first time in potato. Our findings indicate the existence of extensive and complex small RNA population in this crop and suggest their important role in pathways involved in diverse biological processes, including tuber development.

Global insights into high temperature and drought stress regulated genes by RNA-Seq in economically important oilseed crop Brassica juncea

BackgroundBrassica juncea var. Varuna is an economically important oilseed crop of family Brassicaceae which is vulnerable to abiotic stresses at specific stages in its life cycle. Till date no attempts have been made to elucidate genome-wide changes in its transcriptome against high temperature or drought stress. To gain global insights into genes, transcription factors and kinases regulated by these stresses and to explore information on coding transcripts that are associated with traits of agronomic importance, we utilized a combinatorial approach of next generation sequencing and de-novo assembly to discover B. juncea transcriptome associated with high temperature and drought stresses.ResultsWe constructed and sequenced three transcriptome libraries namely Brassica control (BC), Brassica high temperature stress (BHS) and Brassica drought stress (BDS). More than 180 million purity filtered reads were generated which were processed through quality parameters and high quality reads were assembled de-novo using SOAPdenovo assembler. A total of 77750 unique transcripts were identified out of which 69,245 (89%) were annotated with high confidence. We established a subset of 19110 transcripts, which were differentially regulated by either high temperature and/or drought stress. Furthermore, 886 and 2834 transcripts that code for transcription factors and kinases, respectively, were also identified. Many of these were responsive to high temperature, drought or both stresses. Maximum number of up-regulated transcription factors in high temperature and drought stress belonged to heat shock factors (HSFs) and dehydration responsive element-binding (DREB) families, respectively. We also identified 239 metabolic pathways, which were perturbed during high temperature and drought treatments. Analysis of gene ontologies associated with differentially regulated genes forecasted their involvement in diverse biological processes.ConclusionsOur study provides first comprehensive discovery of B. juncea transcriptome under high temperature and drought stress conditions. Transcriptome resource generated in this study will enhance our understanding on the molecular mechanisms involved in defining the response of B. juncea against two important abiotic stresses. Furthermore this information would benefit designing of efficient crop improvement strategies for tolerance against conditions of high temperature regimes and water scarcity.

A Genome-Wide Perspective of miRNAome in Response to High Temperature, Salinity and Drought Stresses in Brassica juncea (Czern) L

Micro RNAs (miRNAs) are involved in diverse biological processes including adaptive response towards abiotic stresses. To unravel small RNAs and more specifically miRNAs that can potentially regulate determinants of abiotic stress tolerance, next generation sequencing of B. juncea seedlings subjected to high temperature, high salt and drought conditions was carried out. With the help of UEA sRNA workbench software package, 51 conserved miRNAs belonging to 30 miRNA families were identified. As there was limited genomic information available for B. juncea, we generated and assembled its genome sequence at a low coverage. Using the generated sequence and other publically available Brassica genomic/transcriptomic resources as mapping reference, 126 novel (not reported in any plant species) were discovered for the first time in B. juncea. Further analysis also revealed existence of 32 and 37 star sequences for conserved and novel miRNAs, respectively. The expression of selected conserved and novel miRNAs under conditions of different abiotic stresses was revalidated through universal TaqMan based real time PCR. Putative targets of identified conserved and novel miRNAs were predicted in B. rapa to gain insights into functional roles manifested by B. juncea miRNAs. Furthermore, SPL2-like, ARF17-like and a NAC domain containing protein were experimentally validated as targets of miR156, miR160 and miR164 respectively. Investigation of gene ontologies linked with targets of known and novel miRNAs forecasted their involvement in various biological functions.

A two gene – two promoter system for enhanced expression of a restorer gene (barstar) and development of improved fertility restorer lines for hybrid seed production in crop plants

We report in this study, strategies for enhancing and extending tissue-specific expression of a restorer gene (barstar) and their use in the development of improved fertility restorer lines for transgenic male sterile (barnase) lines for hybrid seed production. In the first strategy, a chimeric promoter was developed by combining regulatory elements of two overlapping tapetum-specific promoters (A9 and TA29) by placing a 275bp fragment of the A9 promoter (minus the TATA box) upstream to a functional 330bp fragment of the TA29 promoter. Analysis of chimeric promoter activity using the gusA gene showed that it enhanced transgene expression levels in anthers of transgenic Brassica juncea (Indian oilseed mustard) plants by two- to three-fold over those obtained with either promoter alone but retained the temporal expression profile of the TA29 promoter. In the second strategy, the A9 and TA29 promoters were used to independently express two copies of the wild type barstar gene sequence within the same T-DNA borders. Alternatively, one of the promoters was used to express the wild type sequence and the other was used to drive expression of a codon-modified version of the barstar gene designed for enhanced and stable expression in dicot plants. The efficiencies of these constructs for fertility restoration was determined by retransforming a male sterile TA29-barnase line of B. juncea and analyzing the frequency and pollen viability of male fertile (restored) plants thus obtained. A significantly higher frequency of restored plants was obtained with constructs containing two transcription units of the barstar gene (92.9% for constructs with two copies of the wild type gene and 89.8% for constructs using the wild type and codon-modified sequences) as compared to that obtained with the chimeric promoter-barstar construct (77.8%) or the TA29-barstar construct (65.6%). Pollen viability assays on restored plants showed that the desirable extent of restoration is obtained in a higher proportion of male fertile plants generated using the construct containing a combination of the wild type and modified barstar genes. These observations were corroborated by genetic analysis of F1 progeny of crosses between various barstar lines and four different male sterile barnase lines, three of which could not be restored in earlier studies using TA29-barstar constructs. The simultaneous use of the wild type and modified barstar genes with the TA29 and A9 promoters was found to be the most effective method for fertility restoration with efficient restorers being identified for all the male sterile barnase lines tested, including those which could not be restored earlier. The two gene-two promoter strategy could be successfully deployed for enhancing and extending the tissue-specific as well as constitutive expression of other agronomically important transgenes.

Retransformation of a male sterile barnase line with the barstar gene as an efficient alternative method to identify male sterile–restorer combinations for heterosis breeding

We report in this study, an improved method for identifying male sterile–restorer combinations using the barnase–barstar system of pollination control for heterosis breeding in crop plants, as an alternative to the conventional line × tester cross method. In this strategy, a transgenic male sterile barnase line was retransformed with appropriate barstar constructs. Double transformants carrying both the barnase and barstar genes were identified and screened for their male fertility status. Using this strategy, 66–90% of fertile retransformants (restored events) were obtained in Brassica juncea using two different barstar constructs. Restored events were analysed for their pollen viability and copy number of the barstar gene. Around 90% of the restored events showed high pollen viability and ∼30% contained single copy integrations of the barstar gene. These observations were significantly different from those made in our earlier studies using line (barnase) × tester (barstar) crosses, wherein only two viable male sterile–restorer combinations were identified by screening 88 different cross-combinations. The retransformation strategy not only generated several independent restorers for a given male sterile line from a single transformation experiment but also identified potential restorers in the T0 generation itself leading to significant savings in time, cost and labour. Single copy restored plants with high pollen viability were selfed to segregate male sterile (barnase) and restorer (barstar) lines in the T1 progeny which could subsequently be diversified into appropriate combiners for heterosis breeding. This strategy will be particularly useful for crop plants where poor transformation frequencies and/or lengthy transformation protocols are a major limitation.

Assessment of Genetic Diversity and Population Structure in a Global Reference Collection of 531 Accessions of Carthamus tinctorius L. (Safflower) Using AFLP Markers

Carthamus tinctorius L. (safflower) is an important oilseed crop that is cultivated in several countries. The present study investigates the genetic diversity and population structure of 531 safflower accessions from 43 countries representing all safflower growing regions of the world. Diversity analysis was performed using ten informative EcoRI/MseI amplified fragment length polymorphism primer pairs that were identified by screening 150 primer combinations. The selected primer pairs generated 381 fragments of which 157 were polymorphic among the analyzed accessions. The genetic diversity indices obtained for the entire collection (I = 0.4536, H = 0.2955) indicated high levels of molecular variability. The distance-based, neighbor-joining method classified the accessions into six clusters with internal subgroupings that were in consonance with 19 clusters obtained using Bayesian model-based BAPS analysis. Clusters obtained through STRUCTURE analysis (at K = 4) could not be correlated with their geographically diverse origins, while BAPS analysis (at K = 19) revealed geographical delineation with low admixture levels among most of the studied accessions. Accessions from Far East and Egypt clustered in distinct groups, indicating conserved nature of their gene pools. The Near East and Iran–Afghanistan regions were collectively found to harbor maximum diversity in accordance with earlier reports. Accessions from the Indian subcontinent showed substantial diversity that was previously undetected. The American accessions showed low molecular variability in contrast to earlier studies. Genetic sub-structuring within gene pools and inter-relationships between accessions belonging to different regional pools was also observed. To the best of our knowledge, this is the first comprehensive study of existing genetic variability in a large collection of safflower germplasm with a global distribution, which provides a more accurate representation of genetic structuring in the crop. This information will facilitate selection of elite genotypes for broadening the genetic base of various breeding programs in safflower.

Genome-wide identification and characterization of miRNAome from tomato (Solanum lycopersicum) roots and root-knot nematode (Meloidogyne incognita) during susceptible interaction

Root-knot nematodes (RKNs, Meloidogyne spp.) are the most damaging plant parasites causing severe losses to crop production. The present study reports genome-wide identification and characterization of both tomato and RKN miRNAs simultaneously from RKN-infected susceptible tomato roots using high-throughput sequencing technique. RNAseq data from 11 small RNA libraries derived from 5 disease development stages identified 281 novel miRNAs of tomato in addition to 52 conserved and 4 variants of conserved miRNAs. Additionally, the same set of RNAseq data identified 38 conserved and 290 novel RKN miRNAs. Both tomato and RKN miRNAs showed differential expression at 5 stages of disease development based on digital expression profiles. In tomato, further validation through qRT-PCR confirmed that majority of miRNAs were significantly upregulated during susceptible response whereas downregulated during resistance response. The predicted targets of 8 conserved and 1 novel miRNAs were validated through 5’RLM-RACE. A negative correlation between expression profiles of a few conserved miRNAs (miR156, miR159, miR164 and miR396) and their targets (SBP, GAMYB-like, NAC and GRF1 transcription factor) was confirmed. A novel Sly_miRNA996 also showed a negative correlation with its target MYB-like transcription factor. These results indicate that the conserved and novel tomato miRNAs are involved in regulating developmental changes in host root during RKN infection. In RKN, the targets of conserved miRNAs were also predicted and a few of their predicted target genes are known to be involved in nematode parasitism. Further, the potential roles of both tomato and RKN miRNAs have been discussed.

De novo transcriptome profiling of cold-stressed siliques during pod filling stages in Indian mustard (Brassica juncea L.)

Low temperature is a major abiotic stress that impedes plant growth and development. Brassica juncea is an economically important oil seed crop and is sensitive to freezing stress during pod filling subsequently leading to abortion of seeds. To understand the cold stress mediated global perturbations in gene expression, whole transcriptome of B. juncea siliques that were exposed to sub-optimal temperature was sequenced. Manually self-pollinated siliques at different stages of development were subjected to either short (6 h) or long (12 h) durations of chilling stress followed by construction of RNA-seq libraries and deep sequencing using Illuminas NGS platform. De-novo assembly of B. juncea transcriptome resulted in 133,641 transcripts, whose combined length was 117 Mb and N50 value was 1428 bp. We identified 13,342 differentially regulated transcripts by pair-wise comparison of 18 transcriptome libraries. Hierarchical clustering along with Spearman correlation analysis identified that the differentially expressed genes segregated in two major clusters representing early (5–15 DAP) and late stages (20–30 DAP) of silique development. Further analysis led to the discovery of sub-clusters having similar patterns of gene expression. Two of the sub-clusters (one each from the early and late stages) comprised of genes that were inducible by both the durations of cold stress. Comparison of transcripts from these clusters led to identification of 283 transcripts that were commonly induced by cold stress, and were referred to as “core cold-inducible” transcripts. Additionally, we found that 689 and 100 transcripts were specifically up-regulated by cold stress in early and late stages, respectively. We further explored the expression patterns of gene families encoding for transcription factors (TFs), transcription regulators (TRs) and kinases, and found that cold stress induced protein kinases only during early silique development. We validated the digital gene expression profiles of selected transcripts by qPCR and found a high degree of concordance between the two analyses. To our knowledge this is the first report of transcriptome sequencing of cold-stressed B. juncea siliques. The data generated in this study would be a valuable resource for not only understanding the cold stress signaling pathway but also for introducing cold hardiness in B. juncea.