Structural and functional analysis of CCT family genes in pigeonpea.

BACKGROUND: Pigeonpea (Cajanus cajan L.) is a photoperiod-sensitive short-day plant. Understanding the flowering-related genes is critical to developing photoperiod insensitive cultivars. METHODS: The CCT family genes were identified using 'CCT DOMAIN PROTEIN' as a keyword and localized on the chromosomes using the BLAST search option available at the LIS database. The centromeric positions were identified through BLAST search using the centromeric repeat sequence of C. cajan as a query against the chromosome-wise FASTA files downloaded from the NCBI database. The CCT family genes were classified based on additional domains and/or CCT domains. The orthologous and phylogenetic relationships were inferred using the OrthoFinder and MEGA 10.1 software, respectively. The CCT family genes' expression level in photoperiod-sensitive and insensitive genotypes was compared using RNA-seq data and qRT-PCR analysis. RESULTS: We identified 33 CCT family genes in C. cajan distributed on ten chromosomes and nine genomic scaffolds. They were classified into CMF-type, COL-type, PRR-type, and GTCC- type. The CCT family genes of legumes exhibited an extensive orthologous relationship. Glycine max showed the maximum similarity of CCT family genes with C. cajan. The expression analysis of CCT family genes using photoperiod insensitive (ICP20338) and photoperiod sensitive (MAL3) genotypes of C. cajan demonstrated that CcCCT4 and CcCCT23 are the active CONSTANS in ICP20338. In contrast, only CcCCT23 is active in MAL3. CONCLUSION: The CCT family genes in C. cajan vary considerably in structure and domain types. They are maximally similar to soybean's CCT family genes. The differential photoperiod response of pigeonpea genotypes, ICP20338 and MAL3, is possibly due to the difference in the number and types of active CONSTANS in them.

CRISPR-Cas9 directed genome engineering for enhancing salt stress tolerance in rice.

Crop productivity in rice is harshly limited due to high concentration of salt in the soil. To understand the intricacies of the mechanism it is important to unravel the key pathways operating inside the plant cell. Emerging state-of-the art technologies have provided the tools to discover the key components inside the plant cell for salt tolerance. Among the molecular entities, transcription factors and/or other important components of sensing and signaling cascades have been the attractive targets and the role of NHX and SOS1 transporters amply described. Not only marker assisted programs but also transgenic approaches by using reverse genetic strategies (knockout or knockdown) or overexpression have been extensively used to engineer rice crop. CRISPR/Cas is an attractive paradigm and provides the feasibility for manipulating several genes simultaneously. Here, in this review we highlight some of the molecular entities that could be potentially targeted for generating rice amenable to sustain growth under high salinity conditions by employing CRISPR/Cas. We also try to address key questions for rice salt stress tolerance other than what is already known.

Allele mining for a drought responsive gene DRO1 determining root growth angle in donors of drought tolerance in rice (Oryza sativa L.).

Deeper Rooting 1 (DRO1) gene identified from a major QTL on chromosome 9 increases the root growth angle (RGA) and thus facilitates survival under drought and hence is an excellent candidate for rice improvement. Twenty-four major Indian upland and lowland genotypes including the 'yield under drought' (DTY) QTL donors were subjected to allele mining of DRO1 (3058 bp) using four pairs of overlapping primers. A total of 216 and 52 SNPs were identified across all genotypes in the gene and coding region (756 bp) respectively with transversions 3.6 fold more common than transitions in the gene and 2.5 times in the CDS. In 251 amino acid long protein, substitutions were found in 19 positions, wherein change in position 92 was the most frequent. Based on allele mining, the 24 genotypes can be classified into 16 primary structure variants ranging from complete functional allele (Satti, IR36 and DTY 3.1 donor, IR81896-B-B-195) to truncated non-functional alleles in PMK2, IR64, IR20 and Swarna. All the DTY donors, other than IR81896-B-B-195, and most of the upland drought tolerant cultivars (Nagina 22, Vandana and Dhagaddeshi) had accumulated 6-19 SNPs and 4-8 amino acid substitutions resulting in substantial differences in their protein structure. The expression analysis revealed that all the genotypes showed upregulation under drought stress though the degree of upregulation varied among genotypes. The information on structural variations in DRO1 gene will be very useful for the breeders, especially in the light of recent breeding programmes on improving drought tolerance using several DTY donors and upland accessions. SUPPLEMENTARY INFORMATION: The online version of this article (10.1007/s12298-021-00950-2).

