Global Methylome and gene expression analysis during early Peanut pod development

Genome-wide DNA methylation dynamics and RNA-seq analysis during grape (cv. 'Cabernet Franc') skin coloration

This study generated whole genome DNA methylation maps to characterize DNA methylomes of grape (cv. 'Cabernet Franc') skins and examine their functional significance during grape skin coloration. We sampled grape skin tissues at three key stages (the early stage of grape berry swelling, the late stage of grape berry swelling and the veraison) during which the color of grape berries changed from green to red. DNA methylation levels of grape skins at the three stages were higher in transposable element regions than in the genic regions, and the CG and CHG DNA methylation levels of the genic region were higher than the CHH DNA methylation levels. We identified differentially methylated regions (DMRs) in S2_vs_S1 and S3_vs_S1. The results indicated that DMRs predominantly occurred within the CHH context during grape skin coloration. Many gene ontology (GO)-enriched DMR-related genes were involved in "nucleotide binding," "catalytic activity" and "ribonucleotide binding" terms; however, many KEGG-enriched DMR-related genes were involved in the "flavonoid biosynthesis" pathway. Our results could provide an important foundation for future research on the development mechanism of grape berries.

Dynamic DNA methylation modification in peanut seed development

Cytosine methylation is an important epigenetic modification involved in regulation of plant development. However, the epigenetic mechanisms governing peanut seed development remain unclear. Herein, we generated DNA methylation profiles of developmental seeds of peanut H2014 and its smaller seed mutant H1314 at 15 and 60 days after pegging (DAP, S1, S4). Accompanying seed development, globally elevated methylation was observed in both lines. The mutant had a higher methylation level of 31.1% than wild type at S4, and 27.1-35.9% of the differentially methylated regions (DMRs) between the two lines were distributed in promoter or genic regions at both stages. Integrated methylome and transcriptome analysis revealed important methylation variations closely associated with seed development. Furthermore, some genes showed significantly negative correlation of expression with the methylation level within promoter or gene body. The results provide insights into the roles of DNA methylation in peanut seed development.

Legume-wide comparative analysis of pod shatter locus PDH1 reveals phaseoloid specificity, high cowpea expression, and stress responsive genomic context

Pod dehiscence is a major source of yield loss in legumes, which is exacerbated by aridity. Disruptive mutations in "Pod indehiscent 1" (PDH1), a pod sclerenchyma-specific lignin biosynthesis gene, has been linked to significant reductions in dehiscence in several legume species. We compared syntenic PDH1 regions across 12 legumes and two outgroups to uncover key historical evolutionary trends at this important locus. Our results clarified the extent to which PDH1 orthologs are present in legumes, showing the typical genomic context surrounding PDH1 has only arisen relatively recently in certain phaseoloid species (Vigna, Phaseolus, Glycine). The notable absence of PDH1 in Cajanus cajan may be a major contributor to its indehiscent phenotype compared with other phaseoloids. In addition, we identified a novel PDH1 ortholog in Vigna angularis and detected remarkable increases in PDH1 transcript abundance during Vigna unguiculata pod development. Investigation of the shared genomic context of PDH1 revealed it lies in a hotspot of transcription factors and signaling gene families that respond to abscisic acid and drought stress, which we hypothesize may be an additional factor influencing expression of PDH1 under specific environmental conditions. Our findings provide key insights into the evolutionary history of PDH1 and lay the foundation for optimizing the pod dehiscence role of PDH1 in major and understudied legume species.

Genome-wide characterization of phospholipase D family genes in allotetraploid peanut and its diploid progenitors revealed their crucial roles in growth and abiotic stress responses

Abiotic stresses such as cold, drought and salinity are the key environmental factors that limit the yield and quality of oil crop peanut. Phospholipase Ds (PLDs) are crucial hydrolyzing enzymes involved in lipid mediated signaling and have valuable functions in plant growth, development and stress tolerance. Here, 22, 22 and 46 PLD genes were identified in Arachis duranensis, Arachis ipaensis and Arachis hypogaea, respectively, and divided into α, β, γ, δ, ε, ζ and φ isoforms. Phylogenetic relationships, structural domains and molecular evolution proved the conservation of PLDs between allotetraploid peanut and its diploid progenitors. Almost each A. hypogaea PLD except for AhPLDα6B had a corresponding homolog in A. duranensis and A. ipaensis genomes. The expansion of Arachis PLD gene families were mainly attributed to segmental and tandem duplications under strong purifying selection. Functionally, the most proteins interacting with AhPLDs were crucial components of lipid metabolic pathways, in which ahy-miR3510, ahy-miR3513-3p and ahy-miR3516 might be hub regulators. Furthermore, plenty of cis-regulatory elements involved in plant growth and development, hormones and stress responses were identified. The tissue-specific transcription profiling revealed the broad and unique expression patterns of AhPLDs in various developmental stages. The qRT-PCR analysis indicated that most AhPLDs could be induced by specific or multiple abiotic stresses. Especially, AhPLDα3A, AhPLDα5A, AhPLDβ1A, AhPLDβ2A and AhPLDδ4A were highly up-regulated under all three abiotic stresses, whereas AhPLDα9A was neither expressed in 22 peanut tissues nor induced by any abiotic stresses. This genome-wide study provides a systematic analysis of the Arachis PLD gene families and valuable information for further functional study of candidate AhPLDs in peanut growth and abiotic stress responses.

