Naveen C. Bisht

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.

Evaluation of Candidate Reference Genes for Gene Expression Normalization in Brassica juncea Using Real Time Quantitative RT-PCR

The real time quantitative reverse transcription PCR (qRT-PCR) is becoming increasingly important to gain insight into function of genes. Given the increased sensitivity, ease and reproducibility of qRT-PCR, the requirement of suitable reference genes for normalization has become important and stringent. It is now known that the expression of internal control genes in living organism vary considerably during developmental stages and under different experimental conditions. For economically important Brassica crops, only a couple of reference genes are reported till date. In this study, expression stability of 12 candidate reference genes including ACT2, ELFA, GAPDH, TUA, UBQ9 (traditional housekeeping genes), ACP, CAC, SNF, TIPS-41, TMD, TSB and ZNF (new candidate reference genes), in a diverse set of 49 tissue samples representing different developmental stages, stress and hormone treated conditions and cultivars of Brassica juncea has been validated. For the normalization of vegetative stages the ELFA, ACT2, CAC and TIPS-41 combination would be appropriate whereas TIPS-41 along with CAC would be suitable for normalization of reproductive stages. A combination of GAPDH, TUA, TIPS-41 and CAC were identified as the most suitable reference genes for total developmental stages. In various stress and hormone treated samples, UBQ9 and TIPS-41 had the most stable expression. Across five cultivars of B. juncea, the expression of CAC and TIPS-41 did not vary significantly and were identified as the most stably expressed reference genes. This study provides comprehensive information that the new reference genes selected herein performed better than the traditional housekeeping genes. The selection of most suitable reference genes depends on the experimental conditions, and is tissue and cultivar-specific. Further, to attain accuracy in the results more than one reference genes are necessary for normalization.

Whole genome co-expression analysis of soybean cytochrome P450 genes identifies nodulation-specific P450 monooxygenases

BackgroundCytochrome P450 monooxygenases (P450s) catalyze oxidation of various substrates using oxygen and NAD(P)H. Plant P450s are involved in the biosynthesis of primary and secondary metabolites performing diverse biological functions. The recent availability of the soybean genome sequence allows us to identify and analyze soybean putative P450s at a genome scale. Co-expression analysis using an available soybean microarray and Illumina sequencing data provides clues for functional annotation of these enzymes. This approach is based on the assumption that genes that have similar expression patterns across a set of conditions may have a functional relationship.ResultsWe have identified a total number of 332 full-length P450 genes and 378 pseudogenes from the soybean genome. From the full-length sequences, 195 genes belong to A-type, which could be further divided into 20 families. The remaining 137 genes belong to non-A type P450s and are classified into 28 families. A total of 178 probe sets were found to correspond to P450 genes on the Affymetrix soybean array. Out of these probe sets, 108 represented single genes. Using the 28 publicly available microarray libraries that contain organ-specific information, some tissue-specific P450s were identified. Similarly, stress responsive soybean P450s were retrieved from 99 microarray soybean libraries. We also utilized Illumina transcriptome sequencing technology to analyze the expressions of all 332 soybean P450 genes. This dataset contains total RNAs isolated from nodules, roots, root tips, leaves, flowers, green pods, apical meristem, mock-inoculated and Bradyrhizobium japonicum-infected root hair cells. The tissue-specific expression patterns of these P450 genes were analyzed and the expression of a representative set of genes were confirmed by qRT-PCR. We performed the co-expression analysis on many of the 108 P450 genes on the Affymetrix arrays. First we confirmed that CYP93C5 (an isoflavone synthase gene) is co-expressed with several genes encoding isoflavonoid-related metabolic enzymes. We then focused on nodulation-induced P450s and found that CYP728H1 was co-expressed with the genes involved in phenylpropanoid metabolism. Similarly, CYP736A34 was highly co-expressed with lipoxygenase, lectin and CYP83D1, all of which are involved in root and nodule development.ConclusionsThe genome scale analysis of P450s in soybean reveals many unique features of these important enzymes in this crop although the functions of most of them are largely unknown. Gene co-expression analysis proves to be a useful tool to infer the function of uncharacterized genes. Our work presented here could provide important leads toward functional genomics studies of soybean P450s and their regulatory network through the integration of reverse genetics, biochemistry, and metabolic profiling tools. The identification of nodule-specific P450s and their further exploitation may help us to better understand the intriguing process of soybean and rhizobium interaction.

Fine mapping of loci involved with glucosinolate biosynthesis in oilseed mustard (Brassica juncea) using genomic information from allied species

