Ajay Kumar Parida (1963-2022), an eminent plant biotechnologist with a passion for mangrove biology.

We remember Dr Ajay Parida, a leading plant biotechnologist, whose premature passing has deprived the Indian plant science community of a committed scientist and an able administrator. Born on 12 December 1963 in Bhagabanpur, Cuttack District (now Jajpur district), Odisha, he passed away in Guwahati on 19 July 2022. A collegial scientist, his down-to-earth and approachable nature, as well as his resourcefulness were instrumental in advancing the cause of Indian science and harnessing frontier biotechnological tools as vehicles of social consciousness. His expertise in quantitative DNA variation and molecular marker analysis, paved the way for subsequent research on mangrove molecular diversity at the M. S. Swaminathan Research Foundation (MSSRF), Chennai. His contributions to mangrove biology, genetics and genomics as well as extremophile plant species in the Indian context over two decades are a benchmark in his field. He also provided commendable leadership in his capacity as Director, Institute of Life Sciences (ILS), Bhubaneshwar during the COVID-19 pandemic.

De novo transcriptome assembly and analysis of gene expression in different tissues of moth bean (Vigna aconitifolia) (Jacq.) Marechal.

BACKGROUND: The underutilized species Vigna aconitifolia (Moth Bean) is an important legume crop cultivated in semi-arid conditions and is valued for its seeds for their high protein content. It is also a popular green manure cover crop that offers many agronomic benefits including nitrogen fixation and soil nutrients. Despite its economic potential, genomic resources for this crop are scarce and there is limited knowledge on the developmental process of this plant at a molecular level. In the present communication, we have studied the molecular mechanisms that regulate plant development in V. aconitifolia, with a special focus on flower and seed development. We believe that this study will greatly enrich the genomic resources for this plant in form of differentially expressed genes, transcription factors, and genic molecular markers. RESULTS: We have performed the de novo transcriptome assembly using six types of tissues from various developmental stages of Vigna aconitifolia (var. RMO-435), namely, leaves, roots, flowers, pods, and seed tissue in the early and late stages of development, using the Illumina NextSeq platform. We assembled the transcriptome to get 150938 unigenes with an average length of 937.78 bp. About 79.9% of these unigenes were annotated in public databases and 12839 of those unigenes showed a significant match in the KEGG database. Most of the unigenes displayed significant differential expression in the late stages of seed development as compared with leaves. We annotated 74082 unigenes as transcription factors and identified 12096 simple sequence repeats (SSRs) in the genic regions of V.aconitifolia. Digital expression analysis revealed specific gene activities in different tissues which were validated using Real-time PCR analysis. CONCLUSIONS: The Vigna aconitifolia transcriptomic resources generated in this study provide foundational resources for gene discovery with respect to various developmental stages. This study provides the first comprehensive analysis revealing the genes involved in molecular as well as metabolic pathways that regulate seed development and may be responsible for the unique nutritive values of moth bean seeds. Hence, this study would serve as a foundation for characterization of candidate genes which would not only provide novel insights into understanding seed development but also provide resources for improved moth bean and related species genetic enhancement.

Identification and molecular characterization of miRNAs and their target genes associated with seed development through small RNA sequencing in chickpea.

Multiple studies have attempted to dissect the molecular mechanism underlying seed development in chickpea (Cicer arietinum L.). These studies highlight the need to focus on the role of miRNAs in regulating storage protein accumulation in seeds. Therefore, a total of 8,856,691 short-read sequences were generated from a small RNA library of developing chickpea seeds and were analyzed using miRDeep-P to identify 74 known and 26 novel miRNA sequences. Known miRNAs were classified into 22 miRNA families with miRNA156 family being most abundant. Of the 26 putative novel miRNAs identified, only 22 could be experimentally validated using stem loop end point PCR. Differential expression analyses led to the identification of known as well as novel miRNAs that could regulate various stages of chickpea seed development. In silico target prediction revealed several important target genes and transcription factors like SPL, mediator of RNA Polymerase II transcription subunit 12, aspartic proteinase and NACs, which were further validated by real-time PCR analysis. A comparative expression analysis in chickpea genotypes with contrasting seed protein content revealed one known (Car-miR156h) and two novel miRNA (Car-novmiR7 and Car-novmiR23) candidates to be highly expressed in the LPC (low protein content) chickpea genotypes, targets of which are known to regulate seed storage protein accumulation. Therefore, this study provides a useful resource in the form of miRNA and their targets which can be further utilized to understand and manipulate various regulatory mechanisms involved in seed development with the overall aim of improving yield and nutrition attributes in chickpea.