The pattern of DNA methylation alteration, and its association with the expression changes of non-coding RNAs and mRNAs in Moso bamboo under abiotic stress

Epigenetic changes play an important role in plant growth and development and in stress response. However, DNA methylation pattern and its relationship with the expression changes of non-coding RNAs and mRNAs of Moso bamboo in response to abiotic stress is still largely unknown. In this work, we used whole-genome bisulfite sequencing in combination with whole-transcriptome sequencing to analyze the DNA methylation and transcription patterns of mRNAs and non-coding RNAs in Moso bamboo under abiotic stresses such as cold, heat, ultraviolet (UV) and salinity. We found that CHH methylation in the promoter region was positively correlated with gene expression, while CHG and CHH methylations in the gene body regions were negatively associated with gene expression. Moreover, CG and CHG methylations in the promoter regions were negatively correlated with the transcript abundance of long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs). Similarly, the methylation levels of three contexts in the genic regions were negatively correlated with the transcript abundance of lncRNAs and miRNAs but positively correlated with that of circRNAs. In addition, we suggested that the reduction of 21-nt and 24-nt small interfering RNA (siRNA) expression tended to increase methylation levels in the genic regions. We found that stress-responsive genes such as CRPK1, HSFB2A and CIPK were differentially methylated and expressed. Our results also proposed that DNA methylation may regulate the expression of the transcription factors (TFs) and plant hormone signalling genes such as IAA9, MYC2 and ERF110 in response to abiotic stress. This study firstly reports the abiotic stress-responsive DNA methylation pattern and its involvement of expression of coding RNAs and non-coding RNAs in Moso bamboo. The results expand the knowledge of epigenetic mechanisms in Moso bamboo under abiotic stress and support in-depth deciphering of the function of specific non-coding RNAs in future studies.

Genome-wide analyses of the Nodulin-like gene family in bread wheat revealed its potential roles during arbuscular mycorrhizal symbiosis

Nodulin-like (NL) genes are involved in transporting of various substances and may play key roles during the establishment of symbiosis in legumes plants. However, basic biological information of NL genes in the wheat genome is still largely unknown. Here, we identified and characterized NL genes in wheat via integrating genomic information, collinearity analysis, co-expression network analysis (WGCNA) and transcriptome analysis. In addition, we analyzed the polymorphisms and the roles of NL genes during arbuscular mycorrhizal (AM) symbiosis using a large wheat panel consists of 259 wheat genotypes. We identified 181 NL genes in the wheat genome, which were classified into SWEET, Early Nodulin-Like (ENODL), Major Facilitator Superfamily-Nodulin (MFS), Vacuolar Iron Transporter (VIT) and Early nodulin 93 (ENOD93) subfamily. The expansion of NL genes was mainly driven by segmental duplication. The bHLH genes are potential unrecognized transcription factors regulating NL genes. Moreover, two NL genes were more sensitive than other NL genes to AM colonization. The polymorphisms of NL genes are mainly due to random drift, and the natural mutation of NL genes led to significant differences in the mycorrhizal dependence of wheat in phosphorus uptake. The results concluded that NL genes potentially play important roles during AM symbiosis with wheat.

Deep transcriptomic study reveals the role of cell wall biosynthesis and organization networks in the developing shell of peanut pod

BACKGROUND

Peanut (Arachis hypogaea L.) belongs to an exceptional group of legume plants, wherein the flowers are produced aerially, but the pods develop under the ground. In such a unique environment, the pod's outer shell plays a vital role as a barrier against mechanical damage and soilborne pathogens. Recent studies have reported the uniqueness and importance of gene expression patterns that accompany peanut pods' biogenesis. These studies focused on biogenesis and pod development during the early stages, but the late developmental stages and disease resistance aspects still have gaps. To extend this information, we analyzed the transcriptome generated from four pod developmental stages of two genotypes, Hanoch (Virginia-type) and IGC53 (Peruvian-type), which differs significantly in their pod shell characteristics and pathogen resistance.