Fine mapping of six seed glucosinolate QTL (J2Gsl1, J3Gsl2, J9Gsl3, J16Gsl4, J17Gsl5 and J3Gsl6) (Ramchiary et al. in Theor Appl Genet 116:77–85, 2007a) was undertaken by the candidate gene approach. Based on the DNA sequences from Arabidopsis and Brassica oleracea for the different genes involved in the aliphatic glucosinolate biosynthesis, candidate genes were amplified and sequenced from high to low glucosinolate Brassica juncea lines Varuna and Heera, respectively. Of the 20 paralogues identified, 17 paralogues belonging to six gene families were mapped to 12 of the 18 linkage groups of B. juncea genome. Co-mapping of candidate genes with glucosinolate QTL revealed that the candidate gene BjuA.GSL-ELONG.a mapped to the QTL interval of J2Gsl1, BjuA.GSL-ELONG.c, BjuA.GSL-ELONG.d and BjuA.Myb28.a mapped to the QTL interval of J3Gsl2, BjuA.GSL-ALK.a mapped to the QTL interval of J3Gsl6 and BjuB.Myb28.a mapped to the QTL interval of J17Gsl5. The QTL J9Gsl3 and J16Gsl4 did not correspond to any of the mapped candidate genes. The functionality and contribution of different candidate genes/QTL was assessed by allelic variation study using phenotypic data of 785 BC4DH lines. It was observed that BjuA.Myb28.a and J9Gsl3 contributed significantly to the base level glucosinolate production while J16Gsl4, probably GSL-PRO, BjuA.GSL-ELONG.a and BjuA.GSL-ELONG.c contributed to the C3, C4 and C5 elongation pathways, respectively. Three A genome QTL: J2Gsl1harbouring BjuA.GSL-ELONG.a,J3Gsl2 harbouring both BjuA.GSL-ELONG.c and BjuA.Myb28.a and J9Gsl3, possibly the ‘Bronowski genes’, were identified as most important loci for breeding low glucosinolate B. juncea. We observed two-step genetic control of seed glucosinolate in B. juncea mainly effected by these three A genome QTL. This study, therefore, provides clues to the genetic mechanism of ‘Bronowski genes’ controlling the glucosinolate trait and also provides efficient markers for marker-assisted introgression of low glucosinolate trait in B. juncea.

QTL analysis reveals context-dependent loci for seed glucosinolate trait in the oilseed Brassica juncea: importance of recurrent selection backcross scheme for the identification of ‘true’ QTL

Seed glucosinolate content in Brassica juncea is a complex quantitative trait. A recurrent selection backcross (RSB) method with a doubled haploid (DH) generation interspersing backcross generations was used for the introgression of low glucosinolate alleles from an east European gene pool B. juncea line, Heera into an Indian gene pool variety, Varuna. Phenotypic comparisons among the DH populations derived from early to advanced backcrosses revealed a shift in the mean values for various glucosinolates with the advancement of backcrossing, indicating a change in the selective values of the alleles with change in the genetic background due to the existence of epistasis and context dependencies. QTL mapping for various seed glucosinolates from early (F1DH) and advanced generation (BC4DH) populations confirmed the presence of epistasis and context dependency. The common QTL detected in both F1DH and BC4DH changed their R2 values from the former to the later generation. Some of the QTL detected in the F1DH became irrelevant in the BC4DH population. Further, new QTL were detected in the BC4DH population for various glucosinolates. A validation study on a population of low glucosinolate DH lines derived from all the backcross generations of the RSB breeding programme revealed that the QTL detected in BC4DH were the ‘true’ QTL. Using glucosinolate as an example, the study provides strong evidence for the importance of the RSB method for the identification of the ‘true’ QTL which would be significant for marker assisted introgression of a complex quantitative trait whose expression is influenced by epistatic interactions.

Four genes encoding MYB28, a major transcriptional regulator of the aliphatic glucosinolate pathway, are differentially expressed in the allopolyploid Brassica juncea

Glucosinolates are Capparales-specific secondary metabolites that have immense potential in human health and agriculture. Unlike Arabidopsis thaliana, our knowledge about glucosinolate regulators in the Brassica crops is sparse. In the current study, four MYB28 homologues were identified (BjuMYB28-1,-2,-3,-4) from the polyploid Brassica juncea, and the effects of allopolyploidization on the divergence of gene sequence, structure, function, and expression were assessed. The deduced protein sequences of the four BjuMYB28 genes showed 76.1–83.1% identity with the Arabidopsis MYB28. Phylogenetic analysis revealed that the four BjuMYB28 proteins have evolved via the hybridization and duplication processes forming the B. juncea genome (AABB) from B. rapa (AA) and B. nigra (BB), while retaining high levels of sequence conservation. Mutant complementation and over-expression studies in A. thaliana showed that all four BjuMYB28 genes encode functional MYB28 proteins and resulted in similar aliphatic glucosinolate composition and content. Detailed expression analysis using qRT-PCR assays and promoter-GUS lines revealed that the BjuMYB28 genes have both tissue- and cell-specific expression partitioning in B. juncea. The two B-genome origin BjuMYB28 genes had more abundant transcripts during the early stages of plant development than the A-genome origin genes. However, with the onset of the reproductive phase, expression levels of all four BjuMYB28 increased significantly, which may be necessary for producing and maintaining high amounts of aliphatic glucosinolates during the later stages of plant development. Taken together, our results suggest that the four MYB28 genes are differentially expressed and regulated in B. juncea to play discrete though overlapping roles in controlling aliphatic glucosinolate biosynthesis.

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.

