Genome-wide DNA methylation and their transgenerational pattern differ in Arabidopsis thaliana populations originated along the elevation of West Himalaya.

Methylation at 5' cytosine of DNA molecule is an important epigenetic mark. It is known to play critical role in adaptation of organisms under different biotic and abiotic stressors via modulating gene expression and/or chromatin architecture. Plant populations evolved under variable climatic conditions may have evolved different epigenetic marks including DNA methylation. Here we, describe the genome-wide DNA methylation pattern under native field, F1 and F6 generation followed by their association with phenotypes, climate and global gene expression in the three Arabidopsis thaliana populations originated at different elevation ranges of Indian West Himalaya. We show that the global methyl cytosine (mC) content is more or less similar in the three populations but differ in their distribution across genome. There was an increase in differential methylation between the populations as elevation increased. The methylation divergence was the highest between the low and the high elevation populations. The high elevation populations were hypo-methylated than the low elevation population. The methylation in the genes was associated with population specific phenotypes and climate of the region. The genes which were differentially methylated as well as differentially expressed between the low and high elevation populations were mostly related to abiotic stresses. When grown under controlled condition, there was gain of differential methylation over native condition and the maximum percent changes was observed in CHH-sequence context. Further ~ 99.8% methylated cytosines were stably passed on from F1 to F6 generation. Overall, our data suggest that high elevation population is epigenetically more plastic under changing environmental condition.Background Arabidopsis thaliana is the model plant species and has been extensively studied to understand plants life processes. There are numerous reports on its origin, demography, evolution, epigenomes and adaptation etc. however, Indian populations of Arabidopsis thaliana evolved along wide elevation ranging from ~ 700 m amsl to ~ 3400 m amsl not explored yet. Here we, describe the genome-wide DNA methylation pattern under native field, F1 and F6 generation followed by their association with phenotypes, climate and global gene expression in the three Arabidopsis thaliana populations originated at different elevation ranges of Indian West Himalaya.Results In our study we found that total mCs percent was more or less similar in the three populations but differ in their distribution across genome. The proportion of CG-mCs was the highest, followed by CHH-mCs and CHG-mCs in all the three populations. Under native field condition the methylation divergence was more prominent between low and high elevation populations and the high elevation populations were hypo-methylated than the low elevation population. The methylation in the genes was linked to population-specific phenotypes and the regional climate. The genes that showed differential methylation and expression between low and high elevation populations were primarily associated with abiotic stress responses. When grown under controlled condition, there was gain of differential methylation compared to the native condition and the maximum percent changes was observed in CHH-sequence context. Further 99.8% methylated cytosines were stably passed on from F1 to F6 generation.Conclusions The populations of A. thaliana adapted at different climatic conditions were significantly differentially methylated both under native and controlled condition. However, the magnitude and extent of gain or loss of methylation were most significant between the low and the high elevation populations. Overall, our data suggest that high elevation population is epigenetically more plastic under changing environmental condition.

Tubby-like proteins (TLPs) transcription factor in different regulatory mechanism in plants: a review.

Tubby-like proteins (TLPs) transcription factors are found in single-celled to multi-cellular eukaryotes in the form of large multigene families. TLPs are identified through a specific signature of carboxyl terminal tubby domain, required for plasma membrane tethering and amino terminal F-box domain communicate as functional SCF-type E3 ligases. The comprehensive distribution of TLP gene family members in diverse species indicates some conserved functions of TLPs in multicellular organisms. Plant TLPs have higher gene members than animals and these members reported important role in multiple physiological and developmental processes and various environmental stress responses. Although the TLPs are suggested to be a putative transcription factors but their functional mechanism is not much clear. This review provides significant recent updates on TLP-mediated regulation with an insight into its functional roles, origin and evolution and also phytohormones related regulation to combat with various stresses and its involvement in adaptive stress response in crop plants.

