Kishor Gaikwad

Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea [Cajanus cajan (L.) Millspaugh]

BackgroundPigeonpea [Cajanus cajan (L.) Millspaugh], one of the most important food legumes of semi-arid tropical and subtropical regions, has limited genomic resources, particularly expressed sequence based (genic) markers. We report a comprehensive set of validated genic simple sequence repeat (SSR) markers using deep transcriptome sequencing, and its application in genetic diversity analysis and mapping.ResultsIn this study, 43,324 transcriptome shotgun assembly unigene contigs were assembled from 1.696 million 454 GS-FLX sequence reads of separate pooled cDNA libraries prepared from leaf, root, stem and immature seed of two pigeonpea varieties, Asha and UPAS 120. A total of 3,771 genic-SSR loci, excluding homopolymeric and compound repeats, were identified; of which 2,877 PCR primer pairs were designed for marker development. Dinucleotide was the most common repeat motif with a frequency of 60.41%, followed by tri- (34.52%), hexa- (2.62%), tetra- (1.67%) and pentanucleotide (0.76%) repeat motifs. Primers were synthesized and tested for 772 of these loci with repeat lengths of ≥18 bp. Of these, 550 markers were validated for consistent amplification in eight diverse pigeonpea varieties; 71 were found to be polymorphic on agarose gel electrophoresis. Genetic diversity analysis was done on 22 pigeonpea varieties and eight wild species using 20 highly polymorphic genic-SSR markers. The number of alleles at these loci ranged from 4-10 and the polymorphism information content values ranged from 0.46 to 0.72. Neighbor-joining dendrogram showed distinct separation of the different groups of pigeonpea cultivars and wild species. Deep transcriptome sequencing of the two parental lines helped in silico identification of polymorphic genic-SSR loci to facilitate the rapid development of an intra-species reference genetic map, a subset of which was validated for expected allelic segregation in the reference mapping population.ConclusionWe developed 550 validated genic-SSR markers in pigeonpea using deep transcriptome sequencing. From these, 20 highly polymorphic markers were used to evaluate the genetic relationship among species of the genus Cajanus. A comprehensive set of genic-SSR markers was developed as an important genomic resource for diversity analysis and genetic mapping in pigeonpea.

The sequence of rice chromosomes 11 and 12, rich in disease resistance genes and recent gene duplications

Rice is an important staple food and, with the smallest cereal genome, serves as a reference species for studies on the evolution of cereals and other grasses. Therefore, decoding its entire genome will be a prerequisite for applied and basic research on this species and all other cereals. We have determined and analyzed the complete sequences of two of its chromosomes, 11 and 12, which total 55.9 Mb (14.3% of the entire genome length), based on a set of overlapping clones. A total of 5,993 non-transposable element related genes are present on these chromosomes. Among them are 289 disease resistance-like and 28 defense-response genes, a higher proportion of these categories than on any other rice chromosome. A three-Mb segment on both chromosomes resulted from a duplication 7.7 million years ago (mya), the most recent large-scale duplication in the rice genome. Paralogous gene copies within this segmental duplication can be aligned with genomic assemblies from sorghum and maize. Although these gene copies are preserved on both chromosomes, their expression patterns have diverged. When the gene order of rice chromosomes 11 and 12 was compared to wheat gene loci, significant synteny between these orthologous regions was detected, illustrating the presence of conserved genes alternating with recently evolved genes. Because the resistance and defense response genes, enriched on these chromosomes relative to the whole genome, also occur in clusters, they provide a preferred target for breeding durable disease resistance in rice and the isolation of their allelic variants. The recent duplication of a large chromosomal segment coupled with the high density of disease resistance gene clusters makes this the most recently evolved part of the rice genome. Based on syntenic alignments of these chromosomes, rice chromosome 11 and 12 do not appear to have resulted from a single whole-genome duplication event as previously suggested.BackgroundRice is an important staple food and, with the smallest cereal genome, serves as a reference species for studies on the evolution of cereals and other grasses. Therefore, decoding its entire genome will be a prerequisite for applied and basic research on this species and all other cereals.ResultsWe have determined and analyzed the complete sequences of two of its chromosomes, 11 and 12, which total 55.9 Mb (14.3% of the entire genome length), based on a set of overlapping clones. A total of 5,993 non-transposable element related genes are present on these chromosomes. Among them are 289 disease resistance-like and 28 defense-response genes, a higher proportion of these categories than on any other rice chromosome. A three-Mb segment on both chromosomes resulted from a duplication 7.7 million years ago (mya), the most recent large-scale duplication in the rice genome. Paralogous gene copies within this segmental duplication can be aligned with genomic assemblies from sorghum and maize. Although these gene copies are preserved on both chromosomes, their expression patterns have diverged. When the gene order of rice chromosomes 11 and 12 was compared to wheat gene loci, significant synteny between these orthologous regions was detected, illustrating the presence of conserved genes alternating with recently evolved genes.ConclusionBecause the resistance and defense response genes, enriched on these chromosomes relative to the whole genome, also occur in clusters, they provide a preferred target for breeding durable disease resistance in rice and the isolation of their allelic variants. The recent duplication of a large chromosomal segment coupled with the high density of disease resistance gene clusters makes this the most recently evolved part of the rice genome. Based on syntenic alignments of these chromosomes, rice chromosome 11 and 12 do not appear to have resulted from a single whole-genome duplication event as previously suggested.