A high-density SNP-based linkage map using genotyping-by-sequencing and its utilization for improved genome assembly of chickpea (Cicer arietinum L.).

Genotyping-by-sequencing (GBS) allows rapid identification of markers for use in development of linkage maps, which expedite efficient breeding programs. In the present study, we have utilized GBS approach to identify and genotype single-nucleotide polymorphism (SNP) markers in an inter-specific RIL population of Cicer arietinum L. X C. reticulatum. A total of 141,639 raw SNPs were identified using the TASSEL-GBS pipeline. After stringent filtering, 8208 candidate SNPs were identified of which ~ 37% were localized in the intragenic regions followed by genic regions (~ 30%) and intergenic regions (~ 27%). We then utilized 6920 stringent selected SNPs from present study and 6714 SNPs and microsatellite markers available from previous studies for construction of linkage map. The resulting high-density linkage map comprising of eight linkage groups contained 13,590 markers which spanned 1299.14 cM of map length with an average marker density of 0.095 cM. Further, the derived linkage map was used to improve the available assembly of desi chickpea genome by anchoring 443 previously unplaced scaffolds onto eight linkage groups. The present efforts have refined anchoring of the desi chickpea genome assembly to 55.57% of the ~ 520 Mb of assembled desi genome. To the best of our knowledge, the linkage map generated in the present study represents one of the most dense linkage map developed for the crop till date. It will serve as a valuable resource for fine mapping and positional cloning of important quantitative trait loci (QTLs) associated with agronomical traits and also for anchoring and ordering of future genome sequence assemblies.

Transcriptomics reveal useful resources for examining fruit development and variation in fruit size in Coccinia grandis.

Introduction: The Cucurbitaceae family comprises many agronomically important members, that bear nutritious fruits and vegetables of great economic importance. Coccinia grandis, commonly known as Ivy gourd, belongs to this family and is widely consumed as a vegetable. Members of this family are known to display an impressive range of variation in fruit morphology. Although there have been studies on flower development in Ivy gourd, fruit development remains unexplored in this crop. Methods: In this study, comparative transcriptomics of two Ivy gourd cultivars namely "Arka Neelachal Kunkhi" (larger fruit size) and "Arka Neelachal Sabuja" (smaller fruit size) differing in their average fruit size was performed. A de novo transcriptome assembly for Ivy gourd was developed by collecting fruits at different stages of development (5, 10, 15, and 20 days after anthesis i.e. DAA) from these two varieties. The transcriptome was analyzed to identify differentially expressed genes, transcription factors, and molecular markers. Results: The transcriptome of Ivy gourd consisted of 155205 unigenes having an average contig size of 1472bp. Unigenes were annotated on publicly available databases to categorize them into different biological functions. Out of these, 7635 unigenes were classified into 38 transcription factor (TF) families, of which Trihelix TFs were most abundant. A total of 11,165 unigenes were found to be differentially expressed in both the varieties and the in silico expression results were validated through real-time PCR. Also, 98768 simple sequence repeats (SSRs) were identified in the transcriptome of Ivy gourd. Discussion: This study has identified a number of genes, including transcription factors, that could play a crucial role in the determination of fruit shape and size in Ivy gourd. The presence of polymorphic SSRs indicated a possibility for marker-assisted selection for crop breeding in Ivy gourd. The information obtained can help select candidate genes that may be implicated in regulating fruit development and size in other fruit crops.