Identification of miRNA and their target genes in Cestrum nocturnum L. and Cestrum diurnum L. in stress responses.

MicroRNAs (miRNAs) are small, highly conserved non-coding RNA molecules and products of primary miRNAs that regulate the target gene expression. Homology-based approaches were employed to identify miRNAs and their targets in Cestrum nocturnum L. and Cestrum diurnum L. A total of 32 and 12 miRNA candidates were identified in C. nocturnum and C. diurnum. These miRNAs belong to 26 and 10 miRNA families and regulate 1024 and 1007 target genes in C. nocturnum, and C. diurnum, respectively. The functional roles of these miRNAs have not been earlier elucidated in Cestrum. MiR815a, miR849, miR1089 and miR172 have a strong propensity to target genes controlling phytochrome-interacting factor 1 (PIF1), ubiquitin-specific protease 12 (UBP12), leucine-rich repeat (LRR) protein kinase and GAI, RGA, SCR (GRAS) family transcription factor in C. nocturnum. While miR5205a, miR1436 and miR530 regulate PATATIN-like protein 6 (PLP6), PHD finger transcription factor and myb domain protein 48 (MYB48) in C. diurnum. Overall, these miRNAs have regulatory responses in biotic and abiotic stresses in both plant species. Eight putative miRNAs and their target genes were selected for qRT-PCR validation. The validated results suggested the importance of miR815a, miR849, miR5205a, miR1089, miR172, miR1436, and miR530 in exerting control over stress responses in C. nocturnum and C. diurnum. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-022-01127-1.

Genome-wide identification, phylogeny, and expression analysis of the bHLH gene family in tobacco (Nicotiana tabacum).

The basic helix-loop-helix (bHLH) is the second-largest TF family in plants that play important roles in plant growth, development, and stress responses. In this study, a total of 100 bHLHs were identified using Hidden Markov Model profiles in the Nicotiana tabacum genome, clustered into 15 major groups (I-XV) based on their conserved domains and phylogenetic relationships. Group VIII genes were found to be the most abundant, with 27 NtbHLH members. The expansion of NtbHLHs in the genome was due to segmental and tandem duplication. The purifying selection was found to have an important role in the evolution of NtHLHs. Subsequent qRT-PCR validation of five selected genes from transcriptome data revealed that NtbHLH3.1, NtbHLH3.2, NtbHLH24, NtbHLH50, and NtbHLH59.2 have higher expressions at 12 and 24 h in comparison to 0 h (control) of chilling stress. The validated results demonstrated that NtbHLH3.2 and NtbHLH24 genes have 3 and fivefold higher expression at 12 h and 2 and threefold higher expression at 24 h than control plant, shows high sensitivity towards chilling stress. Moreover, the co-expression of positively correlated genes of NtbHLH3.2 and NtbHLH24 confirmed their functional significance in chilling stress response. Therefore, suggesting the importance of NtbHLH3.2 and NtbHLH24 genes in exerting control over the chilling stress responses in tobacco. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-021-01042-x.

Chemometric Based Identification and Validation of Specific Chemical Markers for Geographical, Seasonal and Gender Variations in Tinospora cordifolia Stem using HPLC-ESI-QTOF-MS Analysis.