Combining QTL mapping and transcriptome profiling of bulked RILs for identification of functional polymorphism for salt tolerance genes in rice (Oryza sativa L.).

Identification of genes for quantitative traits is difficult using any single approach due to complex inheritance of the traits and limited resolving power of the individual techniques. Here a combination of genetic mapping and bulked transcriptome profiling was used to narrow down the number of differentially expressed salt-responsive genes in rice in order to identify functional polymorphism of genes underlying the quantitative trait loci (QTL). A population of recombinant inbred lines (RILs) derived from cross between salt-tolerant variety CSR 27 and salt-sensitive variety MI 48 was used to map QTL for salt ion concentrations in different tissues and salt stress susceptibility index (SSI) for spikelet fertility, grain weight, and grain yield. Eight significant QTL intervals were mapped on chromosomes 1, 8, and 12 for the salt ion concentrations and a QTL controlling SSI for spikelet fertility was co-located in one of these intervals on chromosome 8. However, there were total 2,681 genes in these QTL intervals, making it difficult to pinpoint the genes responsible for the functional differences for the traits. Similarly, transcriptome profiling of the seedlings of tolerant and sensitive parents grown under control and salt-stress conditions showed 798 and 2,407 differentially expressed gene probes, respectively. By analyzing pools of RNA extracted from ten each of extremely tolerant and extremely sensitive RILs to normalize the background noise, the number of differentially expressed genes under salt stress was drastically reduced to 30 only. Two of these genes, an integral transmembrane protein DUF6 and a cation chloride cotransporter, were not only co-located in the QTL intervals but also showed the expected distortion of allele frequencies in the extreme tolerant and sensitive RILs, and therefore are suitable for future validation studies and development of functional markers for salt tolerance in rice to facilitate marker-assisted breeding.

A Moricandia arvensis– based cytoplasmic male sterility and fertility restoration system in Brassica juncea