RESULTS

The transcriptome study revealed a significant reprogramming of the number and nature of differentially expressed (DE) genes during shell development. Generally, the numbers of DE genes were higher in IGC53 than in Hanoch, and the R5-R6 transition was the most dynamic in terms of transcriptomic changes. Genes related to cell wall biosynthesis, modification and transcription factors (TFs) dominated these changes therefore, we focused on their differential, temporal and spatial expression patterns. Analysis of the cellulose synthase superfamily identified specific Cellulose synthase (CesAs) and Cellulose synthase-like (Csl) genes and their coordinated interplay with other cell wall-related genes during the peanut shell development was demonstrated. TFs were also identified as being involved in the shell development process, and their pattern of expression differed in the two peanut genotypes. The shell component analysis showed that overall crude fiber, cellulose, lignin, hemicelluloses and dry matter increased with shell development, whereas K, N, protein, and ash content decreased. Genotype IGC53 contained a higher level of crude fiber, cellulose, NDF, ADF, K, ash, and dry matter percentage, while Hanoch had higher protein and nitrogen content.

CONCLUSIONS

The comparative transcriptome analysis identified differentially expressed genes, enriched processes, and molecular processes like cell wall biosynthesis/modifications, carbohydrate metabolic process, signaling, transcription factors, transport, stress, and lignin biosynthesis during the peanut shell development between two contrasting genotypes. TFs and other genes like chitinases were also enriched in peanut shells known for pathogen resistance against soilborne major pathogens causing pod wart disease and pod damages. This study will shed new light on the biological processes involved with underground pod development in an important legume crop.

A novel DNA biosensor for the ultrasensitive detection of DNA methyltransferase activity based on a high-density "hot spot" SERS substrate and rolling circle amplification strategy

Herein, we proposed a novel biosensor based on a high-density "hot spot" Au@SiO₂ array substrate and rolling circle amplification (RCA) strategy for the ultrasensitive detection of CpG methyltransferase (M.SssI) activity. In the presence of M.SssI, the RCA process can be triggered, causing the augmentation of the single-stranded DNA (ssDNA) at the tail of the double-stranded DNA (dsDNA), and the ssDNA can be hybridized with numerous DNA probes labeled with Raman reporters in the next steps. Afterwards, the resultant ssDNA can be modified to the Au@SiO₂ array substrate with the SERS enhancement factor of 7.49 × 10⁶. The substrate was synthesized by using a monolayer SiO₂ array to pick up the Au nanoparticle (AuNP) array and finite-difference time-domain (FDTD) simulation showed its excellent SERS effect. Particularly, the developed biosensor displayed a significant sensitivity with a broad detection range covering from 0.005 to 50 U mL^-1, and the limits of detection (LODs) in PBS buffer and human serum were 2.37 × 10^-4 U mL^-1 and 2.51 × 10^-4 U mL^-1, respectively. Finally, in order to verify the feasibility of its clinical application, the serum samples of healthy subjects and breast cancer, prostate cancer, gastric cancer and cervical cancer patients were analyzed, and the reliability of the results was also confirmed by western blot (WB) experiments. Taking advantage of these merits, the proposed biosensor can be a very promising alternative tool for the detection of M.SssI activity, which is of vital importance in the early detection and prevention of tumors.

Comprehensive Transcriptome Analyses Reveal Candidate Genes for Variation in Seed Size/Weight During Peanut (Arachis hypogaea L.) Domestication

Seed size/weight, a key domestication trait, is also an important selection target during peanut breeding. However, the mechanisms that regulate peanut seed development are unknown. We re-sequenced 12 RNA samples from developing seeds of two cultivated peanut accessions (Lines 8106 and 8107) and wild Arachis monticola at 15, 30, 45, and 60 days past flowering (DPF). Transcriptome analyses showed that ∼36,000 gene loci were expressed in each of the 12 RNA samples, with nearly half exhibiting moderate (2 ≤ FPKM < 10) expression levels. Of these genes, 12.2% (4,523) were specifically expressed during seed development, mainly at 15 DPF. Also, ∼12,000 genes showed significant differential expression at 30, 45, and/or 60 DPF within each of the three peanut accessions, accounting for 31.8-34.1% of the total expressed genes. Using a method that combined comprehensive transcriptome analysis and previously mapped QTLs, we identified several candidate genes that encode transcription factor TGA7, topless-related protein 2, IAA-amino acid hydrolase ILR1-like 5, and putative pentatricopeptide repeat-containing (PPR) protein. Based on sequence variations identified in these genes, SNP markers were developed and used to genotype both 30 peanut landraces and a genetic segregated population, implying that EVM0025654 encoding a PPR protein may be associated with the increased seed size/weight of the cultivated accessions in comparison with the allotetraploid wild peanut. Our results provide additional knowledge for the identification and functional research into candidate genes responsible for the seed size/weight phenotype in peanut.