INTRODUCTION: Tinospora cordifolia is a widely distributed medicinal plant used in various traditional and commercial Ayurvedic formulations. Due to the wide use of this plant it is important to know the extent of variability in the metabolite profile resulting from geographical location, season and gender. OBJECTIVE: To develop a statistical approach based on phytochemical markers for confident prediction of variations in metabolic profile and cytotoxicity due to geographical, seasonal and gender difference in T. cordifolia stem. METHODS: A HPLC-ESI-QTOF-MS method was used for the metabolite profiling of T. cordifolia stem. The data were analysed using chemometric methods including Student's t-test, ANOVA, FA/PCA and ROC curve analysis and validated for the identification of chemical variations. The bioactivity of selected samples was also tested using a cell cytotoxicity assay to assess the functional aspect of the phytochemical variability. RESULTS: The chemometric approach applied here identified marker ions for geographical locations (m/z 294.1139 and 445.2136), seasons (m/z 344.1482, 359.1501, and 373.1305) and gender (m/z 257.1380) with 100% statistical sensitivity and specificity. An in vitro cytotoxicity evaluation revealed that male T. cordifolia stem was the most effective in inhibiting the growth of cancerous cell lines. CONCLUSIONS: The developed and validated chemometric approach identified the analytical markers for phytochemical variations in unknown T. cordifolia stem samples from male or female plants and samples collected from different geographical locations and seasons. The results are supported by comparative cytotoxic activity data. Copyright © 2017 John Wiley & Sons, Ltd.

Global nucleosome positioning regulates salicylic acid mediated transcription in Arabidopsis thaliana.

BACKGROUND: The nucleosome positioning regulates the gene expression and many other DNA-related processes in eukaryotes. Genome-wide mapping of nucleosome positions and correlation of genome-wide nucleosomal remodeling with the changes in the gene expression can help us understanding gene regulation on genome level. RESULTS: In the present study, we correlate the gene expression and the genomic nucleosomal remodeling in response to salicylic acid (SA) treatment in A. thaliana. We have mapped genome-wide nucleosomes by performing tiling microarray using 146 bp mononucleosomal template DNA. The average nucleosomal coverage is approximately 346 bp per nucleosome both under the control and the SA-treated conditions. The nucleosomal coverage is more in the coding region than in the 5' regulatory regions. We observe approximately 50% nucleosomal remodeling on SA treatment where significant nucleosomal depletion and nucleosomal enrichment around the transcription start site (TSS) occur in SA induced genes and SA repressed genes respectively in response to SA treatment. Especially in the case of the SA-induced group, the nucleosomal remodeling over the minimal promoter in response to SA is especially significant in the Non-expresser of PR1 (NPR1)-dependent genes. A detailed investigation of npr1-1 mutant confirms a distinct role of NPR1 in the nucleosome remodeling over the core promoter. We have also identified several motifs for various hormonal responses; including ABRE elements in the remodeled nucleosomal regions around the promoter region in the SA regulated genes. We have further identified that the W-box and TGACG/C motif, reported to play an important role in SA-mediated induction, are enriched in nucleosome free regions (NFRs) of the promoter region of the SA induced genes. CONCLUSIONS: This is the first study reporting genome-wide effects of SA treatment on the chromatin architecture of A. thaliana. It also reports significant role of NPR1 in genome-wide nucleosomal remodeling in response to SA.

Heavy metals induce oxidative stress and genome-wide modulation in transcriptome of rice root.

Industrial growth, ecological disturbances and agricultural practices have contaminated the soil and water with many harmful compounds, including heavy metals. These heavy metals affect growth and development of plants as well as cause severe human health hazards through food chain contamination. In past, studies have been made to identify biochemical and molecular networks associated with heavy metal toxicity and uptake in plants. Studies suggested that most of the physiological and molecular processes affected by different heavy metals are similar to those affected by other abiotic stresses. To identify common and unique responses by different metals, we have studied biochemical and genome-wide modulation in transcriptome of rice (IR-64 cultivar) root after exposure to cadmium (Cd), arsenate [As(V)], lead (Pb) and chromium [Cr(VI)] in hydroponic condition. We observed that root tissue shows variable responses for antioxidant enzyme system for different heavy metals. Genome-wide expression analysis suggests variable number of genes differentially expressed in root in response to As(V), Cd, Pb and Cr(VI) stresses. In addition to unique genes, each heavy metal modulated expression of a large number of common genes. Study also identified cis-acting regions of the promoters which can be determinants for the modulated expression of the genes in response to different heavy metals. Our study advances understanding related to various processes and networks which might be responsible for heavy metal stresses, accumulation and detoxification.