Abstract A cytoplasmic male-sterility system has been developed in mustard (Brassica juncea) following repeated backcrossings of the somatic hybrid Moricandia arvensis (2n=28, MM)+B. juncea (2n=36, AABB), carrying mitochondria and chloroplasts from M. arvensis, to Brassica juncea. Cytoplasmic male-sterile (CMS) plants are similar to normal B. juncea; however, the leaves exhibit severe chlorosis resulting in delayed flowering. Flowers are normal with slender, non-dehiscent anthers and excellent nectaries. CMS plants show regular meiosis with pollen degeneration occurring during microsporogenesis. Female fertility was normal. Genetic information for fertility restoration was introgressed following the development of a M. arvensis monosomic addition line on CMS B. juncea. The additional chromosome paired allosyndetically with one of the B. juncea bivalents and allowed introgression. The putative restorer plant also exhibited severe chlorosis similar to CMS plants but possessed 89% and 73% pollen and seed fertility, respectively, which subsequently increased to 96% and 87% in the selfed progeny. The progeny of the cross of CMS line with the restorer line MJR-15, segregated into 1 fertile : 1 sterile. The CMS (Moricandia) B. juncea, the restorer (MJR-15), and fertility restored F1 plants possess similar cytoplasmic organellar genomes as revealed by ‘Southern’ analysis.

A Comprehensive Transcriptome Assembly of Pigeonpea (Cajanus cajan L.) using Sanger and Second-Generation Sequencing Platforms

A comprehensive transcriptome assembly for pigeonpea has been developed by analyzing 128.9 million short Illumina GA IIx single end reads, 2.19 million single end FLX/454 reads, and 18 353 Sanger expressed sequenced tags from more than 16 genotypes. The resultant transcriptome assembly, referred to as CcTA v2, comprised 21 434 transcript assembly contigs (TACs) with an N50 of 1510 bp, the largest one being ∼8 kb. Of the 21 434 TACs, 16 622 (77.5%) could be mapped on to the soybean genome build 1.0.9 under fairly stringent alignment parameters. Based on knowledge of intron junctions, 10 009 primer pairs were designed from 5033 TACs for amplifying intron spanning regions (ISRs). By using in silico mapping of BAC-end-derived SSR loci of pigeonpea on the soybean genome as a reference, putative mapping positions at the chromosome level were predicted for 6284 ISR markers, covering all 11 pigeonpea chromosomes. A subset of 128 ISR markers were analyzed on a set of eight genotypes. While 116 markers were validated, 70 markers showed one to three alleles, with an average of 0.16 polymorphism information content (PIC) value. In summary, the CcTA v2 transcript assembly and ISR markers will serve as a useful resource to accelerate genetic research and breeding applications in pigeonpea.

Single-copy genes define a conserved order between rice and wheat for understanding differences caused by duplication, deletion, and transposition of genes

The high-quality rice genome sequence is serving as a reference for comparative genome analysis in crop plants, especially cereals. However, early comparisons with bread wheat showed complex patterns of conserved synteny (gene content) and colinearity (gene order). Here, we show the presence of ancient duplicated segments in the progenitor of wheat, which were first identified in the rice genome. We also show that single-copy (SC) rice genes, those representing unique matches with wheat expressed sequence tag (EST) unigene contigs in the whole rice genome, show more than twice the proportion of genes mapping to syntenic wheat chromosome as compared to the multicopy (MC) or duplicated rice genes. While 58.7% of the 1,244 mapped SC rice genes were located in single syntenic wheat chromosome groups, the remaining 41.3% were distributed randomly to the other six non-syntenic wheat groups. This could only be explained by a background dispersal of genes in the genome through transposition or other unknown mechanism. The breakdown of rice–wheat synteny due to such transpositions was much greater near the wheat centromeres. Furthermore, the SC rice genes revealed a conserved primordial gene order that gives clues to the origin of rice and wheat chromosomes from a common ancestor through polyploidy, aneuploidy, centromeric fusions, and translocations. Apart from the bin-mapped wheat EST contigs, we also compared 56,298 predicted rice genes with 39,813 wheat EST contigs assembled from 409,765 EST sequences and identified 7,241 SC rice gene homologs of wheat. Based on the conserved colinearity of 1,063 mapped SC rice genes across the bins of individual wheat chromosomes, we predicted the wheat bin location of 6,178 unmapped SC rice gene homologs and validated the location of 213 of these in the telomeric bins of 21 wheat chromosomes with 35.4% initial success. This opens up the possibility of directed mapping of a large number of conserved SC rice gene homologs in wheat. Overall, only 46.4% of these SC genes code for proteins with known functional domains; the remaining 53.6% have unknown function, and hence, represent an important, but yet, under explored category of genes.