Alternative splicing profiling provides insights into the molecular mechanisms of peanut peg development

BACKGROUND

The cultivated peanut (Arachis hypogaea) is one of the most important oilseed crops worldwide, and the generation of pegs and formation of subterranean pods are essential processes in peanut reproductive development. However, little information has been reported about alternative splicing (AS) in peanut peg formation and development.

RESULTS

Herein, we presented a comprehensive full-length (FL) transcriptome profiling of AS isoforms during peanut peg and early pod development. We identified 1448, 1102, 832, and 902 specific spliced transcripts in aerial pegs, subterranean pegs, subterranean unswollen pegs, and early swelling pods, respectively. A total of 184 spliced transcripts related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins were identified as possibly involved in peanut peg development. For aerial pegs, spliced transcripts we got were mainly involved in gravity stimulation and cell wall morphogenetic processes. The genes undergoing AS in subterranean peg were possibly involved in gravity stimulation, cell wall morphogenetic processes, and abiotic response. For subterranean unswollen pegs, spliced transcripts were predominantly related to the embryo development and root formation. The genes undergoing splice in early swelling pods were mainly related to ovule development, root hair cells enlargement, root apex division, and seed germination.

CONCLUSION

This study provides evidence that multiple genes are related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins undergoing AS express development-specific spliced isoforms or exhibit an obvious isoform switch during the peanut peg development. AS isoforms in subterranean pegs and pods provides valuable sources to further understand post-transcriptional regulatory mechanisms of AS in the generation of pegs and formation of subterranean pods.

Comprehensive identification and analysis of DELLA genes throughout the plant kingdom

BACKGROUND

DELLAs play key roles in plant gibberellin signaling pathways and are generally important in plant development and growth. However, DELLAs in many plant taxa have not yet been systematically analyzed.

RESULTS

In our study, we searched for DELLA genes across 58 green plant genomes and found 181 DELLAs. Structure analysis showed some DELLA domains do not contain "D-E-L-L-A" sequences and instead contain similar domains, including DGLLA and DSLLH domains. "VHYNP" motifs in plant DELLAs comprise 23 types of sequences, while some DELLAs did not contain GRAS domains. In grape, we found that the DELLA protein GSVIVT01015465001 contains an F-box domain, while apple DELLA proteins MDP0000220512 and MDP0000403162 contain a WW domain and a BCIP domain, respectively. These DELLAs can be divided into 22 homologous groups and 17 orthologous groups, and 35 paralogous genes were identified. In total, 35 positively selected genes (PSGs) and 121 negatively selected genes (NSGs) were found among DELLAs based on selective pressure analysis, with an average K_s of NSGs that was significantly higher than that of PSGs (P < 0.05). Among the paralogous groups, CBI and Fop were significantly positively correlated with GC, GC1, GC2, GC12, and GC3, while CAI was significantly positively correlated with GC, GC1, GC12, and GC. The paralogous groups with ω values exceeding 1 had significantly higher K_a values. We also found some paralogous groups with ω values exceeding 1 that differed in their motifs.

CONCLUSIONS

This study provides helpful insights into the evolution of DELLA genes and offers exciting opportunities for the investigation of DELLA functions in different plants.

Genetic mechanisms of salt stress responses in halophytes

Abiotic stress is a major threat to plant growth and development, resulting in extensive crop loss worldwide. Plants react to abiotic stresses through physiological, biochemical, molecular, and genetic adaptations that promote survival. Exploring the molecular mechanisms involved in abiotic stress responses across various plant species is essential for improving crop yields in unfavorable environments. Halophytes are characterized as plants that survive to reproduce in soils containing high salt concentrations, and thus act as an ideal model to comprehend complicated genetic and physiological mechanisms of salinity stress tolerance. Plant ecologists classify halophytes into three main groups: euhalophytes, recretohalophytes, and pseudo-halophytes. Recent genetic and molecular research has showed complicated regulatory networks by which halophytes coordinate stress adaptation and tolerance. Furthermore, investigation of natural variations in these stress responses has supplied new perspectives on the evolution of mechanisms that regulate tolerance and adaptation. This review discusses the current understanding of the genetic mechanisms that contribute to salt-stress tolerance among different classes of halophytes.