Unveiling the molecular mechanisms of arsenic tolerance and resilience in the primitive bryophyte Marchantia polymorpha L.

This study addresses the pressing issue of high arsenic (As) contaminations, which poses a severe threat to various life forms in our ecosystem. Despite this prevailing concern, all organisms have developed some techniques to mitigate the toxic effects of As. Certain plants, such as bryophytes, the earliest land plants, exhibit remarkable tolerance to wide range of harsh environmental conditions, due to their inherent competence. In this study, bryophytes collected from West Bengal, India, across varying contamination levels were investigated for their As tolerance capabilities. Assessment of As accumulation potential and antioxidant defense efficiency, including SOD, CAT, APX, GPX etc. revealed Marchantia polymorpha as the most tolerant species. It exhibited highest As accumulation, antioxidative proficiency, and minimal damage. Transcriptomic analysis of M. polymorpha exposed to 40 µM As(III) for 24 and 48 h identified several early responsive differentially expressing genes (DEGs) associated with As tolerance. These includes GSTs, GRXs, Hsp20s, SULTR1;2, ABCC2 etc., indicating a mechanism involving vacuolar sequestration. Interestingly, one As(III) efflux-transporter ACR3, an extrusion pump, known to combat As toxicity was found to be differentially expressed compared to control. The SEM-EDX analysis, further elucidated the operation of As extrusion mechanism, which contributes added As resilience in M. polymorpha. Yeast complementation assay using Δacr3 yeast cells, showed increased tolerance towards As(III), compared to the mutant cells, indicating As tolerant phenotype. Overall, these findings significantly enhance our understanding of As tolerance mechanisms in bryophytes. This can pave the way for the development of genetically engineered plants with heightened As tolerance and the creation of improved plant varieties.

Genome-wide identification and expression analysis of the HD2 protein family and its response to drought and salt stress in Gossypium species.

Histone deacetylase 2 (HD2) proteins play an important role in the regulation of gene expression. This helps with the growth and development of plants and also plays a crucial role in responses to biotic and abiotic stress es. HD2s comprise a C2H2-type Zn2+ finger at their C-terminal and an HD2 label, deacetylation and phosphorylation sites, and NLS motifs at their N-terminal. In this study, a total of 27 HD2 members were identified, using Hidden Markov model profiles, in two diploid cotton genomes (Gossypium raimondii and Gossypium arboretum) and two tetraploid cotton genomes (Gossypium hirsutum and Gossypium barbadense). These cotton HD2 members were classified into 10 major phylogenetic groups (I-X), of which group III was found to be the largest with 13 cotton HD2 members. An evolutionary investigation showed that the expansion of HD2 members primarily occurred as a result of segmental duplication in paralogous gene pairs. Further qRT-PCR validation of nine putative genes using RNA-Seq data suggested that GhHDT3D.2 exhibits significantly higher levels of expression at 12h, 24h, 48h, and 72h of exposure to both drought and salt stress conditions compared to a control measure at 0h. Furthermore, gene ontology, pathways, and co-expression network study of GhHDT3D.2 gene affirmed their significance in drought and salt stress responses.

Transcriptome Meta-Analysis Associated Targeting Hub Genes and Pathways of Drought and Salt Stress Responses in Cotton (Gossypium hirsutum): A Network Biology Approach.