De novo Whole-Genome Assembly of Moringa oleifera Helps Identify Genes Regulating Drought Stress Tolerance.

Abiotic stresses, especially drought stress, are responsible for heavy losses in productivity, which in turn poses an imminent threat for future food security. Understanding plants' response to abiotic stress at the molecular level is crucially important for mitigating the impacts of climate change. Moringa oleifera is an important multipurpose plant with medicinal and nutritional properties and with an ability to grow in low water conditions, which makes the species an ideal candidate to study the regulatory mechanisms that modulate drought tolerance and its possible use in agroforestry system. In the present communication, we report whole-genome sequencing (WGS) of this species and assemble about 90% of the genome of M. oleifera var. Bhagya into 915 contigs with a N50 value of 4.7 Mb and predicted 32,062 putative protein-coding genes. After annotating the genome, we have chosen to study the heat shock transcription factor (HSF) family of genes to analyze their role in drought tolerance in M. oleifera. We predicted a total of 21 HSFs in the M. oleifera genome and carried out phylogenetic analyses, motif identification, analysis of gene duplication events, and differential expression of the HSF-coding genes in M. oleifera. Our analysis reveals that members of the HSF family have an important role in the plant's response to abiotic stress and are viable candidates for further characterization.

De-novo transcriptome analysis unveils differentially expressed genes regulating drought and salt stress response in Panicum sumatrense.

Screening the transcriptome of drought tolerant variety of little millet (Panicum sumatrense), a marginally cultivated, nutritionally rich, susbsistent crop, can identify genes responsible for its hardiness and enable identification of new sources of genetic variation which can be used for crop improvement. RNA-Seq generated ~ 230 million reads from control and treated tissues, which were assembled into 86,614 unigenes. In silico differential gene expression analysis created an overview of patterns of gene expression during exposure to drought and salt stress. Separate gene expression profiles for leaf and root tissue revealed the differences in regulatory mechanisms operating in these tissues during exposure to abiotic stress. Several transcription factors were identified and studied for differential expression. 61 differentially expressed genes were found to be common to both tissues under drought and salinity stress and were further validated using qRT-PCR. Transcriptome of P. sumatrense was also used to mine for genic SSR markers relevant to abiotic stress tolerance. This study is first report on a detailed analysis of molecular mechanisms of drought and salinity stress tolerance in a little millet variety. Resources generated in this study can be used as potential candidates for further characterization and to improve abiotic stress tolerance in food crops.

De novo transcriptome assembly and analysis of Phragmites karka, an invasive halophyte, to study the mechanism of salinity stress tolerance.

With the rapidly deteriorating environmental conditions, the development of stress tolerant plants has become a priority for sustaining agricultural productivity. Therefore, studying the process of stress tolerance in naturally tolerant species hold significant promise. Phragmites karka is an invasive plant species found abundantly in tropical and sub tropical regions, fresh water regions and brackish marshy areas, such as river banks and lake shores. The plant possesses the ability to adapt and survive under conditions of high salinity. We subjected P. karka seedlings to salt stress and carried out whole transcriptome profiling of leaf and root tissues. Assessing the global transcriptome changes under salt stress resulted in the identification of several genes that are differentially regulated under stress conditions in root and leaf tissue. A total of 161,403 unigenes were assembled and used as a reference for digital gene expression analysis. A number of key metabolic pathways were found to be over-represented. Digital gene expression analysis was validated using qRT-PCR. In addition, a number of different transcription factor families including WRKY, MYB, CCCH, NAC etc. were differentially expressed under salinity stress. Our data will facilitate further characterisation of genes involved in salinity stress tolerance in P. karka. The DEGs from our results are potential candidates for understanding and engineering abiotic stress tolerance in plants.

Genome-wide analysis of the CCCH zinc finger family identifies tissue specific and stress responsive candidates in chickpea (Cicer arietinum L.).