Allele-specific analysis of single parent backcross population identifies HOX10 transcription factor as a candidate gene regulating rice root growth.

Understanding the molecular and physiological mechanisms of trait diversity is crucial for crop improvement to achieve drought adaptation. Root traits such as high biomass and/or deep rootedness are undoubtedly important drought adaptive traits. The major aim of this investigation was to functionally characterize a set of ethyl methane sulfonate-induced rice mutants for root traits. We report the identification of a high-root biomass mutant through a novel screening strategy for yield and Δ13 C measurements. The high-root mutant (392-9-1) thus identified, had a 66% higher root biomass compared to wild-type (Nagina-22). Better maintenance of leaf turgor and carbon assimilation rates resulted in lower drought susceptibility index in 392-9-1. Targeted resequencing revealed three non-synonymous single nucleotide variations in 392-9-1 for the genes HOX10, CITRATE SYNTHASE and ZEAXANTHIN EPOXIDASE. Segregation pattern of phenotype and mutant alleles in a single parent backcross F2 population revealed a typical 3:1 segregation for each of the mutant alleles. The number of F2 progeny with root biomass equal to or greater than that of 392-9-1 represented approximately one-third of the population indicating a major role played by HOX10 gene in regulating root growth in rice. Allele-specific Sanger sequencing in contrasting F2 progenies confirmed the co-segregation of HOX10 allele with the root biomass. The non-synonymous mutations in the other two genes did not reveal any specific pattern of co-segregation with root phenotype, indicating a strong role of HOX10, an upstream transcription factor, in regulating root biomass in rice.

Genome-wide discovery of tissue-specific miRNAs in clusterbean (Cyamopsis tetragonoloba) indicates their association with galactomannan biosynthesis.

Owing to the presence of 80% soluble dietary fibre, high protein content and high value gum, clusterbean (Cyamopsis tetragonoloba) has recently emerged as an economically important legume. The developing clusterbean seeds accumulate 90% galactomannans in the endosperm and, therefore, can be used as a model crop to understand galactomannan biosynthesis and its regulation. miRNAs are tiny master regulators of their corresponding target genes, resulting in variations in the amounts of their metabolic end products. To understand the role of these regulators in galactomannan biosynthesis regulation, small RNA libraries were prepared and sequenced from five tissues of clusterbean genotype RGC-936, and miRanalyzer and DSAP programs were used to identify conserved miRNAs and novel small RNAs. A total of 187 known and 171 novel miRNAs were found to be differentially expressed, of which 10 miRNAs were validated. A complicated network topology and 35% sharing of the target mRNAs between known and novel miRNAs suggest random evolution of novel miRNAs. The gene ontology (GO) annotation of potential target genes revealed the genes coding for signalling and carbohydrate metabolism (50.10%), kinases and other enzymes (20.75%), transcription factors (10.20%), transporters (8.35%) and other targets (10.6%). Two novel unigenes were annotated as ManS (mannosyltransferase/mannan synthase) and UGE (UDP- D-glucose 4-epimerase) and validated as targets for three novel miRNAs, that is Ct-miR3130, Ct-miR3135 and Ct-miR3157. Our findings reveal that these novel miRNAs could play an important role in the regulation of the galactomannan pathway in C. tetragonoloba and possibly other galactomannan-producing species.

Genetics, Physiological Mechanisms and Breeding of Flood-Tolerant Rice (Oryza sativa L.).