Molecular mapping of QTLs for plant type and earliness traits in pigeonpea ( Cajanus cajan L. Millsp.)

BackgroundPigeonpea is an important grain legume of the semi-arid tropics and sub-tropical regions where it plays a crucial role in the food and nutritional security of the people. The average productivity of pigeonpea has remained very low and stagnant for over five decades due to lack of genomic information and intensive breeding efforts. Previous SSR-based linkage maps of pigeonpea used inter-specific crosses due to low inter-varietal polymorphism. Here our aim was to construct a high density intra-specific linkage map using genic-SNP markers for mapping of major quantitative trait loci (QTLs) for key agronomic traits, including plant height, number of primary and secondary branches, number of pods, days to flowering and days to maturity in pigeonpea.ResultsA population of 186 F2:3 lines derived from an intra-specific cross between inbred lines ‘Pusa Dwarf’ and ‘HDM04-1’ was used to construct a dense molecular linkage map of 296 genic SNP and SSR markers covering a total adjusted map length of 1520.22 cM for the 11 chromosomes of the pigeonpea genome. This is the first dense intra-specific linkage map of pigeonpea with the highest genome length coverage. Phenotypic data from the F2:3 families were used to identify thirteen QTLs for the six agronomic traits. The proportion of phenotypic variance explained by the individual QTLs ranged from 3.18% to 51.4%. Ten of these QTLs were clustered in just two genomic regions, indicating pleiotropic effects or close genetic linkage. In addition to the main effects, significant epistatic interaction effects were detected between the QTLs for number of pods per plant.ConclusionsA large amount of information on transcript sequences, SSR markers and draft genome sequence is now available for pigeonpea. However, there is need to develop high density linkage maps and identify genes/QTLs for important agronomic traits for practical breeding applications. This is the first report on identification of QTLs for plant type and maturity traits in pigeonpea. The QTLs identified in this study provide a strong foundation for further validation and fine mapping for utilization in the pigeonpea improvement.

Functionally relevant microsatellites in sugarcane unigenes

BackgroundUnigene sequences constitute a rich source of functionally relevant microsatellites. The present study was undertaken to mine the microsatellites in the available unigene sequences of sugarcane for understanding their constitution in the expressed genic component of its complex polyploid/aneuploid genome, assessing their functional significance in silico, determining the extent of allelic diversity at the microsatellite loci and for evaluating their utility in large-scale genotyping applications in sugarcane.ResultsThe average frequency of perfect microsatellite was 1/10.9 kb, while it was 1/44.3 kb for the long and hypervariable class I repeats. GC-rich trinucleotides coding for alanine and the GA-rich dinucleotides were the most abundant microsatellite classes. Out of 15,594 unigenes mined in the study, 767 contained microsatellite repeats and for 672 of these putative functions were determined in silico. The microsatellite repeats were found in the functional domains of proteins encoded by 364 unigenes. Its significance was assessed by establishing the structure-function relationship for the beta-amylase and protein kinase encoding unigenes having repeats in the catalytic domains. A total of 726 allelic variants (7.42 alleles per locus) with different repeat lengths were captured precisely for a set of 47 fluorescent dye labeled primers in 36 sugarcane genotypes and five cereal species using the automated fragment analysis system, which suggested the utility of designed primers for rapid, large-scale and high-throughput genotyping applications in sugarcane. Pair-wise similarity ranging from 0.33 to 0.84 with an average of 0.40 revealed a broad genetic base of the Indian varieties in respect of functionally relevant regions of the large and complex sugarcane genome.ConclusionMicrosatellite repeats were present in 4.92% of sugarcane unigenes, for most (87.6%) of which functions were determined in silico. High level of allelic diversity in repeats including those present in the functional domains of proteins encoded by the unigenes demonstrated their use in assay of useful variation in the genic component of complex polyploid sugarcane genome.