Abiotic stress tolerance is an intricate feature controlled through several genes and networks in the plant system. In abiotic stress, salt, and drought are well known to limit cotton productivity. Transcriptomics meta-analysis has arisen as a robust method to unravel the stress-responsive molecular network in crops. In order to understand drought and salt stress tolerance mechanisms, a meta-analysis of transcriptome studies is crucial. To confront these issues, here, we have given details of genes and networks associated with significant differential expression in response to salt and drought stress. The key regulatory hub genes of drought and salt stress conditions have notable associations with functional drought and salt stress-responsive (DSSR) genes. In the network study, nodulation signaling pathways 2 (NSP2), Dehydration-responsive element1 D (DRE1D), ethylene response factor (ERF61), cycling DOF factor 1 (CDF1), and tubby like protein 3 (TLP3) genes in drought and tubby like protein 1 (TLP1), thaumatin-like proteins (TLP), ethylene-responsive transcription factor ERF109 (EF109), ETS-Related transcription Factor (ELF4), and Arabidopsis thaliana homeodomain leucine-zipper gene (ATHB7) genes in salt showed the significant putative functions and pathways related to providing tolerance against drought and salt stress conditions along with the significant expression values. These outcomes provide potential candidate genes for further in-depth functional studies in cotton, which could be useful for the selection of an improved genotype of Gossypium hirsutum against drought and salt stress conditions.

Genome-Wide Identification and Evolutionary Analysis of Gossypium Tubby-Like Protein (TLP) Gene Family and Expression Analyses During Salt and Drought Stress.

Tubby-like proteins (TLPs) possess a highly conserved closed ß barrel tubby domain at C-terminal and N-terminal F-box. The role of TLP gene family members has been widely discussed in numerous organisms; however, the detailed genome-wide study of this gene family in Gossypium species has not been reported till date. Here, we systematically identified 105 TLP gene family members in cotton (Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense) genomes and classified them into eight phylogenetic groups. Cotton TLP12 gene family members clustered into two groups, 4 and 8. They experienced higher evolutionary pressure in comparison to others, indicating the faster evolution in both diploid as well as in tetraploid cotton. Cotton TLP gene family members expanded mainly due to segmental duplication, while only one pair of tandem duplication was found in cotton TLPs paralogous gene pairs. Subsequent qRT-PCR validation of seven putative key candidate genes of GhTLPs indicated that GhTLP11A and GhTLP12A.1 genes were highly sensitive to salt and drought stress. The co-expression network, pathways, and cis-regulatory elements of GhTLP11A and GhTLP12A.1 genes confirmed their functional importance in salt and drought stress responses. This study proposes the significance of GhTLP11A and GhTLP12A.1 genes in exerting control over salt and drought stress responses in G. hirsutum and also provides a reference for future research, elaborating the biological roles of G. hirsutum TLPs in both stress responses.

Transcriptome-Guided Insights Into Plastic Degradation by the Marine Bacterium.

Polyethylene terephthalate (PET) is a common single-use plastic that accumulated in the environment because of its non-degradable characteristics. In recent years, microbes from different environments were found to degrade plastics and suggested their capability to degrade plastics under varying environmental conditions. However, complete degradation of plastics is still a void for large-scale implications using microbes because of the lack of knowledge about genes and pathways intricate in the biodegradation process. In the present study, the growth and adherence of marine Bacillus species AIIW2 on PET surface instigating structural deterioration were confirmed through weight loss and hydrophobicity reduction, as well as analyzing the change in bond indexes. The genome-wide comparative transcriptomic analysis of strain AIIW2 was completed to reveal the genes during PET utilization. The expression level of mRNA in the strain AIIW2 was indexed based on the log-fold change between the presence and absence of PET in the culture medium. The genes represent carbon metabolism, and the cell transport system was up-regulated in cells growing with PET, whereas sporulation genes expressed highly in the absence of PET. This indicates that the strain AIIW2 hydrolyzes PET and assimilated via cellular carbon metabolism. A protein-protein interaction network was built to obtain the interaction between genes during PET utilization. The genes traced to degrade PET were confirmed by detecting the hydrolytic product of PET, and genes were cloned to improve PET utilization by microbial system as an eco-friendly solution.

Genome-wide analysis of rice dehydrin gene family: Its evolutionary conservedness and expression pattern in response to PEG induced dehydration stress.