Flooding of rice fields is a serious problem in the river basins of South and South-East Asia where about 15 Mha of lowland rice cultivation is regularly affected. Flooding creates hypoxic conditions resulting in poor germination and seedling establishment. Flash flooding, where rice plants are completely submerged for 10-15 d during their vegetative stage, causes huge losses. Water stagnation for weeks to months also leads to substantial yield losses when large parts of rice aerial tissues are inundated. The low-yielding traditional varieties and landraces of rice adapted to these flooding conditions have been replaced by flood-sensitive high-yielding rice varieties. The 'FR13A' rice variety and the Submergence 1A (SUB1A) gene were identified for flash flooding and subsequently introgressed to high-yielding rice varieties. The challenge is to find superior alleles of the SUB1A gene, or even new genes that may confer greater tolerance to submergence. Similarly, genes have been identified in tolerant landraces of rice for their ability to survive by rapid stem elongation (SNORKEL1 and SNORKEL2) during deep-water flooding, and for anaerobic germination ability (TPP7). Research on rice genotypes and novel genes that are tolerant to prolonged water stagnation is in progress. These studies will greatly assist in devising more efficient and precise molecular breeding strategies for developing climate-resilient high-yielding rice varieties for flood-prone regions. Here we review the state of our knowledge of flooding tolerance in rice and its application in varietal improvement.

De novo transcriptome sequencing facilitates genomic resource generation in Tinospora cordifolia.

Tinospora cordifolia is known for its medicinal properties owing to the presence of useful constituents such as terpenes, glycosides, steroids, alkaloids, and flavonoids belonging to secondary metabolism origin. However, there is little information available pertaining to critical genomic elements (ESTs, molecular markers) necessary for judicious exploitation of its germplasm. We employed 454 GS-FLX pyrosequencing of entire transcripts and altogether ∼25 K assembled transcripts or Expressed sequence tags (ESTs) were identified. As the interest in T. cordifolia is primarily due to its secondary metabolite constituents, the ESTs pertaining to terpenoids biosynthetic pathway were identified in the present study. Additionally, several ESTs were assigned to different transcription factor families. To validate our transcripts dataset, the novel EST-SSR markers were generated to assess the genetic diversity among germplasm of T. cordifolia. These EST-SSR markers were found to be polymorphic and the dendrogram based on dice similarity index revealed three distinct clustering of accessions. The present study demonstrates effectiveness in using both NEWBLER and MIRA sequence read assembler software for enriching transcript-dataset and thus enables better exploitation of EST resources for mining candidate genes and designing molecular markers.

Single nucleotide polymorphism in sugar pathway and disease resistance genes in sugarcane.

KEY MESSAGE: Single nucleotide polymorphism in sugar pathway and disease resistance genes showing genetic association with sugar content and red rot resistance would be useful in marker-assisted genetic improvement of sugarcane. Validation and genotyping of potential sequence variants in candidate genes are necessary to understand their functional significance and trait association potential. We discovered, characterized, validated and genotyped SNPs and InDels in sugar pathway and disease resistance genes of Saccharum complex and sugarcane varieties using amplicon sequencing and CAPS assays. The SNPs were abundant in the non-coding 3'UTRs than 5'UTRs and coding sequences depicting a strong bias toward C to T transition substitutions than transversions. Sequencing of cloned amplicons validated 61.6 and 45.2 % SNPs detected in silico in 21 sugar pathway and 16 disease resistance genes, respectively. Sixteen SNPs in four sugar pathway genes and 10 SNPs in nine disease resistance genes were validated through cost-effective CAPS assay. Functional and adaptive significance of SNP and protein haplotypes identified in sugar pathway and disease resistance genes was assessed by correlating their allelic variation with missense amino acid substitutions in the functional domains, alteration in protein structure models and possible modulation of catalytic enzyme activity in contrasting high and low sugar and moderately red rot resistant and highly susceptible sugarcane genotypes. A strong genetic association of five SNPs in the sugar pathway and disease resistance genes, and an InDel marker in the promoter sequence of sucrose synthase-2 gene, with sugar content and red rot resistance, was evident. The functionally relevant SNPs and InDels, detected and validated in sugar pathway and disease resistance genes, and genic CAPS markers designed, would be of immense use in marker-assisted genetic improvement of sugarcane for sugar content and disease resistance.