Mapping of QTLs Controlling Na+, K+ and CI− Ion Concentrations in Salt Tolerant Indica Rice Variety CSR27

Soil salinity and sodicity are major constraints to rice production in about twenty per cent of the irrigated crop land. Inbuilt genetic tolerance to salinity is the most economical and environmentally sustainable way to solve this problem. A mapping population of 200 F2 plants and their corresponding F3 families, derived from a cross between a salt tolerant indica rice variety CSR27 and a salt sensitive variety MI48 were used to map OTLs for salt tolerance. Seventeen different parameters, including seedling salt injury score, Na+, K+, CI− concentrations and Na+/K+ ratio in leaf and stem tissues at vegetative and reproductive stages were mapped. A framework linkage map was constructed using 79 SSR and EST markers distributed over the twelve rice chromosomes at an average interval of 20.7cM and total map distance of 1634.5 cM. Twenty five major OTLs, each explaining more than ten per cent of the trait phenotypic variance, were mapped on chromosomes 1, 2, 3 and 8. These included one OTL for seedling salt injury score, nine for Na+ concentration, three for K+ concentration and four for Cl− concentration in leaf and stem tissues at vegetative and reproductive stages. The Na+/K+ ratio, an important ion balancing parameter for the salt tolerance, was controlled by eight OTLs explaining phenotypic variance in the range of 42.88–52.63%. Four OTL intervals were robust with major effect and having OTLs for multiple salt tolerance parameters that might be governed by common or tightly linked genes. One major OTL for multiple salt tolerance parameters on chromosome 8 and three major OTLs for CI− ion concentration are novel for this study. The OTLs identified here will serve as a base for fine mapping, gene tagging and marker assisted selection for salt tolerance in rice.

Molecular mapping of rice blast resistance gene Pi-kh in the rice variety Tetep.

We have used rice line Tetep as a resistant donor with the aim of mapping a durable blast resistance gene Pi-kh using RAPD and AFLP techniques in conjunction with bulk segregant analysis. An F2 mapping population consisting of 205 plants was generated by crossing Tetep with HP2216, a highly susceptible cultivar. Inoculation with specific isolate (PLP-1) of Magnaporthe grisea at seeding stage showed that the Pi-kh gene inherited as a single dominant gene in F2 population. RAPD analysis was performed with 240 primers to detect polymorphism between resistant and susceptible parents. Of these, 48 primers produced polymorphic banding pattern between resistant and susceptible parents. Bulk segregant analysis was performed with 48 primers of which 5 showed polymorphism between resistant and susceptible bulks. A 700 bp DNA band was obtained in resistant F2 plants with primer 5-129 indicating its linkage to the resistance gene. Out of 64 AFLP primer combinations used for polymorphism survey between HP 2216 and Tetep, 11 AFLP primer combinations were able to distinguish the resistant and susceptible bulks. An AFLP band of 75 bp obtained with primer combination, E-TAlM-CTC co-segregated with the resistance gene. The RAPD marker 5-129700 and AFLP75 were placed on the linkage map at a distance of 2.1 eM and 15.1 eM flanking to Pi-khgene, respectively. The RAPD band closely linked to Pi-kh gene was sequenced and used for the development of CAPs markers which also co-segregated with resistant phenotype in the mapping population. On sequence analysis and homology search of RAPD fragment with whole rice genome sequence database and the information available on physical, genetic and sequence maps of rice, the co-segregating CAPs marker was placed at long arm of rice chromosome 11. CAPs marker developed in this study showed polymorphism in different rice cultivars grown in North-Western Himalayan region and is being used for the pyramiding of Pi-kh gene along with other blast resistance genes using marker-assisted selection.