Toward the understanding of the role of CDC48, a major component of the protein quality control, in plant immunity

Integrated comparative physiological and transcriptomic analyses of Elymus sibiricus L. reveal the similarities and differences in the molecular mechanisms in response to drought and cold stress

Drought and cold crucially affect plant growth and distribution. Plants have evolved complex molecular mechanisms to adapt to such adverse environmental conditions. This study examines two Elymus sibiricus (Es) germplasms differing in resilience to these stresses. Analyzing physiological responses and gene expression changes under drought and cold, it reveals the similarities and differences in their molecular mechanisms that underlie these responses. The results indicate that both drought stress and cold stress severely damage the integrity of the cell membrane in Es. Notably, under cold stress, the accumulation of osmotic regulation substances in Es is more significant, which may be related to the regulation of carbohydrate metabolism (CM)-related genes in cold environments. Furthermore, the response to oxidative stress triggered by cold stress in Es is partially inhibited. The enrichment analysis showed that the DEGs responsive to drought stress in Es were mainly related to the pathway of photosynthesis, whereas the DEGs responsive to cold stress were more associated with the protein processing in endoplasmic reticulum (PPER), highlighting distinct molecular responses. In addition, we discovered that the abscisic acid (ABA) signaling transduction plays a dominant role in mediating the drought resistance mechanism of Es. We have identified 86 key candidate genes related to photosynthesis, Phst, CM, and PPER, including 5 genes that can respond to both drought and cold stress. This study provides a foundation for the molecular mechanisms underlying cold and drought resistance in Es, with insight into its future genetic improvement for stress resistance.

Cullin-Conciliated Regulation of Plant Immune Responses: Implications for Sustainable Crop Protection

Cullins are crucial components of the ubiquitin-proteasome system, playing pivotal roles in the regulation of protein metabolism. This review provides insight into the wide-ranging functions of cullins, particularly focusing on their impact on plant growth, development, and environmental stress responses. By modulating cullin-mediated protein mechanisms, researchers can fine-tune hormone-signaling networks to improve various agronomic traits, including plant architecture, flowering time, fruit development, and nutrient uptake. Furthermore, the targeted manipulation of cullins that are involved in hormone-signaling pathways, e.g., cytokinin, auxin, gibberellin, abscisic acids, and ethylene, can boost crop growth and development while increasing yield and enhancing stress tolerance. Furthermore, cullins also play important roles in plant defense mechanisms through regulating the defense-associated protein metabolism, thus boosting resistance to pathogens and pests. Additionally, this review highlights the potential of integrating cullin-based strategies with advanced biological tools, such as CRISPR/Cas9-mediated genome editing, genetic engineering, marker-associated selections, gene overexpression, and gene knockout, to achieve precise modifications for crop improvement and sustainable agriculture, with the promise of creating resilient, high-yielding, and environmentally friendly crop varieties.

P-NADs: PUX-based NAnobody degraders for ubiquitin-independent degradation of target proteins

Targeted protein degradation (TPD) allows cells to maintain a functional proteome and to rapidly adapt to changing conditions. Methods that repurpose TPD for the deactivation of specific proteins have demonstrated significant potential in therapeutic and research applications. Most of these methods are based on proteolysis targeting chimaeras (PROTACs) which link the protein target to an E3 ubiquitin ligase, resulting in the ubiquitin-based degradation of the target protein. In this study, we introduce a method for ubiquitin-independent TPD based on nanobody-conjugated plant ubiquitin regulatory X domain-containing (PUX) adaptor proteins. We show that the PUX-based NAnobody Degraders (P-NADs) can unfold a target protein through the Arabidopsis and human orthologues of the CDC48 unfoldase without the need for ubiquitination or initiating motifs. We demonstrate that P-NAD plasmids can be transfected into a human cell line, where the produced P-NADs use the endogenous CDC48 machinery for ubiquitin-independent TPD of a 143 kDa multidomain protein. Thus, P-NADs pave the road for ubiquitin-independent therapeutic TPD approaches. In addition, the modular P-NAD design combined with in vitro and cellular assays provide a versatile platform for elucidating functional aspects of CDC48-based TPD in plants and animals.