Comparative transcriptome analysis of two contrasting genotypes provides new insights into the drought response mechanism in pigeon pea (Cajanus cajan L. Millsp.).

BACKGROUND: Despite plant's ability to adapt and withstand challenging environments, drought poses a severe threat to their growth and development. Although pigeon pea is already quite resistant to drought, the prolonged dehydration induced by the aberrant climate poses a serious threat to their survival and productivity. OBJECTIVE: Comparative physiological and transcriptome analyses of drought-tolerant (CO5) and drought-sensitive (CO1) pigeon pea genotypes subjected to drought stress were carried out in order to understand the molecular basis of drought tolerance in pigeon pea. METHODS: The transcriptomic analysis allowed us to examine how drought affects the gene expression of C. cajan. Using bioinformatics tools, the unigenes were de novo assembled, annotated, and functionally evaluated. Additionally, a homology-based sequence search against the droughtDB database was performed to identify the orthologs of the DEGs. RESULTS: 1102 potential drought-responsive genes were found to be differentially expressed genes (DEGs) between drought-tolerant and drought-sensitive genotypes. These included Abscisic acid insensitive 5 (ABI5), Nuclear transcription factor Y subunit A-7 (NF-YA7), WD40 repeat-containing protein 55 (WDR55), Anthocyanidin reductase (ANR) and Zinc-finger homeodomain protein 6 (ZF-HD6) and were highly expressed in the tolerant genotype. Further, GO analysis revealed that the most enriched classes belonged to biosynthetic and metabolic processes in the biological process category, binding and catalytic activity in the molecular function category and nucleus and protein-containing complex in the cellular component category. Results of KEGG pathway analysis revealed that the DEGs were significantly abundant in signalling pathways such as plant hormone signal transduction and MAPK signalling pathways. Consequently, in our investigation, we have identified and validated by qPCR a group of genes involved in signal reception and propagation, stress-specific TFs, and basal regulatory genes associated with drought response. CONCLUSION: In conclusion, our comprehensive transcriptome dataset enabled the discovery of candidate genes connected to pathways involved in pigeon pea drought response. Our research uncovered a number of unidentified genes and transcription factors that could be used to understand and improve susceptibility to drought.

Phylogenomic Analysis of micro-RNA Involved in Juvenile to Flowering-Stage Transition in Photophilic Rice and Its Sister Species.

Vegetative to reproductive phase transition in phototropic plants is an important developmental process and is sequentially mediated by the expression of micro-RNA MIR172. To obtain insight into the evolution, adaptation, and function of MIR172 in photophilic rice and its wild relatives, we analyzed the genescape of a 100 kb segment harboring MIR172 homologs from 11 genomes. The expression analysis of MIR172 revealed its incremental accumulation from the 2-leaf to 10-leaf stage, with maximum expression coinciding with the flag-leaf stage in rice. Nonetheless, the microsynteny analysis of MIR172s revealed collinearity within the genus Oryza, but a loss of synteny was observed in (i) MIR172A in O. barthii (AA) and O. glaberima (AA); (ii) MIR172B in O. brachyantha (FF); and (iii) MIR172C in O. punctata (BB). Phylogenetic analysis of precursor sequences/region of MIR172 revealed a distinct tri-modal clade of evolution. The genomic information generated in this investigation through comparative analysis of MIRNA, suggests mature MIR172s to have evolved in a disruptive and conservative mode amongst all Oryza species with a common origin of descent. Further, the phylogenomic delineation provided an insight into the adaptation and molecular evolution of MIR172 to changing environmental conditions (biotic and abiotic) of phototropic rice through natural selection and the opportunity to harness untapped genomic regions from rice wild relatives (RWR).