The CCCH zinc finger is a group of proteins characterised by a typical motif consisting of three cysteine residues and one histidine residue. These proteins have been reported to play important roles in regulation of plant growth, developmental processes and environmental responses. In the present study, genome wide analysis of the CCCH zinc finger gene family was carried out in the available chickpea genome. Various bioinformatics tools were employed to predict 58 CCCH zinc finger genes in chickpea (designated CarC3H1-58), which were analysed for their physio-chemical properties. Phylogenetic analysis classified the proteins into 12 groups in which members of a particular group had similar structural organization. Further, the numbers as well as the types of CCCH motifs present in the CarC3H proteins were compared with those from Arabidopsis and Medicago truncatula. Synteny analysis revealed valuable information regarding the evolution of this gene family. Tandem and segmental duplication events were identified and their Ka/Ks values revealed that the CarC3H gene family in chickpea had undergone purifying selection. Digital, as well as real time qRT-PCR expression analysis was performed which helped in identification of several CarC3H members that expressed preferentially in specific chickpea tissues as well as during abiotic stresses (desiccation, cold, salinity). Moreover, molecular characterization of an important member CarC3H45 was carried out. This study provides comprehensive genomic information about the important CCCH zinc finger gene family in chickpea. The identified tissue specific and abiotic stress specific CCCH genes could be potential candidates for further characterization to delineate their functional roles in development and stress.

Dissecting the Root Nodule Transcriptome of Chickpea (Cicer arietinum L.).

A hallmark trait of chickpea (Cicer arietinum L.), like other legumes, is the capability to convert atmospheric nitrogen (N2) into ammonia (NH3) in symbiotic association with Mesorhizobium ciceri. However, the complexity of molecular networks associated with the dynamics of nodule development in chickpea need to be analyzed in depth. Hence, in order to gain insights into the chickpea nodule development, the transcriptomes of nodules at early, middle and late stages of development were sequenced using the Roche 454 platform. This generated 490.84 Mb sequence data comprising 1,360,251 reads which were assembled into 83,405 unigenes. Transcripts were annotated using Gene Ontology (GO), Cluster of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways analysis. Differential expression analysis revealed that a total of 3760 transcripts were differentially expressed in at least one of three stages, whereas 935, 117 and 2707 transcripts were found to be differentially expressed in the early, middle and late stages of nodule development respectively. MapMan analysis revealed enrichment of metabolic pathways such as transport, protein synthesis, signaling and carbohydrate metabolism during root nodulation. Transcription factors were predicted and analyzed for their differential expression during nodule development. Putative nodule specific transcripts were identified and enriched for GO categories using BiNGO which revealed many categories to be enriched during nodule development, including transcription regulators and transporters. Further, the assembled transcriptome was also used to mine for genic SSR markers. In conclusion, this study will help in enriching the transcriptomic resources implicated in understanding of root nodulation events in chickpea.

High density linkage mapping of genomic and transcriptomic SNPs for synteny analysis and anchoring the genome sequence of chickpea.

This study presents genome-wide discovery of SNPs through next generation sequencing of the genome of Cicer reticulatum. Mapping of the C. reticulatum sequenced reads onto the draft genome assembly of C. arietinum (desi chickpea) resulted in identification of 842,104 genomic SNPs which were utilized along with an additional 36,446 genic SNPs identified from transcriptome sequences of the aforementioned varieties. Two new chickpea Oligo Pool All (OPAs) each having 3,072 SNPs were designed and utilized for SNP genotyping of 129 Recombinant Inbred Lines (RILs). Using Illumina GoldenGate Technology genotyping data of 5,041 SNPs were generated and combined with the 1,673 marker data from previously published studies, to generate a high resolution linkage map. The map comprised of 6698 markers distributed on eight linkage groups spanning 1083.93 cM with an average inter-marker distance of 0.16 cM. Utility of the present map was demonstrated for improving the anchoring of the earlier reported draft genome sequence of desi chickpea by ~30% and that of kabuli chickpea by 18%. The genetic map reported in this study represents the most dense linkage map of chickpea , with the potential to facilitate efficient anchoring of the draft genome sequences of desi as well as kabuli chickpea varieties.

Global transcriptome analysis of developing chickpea (Cicer arietinum L.) seeds.