CDC48 in plants and its emerging function in plant immunity

Protein homeostasis, namely the balance between protein synthesis and degradation, must be finely controlled to ensure cell survival, notably through the ubiquitin-proteasome system (UPS). In all species, including plants, homeostasis is disrupted by biotic and abiotic stresses. A key player in the maintenance of protein balance, the protein CDC48, shows emerging functions in plants, particularly in response to biotic stress. In this review on CDC48 in plants, we detail its highly conserved structure, describe a gene expansion that is only present in Viridiplantae, discuss its various functions and regulations, and finally highlight its recruitment, still not clear, during the plant immune response.

Integrated analysis of transcriptome and small RNAome reveals regulatory network of rapid and long-term response to heat stress in Rhododendron moulmainense

MAIN CONCLUSION

The post-transcriptional gene regulatory pathway and small RNA pathway play important roles in regulating the rapid and long-term response of Rhododendron moulmainense to high-temperature stress. The Rhododendron plays an important role in maintaining ecological balance. However, it is difficult to domesticate for use in urban ecosystems due to their strict optimum growth temperature condition, and its evolution and adaptation are little known. Here, we combined transcriptome and small RNAome to reveal the rapid response and long-term adaptability regulation strategies in Rhododendron moulmainense under high-temperature stress. The post-transcriptional gene regulatory pathway plays important roles in stress response, in which the protein folding pathway is rapidly induced at 4 h after heat stress, and alternative splicing plays an important role in regulating gene expression at 7 days after heat stress. The chloroplasts oxidative damage is the main factor inhibiting photosynthesis efficiency. Through WGCNA analysis, we identified gene association patterns and potential key regulatory genes responsible for maintaining the ROS steady-state under heat stress. Finally, we found that the sRNA synthesis pathway is induced under heat stress. Combined with small RNAome, we found that more miRNAs are significantly changed under long-term heat stress. Furthermore, MYBs might play a central role in target gene interaction network of differentially expressed miRNAs in R. moulmainense under heat stress. MYBs are closely related to ABA, consistently, ABA synthesis and signaling pathways are significantly inhibited, and the change in stomatal aperture is not obvious under heat stress. Taken together, we gained valuable insights into the transplantation and long-term conservation domestication of Rhododendron, and provide genetic resources for genetic modification and molecular breeding to improve heat resistance in Rhododendron.

A fungal core effector exploits the OsPUX8B.2-OsCDC48-6 module to suppress plant immunity

Proteins containing a ubiquitin regulatory X (UBX) domain are cofactors of Cell Division Cycle 48 (CDC48) and function in protein quality control. However, whether and how UBX-containing proteins participate in host-microbe interactions remain unclear. Here we show that MoNLE1, an effector from the fungal pathogen Magnaporthe oryzae, is a core virulence factor that suppresses rice immunity by specifically interfering with OsPUX8B.2. The UBX domain of OsPUX8B.2 is required for its binding to OsATG8 and OsCDC48-6 and controls its 26 S proteasome-dependent stability. OsPUX8B.2 and OsCDC48-6 positively regulate plant immunity against blast fungus, while the high-temperature tolerance heat-shock protein OsBHT, a putative cytoplasmic substrate of OsPUX8B.2-OsCDC48-6, negatively regulates defense against blast infection. MoNLE1 promotes the nuclear migration and degradation of OsPUX8B.2 and disturbs its association with OsBHT. Given the high conservation of MoNLE1 among fungal isolates, plants with broad and durable blast resistance might be generated by engineering intracellular proteins resistant to MoNLE1.

Genomic and pangenomic analyses provide insights into the population history and genomic diversification of bottle gourd