Drought and Oxidative Stress in Flax (Linum usitatissimum L.) Entails Harnessing Non-Canonical Reference Gene for Precise Quantification of qRT-PCR Gene Expression.

Flax (Linum usitatissimum L.) is a self-pollinating, annual, diploid crop grown for multi-utility purposes for its quality oil, shining bast fiber, and industrial solvent. Being a cool (Rabi) season crop, it is affected by unprecedented climatic changes such as high temperature, drought, and associated oxidative stress that, globally, impede its growth, production, and productivity. To precisely assess the imperative changes that are inflicted by drought and associated oxidative stress, gene expression profiling of predominant drought-responsive genes (AREB, DREB/CBF, and ARR) was carried out by qRT-PCR. Nevertheless, for normalization/quantification of data obtained from qRT-PCR results, a stable reference gene is mandatory. Here, we evaluated a panel of four reference genes (Actin, EF1a, ETIF5A, and UBQ) and assessed their suitability as stable reference genes for the normalization of gene expression data obtained during drought-induced oxidative stress in flax. Taking together, from the canonical expression of the proposed reference genes in three different genotypes, we report that EF1a as a stand-alone and EF1a and ETIF5A in tandem are suitable reference genes to be used for the real-time visualization of cellular impact of drought and oxidative stress on flax.

SNPs in stress-responsive rice genes: validation, genotyping, functional relevance and population structure.

BACKGROUND: Single nucleotide polymorphism (SNP) validation and large-scale genotyping are required to maximize the use of DNA sequence variation and determine the functional relevance of candidate genes for complex stress tolerance traits through genetic association in rice. We used the bead array platform-based Illumina GoldenGate assay to validate and genotype SNPs in a select set of stress-responsive genes to understand their functional relevance and study the population structure in rice. RESULTS: Of the 384 putative SNPs assayed, we successfully validated and genotyped 362 (94.3%). Of these 325 (84.6%) showed polymorphism among the 91 rice genotypes examined. Physical distribution, degree of allele sharing, admixtures and introgression, and amino acid replacement of SNPs in 263 abiotic and 62 biotic stress-responsive genes provided clues for identification and targeted mapping of trait-associated genomic regions. We assessed the functional and adaptive significance of validated SNPs in a set of contrasting drought tolerant upland and sensitive lowland rice genotypes by correlating their allelic variation with amino acid sequence alterations in catalytic domains and three-dimensional secondary protein structure encoded by stress-responsive genes. We found a strong genetic association among SNPs in the nine stress-responsive genes with upland and lowland ecological adaptation. Higher nucleotide diversity was observed in indica accessions compared with other rice sub-populations based on different population genetic parameters. The inferred ancestry of 16% among rice genotypes was derived from admixed populations with the maximum between upland aus and wild Oryza species. CONCLUSIONS: SNPs validated in biotic and abiotic stress-responsive rice genes can be used in association analyses to identify candidate genes and develop functional markers for stress tolerance in rice.

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

BACKGROUND: Pigeonpea 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. RESULTS: A 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. CONCLUSIONS: A 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.

The first draft of the pigeonpea genome sequence.