Abiotic stresses adversely affect cellular homeostasis, impairing overall growth and development of plants. These initial stress signals activate downstream signalling processes, which, subsequently, activate stress-responsive mechanisms to re-establish homeostasis. Dehydrins (DHNs) play an important role in combating dehydration stress. Rice (Oryza sativa L.), which is a paddy crop, is susceptible to drought stress. As drought survival in rice might be viewed as a trait with strong evolutionary selection pressure, we observed DHNs in the light of domestication during the course of evolution. Overall, 65 DHNs were identified by a genome-wide survey of 11 rice species, and 3 DHNs were found to be highly conserved. The correlation of a conserved pattern of DHNs with domestication and diversification of wild to cultivated rice was validated by synonymous substitution rates, indicating that Oryza rufipogon and Oryza sativa ssp. japonica follow an adaptive evolutionary pattern; whereas Oryza nivara and Oryza sativa ssp. indica demonstrate a conserved evolutionary pattern. A comprehensive analysis of tissue-specific expression of DHN genes in japonica and their expression profiles in normal and PEG (poly ethylene glycol)-induced dehydration stress exhibited a spatiotemporal expression pattern. Their interaction network reflects the cross-talk between gene expression and the physiological processes mediating adaptation to dehydration stress. The results obtained strongly indicated the importance of DHNs, as they are conserved during the course of domestication.

Genome-wide analysis of the AP2/ERF family in Musa species reveals divergence and neofunctionalisation during evolution.

AP2/ERF domain containing transcription factor super family is one of the important regulators in the plant kingdom. The involvement of AP2/ERF family members has been elucidated in various processes associated with plant growth, development as well as in response to hormones, biotic and abiotic stresses. In this study, we carried out genome-wide analysis to identify members of AP2/ERF family in Musa acuminata (A genome) and Musa balbisiana (B genome) and changes leading to neofunctionalisation of genes. Analysis identified 265 and 318 AP2/ERF encoding genes in M. acuminata and M. balbisiana respectively which were further classified into ERF, DREB, AP2, RAV and Soloist groups. Comparative analysis indicated that AP2/ERF family has undergone duplication, loss and divergence during evolution and speciation of the Musa A and B genomes. We identified nine genes which are up-regulated during fruit ripening and might be components of the regulatory machinery operating during ethylene-dependent ripening in banana. Tissue-specific expression analysis of the genes suggests that different regulatory mechanisms might be involved in peel and pulp ripening process through recruiting specific ERFs in these tissues. Analysis also suggests that MaRAV-6 and MaERF026 have structurally diverged from their M. balbisiana counterparts and have attained new functions during ripening.

Decoding regulatory landscape of somatic embryogenesis reveals differential regulatory networks between japonica and indica rice subspecies.

Somatic embryogenesis is a unique process in plants and has considerable interest for biotechnological application. Compare to japonica, indica rice has been less responsive to in vitro culture. We used Illumina Hiseq 2000 sequencing platform for comparative transcriptome analysis between two rice subspecies at six different developmental stages combined with a tag-based digital gene expression profiling. Global gene expression among different samples showed greater complexity in japonica rice compared to indica which may be due to polyphyletic origin of two rice subspecies. Expression pattern in initial stage indicate major differences in proembryogenic callus induction phase that may serve as key regulator to observe differences between both subspecies. Our data suggests that phytohormone signaling pathways consist of elaborate networks with frequent crosstalk, thereby allowing plants to regulate somatic embryogenesis pathway. However, this crosstalk varies between the two rice subspecies. Down regulation of positive regulators of meristem development (i.e. KNOX, OsARF5) and up regulation of its counterparts (OsRRs, MYB, GA20ox1/GA3ox2) in japonica may be responsible for its better regeneration and differentiation of somatic embryos. Comprehensive gene expression information in the present experiment may also facilitate to understand the monocot specific meristem regulation for dedifferentiation of somatic cell to embryogenic cells.