Bottle gourd (Lagenaria siceraria (Mol.) Strandl.) is an economically important vegetable crop and one of the earliest domesticated crops. However, the population history and genomic diversification of bottle gourd have not been extensively studied. We generated a comprehensive bottle gourd genome variation map from genome sequences of 197 world-wide representative accessions, which enables a genome-wide association study for identifying genomic loci associated with resistance to zucchini yellow mosaic virus, and constructed a bottle gourd pangenome that harbors 1534 protein-coding genes absent in the reference genome. Demographic analyses uncover that domesticated bottle gourd originated in Southern Africa c. 12 000 yr ago, and subsequently radiated to the New World via the Atlantic drift and to Eurasia through the efforts of early farmers in the initial Holocene. The identified highly differentiated genomic regions among different bottle gourd populations harbor many genes contributing to their local adaptations such as those related to disease resistance and stress tolerance. Presence/absence variation analysis of genes in the pangenome reveals numerous genes including those involved in abiotic/biotic stress responses that have been under selection during the world-wide expansion of bottle gourds. The bottle gourd variation map and pangenome provide valuable resources for future functional studies and genomics-assisted breeding.

Zeocin-induced DNA damage response in barley and its dependence on ATR

DNA damage response (DDR) is an essential mechanism by which living organisms maintain their genomic stability. In plants, DDR is important also for normal growth and yield. Here, we explored the DDR of a temperate model crop barley (Hordeum vulgare) at the phenotypic, physiological, and transcriptomic levels. By a series of in vitro DNA damage assays using the DNA strand break (DNA-SB) inducing agent zeocin, we showed reduced root growth and expansion of the differentiated zone to the root tip. Genome-wide transcriptional profiling of barley wild-type and plants mutated in DDR signaling kinase ATAXIA TELANGIECTASIA MUTATED AND RAD3-RELATED (hvatr.g) revealed zeocin-dependent, ATR-dependent, and zeocin-dependent/ATR-independent transcriptional responses. Transcriptional changes were scored also using the newly developed catalog of 421 barley DDR genes with the phylogenetically-resolved relationships of barley SUPRESSOR OF GAMMA 1 (SOG1) and SOG1-LIKE (SGL) genes. Zeocin caused up-regulation of specific DDR factors and down-regulation of cell cycle and histone genes, mostly in an ATR-independent manner. The ATR dependency was obvious for some factors associated with DDR during DNA replication and for many genes without an obvious connection to DDR. This provided molecular insight into the response to DNA-SB induction in the large and complex barley genome.

Comprehensive meta-analysis of QTL and gene expression studies identify candidate genes associated with Aspergillus flavus resistance in maize

Aflatoxin (AF) contamination, caused by Aspergillus flavus, compromises the food safety and marketability of commodities, such as maize, cotton, peanuts, and tree nuts. Multigenic inheritance of AF resistance impedes conventional introgression of resistance traits into high-yielding commercial maize varieties. Several AF resistance-associated quantitative trait loci (QTLs) and markers have been reported from multiple biparental mapping and genome-wide association studies (GWAS) in maize. However, QTLs with large confidence intervals (CI) explaining inconsistent phenotypic variance limit their use in marker-assisted selection. Meta-analysis of published QTLs can identify significant meta-QTLs (MQTLs) with a narrower CI for reliable identification of genes and linked markers for AF resistance. Using 276 out of 356 reported QTLs controlling resistance to A. flavus infection and AF contamination in maize, we identified 58 MQTLs on all 10 chromosomes with a 66.5% reduction in the average CI. Similarly, a meta-analysis of maize genes differentially expressed in response to (a)biotic stresses from the to-date published literature identified 591 genes putatively responding to only A. flavus infection, of which 14 were significantly differentially expressed (-1.0 ≤ Log2Fc ≥ 1.0; p ≤ 0.05). Eight MQTLs were validated by their colocalization with 14 A. flavus resistance-associated SNPs identified from GWAS in maize. A total of 15 genes were physically close between the MQTL intervals and SNPs. Assessment of 12 MQTL-linked SSR markers identified three markers that could discriminate 14 and eight cultivars with resistance and susceptible responses, respectively. A comprehensive meta-analysis of QTLs and differentially expressed genes led to the identification of genes and makers for their potential application in marker-assisted breeding of A. flavus-resistant maize varieties.

Expanding roles for S-nitrosylation in the regulation of plant immunity