Pigeonpea (Cajanus cajan) is an important grain legume of the Indian subcontinent, South-East Asia and East Africa. More than eighty five percent of the world pigeonpea is produced and consumed in India where it is a key crop for food and nutritional security of the people. Here we present the first draft of the genome sequence of a popular pigeonpea variety 'Asha'. The genome was assembled using long sequence reads of 454 GS-FLX sequencing chemistry with mean read lengths of >550 bp and >10-fold genome coverage, resulting in 510,809,477 bp of high quality sequence. Total 47,004 protein coding genes and 12,511 transposable elements related genes were predicted. We identified 1,213 disease resistance/defense response genes and 152 abiotic stress tolerance genes in the pigeonpea genome that make it a hardy crop. In comparison to soybean, pigeonpea has relatively fewer number of genes for lipid biosynthesis and larger number of genes for cellulose synthesis. The sequence contigs were arranged in to 59,681 scaffolds, which were anchored to eleven chromosomes of pigeonpea with 347 genic-SNP markers of an intra-species reference genetic map. Eleven pigeonpea chromosomes showed low but significant synteny with the twenty chromosomes of soybean. The genome sequence was used to identify large number of hypervariable 'Arhar' simple sequence repeat (HASSR) markers, 437 of which were experimentally validated for PCR amplification and high rate of polymorphism among pigeonpea varieties. These markers will be useful for fingerprinting and diversity analysis of pigeonpea germplasm and molecular breeding applications. This is the first plant genome sequence completed entirely through a network of Indian institutions led by the Indian Council of Agricultural Research and provides a valuable resource for the pigeonpea variety improvement.

TinoTranscriptDB: A Database of Transcripts and Microsatellite Markers of Tinospora cordifolia, an Important Medicinal Plant.

Tinospora cordifolia, commonly known as "Giloe" in India, is a shrub belonging to the family Menispermaceae. It is an important medicinal plant known for its antipyretic, anti-inflammatory, antispasmodic, and antidiabetic properties and is used in the treatment of jaundice, gout, and rheumatism. Despite its economic importance, the limited information related to its genomic resources prohibits its judicious exploitation through molecular breeding or biotechnological approaches. In this study, we generated a meta-transcriptome assembly of 43,090 non-redundant transcripts by merging the RNASeq data obtained from Roche 454 GS-FLX, and Illumina platforms, and report the first transcriptome-based database for simple sequence repeats and transcription factors ("TinoTranscriptDB" (Tinospora cordifolia Transcriptome Database)). We annotated 26,716 (62%) of the total transcripts successfully from National Center for Biotechnology Information non-redundant protein (NCBI-NR), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-Prot, and Pfam databases. This database contains information of 2620 perfect simple sequence repeats (P-SSRs) with a relative abundance of 340.12 (loci/Mb), and relative density of 6309.29 (bp/Mb). Excluding mono-nucleotides, the most abundant SSR motifs were tri-nucleotides (54.31%), followed by di-nucleotides (37.51%), tetra-nucleotides (4.54%), penta-nucleotides (3.16%) and hexa-nucleotides (0.45%). Additionally, we also identified 4,311 transcription factors (TFs) and categorized them into 55 sub-families. This database is expected to fill the gap in genomic resource availability in T. cordifolia and thus accelerate molecular breeding and related functional and other applied studies aimed towards genetic improvements of T. cordifolia and related species.

Identification of novel resources for panicle blast resistance from wild rice accessions and mutants of cv. Nagina 22 by syringe inoculation under field conditions.

Panicle blast is the most severe type of rice blast disease. Screening of rice genotypes for panicle blast resistance at the field level requires an efficient and robust method of inoculation. Here, we standardized a method that can be utilized for both small- and large-scale screening and assessment of panicle blast infection and disease reaction. The method involves inoculation of Magnaporthe oryzae spore culture in the neck of the rice panicle using a syringe and covering the inoculation site with wet cotton wrapped with aluminum foil to provide the required humidity for spore germination. The method was standardized using panicle blast-resistant cv. Tetep and susceptible cv. HP2216 inoculated with Mo-ni-025 isolate of M. oryzae. The method was evaluated at phenotypic as well as molecular level by expression analysis of disease responsive pathogenesis-related (PR) genes. We found this method simple, robust, reliable, and highly efficient for screening of large germplasm sets of rice for panicle blast. This was validated by screening the wild rice germplasm for panicle blast response in the field using three M. oryzae strains and subsequently with the most virulent strain in 45 EMS-induced mutants of Nagina 22 shortlisted based on field screening in a blast hotspot region. We identified five novel blast disease-resistant wild rice genotypes and 15 Nagina 22 mutants that can be used in breeding programmes.

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