Understanding developmental processes, especially in non-model crop plants, is extremely important in order to unravel unique mechanisms regulating development. Chickpea (C. arietinum L.) seeds are especially valued for their high carbohydrate and protein content. Therefore, in order to elucidate the mechanisms underlying seed development in chickpea, deep sequencing of transcriptomes from four developmental stages was undertaken. In this study, next generation sequencing platform was utilized to sequence the transcriptome of four distinct stages of seed development in chickpea. About 1.3 million reads were generated which were assembled into 51,099 unigenes by merging the de novo and reference assemblies. Functional annotation of the unigenes was carried out using the Uniprot, COG and KEGG databases. RPKM based digital expression analysis revealed specific gene activities at different stages of development which was validated using Real time PCR analysis. More than 90% of the unigenes were found to be expressed in at least one of the four seed tissues. DEGseq was used to determine differentially expressing genes which revealed that only 6.75% of the unigenes were differentially expressed at various stages. Homology based comparison revealed 17.5% of the unigenes to be putatively seed specific. Transcription factors were predicted based on HMM profiles built using TF sequences from five legume plants and analyzed for their differential expression during progression of seed development. Expression analysis of genes involved in biosynthesis of important secondary metabolites suggested that chickpea seeds can serve as a good source of antioxidants. Since transcriptomes are a valuable source of molecular markers like simple sequence repeats (SSRs), about 12,000 SSRs were mined in chickpea seed transcriptome and few of them were validated. In conclusion, this study will serve as a valuable resource for improved chickpea breeding.

Castleman's disease - hyaline vascular type - clinical, cytological and histological features with review of literature.

BACKGROUND AND AIMS: Castleman's disease (CD) is a rare disease of lymph node with two identified forms - the hyaline vascular (HV) type or angiofollicular type and plasma cell (PC) type. It usually presents as localized or systemic lymphadenopathy or even as extra nodal mass and may give rise to several differential diagnoses. Fine needle aspiration cytology (FNAC), as the initial investigation, may be misleading. Excision biopsy usually gives the diagnosis. This study analyzes the clinical, cytological and histological features of CD of HV type. MATERIALS AND METHODS: All the cases of CD reported between 2001 and March 2008 have been reviewed. Cytology and histopathology slides of all cases were reviewed by two reporting pathologists independently and the clinical records were analyzed. RESULT: We had five cases of CD all of which presented as cervical lymphadenopathy of 1 to 5 cms in young patients. Patients were mostly asymptomatic with all relevant investigations within normal limits. The cytology findings mostly showed a predominance of small lymphocytes with capillary fragments. On excision biopsy, all cases were diagnosed as CD of HV type. CONCLUSION: CD should be kept as a differential of lymphadenopathy, especially in asymptomatic and young patients.

Chromoblastomycosis in Nepal: a study of 13 cases.

BACKGROUND: Chromoblastomycosis is a chronic fungal infection caused by several pigmented fungi commonly seen in tropical and subtropical climates. AIM: To evaluate the epidemiologic, clinical and pathological characteristics of chromoblastomycosis in our patients. METHODS: This retrospective and prospective study was conducted at the Manipal Teaching Hospital, Pokhara, Nepal. Clinical features and histopathology of all the cases diagnosed as chromoblastomycosis during the last eight years were studied. RESULTS: A total of 13 cases of chromoblastomycosis were diagnosed during the period of 1999-2006. The disease was seen predominantly in middle-aged male farmers and those from rural areas. The lesions commonly involved the lower extremity and were single or multiple in number. They clinically presented as verrucous or nodular growths. Out of these 13 cases, three were diagnosed clinically as squamous cell carcinoma and one as psoriasis. The histopathological features included sclerotic bodies in 12 cases (92%), microabscess formation in 10 cases (76.9%), pseudoepitheliomatous hyperplasia in nine cases (69.2%) and granuloma in eight cases (61.5%). CONCLUSION: Farming is the commonest occupation in patients with chromoblastomycosis. Early histological diagnosis helps in effective management of the condition.