Following pathogen recognition, plant cells produce a nitrosative burst resulting in a striking increase in nitric oxide (NO), altering the redox state of the cell, which subsequently helps orchestrate a plethora of immune responses. NO is a potent redox cue, efficiently relayed between proteins through its co-valent attachment to highly specific, powerfully reactive protein cysteine (Cys) thiols, resulting in formation of protein S-nitrosothiols (SNOs). This process, known as S-nitrosylation, can modulate the function of target proteins, enabling responsiveness to cellular redox changes. Key targets of S-nitrosylation control the production of reactive oxygen species (ROS), the transcription of immune-response genes, the triggering of the hypersensitive response (HR) and the establishment of systemic acquired resistance (SAR). Here, we bring together recent advances in the control of plant immunity by S-nitrosylation, furthering our appreciation of how changes in cellular redox status reprogramme plant immune function.

Activation of CDC48 and acetyltransferase encoding genes contributes to enhanced abiotic stress tolerance and improved productivity traits in rice

World-wide crop productivity is highly impacted by various extreme environmental conditions. In the present investigation, activation tagged (AT) line A10-Ds-RFP6 of rice endowed with improved agronomic attributes was tested for its tolerance ability against drought and salinity stress conditions as well as identification of genes associated with these traits. Under both drought and salinity stress conditions, A10-Ds-RFP6 line exhibited increased seed germination rates and improved plant growth characteristics at seedling, vegetative and reproductive stages as compared to wild-type (WT) plants. Moreover, A10-Ds-RFP6 revealed effective antioxidant systems resulting in decreased accumulation of reactive oxygen species and delayed stress symptoms compared to WT plants. Reduced accumulation of malondialdehyde with concomitant increase in proline and soluble sugars in A10-Ds-RFP6 line further endorse its improved stress tolerance levels. Furthermore, A10-Ds-RFP6 disclosed enhanced plant water content, photosynthetic efficiency, stomatal conductance, water use efficiency and maximum quantum yield compared to WT plants. TAIL and qRT-PCR analyses of AT rice line revealed the integration site of Ds element in the genome and increased expression levels of CDC48 and acetyltransferase genes involved in various aspects of plant development and stress tolerance. As such, the promising AT line plausibly serve as a rare genetic resource for fortifying stress tolerance and productivity traits of elite rice cultivars.

The chaperone-like protein Cdc48 regulates ubiquitin-proteasome system in plants

The degradation of misfolded proteins is mainly mediated by the ubiquitin-proteasome system (UPS). UPS can be assisted by the protein Cdc48 but the relationship between UPS and Cdc48 in plants has been poorly investigated. Here, we analysed the regulation of UPS by Cdc48 in tobacco thanks to two independent cell lines overexpressing Cdc48 constitutively and plant leaves overexpressing Cdc48 transiently. In the cell lines, the accumulation of ubiquitinated proteins was affected both quantitatively and qualitatively and the number of proteasomal subunits was modified, while proteolytic activities were unchanged. Similarly, the over-expression of Cdc48 in planta impacted the accumulation of ubiquitinated proteins. A similar process occurred in leaves overexpressing transiently Rpn3, a proteasome subunit. Cdc48 being involved in plant immunity, its regulation of UPS was also investigated in response to cryptogein, an elicitor of immune responses. In the cell lines stably overexpressing Cdc48 and in leaves transiently overexpressing Cdc48 and/or Rpn3, cryptogein triggered a premature cell death while no increase of the proteasomal activity occurred. Overall, this study highlights a role for Cdc48 in ubiquitin homeostasis and confirms its involvement, as well as that of Rpn3, in the processes underlying the hypersensitive response.

Versatile control of the CDC48 segregase by the plant UBX-containing (PUX) proteins

In plants, AAA-adenosine triphosphatase (ATPase) Cell Division Control Protein 48 (CDC48) uses the force generated through ATP hydrolysis to pull, extract, and unfold ubiquitylated or sumoylated proteins from the membrane, chromatin, or protein complexes. The resulting changes in protein or RNA content are an important means for plants to control protein homeostasis and thereby adapt to shifting environmental conditions. The activity and targeting of CDC48 are controlled by adaptor proteins, of which the plant ubiquitin regulatory X (UBX) domain-containing (PUX) proteins constitute the largest family. Emerging knowledge on the structure and function of PUX proteins highlights that these proteins are versatile factors for plant homeostasis and adaptation that might inspire biotechnological applications.