Insights into structural and functional diversity of Dof (DNA binding with one finger) transcription factor.

AbstractMain conclusionThe structural, functional and in-silico studies of Dof transcription factor attempted so far reveals immense opportunity to analyze the plant genomes in terms of number of Dof genes and discuss in light of the evolution. The multiple functions of Dof genes needs to explored for crop improvement. Transcription factors play a very vital role in gene regulation at transcriptional level and are being extensively studied across phylas. In recent years, sequencing of plant genomes has led to genome-wide identification and characterizations of diverse types of plant-specific transcription factor gene family providing key insights into their structural and functional diversity. The DNA binding with one finger (Dof), a class belonging to C2H2-type zinc finger family proteins, is a plant-specific transcription factor having multiple roles such as seed maturation and germination, phytohormone and light-mediated regulation and plant responses to biotic and abiotic stresses. Dof proteins are present across plant lineage, from green algae to higher angiosperm, and represent a unique class of transcription factor having bifunctional binding activities, with both DNA and proteins, to regulate the complex transcriptional machinery in plant cells. The structural and functional diversity of the Dof transcription factor family along with the bioinformatics analysis highlighting the phylogeny of Dof families is reviewed in light of its importance in plant biotechnology for crop improvement.

Class-Specific Evolution and Transcriptional Differentiation of 14-3-3 Family Members in Mesohexaploid Brassica rapa

14-3-3s are highly conserved, multigene family proteins that have been implicated in modulating various biological processes. The presence of inherent polyploidy and genome complexity has limited the identification and characterization of 14-3-3 proteins from globally important Brassica crops. Through data mining of Brassica rapa, the model Brassica genome, we identified 21 members encoding 14-3-3 proteins namely, BraA.GRF14.a to BraA.GRF14.u. Phylogenetic analysis indicated that B. rapa contains both ε (epsilon) and non-ε 14-3-3 isoforms, having distinct intron-exon structural organization patterns. The non-ε isoforms showed lower divergence rate (Ks < 0.45) compared to ε protein isoforms (Ks > 0.48), suggesting class-specific divergence pattern. Synteny analysis revealed that mesohexaploid B. rapa genome has retained 1–5 orthologs of each Arabidopsis 14-3-3 gene, interspersed across its three fragmented sub-genomes. qRT-PCR analysis showed that 14 of the 21 BraA.GRF14 were expressed, wherein a higher abundance of non-ε transcripts was observed compared to the ε genes, indicating class-specific transcriptional bias. The BraA.GRF14 genes showed distinct expression pattern during plant developmental stages and in response to abiotic stress, phytohormone treatments, and nutrient deprivation conditions. Together, the distinct expression pattern and differential regulation of BraA.GRF14 genes indicated the occurrence of functional divergence of B. rapa 14-3-3 proteins during plant development and stress responses.

Two Chimeric Regulators of G-protein Signaling (RGS) Proteins Differentially Modulate Soybean Heterotrimeric G-protein Cycle

Background: The soybean genome encodes the most expanded plant heterotrimeric G-protein network reported to date. Results: Each Gα has distinct biochemical properties, and the RGS proteins have different GTPase-activating effects on each Gα. Conclusion: The core G-protein components, their interactions, and biochemical properties are conserved across phyla, but important mechanistic differences exist. Significance: This study provides insight into the complexity of plant G-protein networks. Heterotrimeric G-proteins and the regulator of G-protein signaling (RGS) proteins, which accelerate the inherent GTPase activity of Gα proteins, are common in animals and encoded by large gene families; however, in plants G-protein signaling is thought to be more limited in scope. For example, Arabidopsis thaliana contains one Gα, one Gβ, three Gγ, and one RGS protein. Recent examination of the Glycine max (soybean) genome reveals a larger set of G-protein-related genes and raises the possibility of more intricate G-protein networks than previously observed in plants. Stopped-flow analysis of GTP-binding and GDP/GTP exchange for the four soybean Gα proteins (GmGα1–4) reveals differences in their kinetic properties. The soybean genome encodes two chimeric RGS proteins with an N-terminal seven transmembrane domain and a C-terminal RGS box. Both GmRGS interact with each of the four GmGα and regulate their GTPase activity. The GTPase-accelerating activities of GmRGS1 and -2 differ for each GmGα, suggesting more than one possible rate of the G-protein cycle initiated by each of the Gα proteins. The differential effects of GmRGS1 and GmRGS2 on GmGα1–4 result from a single valine versus alanine difference. The emerging picture suggests complex regulation of the G-protein cycle in soybean and in other plants with expanded G-protein networks.