BACKGROUND: Pigeonpea [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. RESULTS: In 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. CONCLUSION: We 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.

Analysis of BAC-end sequences (BESs) and development of BES-SSR markers for genetic mapping and hybrid purity assessment in pigeonpea (Cajanus spp.).

BACKGROUND: Pigeonpea [Cajanus cajan (L.) Millsp.] is an important legume crop of rainfed agriculture. Despite of concerted research efforts directed to pigeonpea improvement, stagnated productivity of pigeonpea during last several decades may be accounted to prevalence of various biotic and abiotic constraints and the situation is exacerbated by availability of inadequate genomic resources to undertake any molecular breeding programme for accelerated crop improvement. With the objective of enhancing genomic resources for pigeonpea, this study reports for the first time, large scale development of SSR markers from BAC-end sequences and their subsequent use for genetic mapping and hybridity testing in pigeonpea. RESULTS: A set of 88,860 BAC (bacterial artificial chromosome)-end sequences (BESs) were generated after constructing two BAC libraries by using HindIII (34,560 clones) and BamHI (34,560 clones) restriction enzymes. Clustering based on sequence identity of BESs yielded a set of >52K non-redundant sequences, comprising 35 Mbp or >4% of the pigeonpea genome. These sequences were analyzed to develop annotation lists and subdivide the BESs into genome fractions (e.g., genes, retroelements, transpons and non-annotated sequences). Parallel analysis of BESs for microsatellites or simple sequence repeats (SSRs) identified 18,149 SSRs, from which a set of 6,212 SSRs were selected for further analysis. A total of 3,072 novel SSR primer pairs were synthesized and tested for length polymorphism on a set of 22 parental genotypes of 13 mapping populations segregating for traits of interest. In total, we identified 842 polymorphic SSR markers that will have utility in pigeonpea improvement. Based on these markers, the first SSR-based genetic map comprising of 239 loci was developed for this previously uncharacterized genome. Utility of developed SSR markers was also demonstrated by identifying a set of 42 markers each for two hybrids (ICPH 2671 and ICPH 2438) for genetic purity assessment in commercial hybrid breeding programme. CONCLUSION: In summary, while BAC libraries and BESs should be useful for genomics studies, BES-SSR markers, and the genetic map should be very useful for linking the genetic map with a future physical map as well as for molecular breeding in pigeonpea.

Analysis of genetic diversity in earthworms using DNA markers.

Earthworms are one of the most important and beneficial macrofauna, and are used extensively in organic farming. Earthworms mediate soil biological regulation systems, and produce biogenic structures. They help to maintain soil structure, water infiltration, and regulate the availability of nutrients assimilated by plants. The objectives of this study were to perform morphological and molecular characterizations of 24 earthworm individuals collected from geographically diverse locations to assess the level of genetic variation. For molecular analysis, the effectiveness of RAPD, ISSR, and Universal rice primers (URPs) markers was investigated to identify polymorphism among 24 isolates of earthworms. A total of 62 molecular markers were used for amplification of genomic DNA of earthworms. Of these, 10 RAPD, 10 ISSR, and 10 URPs markers were used for characterization, which showed 95.7%, 96.7% and 98.3% polymorphism, respectively. The dendrogram, generated from the DNA markers by the unweighted pair group method using arithmetic averages, grouped all the isolates into two main clusters. All Eisenia fetida isolates were clustered in group A, whereas group B included three isolates belonging to Eudrilus eugeniae. Molecular markers allowed a rapid assessment of genetic variation among these closely related isolates of earthworms. These results suggest that molecular markers are a good choice for diversity analysis of earthworm individuals.

Functionally relevant microsatellites in sugarcane unigenes.

BACKGROUND: Unigene 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. RESULTS: The 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. CONCLUSION: Microsatellite 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.