Endoplasmic reticulum-related E3 ubiquitin ligases: Key regulators of plant growth and stress responses

Accumulating evidence has revealed that the ubiquitin proteasome system plays fundamental roles in the regulation of diverse cellular activities in eukaryotes. The ubiquitin protein ligases (E3s) are central to the proteasome system because of their ability to determine its substrate specificity. Several studies have demonstrated the essential role of a group of ER (endoplasmic reticulum)-localized E3s in the positive or negative regulation of cell homeostasis. Most ER-related E3s are conserved between plants and mammals, and a few plant-specific components have been reported. In this review, we summarize the functions of ER-related E3s in plant growth, ER-associated protein degradation and ER-phagy, abiotic and biotic stress responses, and hormone signaling. Furthermore, we highlight several questions that remain to be addressed and suggest directions for further research on ER-related E3 ubiquitin ligases.

Cross Kingdom Immunity: The Role of Immune Receptors and Downstream Signaling in Animal and Plant Cell Death

Both plants and animals are endowed with sophisticated innate immune systems to combat microbial attack. In these multicellular eukaryotes, innate immunity implies the presence of cell surface receptors and intracellular receptors able to detect danger signal referred as damage-associated molecular patterns (DAMPs) and pathogen-associated molecular patterns (PAMPs). Membrane-associated pattern recognition receptors (PRRs), such as Toll-like receptors (TLRs), C-type lectin receptors (CLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are employed by these organisms for sensing different invasion patterns before triggering antimicrobial defenses that can be associated with a form of regulated cell death. Intracellularly, animals nucleotide-binding and oligomerization domain (NOD)-like receptors or plants nucleotide-binding domain (NBD)-containing leucine rich repeats (NLRs) immune receptors likely detect effectors injected into the host cell by the pathogen to hijack the immune signaling cascade. Interestingly, during the co-evolution between the hosts and their invaders, key cross-kingdom cell death-signaling macromolecular NLR-complexes have been selected, such as the inflammasome in mammals and the recently discovered resistosome in plants. In both cases, a regulated cell death located at the site of infection constitutes a very effective mean for blocking the pathogen spread and protecting the whole organism from invasion. This review aims to describe the immune mechanisms in animals and plants, mainly focusing on cell death signaling pathways, in order to highlight recent advances that could be used on one side or the other to identify the missing signaling elements between the perception of the invasion pattern by immune receptors, the induction of defenses or the transmission of danger signals to other cells. Although knowledge of plant immunity is less advanced, these organisms have certain advantages allowing easier identification of signaling events, regulators and executors of cell death, which could then be exploited directly for crop protection purposes or by analogy for medical research.

SCFSNIPER7 controls protein turnover of unfoldase CDC48A to promote plant immunity

The unfoldase CDC48 (Cell Division Cycle 48) is highly conserved in eukaryotes, serving as an AAA + ATPase to extract ubiquitinated proteins from large protein complexes and membranes. Although its biochemical properties have been studied extensively in yeast and animal systems, the biological roles and regulations of the plant CDC48s have been explored only recently. Here we describe the identification of a novel E3 ligase from the SNIPER (snc1-influencing plant E3 ligase reverse genetic) screen, which contributes to plant defense regulation by targeting CDC48A for degradation. SNIPER7 encodes an F-box protein and its overexpression leads to autoimmunity. We identified CDC48s as interactors of SNIPER7 through immunoprecipitation followed by mass spectrometry proteomic analysis. SNIPER7 overexpression lines phenocopy the autoimmune mutant Atcdc48a-4. Furthermore, CDC48A protein levels are reduced or stabilized when SNIPER7 is overexpressed or inhibited, respectively, suggesting that CDC48A is the ubiquitination substrate of SCF^SNIPER7 . Taken together, this study reveals a new mechanism where a SCF^SNIPER7 complex regulates CDC48 unfoldase levels and modulates immune output.

Protein Quality Control in Plant Organelles: Current Progress and Future Perspectives

The endoplasmic reticulum, chloroplasts, and mitochondria are major plant organelles for protein synthesis, photosynthesis, metabolism, and energy production. Protein homeostasis in these organelles, maintained by a balance between protein synthesis and degradation, is essential for cell functions during plant growth, development, and stress resistance. Nucleus-encoded chloroplast- and mitochondrion-targeted proteins and ER-resident proteins are imported from the cytosol and undergo modification and maturation within their respective organelles. Protein folding is an error-prone process that is influenced by both developmental signals and environmental cues; a number of mechanisms have evolved to ensure efficient import and proper folding and maturation of proteins in plant organelles. Misfolded or damaged proteins with nonnative conformations are subject to degradation via complementary or competing pathways: intraorganelle proteases, the organelle-associated ubiquitin-proteasome system, and the selective autophagy of partial or entire organelles. When proteins in nonnative conformations accumulate, the organelle-specific unfolded protein response operates to restore protein homeostasis by reducing protein folding demand, increasing protein folding capacity, and enhancing components involved in proteasome-associated protein degradation and autophagy. This review summarizes recent progress on the understanding of protein quality control in the ER, chloroplasts, and mitochondria in plants, with a focus on common mechanisms shared by these organelles during protein homeostasis.

Endoplasmic reticulum-associated degradation mediated by MoHrd1 and MoDer1 is pivotal for appressorium development and pathogenicity of Magnaporthe oryzae

Most secretory proteins are folded and modified in the endoplasmic reticulum (ER); however, protein folding is error-prone, resulting in toxic protein aggregation and cause ER stress. Irreversibly misfolded proteins are subjected to ER-associated degradation (ERAD), modified by ubiquitination, and degraded by the 26S proteasome. The yeast ERAD ubiquitin ligase Hrd1p and multispanning membrane protein Der1p are involved in ubiquitination and transportation of the folding-defective proteins. Here, we performed functional characterization of MoHrd1 and MoDer1 and revealed that both of them are localized to the ER and are pivotal for ERAD substrate degradation and the ER stress response. MoHrd1 and MoDer1 are involved in hyphal growth, asexual reproduction, infection-related morphogenesis, protein secretion and pathogenicity of M. oryzae. Importantly, MoHrd1 and MoDer1 mediated conidial autophagic cell death and subsequent septin ring assembly at the appressorium pore, leading to abnormal appressorium development and loss of pathogenicity. In addition, deletion of MoHrd1 and MoDer1 activated the basal unfolded protein response (UPR) and autophagy, suggesting that crosstalk between ERAD and two other closely related mechanisms in ER quality control system (UPR and autophagy) governs the ER stress response. Our study indicates the importance of ERAD function in fungal development and pathogenesis of M. oryzae.

Comparative Transcriptome Analysis of Two Cucumber Cultivars with Different Sensitivity to Cucumber Mosaic Virus Infection

Cucumber mosaic virus (CMV), with extremely broad host range including both monocots and dicots around the world, belongs to most important viral crop threats. Either natural or genetically constructed sources of resistance are being intensively investigated; for this purpose, exhaustive knowledge of molecular virus-host interaction during compatible and incompatible infection is required. New technologies and computer-based “omics” on various levels contribute markedly to this topic. In this work, two cucumber cultivars with different response to CMV challenge were tested, i.e., sensitive cv. Vanda and resistant cv. Heliana. The transcriptomes were prepared from both cultivars at 18 days after CMV or mock inoculation. Subsequently, four independent comparative analyses of obtained data were performed, viz. mock- and CMV-inoculated samples within each cultivar, samples from mock-inoculated cultivars to each other and samples from virus-inoculated cultivars to each other. A detailed picture of CMV-influenced genes, as well as constitutive differences in cultivar-specific gene expression was obtained. The compatible CMV infection of cv. Vanda caused downregulation of genes involved in photosynthesis, and induction of genes connected with protein production and modification, as well as components of signaling pathways. CMV challenge caused practically no change in the transcription profile of the cv. Heliana. The main differences between constitutive transcription activity of the two cultivars relied in the expression of genes responsible for methylation, phosphorylation, cell wall organization and carbohydrate metabolism (prevailing in cv. Heliana), or chromosome condensation and glucan biosynthesis (prevailing in cv. Vanda). Involvement of several genes in the resistant cucumber phenotype was predicted; this can be after biological confirmation potentially applied in breeding programs for virus-resistant crops.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.